<div dir="ltr">Hi,<div><br></div><div>  in my works, I'm computing a formula (involving an L2 projection) of the form:</div><div><br></div><div>(D - M^{-1} N) n</div><div><br></div><div>where D,M and B are matrices and n is a vector.</div><div><br></div><div>Operator (D - M^{-1} N) appears to be the Schur complement of the system</div><div>( M   N )</div><div>(  I    D )</div><div>and it is convenient for me to use a Mat of type schurcomplement to compute it.</div><div>However, my matrix D is not a square matrix (this system is not invertible, it is rectangular).</div><div><br></div><div>Assertions in src/ksp/ksp/utils/schurm.c prohibit the use of schurcomplement for rectangular system, but I can't figure where it is needed by schurcomplement itself.</div><div>Could it be possible to remove those two assertions (see my patch)?</div><div><br></div><div><br></div><div>Patrick Lacasse</div></div>