Hello,
I’ve been trying to solve a sparse matrix following the example given in “pardiso_sym.c”. However, I’ve noticed that if you back substitute the solution obtained from that example, the RHS are not found again. i.e., if I’m trying to solve A*x = b, I know A and b, pardiso solves for x, I plug x back into the equation and should obtain b, however, that is not the case!
Can anyone help me with this? Or has encountered this problem before?
Kind regards,
Ben Jakobs
pardiso_sym.c (7.6 KB)