determinant()

template<class RingType>

void DMatLinAlg< RingType >::determinant

(

ElementType &

result

)

Output: result, the determinant of A.

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

{
  const Mat& LU = mLUObject.LUinPlace();
 
  // This is just the product of the diagonal entries of mLU.
  assert(LU.numRows() == LU.numColumns());
 
  if (mLUObject.signOfPermutation())
    ring().set_from_long(result, 1);
  else
    ring().set_from_long(result, -1);
  for (size_t i = 0; i < LU.numRows(); i++)
    ring().mult(result, result, LU.entry(i, i));
}

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

Referenced by SLEvaluatorConcrete< RT >::computeNextNode(), and MatrixOps::determinant().