Re: [petsc-users] petsc4py help with parallel execution

2018-11-16 Thread Balay, Satish via petsc-users
Likely this build is with gnu compilers. https://software.intel.com/en-us/articles/intel-mkl-link-line-advisor/ The above suggests - the following might work: '--with-blaslapack-lib=-L/opt/intel/mkl/lib/intel64 -lmkl_intel_lp64 -lmkl_gnu_thread -lmkl_core -lgomp -lpthread -lm -ldl' Satish On

Re: [petsc-users] petsc4py help with parallel execution

2018-11-16 Thread Balay, Satish via petsc-users
On Fri, 16 Nov 2018, Ivan Voznyuk via petsc-users wrote: > Hi Satish, > Thanks for your reply. > > Bad news... I tested 2 solutions that you proposed, none has worked. > > 1. --with-blaslapack-dir=/opt/intel/mkl > --with-mkl_pardiso-dir=/opt/intel/mkl installed well, without any problems. > Howe

Re: [petsc-users] petsc4py help with parallel execution

2018-11-16 Thread Balay, Satish via petsc-users
Yes PETSc prefers sequential MKL - as MPI handles parallelism. One way to trick petsc configure to use threaded MKL is to enable pardiso. i.e: --with-blaslapack-dir=/opt/intel/mkl --with-mkl_pardiso-dir=/opt/intel/mkl http://ftp.mcs.anl.gov/pub/petsc/nightlylogs/archive/2018/11/15/configure_mast

Re: [petsc-users] petsc4py help with parallel execution

2018-11-15 Thread Ivan via petsc-users
Matthew, */As I wrote before, its not impossible. You could be directly calling PMI, but I do not think you are doing that./* Could you precise what is PMI? and how can we directly use it? It might be a key to this mystery! */Why do you think its running on 8 processes?/* Well, we base our

Re: [petsc-users] petsc4py help with parallel execution

2018-11-15 Thread Matthew Knepley via petsc-users
On Thu, Nov 15, 2018 at 11:59 AM Ivan Voznyuk wrote: > Hi Matthew, > > Does it mean that by using just command python3 simple_code.py (without > mpiexec) you *cannot* obtain a parallel execution? > As I wrote before, its not impossible. You could be directly calling PMI, but I do not think you a

Re: [petsc-users] petsc4py help with parallel execution

2018-11-15 Thread Ivan Voznyuk via petsc-users
Hi Matthew, Does it mean that by using just command python3 simple_code.py (without mpiexec) you *cannot* obtain a parallel execution? It s been 5 days we are trying to understand with my colleague how he managed to do so. It means that by using simply python3 simple_code.py he gets 8 processors w

[petsc-users] petsc4py help with parallel execution

2018-11-15 Thread Ivan Voznyuk via petsc-users
Dear PETSC community, I have a question regarding the parallel execution of petsc4py. I have a simple code (here attached simple_code.py) which solves a system of linear equations Ax=b using petsc4py. To execute it, I use the command python3 simple_code.py which yields a sequential performance. W