inverse()

bool ffpackInterface::inverse	(	const DMatZZpFFPACK &	A,
		DMatZZpFFPACK &	result_inv )

Definition at line 113 of file dmat.cpp.

{
  assert(mat.numRows() == mat.numColumns());
  result_inv.resize(mat.numRows(), mat.numRows());
 
  assert(result_inv.numRows() == mat.numRows());
  assert(result_inv.numColumns() == mat.numRows());
 
  if (mat.numRows() == 0)
    {
      // 26 April 2014: this branch is needed as FFPACK gives answer of 0 in
      // this case.
      return true;
    }
 
  DMatZZpFFPACK N(mat);
  size_t n = mat.numRows();
  int nullspacedim;
  FFPACK::Invert2(
      mat.ring().field(), n, N.rowMajorArray(), n, result_inv.rowMajorArray(), n, nullspacedim);
  return (nullspacedim == 0);
}

References DMatZZpFFPACK.

Referenced by DMatLinAlg< M2::ARingZZpFFPACK >::inverse().