enable testing of complex numbers for taucs

This commit is contained in:
Gael Guennebaud 2009-05-12 13:43:40 +00:00
parent 159c99a288
commit 877c3c00a2

View File

@ -122,20 +122,18 @@ template<typename Scalar> void sparse_solvers(int rows, int cols)
SparseLLT<SparseSelfAdjointMatrix,Cholmod>(m2).solveInPlace(x); SparseLLT<SparseSelfAdjointMatrix,Cholmod>(m2).solveInPlace(x);
VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LLT: cholmod"); VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LLT: cholmod");
#endif #endif
if (!NumTraits<Scalar>::IsComplex)
{ #ifdef EIGEN_TAUCS_SUPPORT
#ifdef EIGEN_TAUCS_SUPPORT x = b;
x = b; SparseLLT<SparseSelfAdjointMatrix,Taucs>(m2,IncompleteFactorization).solveInPlace(x);
SparseLLT<SparseSelfAdjointMatrix,Taucs>(m2,IncompleteFactorization).solveInPlace(x); VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LLT: taucs (IncompleteFactorization)");
VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LLT: taucs (IncompleteFactorization)"); x = b;
x = b; SparseLLT<SparseSelfAdjointMatrix,Taucs>(m2,SupernodalMultifrontal).solveInPlace(x);
SparseLLT<SparseSelfAdjointMatrix,Taucs>(m2,SupernodalMultifrontal).solveInPlace(x); VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LLT: taucs (SupernodalMultifrontal)");
VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LLT: taucs (SupernodalMultifrontal)"); x = b;
x = b; SparseLLT<SparseSelfAdjointMatrix,Taucs>(m2,SupernodalLeftLooking).solveInPlace(x);
SparseLLT<SparseSelfAdjointMatrix,Taucs>(m2,SupernodalLeftLooking).solveInPlace(x); VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LLT: taucs (SupernodalLeftLooking)");
VERIFY(refX.isApprox(x,test_precision<Scalar>()) && "LLT: taucs (SupernodalLeftLooking)"); #endif
#endif
}
} }
// test LDLT // test LDLT