How to use the pardiso for partial solution with correct iparm(31) value?

Dear All

I am testing the example given in the MKL manual for pardiso with partial solution in Fortran.

As stated for iparm(31), when its value is 3, selected components of the solution vectors will be computed.  However, I tried many times, with different perm values, iparm(31) values(=1, 2 or 3), and paramater locations,  it always give the whole solution.

Since for my real problem(with thousands of DOF), only a few elements in rhs are nonzero and  only a few selected components of the solution are needed. I have to realize the partial solution function to cut the time cost.

I hope you can help me. Thanks a lot.

The attachment is my fortran code.

Jim Cheng