block_index: Find independent blocks of equations.

In lintools: Manipulation of Linear Systems of (in)Equalities

Find independent blocks of equations.

Description

Find independent blocks of equations.

Usage

block_index(A, eps = 1e-08)

Arguments

A	[numeric] Matrix
eps	[numeric] Coefficients with absolute value < eps are treated as zero.

Value

A list containing numeric vectors, each vector indexing an independent block of rows in the system Ax <= b.

Examples

 A <- matrix(c(
   1,0,2,0,0,
   3,0,4,0,0,
   0,5,0,6,7,
   0,8,0,0,9
 ),byrow=TRUE,nrow=4)
 b <- rep(0,4)
 bi <- block_index(A)
 lapply(bi,function(ii) compact(A[ii,,drop=FALSE],b=b[ii])$A)

lintools documentation built on Jan. 17, 2023, 1:06 a.m.