Re: [petsc-users] Solving block systems with some null diagonal blocks

2019-03-26 Thread Manuel Colera Rico via petsc-users
OK, thank you Matt. Manuel --- On 3/25/19 6:27 PM, Matthew Knepley wrote: On Mon, Mar 25, 2019 at 8:07 AM Manuel Colera Rico via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: Hello, I would like to solve a N*N block system (with N>2) in which some of the diagonal

Re: [petsc-users] Bad memory scaling with PETSc 3.10

2019-03-26 Thread Myriam Peyrounette via petsc-users
How can I be sure they are indeed used? Can I print this information in some log file? Thanks in advance Myriam Le 03/25/19 à 18:24, Matthew Knepley a écrit : > On Mon, Mar 25, 2019 at 10:54 AM Myriam Peyrounette via petsc-users > mailto:petsc-users@mcs.anl.gov>> wrote: > > Hi, > >

Re: [petsc-users] Bad memory scaling with PETSc 3.10

2019-03-26 Thread Myriam Peyrounette via petsc-users
Oh you were right, the three options are unsused (-matptap_via scalable, -inner_offdiag_matmatmult_via scalable and -inner_diag_matmatmult_via scalable). Does this mean I am not using the associated PtAP functions? Myriam Le 03/26/19 à 11:10, Dave May a écrit : > > On Tue, 26 Mar 2019 at 09:52,

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Kun Jiao via petsc-users
It is compiling error, error message is: error: identifier "MatCreateMPIAIJMKL" is undefined. From: Mark Adams Sent: Tuesday, March 26, 2019 6:48 AM To: Kun Jiao Cc: petsc-users@mcs.anl.gov Subject: [Ext] Re: [petsc-users] error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

Re: [petsc-users] error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Mark Adams via petsc-users
Please send the output of the error (runtime, compile time, link time?) On Mon, Mar 25, 2019 at 10:50 PM Kun Jiao via petsc-users < petsc-users@mcs.anl.gov> wrote: > Hi Petsc Experts, > > > > Is MatCreateMPIAIJMKL retired in 3.10.4? > > > > I got this error with my code which works fine in 3.8.3

Re: [petsc-users] Bad memory scaling with PETSc 3.10

2019-03-26 Thread Myriam Peyrounette via petsc-users
I checked with -ksp_view (attached) but no prefix is associated with the matrix. Some are associated to the KSP and PC, but none to the Mat. Le 03/26/19 à 11:55, Dave May a écrit : > > > On Tue, 26 Mar 2019 at 10:36, Myriam Peyrounette > mailto:myriam.peyroune...@idris.fr>> wrote: > > Oh you

Re: [petsc-users] Bad memory scaling with PETSc 3.10

2019-03-26 Thread Matthew Knepley via petsc-users
On Tue, Mar 26, 2019 at 9:27 AM Myriam Peyrounette < myriam.peyroune...@idris.fr> wrote: > I checked with -ksp_view (attached) but no prefix is associated with the > matrix. Some are associated to the KSP and PC, but none to the Mat > Another thing that could prevent options being used is that

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Mark Adams via petsc-users
I assume the whole error message will have the line of code. Please send the whole error message and line of offending code if not included. On Tue, Mar 26, 2019 at 10:08 AM Kun Jiao wrote: > It is compiling error, error message is: > > > > error: identifier "MatCreateMPIAIJMKL" is undefined. >

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Kun Jiao via petsc-users
Strange things, when I compile my code in the test dir in PETSC, it works. After I "make install" PETSC, and try to compile my code against the installed PETSC, it doesn't work any more. I guess this is what you means. Is there any way to reenable MatCreateMPIAIJMKL public interface? And, I

Re: [petsc-users] Confusing Schur preconditioner behaviour

2019-03-26 Thread Cotter, Colin J via petsc-users
Hi Dave, Thanks for the tip - you were right, and this works better for higher resolutions now. all the best --Colin From: Dave May Sent: 19 March 2019 11:25:11 To: Cotter, Colin J Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] Confusing Schur

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Mark Adams via petsc-users
So this works with v3.8? I don't see any differences (I see Satish figured this out and has suggestions). You could also work around it with code like this: ierr = MatCreate(PETSC_COMM_WORLD,);CHKERRQ(ierr); ierr = MatSetType(A,MATAIJMKL);CHKERRQ(ierr); ierr =

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Kun Jiao via petsc-users
[kjiao@hyi0016 src/lsqr]% make [ 50%] Building CXX object lsqr/CMakeFiles/p_lsqr.dir/lsqr.cc.o /wgdisk/hy3300/source_code_dev/imaging/kjiao/src/git/src/lsqr/lsqr.cc(318): error: identifier "MatCreateMPIAIJMKL" is undefined ierr =

Re: [petsc-users] Bad memory scaling with PETSc 3.10

2019-03-26 Thread Myriam Peyrounette via petsc-users
*SetFromOptions() was not called indeed... Thanks! The code performance is better now with regard to memory usage! I still have to plot the memory scaling on bigger cases to see if it has the same good behaviour as when using the 3.6 version. I'll let ou know as soon as I have plotted it.

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Balay, Satish via petsc-users
>>> balay@sb /home/balay/petsc (maint=) $ git grep MatCreateMPIAIJMKL maint-3.8 maint-3.8:src/mat/impls/aij/mpi/aijmkl/mpiaijmkl.c: MatCreateMPIAIJMKL - Creates a sparse parallel matrix whose local maint-3.8:src/mat/impls/aij/mpi/aijmkl/mpiaijmkl.c:PetscErrorCode

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Mark Adams via petsc-users
On Tue, Mar 26, 2019 at 3:00 PM Kun Jiao wrote: > Strange things, when I compile my code in the test dir in PETSC, it works. > After I "make install" PETSC, and try to compile my code against the > installed PETSC, it doesn't work any more. > I'm not sure I follow what you are doing exactly but

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Balay, Satish via petsc-users
Please apply the patch I sent earlier and retry. Satish On Tue, 26 Mar 2019, Kun Jiao via petsc-users wrote: > Strange things, when I compile my code in the test dir in PETSC, it works. > After I "make install" PETSC, and try to compile my code against the > installed PETSC, it doesn't work

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Kun Jiao via petsc-users
One strange thing I just found out. Compile a *.c file make it work. mpiicc -o ex5.o -c -fPIC -wd1572 -Ofast -xHost -I /NFS/home/home3/kjiao/software/petsc_3.10.4/include/ -I/wgdisk/hy3300/source_code_dev/imaging/kjiao/software/intel/compilers_and_libraries_2019.2.187/linux/mkl/include

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Kun Jiao via petsc-users
And yes, by applying the patch in the petscmat.h, everything works. Thanks for the help. Regards, Kun Schlumberger-Private -Original Message- From: Balay, Satish Sent: Tuesday, March 26, 2019 3:42 PM To: Kun Jiao Cc: petsc-users Subject: Re: [petsc-users] [Ext] Re: error:

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Mills, Richard Tran via petsc-users
Hi Kun, I'm the author of most of the AIJMKL stuff in PETSc. My apologies for having inadvertently omitted these function prototypes for these interfaces; I'm glad that Satish's patch has fixed this. I want to point out that -- though I can envision some scenarios in which one would want to

Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4

2019-03-26 Thread Kun Jiao via petsc-users
Hi Richard, Understood! Thanks very much for you advice. Regards, Kun Schlumberger-Private From: Mills, Richard Tran Sent: Tuesday, March 26, 2019 8:11 PM To: petsc-users@mcs.anl.gov Cc: Kun Jiao Subject: Re: [petsc-users] [Ext] Re: error: identifier "MatCreateMPIAIJMKL" is undefined in