[petsc-users] mkl cpardiso iparm 31

[Junchao Zhang]

Marius,
  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.
--Junchao Zhang


On Thu, May 7, 2020 at 11:56 PM Marius Buerkle <mbuerkle at web.de> wrote:

> Hi Junchao,
>
> 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.
>
> Best,
> Marius
>
>
>
> Marius,
>   You are right. perm is not referenced. I searched and found this,
>
> 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
> .
>   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.
>
>  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).
>  Do you know good Intel MKL pardiso documentation or examples for me to
> reference?
>
> Thank you.
> --Junchao Zhang
>
> On Thu, May 7, 2020 at 3:14 AM Marius Buerkle <mbuerkle at web.de> wrote:
>
>> Hi,
>>
>> 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?
>>
>> Best,
>> Marius
>>
>>
>>
>>
>> On Fri, May 1, 2020 at 3:33 AM Marius Buerkle <mbuerkle at web.de> wrote:
>>
>>> Hi,
>>>
>>> Is the option "-mat_mkl_cpardiso_31" to calculate Partial solve and
>>> computing selected components of the solution vectors actually
>>> supported by PETSC?
>>>
>> From the code, it seems so.
>>
>>>
>>> Marius
>>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200508/cc9a7693/attachment.html>