complement package:lapack
For an m-by-n-matrix a with m>=n this function computes an
m-by-(m-n)-matrix b such that Matrix.multiply (adjoint b)
a is a zero matrix. The function does not try to compensate a
rank deficiency of a. That is, a|||b has full rank
if and only if a has full rank.
For full-rank matrices you might also call this kernel or
nullspace.
schurComplement :: (Solve typ11, SolveExtra typ11 xl11, SolveExtra typ11 xu11, Floating a) => Quadratic (Square TypeFull measOff vertOff horizOff typ11) ((), xl11, (TypeFull, (), ())) ((), xu11, (TypeFull, (), ())) Filled Filled (sh0 ::+ sh1) a -> Square sh0 a