inverse()

template<class RingType>

bool DMatLinAlg< RingType >::inverse

(

Mat &

X

)

Output: X, a matrix, the inverse of A returns false iff A is (near)singular

Definition at line 564 of file dmat-lu.hpp.

{
  const Mat& LU = mLUObject.LUinPlace();
 
  // Make the identity matrix
  // printf("in DMatLinAlg<RingType>::inverse\n");
 
  Mat id(ring(), LU.numRows(), LU.numRows());
  for (size_t i = 0; i < LU.numRows(); i++)
    ring().set_from_long(id.entry(i, i), 1);
 
  solve(id, X);
  return true;
}

References mLUObject, ring(), and solve().

Referenced by MatrixOps::inverse().