matrixPLU()

template<class RingType>

bool DMatLinAlg< RingType >::matrixPLU	(	std::vector< size_t > &	P,
		Mat &	L,
		Mat &	U )

Output: P = mPerm, L,U; the corresponding parts mLU returns false iff there is an error

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

{
  const Mat& LU = mLUObject.LUinPlace();
 
  //  std::cout << "calling matrixPLU generic\n";  
  setUpperLower(LU, L, U);
  P = mLUObject.permutation();
  return mLUObject.signOfPermutation();
}

References mLUObject, and setUpperLower().

Referenced by MatrixOps::LU().