I am testing an example with sparse right hand side and sparse solution iparm[30 ]=1 in my C++.
What I have noticed that
the solution of (perm = 0, b !=0) and (perm, b=0) are different.
But I expected that the solution should be the same since perm= 0,
Furthermore the solution x is 0th entry of the solution vectro.
However, from the pardiso manul of iparm(31), x should be a random number, since I did't request to caculate it.
The attachment is my C++ code.
Could you tell me why this happened?