https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/PARDISO-on-multiple-systems-with-the-same-sparsity-pattern/m-p/765164#M182 Hi,<DIV>Nonsymmetric permutation and scaling matrices are calculated during reordering step of PARDISO, so it's strongly recommended to calculate it for each stiffness matrix independently.</DIV><DIV>With best regards,</DIV><DIV>Alexander Kalinkin</DIV> Tue, 21 Aug 2012 06:56:26 GMT Alexander_K_Intel2 2012-08-21T06:56:26Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/PARDISO-on-multiple-systems-with-the-same-sparsity-pattern/m-p/765163#M181 <P>I need to use PARDISO on multiple&nbsp;(millions) structurally-symmetric&nbsp;systems with the same sparsity pattern. These systems are relatively small (only a few hundred unknowns per system) but tend to be wildly different in conditioning. I use the default parameters (and mtype=11 because mtype=1 gives bad results) and solve the systems with a single phase=11 run and multiple phase=23 runs. I want to know if the same nonsymmetric permutation and scaling matrices are used for all systems and, if that is the case, whether that will has any impact on numerical accuracy.</P> Sat, 18 Aug 2012 00:02:35 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/PARDISO-on-multiple-systems-with-the-same-sparsity-pattern/m-p/765163#M181 styc 2012-08-18T00:02:35Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/PARDISO-on-multiple-systems-with-the-same-sparsity-pattern/m-p/765164#M182 Hi,<DIV>Nonsymmetric permutation and scaling matrices are calculated during reordering step of PARDISO, so it's strongly recommended to calculate it for each stiffness matrix independently.</DIV><DIV>With best regards,</DIV><DIV>Alexander Kalinkin</DIV> Tue, 21 Aug 2012 06:56:26 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/PARDISO-on-multiple-systems-with-the-same-sparsity-pattern/m-p/765164#M182 Alexander_K_Intel2 2012-08-21T06:56:26Z