Hi,

I've run into the following issue: The SimplicialLLT solver for the tridiagonal matrix given below does not give the right result, even though Eigen::Success is returned. I'm using the 3.1.0-rc1 version.

Code: Select all

A =
1 0 0 0 0 0 0 0 0 0 0
-1 2 -1 0 0 0 0 0 0 0 0
0 -1 2 -1 0 0 0 0 0 0 0
0 0 -1 2 -1 0 0 0 0 0 0
0 0 0 -1 2 -1 0 0 0 0 0
0 0 0 0 -1 2 -1 0 0 0 0
0 0 0 0 0 -1 2 -1 0 0 0
0 0 0 0 0 0 -1 2 -1 0 0
0 0 0 0 0 0 0 -1 2 -1 0
0 0 0 0 0 0 0 0 -1 2 -1
0 0 0 0 0 0 0 0 0 0 1


Roman

That matrix is not symmetric: A(0,1) = 0 and A(1,0) = -1.

You are right... Totally missed that.

Roman