<div dir="ltr"><div><div><div><div><div>Hi,<br><br></div>I need to compute the following in petsc;<br><br></div>P = I -(Z * inv(Z'*A*Z) * Z'<br><br></div>What is the best way to compute this?<br><br></div>I am planing to use MatPtAP in order to compute the expression inside of the inverse. And then use<br>MatMatSolve instead of inv and comtinue with MatMatMult and use MatShift finally.<br><br></div>Any suggestions or comments are welcome..<br></div>