Dear Numpy users,<br><br>I have to solve for Z in the following equation Z^(-1) = A^(-1) - B^(-1),<br>where A and B are covariance matrices with zero determinant.<br><br>I have never used pseudoinverse matrixes, could anybody please point to me any cautions I have to take when solving this equation for Z? The brute force approach linalg.pinv( linalg.pinv(A) - lingal.pinv(B) ) gives me a matrix with all entries equal to &#39;infinity&#39;.<br>

<br>-Jose Borreguero<br>