topic Pardiso parameters : phase, pt maxfct and mnum. Matrix inverse in IntelĀ® oneAPI Math Kernel Library
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849501#M6457
Hi,<BR /><BR />I have a few questions (sorry if some of them look basic). First, is there a way to retrieve the result of the factorisation of phase 2 (the same result one would get from Lapack's posv)? I think there is something to do with parameters pt, maxfct and mnum, but I don't understand exactly how to use them. I suppose I would have to call Pardiso two times : one with phase = 12, and the second with phase = 33 (I want to mimick posv). Is this correct?<BR /><BR />Second, is there a way to use Pardiso to find the inverse of a matrix (like finding A^-1 from Ax = b, if x and b are known)? I was hoping to use it instead of Lapack's potri.<BR /><BR />Thanks and have a nice day!Fri, 09 Apr 2010 15:28:40 GMTxdominique2010-04-09T15:28:40ZPardiso parameters : phase, pt maxfct and mnum. Matrix inverse
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849501#M6457
Hi,<BR /><BR />I have a few questions (sorry if some of them look basic). First, is there a way to retrieve the result of the factorisation of phase 2 (the same result one would get from Lapack's posv)? I think there is something to do with parameters pt, maxfct and mnum, but I don't understand exactly how to use them. I suppose I would have to call Pardiso two times : one with phase = 12, and the second with phase = 33 (I want to mimick posv). Is this correct?<BR /><BR />Second, is there a way to use Pardiso to find the inverse of a matrix (like finding A^-1 from Ax = b, if x and b are known)? I was hoping to use it instead of Lapack's potri.<BR /><BR />Thanks and have a nice day!Fri, 09 Apr 2010 15:28:40 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849501#M6457xdominique2010-04-09T15:28:40ZPardiso parameters : phase, pt maxfct and mnum. Matrix inverse
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849502#M6458
There's really no way currently to access the intermediate factors in PARDISO. This request has been made and I can add this thread to the record for tracking purposes.<BR /><BR />There are no special settings for this if that is what you mean. You could set b to the idenity, I. I'll let others comment on whether this is a good idea or not.<BR /><BR />-ToddSat, 10 Apr 2010 00:17:33 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849502#M6458Todd_R_Intel2010-04-10T00:17:33ZPardiso parameters : phase, pt maxfct and mnum. Matrix inverse
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849503#M6459
<P>Hi,</P><P>Are there plans to release Lapack routines and drivers using sparse storage format? That would be really helpful as it will decrease significantly memory storage size. I'm thinking about computing matrix inverses, as they are computer intensive.</P><P>Thanks and have a nice day!</P>Wed, 14 Apr 2010 15:25:07 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849503#M6459xdominique2010-04-14T15:25:07ZPardiso parameters : phase, pt maxfct and mnum. Matrix inverse
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849504#M6460
then it would be something like SparseLapack :). There no such plan in the nearest versions. What we have for today - Sparse BLAS Level 1 and 2 Routines and Direct and iterative Sparse Solvers for sparse storage format.<DIV>--Gennady</DIV>Wed, 14 Apr 2010 15:59:42 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849504#M6460Gennady_F_Intel2010-04-14T15:59:42ZPardiso parameters : phase, pt maxfct and mnum. Matrix inverse
https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849505#M6461
Sorry for digging through an old post, but I have some news about matrix inverses (for the posterity...). I tried using Pardiso to calculate the inverse, settingthe right-hand-sideas the identity matrix. It works perfectly! The only drawback is that it takes thrice the memory space as you have to allocate space for the right-hand-side matrix and the solution matrix the same size as the input matrix. So for large problems, it is not very useful. But it works. ;)<BR /><BR />And yes, Mr. Gennady, if I may made a request about it, I would really appreciate a Sparse Lapack. (:Mon, 17 May 2010 15:35:01 GMThttps://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/Pardiso-parameters-phase-pt-maxfct-and-mnum-Matrix-inverse/m-p/849505#M6461xdominique2010-05-17T15:35:01Z