Hello everyone

Recently I have started to use Eigen library. I find it very good, but I have found that sparse solvers return not correct solution.
The rank of the linear system of equations is about 1600. I have tried SparseLU, BiCGSTAB, SimplicialCholesky.
I have checked the solution with Blas and Matlab.
Can someone help me with this problem?
How can I attach the files with A and B matrices (A x = B)?

Code: Select all

#include <unsupported/Eigen/SparseExtra>
...
Eigen::saveMarket(A, "filename.mtx", Eigen::Symmetric);
Eigen::saveMarketVector(B, "filename_b.mtx");


You might also try the benchmark routine: http://eigen.tuxfamily.org/dox/group__T ... tml#title3

I am sorry, I have made mistake in my code.