<div dir="ltr">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. <br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">--Junchao Zhang</div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, May 7, 2020 at 11:56 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>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 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"> </div>

<div name="quoted-content">
<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>