[petsc-users] MKL_PARDISO question

[Xiangdong]

Thank you very much for your information. I pulled the master branch but
got the error when configuring it.

When I run configure without mkl_cpardiso
(configure.log_nocpardiso):  ./configure PETSC_ARCH=arch-debug
 --with-debugging=1 --with-mpi-dir=$MPI_ROOT
--with-blaslapack-dir=${MKL_ROOT} , it works fine.

However, when I add mkl_cpardiso (configure.log_withcpardiso): ./configure
PETSC_ARCH=arch-debug  --with-debugging=1 --with-mpi-dir=$MPI_ROOT
 -with-blaslapack-dir=${MKL_ROOT} --with-mkl_cpardiso-dir=${MKL_ROOT} , it
complains about "Could not find a functional BLAS.", but the blas was
provided through mkl as same as previous configuration.

Can you help me on the configuration? Thank you.

Xiangdong

On Wed, Sep 18, 2019 at 2:39 PM Smith, Barry F. <bsmith at mcs.anl.gov> wrote:

>
>
> > On Sep 18, 2019, at 9:15 AM, Xiangdong via petsc-users <
> petsc-users at mcs.anl.gov> wrote:
> >
> > Hello everyone,
> >
> > From here,
> >
> https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MATSOLVERMKL_PARDISO.html
> >
> > It seems thatMKL_PARDISO only works for seqaij. I am curious that
> whether one can use mkl_pardiso in petsc with multi-thread.
>
>    You can use  mkl_pardiso for multi-threaded and  mkl_cpardiso for MPI
> parallelism.
>
>    In both cases you must use the master branch of PETSc (or the next
> release of PETSc) to do this this easily.
>
>    Note that when you use mkl_pardiso with multiple threads the matrix is
> coming from a single MPI process (or the single program if not running with
> MPI). So it is not MPI parallel that matches the rest of the parallelism
> with PETSc. So one much be a little careful: for example if one has 4 cores
> and uses them all with mpiexec -n 4 and then uses mkl_pardiso with 4
> threads (each) then you have 16 threads fighting over 4 cores. So you need
> to select the number of MPI processes and number of threads wisely.
>
> >
> > Is there any reason that MKL_PARDISO is not listed in the linear solver
> table?
> > https://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html
> >
>
>    Just an oversight, thanks for letting us know, I have added it.
>
>
> > Thank you.
> >
> > Best,
> > Xiangdong
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190918/fb10afda/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: configure.log_nocpardiso
Type: application/octet-stream
Size: 1021259 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190918/fb10afda/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: configure.log_withcpardiso
Type: application/octet-stream
Size: 1790213 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190918/fb10afda/attachment-0003.obj>