Hi,I found the part of the pardiso manual explaining "iparm(31)- partial solution for sparse right-hand sides and sparse solution" confusing. You can find the manual here:http://software.intel.com/sites/products/documentation/hpc/mkl/mklman/index.htm1) There are emphasises like "the right hand side must be sparse". I suppose this does not mean structurally sparse because the right hand side is passed into pardiso using dense matrix.2) When iparm(31) = 1, "You can setperm(i)= 1 only if thei-th component of the right hand side is nonzero." Can I trick the solver by setting the i-th component of the right hand side to 1e-20? I don't understand what the logic is for this constraint.3) "To useiparm(31)=2, you must set thei-th component of the right hand side to zero explicitly ifperm(i)is not equal to 1." This is related to question 2. And this statement sounds like instructing me to change my equation to be solved because I "must set thei-th component of the right hand side to zero explicitly"4) "Ifiparm(31)=3, the right hand side can be of any type and you must setperm(i)= 1 to compute thei-th component in the solution vector." What "type" do you mean the 'right hand side' can be?5) "The permutation vectorpermmust be present in all phases of Intel MKL PARDISO software. At the reordering step, the software overwrites the input vectorpermby a permutation vector used by the software at the factorization and solver step. Ifmis the number of components such thatperm(i)= 1, then the lastmcomponents of the output vectorpermare a set of the indicesisatisfying the conditionperm(i)= 1 on input."I have to admit that I am not an experienced user of PARDISO. By reading the above discription, I don't know how to set perm for the partial solve purpose because it will be overwriten by PARDISO at the factorization and solver step.What is the purpose of "Ifmis the number of components such thatperm(i)= 1, then the lastmcomponents of the output vectorpermare a set of the indicesisatisfying the conditionperm(i)= 1 on input"?Thanks in advance for any help in clarification!!Regards,Xin
For more complete information about compiler optimizations, see our Optimization Notice.