Incomplete Cholesky preconditioner... not yet stable
diff --git a/bench/spbench/sp_solver.cpp b/bench/spbench/sp_solver.cpp index e18f2d1..a1f4bac 100644 --- a/bench/spbench/sp_solver.cpp +++ b/bench/spbench/sp_solver.cpp
@@ -13,7 +13,7 @@ #include <Eigen/SuperLUSupport> // #include <unsupported/Eigen/src/IterativeSolvers/Scaling.h> #include <bench/BenchTimer.h> - +#include <unsupported/Eigen/IterativeSolvers> using namespace std; using namespace Eigen; @@ -26,7 +26,8 @@ VectorXd b, x, tmp; BenchTimer timer,totaltime; //SparseLU<SparseMatrix<double, ColMajor> > solver; - SuperLU<SparseMatrix<double, ColMajor> > solver; +// SuperLU<SparseMatrix<double, ColMajor> > solver; + ConjugateGradient<SparseMatrix<double, ColMajor>, Lower,IncompleteCholesky<double,Lower> > solver; ifstream matrix_file; string line; int n;