compositions (version 2.0-0)

gsi.svdsolve: Internal function: Solves singular and non square equation systems

Description

Based on the singular value decomposition, a singular equation system ax=b is solved.

Usage

gsi.svdsolve(a,b,...,cond=1E-10)

Arguments

a

the matrix of ax=b (a.k.a. left-hand side matrix)

b

the vector or matrix b of ax=b (a.k.a right-hand side, independent element)

cond

the smallest-acceptable condition of the matrix. Smaller singular values are truncate

…

additional arguments to svd

Value

The "smallest" vector or matrix solving this system with minimal joint error among all vectors.

Examples

# NOT RUN {
#A <- matrix(c(0,1,0,0,0,0),ncol=2)
#b <- diag(3)
#erg <- gsi.svdsolve(A,b)
#erg
#A %*% erg 
#diag(c(0,1,0))  # richtig
# }