[petsc-users] change matrix

2020-08-31 Thread Kun Jiao via petsc-users
Hi Petsc Experts, Trying to do something like appending some rows (~100 rows) to an already created matrix, but could not find any document about it. Could anyone provide some information about it? Regards, Kun Schlumberger-Private

Re: [petsc-users] [Ext] Re: change matrix

2020-08-31 Thread Kun Jiao via petsc-users
alues() to copy that row into the new matrix. Barry On Aug 31, 2020, at 1:51 PM, Kun Jiao via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: Hi Petsc Experts, Trying to do something like appending some rows (~100 rows) to an already created matrix, but could not find any docume

Re: [petsc-users] [Ext] change matrix

2020-08-31 Thread Kun Jiao via petsc-users
f it is more efficient to reuse the matrix entries than recompute them you can create a larger matrix and then loop over the old matrix calling MatGetRow() and then call MatSetValues() to copy that row into the new matrix. Barry On Aug 31, 2020, at 1:51 PM, Kun Jiao via petsc-users mailto:p

[petsc-users] interchange between two different DMDAs

2019-02-13 Thread Kun Jiao via petsc-users
Hi, I have a problem that need to interchange between two different DMDAs. They have the exactly same size, however they are distributed/partitioned differently across multiple processes. For instance, one DMDA is created with DMDACreated1d, the other is created with DMDACreated3d. How can I in

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

2019-03-25 Thread Kun Jiao via petsc-users
Hi Petsc Experts, Is MatCreateMPIAIJMKL retired in 3.10.4? I got this error with my code which works fine in 3.8.3 version. Regards, Kun Schlumberger-Private

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

2019-03-26 Thread Kun Jiao via petsc-users
undefined in 3.10.4 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 mailto:petsc-users@mcs.anl.gov>> wrote: Hi Petsc Experts, Is MatCreateMPIAIJMKL retired in 3.10.4? I got this error with my code which w

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

2019-03-26 Thread Kun Jiao via petsc-users
ot;MatCreateMPIAIJMKL" is undefined in 3.10.4 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 mailto:petsc-users@mcs.anl.gov>> wrote: Hi Petsc Experts, Is MatCreateMPIAIJMKL retired in 3.10.4? I got

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

2019-03-26 Thread Kun Jiao via petsc-users
ing error, error message is: > > error: identifier "MatCreateMPIAIJMKL" is undefined. > > > > > > From: Mark Adams mailto:mfad...@lbl.gov>> > Sent: Tuesday, March 26, 2019 6:48 AM > To: Kun Jiao mailto:kj...@slb.com>> > Cc: petsc-user

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

2019-03-26 Thread Kun Jiao via petsc-users
t,PetscInt,PetscIn +t,const PetscInt[],PetscInt,const PetscInt[],Mat*); #endif PETSC_EXTERN PetscErrorCode MatCreateSeqSELL(MPI_Comm,PetscInt,PetscInt,PetscInt,const PetscInt[],Mat*); Also note: - this routine is available only when PETSc is built with Intel MKL Satish On Tue, 26 Mar 2019, K

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

2019-03-26 Thread Kun Jiao via petsc-users
: identifier "MatCreateMPIAIJMKL" is undefined in 3.10.4 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,

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

2019-03-26 Thread Kun Jiao via petsc-users
ssing function prototypes. Best regards, Richard On 3/26/19 2:37 PM, Kun Jiao via petsc-users wrote: 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 <mailto

[petsc-users] MatMultTranspose memory usage

2019-07-17 Thread Kun Jiao via petsc-users
Hi, I was using MatMultTranpose and MatMult to solver a linear system. However we found out, MatMultTranpose create an extra memory copy of matrix for its operation. This extra memory copy is not stated everywhere in petsc manual. This basically double my memory requirement to solve my system.

[petsc-users] PetscInitialize take various time to run on different node

2019-08-19 Thread Kun Jiao via petsc-users
Hi, I am running into a problem, different node/process take different time to run petscinitialize. The difference could be 5~10mins. I put some check before petscinitialize, and after petscinitialize. MPI_Init(argc, args); int rank; MPI_Comm_rank(MPI_COMM_WORLD, &rank ); std:

Re: [petsc-users] [Ext] Re: configure error

2020-06-15 Thread Kun Jiao via petsc-users
Thanks. Let me fix that. Try it again. From: Matthew Knepley Sent: Monday, June 15, 2020 12:25 PM To: Kun Jiao Cc: petsc-users@mcs.anl.gov Subject: [Ext] Re: [petsc-users] configure error On Mon, Jun 15, 2020 at 12:13 PM Kun Jiao via petsc-users mailto:petsc-users@mcs.anl.gov>> wro

Re: [petsc-users] [Ext] Re: configure error

2020-06-15 Thread Kun Jiao via petsc-users
/linux/mkl/lib/intel64 Also you are mixing using the Intel MPI and OpenMPI, I am not sure this will work. /opt/openmpi_4.0.3_intel2019.3/lib -lmkl_blacs_intelmpi_lp64 On Jun 15, 2020, at 12:11 PM, Kun Jiao via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: Hi, Lately I run

Re: [petsc-users] [Ext] Re: configure error

2020-06-15 Thread Kun Jiao via petsc-users
Great, it fix everything. Sorry, my bad I did not catch this silly error. Regards, Kun From: Matthew Knepley Sent: Monday, June 15, 2020 12:25 PM To: Kun Jiao Cc: petsc-users@mcs.anl.gov Subject: [Ext] Re: [petsc-users] configure error On Mon, Jun 15, 2020 at 12:13 PM Kun Jiao via petsc

Re: [petsc-users] [Ext] Re: configure error

2020-06-15 Thread Kun Jiao via petsc-users
ith-blaslapack-dir=/mnt/disks/intel-2018-3-222-blade-runtime-env-2018-1-07-08-2018-132838/2019.3/intel/compilers_and_libraries_2019.3.199/linux/mkl Satish On Mon, 15 Jun 2020, Kun Jiao via petsc-users wrote: > Actually, that is a question I should ask. When I use --with-64-bit-indices >