SVD()

template<typename T>

bool MutableMat< T >::SVD ( MutableMatrix * Sigma, MutableMatrix * U, MutableMatrix * Vt, bool use_divide_and_conquer ) const

virtual

Implements MutableMatrix.

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

{
  auto Sigma2 = Sigma->coerce<DMat<HermitianEigenvalueType> >();
  T* U2 = U->coerce<T>();
  T* Vt2 = Vt->coerce<T>();
  if (Sigma2 == 0 || U2 == 0 || Vt2 == 0)
    throw exc::engine_error("expected matrices of the same type");
  int strategy = (use_divide_and_conquer ? 1 : 0);
  return MatrixOps::SVD(mat, *Sigma2, *U2, *Vt2, strategy);
}

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