Dear All
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] = 0, b[0] !=0) and (perm[0], b[0]=0) are different.
But I expected that the solution should be the same since perm[0]= 0,
Furthermore the solution x[0] is 0th entry of the solution vectro.
However, from the pardiso manul of iparm(31), x[0] should be a random number, since I did't request to caculate it.
http://software.intel.com/sites/products/documentation/hpc/mkl/mklman/GU...
The attachment is my C++ code.
Could you tell me why this happened?
Hailong
