<div dir="ltr">Thanks for the update. Let's hope Intel will add it.<br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">--Junchao Zhang</div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, May 31, 2020 at 9:53 PM Marius Buerkle <<a href="mailto:mbuerkle@web.de">mbuerkle@web.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:Verdana;font-size:12px"><div>Sorry for the late reply, but it took some time to figure things out with intel support. It turned out that the feature is at the present only available in the non-mpi MKL pardiso but not in the MKL cluster solver, so the MKL manual is actually wrong here. I filled a feature request at intel premium support but no clue how long this will take (if ever) to add it. in case it will be available in mkl at some point i will come back to you.</div>
<div> </div>
<div>best and thanks</div>
<div>marius</div>
<div>
<div>
<div name="quote" style="margin:10px 5px 5px 10px;padding:10px 0px 10px 10px;border-left:2px solid rgb(195,217,229)">
<div style="margin:0px 0px 10px"><br>
<b>Betreff:</b> Re: Re: Re: [petsc-users] mkl cpardiso iparm 31</div>
<div name="quoted-content">
<div>Marius,
<div> Thanks for the update. Once you get feedback from Intel, please let us know. If Intel supports it, we can add it. I am new to pardiso, but I think it is doable.
<div>
<div>
<div>--Junchao Zhang</div>
</div>
</div>
</div>
</div>
<div class="gmail_quote">
<div class="gmail_attr">On Thu, May 7, 2020 at 11:56 PM Marius Buerkle <<a href="mailto:mbuerkle@web.de" target="_blank">mbuerkle@web.de</a>> wrote:</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div style="font-family:Verdana;font-size:12px">
<div>Hi Junchao,</div>
<div> </div>
<div>I contacted intel support regarding this, they told me that this is a typo in the manual and that iparm[30] is indeed used. However, while it works for the non-MPI (MKL_PARDISO) version, it does not work, or I could not get it working, for the cluster sparse solver (MKL_CPARDISO). I reported this also to intel support but no reply yet. I was also wondering about how perm(N) is distributed and I don't know at the moment.</div>
<div> </div>
<div>Best,</div>
<div>Marius</div>
<div>
<div>
<div style="margin:10px 5px 5px 10px;padding:10px 0px 10px 10px;border-left:2px solid rgb(195,217,229)">
<div style="margin:0px 0px 10px"> </div>
<div>
<div>Marius,
<div> You are right. perm is not referenced. I searched and found this, </div>
<div><a href="https://software.intel.com/content/www/us/en/develop/documentation/mkl-developer-reference-c/top/sparse-solver-routines/parallel-direct-sparse-solver-for-clusters-interface/cluster-sparse-solver.html" target="_blank">https://software.intel.com/content/www/us/en/develop/documentation/mkl-developer-reference-c/top/sparse-solver-routines/parallel-direct-sparse-solver-for-clusters-interface/cluster-sparse-solver.html</a>.</div>
<div> It says "perm Ignored". But from other parts of the document, it seems perm is used. I'm puzzled whether Intel MKL pardiso supports this feature or not.</div>
<div> </div>
I am thinking about adding MatMkl_CPardisoSetPerm(Mat A, IS perm) or MatMkl_CPardisoSetPerm(Mat A, const PetscInt *perm). But I don't know whether perm(N) is distributed or every mpi rank has the same perm(N).<br>
Do you know good Intel MKL pardiso documentation or examples for me to reference?<br>
<br>
Thank you.
<div>
<div>
<div>
<div>
<div>
<div>--Junchao Zhang</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div class="gmail_quote">
<div class="gmail_attr">On Thu, May 7, 2020 at 3:14 AM Marius Buerkle <<a href="mailto:mbuerkle@web.de" target="_blank">mbuerkle@web.de</a>> wrote:</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div style="font-family:Verdana;font-size:12px">
<div>Hi,</div>
<div> </div>
<div>Thanks for the info. But how do I set the values to be calculated. According to the intel parallel sparse cluster solver manual the entries have to be defined in the permutation vector (before each call). However, if I understand what is happening in mkl_cpardiso.c. correctly, perm is set to 0 during the initialization phase and then not referenced anymore. Is this correct? How can I specify the necessary entries in perm?</div>
<div> </div>
<div>Best,</div>
<div>Marius</div>
<div>
<div>
<div style="margin:10px 5px 5px 10px;padding:10px 0px 10px 10px;border-left:2px solid rgb(195,217,229)">
<div style="margin:0px 0px 10px"> </div>
<div>
<div>
<div> </div>
<div class="gmail_quote">
<div class="gmail_attr">On Fri, May 1, 2020 at 3:33 AM Marius Buerkle <<a href="mailto:mbuerkle@web.de" target="_blank">mbuerkle@web.de</a>> wrote:</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div style="font-family:Verdana;font-size:12px">
<div>Hi,</div>
<div> </div>
<div>Is the option "<span style="color:rgb(102,102,102)">-mat_mkl_cpardiso_31" to calculate </span>Partial solve and computing selected components of the solution vectors <span style="color:rgb(102,102,102)">actually supported by PETSC?</span></div>
</div>
</div>
</blockquote>
<div>From the code, it seems so. </div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div style="font-family:Verdana;font-size:12px">
<div> </div>
<div><font color="#666666"><span>Marius</span></font></div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div></div></div>
</blockquote></div>