LU()

template<typename T>

M2_arrayintOrNull MutableMat< T >::LU ( MutableMatrix * L, MutableMatrix * U ) const

virtual

Implements MutableMatrix.

Definition at line 222 of file mutablemat-imp.hpp.

{
  //  std::cout << "MutableMat<T>::LU\n";
  T* L1 = L->coerce<T>();
  T* U1 = U->coerce<T>();
  if (L1 == 0 or U1 == 0)
    throw exc::engine_error("expected matrices of the same ring/type");
  return MatrixOps::LU(mat, *L1, *U1);
}

References MutableMatrix::coerce(), MatrixOps::LU(), mat, MutableMatrix::MutableMatrix(), and T.