! Pardiso needs X(N) even if sol is returned in b
Solution is not returned in the b matrix from your non-symmetric example, is there a variable I have to set to do this
Have you ever considered hooking the PARDISO straight in front of EIGEN and taking out the need for a middle man so to speak
Please try to set parameter iparm = 1 (iparm(6) in Fortran notation), which means 'The solver stores the solution on the right-hand side .'