Re: [petsc-users] MUMPS Error 'INFOG(1)=-3 INFO(2)=3' (SPARSE MATRIX INVERSE)

2023-07-27 Thread Pierre Jolivet
MUMPS errors are documented in section 8 of https://mumps-solver.org/doc/userguide_5.6.1.pdf Thanks, Pierre > On 27 Jul 2023, at 3:50 PM, maitri ksh wrote: > > I am using 'MatMumpsGetInverse()' to get the inverse of a sparse matrix. I am > using parts of ex214.c >

[petsc-users] MUMPS Error 'INFOG(1)=-3 INFO(2)=3' (SPARSE MATRIX INVERSE)

2023-07-27 Thread maitri ksh
I am using 'MatMumpsGetInverse()' to get the inverse of a sparse matrix. I am using parts of ex214.c code to get the inverse, but I get an error that seems to be coming from

Re: [petsc-users] MUMPS error

2021-05-06 Thread Bhamidipati, Vikram
The problem did not go away on a second attempt but was overcome with '--with-make-np=1' flag. Thank you, Vikram -Original Message- From: Satish Balay Sent: Thursday, May 6, 2021 5:03 PM To: Bhamidipati, Vikram Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] MUMPS error

Re: [petsc-users] MUMPS error

2021-05-06 Thread Satish Balay via petsc-users
Can you retry and see if the problem persists? There might be an issue with missing dependencies in mumps parallel build.. If it still fails - you can try the additional configure option: --with-make-np=1 Satish On Thu, 6 May 2021, Bhamidipati, Vikram wrote: > I configured a PetSc build to

Re: [petsc-users] MUMPS Error

2019-01-08 Thread Zhang, Hong via petsc-users
Sal: How large is the matrix, its dimension and number of non-zeros? Can you run it sequentially or with less number of processors? The matrix is 139603x139603 and with 15975113 of non-zero elements. This is quite large for direct solver, but it is very sparse. Have you tried different matrix

Re: [petsc-users] MUMPS Error

2019-01-07 Thread Zhang, Hong via petsc-users
Sal: Recent versions of MUMPS have bugs when using SCALAPACK. You may also try to disable parallelism on the root node by using tried that, does not seem to help... but it is odd. this is not a big system, it is small and made for testing. How large is the matrix, its dimension and number of

Re: [petsc-users] MUMPS Error

2019-01-03 Thread Zhang, Hong via petsc-users
You may try different matrix orderings, or try superlu_dist. Hong On Thu, Jan 3, 2019 at 7:39 AM Matthew Knepley via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: On Thu, Jan 3, 2019 at 5:58 AM Sal Am via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: [0]PETSC ERROR:

Re: [petsc-users] MUMPS Error

2018-12-12 Thread Zhang, Junchao via petsc-users
On Wed, Dec 12, 2018 at 7:14 AM Sal Am via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: Hi I am getting an error using MUMPS. How I run it: bash-4.2$ mpiexec -n 16 ./solveCSys -ksp_type richardson -pc_type lu -pc_factor_mat_solver_type mumps -ksp_max_it 1 -ksp_monitor_true_residual

Re: [petsc-users] MUMPS error reporting in PETSc-3.7

2016-06-06 Thread David Knezevic
On Mon, Jun 6, 2016 at 4:16 PM, Barry Smith wrote: > > David, > >Actually all the routines for accessing the MUMPS information already > exist. I have added this to the manual page for MATSOLVERMUMPS > > Notes: When a MUMPS factorization fails inside a KSP solve,

Re: [petsc-users] MUMPS error reporting in PETSc-3.7

2016-06-06 Thread Barry Smith
David, Actually all the routines for accessing the MUMPS information already exist. I have added this to the manual page for MATSOLVERMUMPS Notes: When a MUMPS factorization fails inside a KSP solve, for example with a KSP_DIVERGED_PCSETUP_FAILED, one can find the MUMPS information

Re: [petsc-users] MUMPS error reporting in PETSc-3.7

2016-06-05 Thread David Knezevic
On Sun, Jun 5, 2016 at 1:42 PM, Barry Smith wrote: > > > On Jun 5, 2016, at 12:12 PM, David Knezevic > wrote: > > > > In PETSc-3.6, when I hit an error using the MUMPS solver, I would get an > error like this: > > > > [0]PETSC ERROR:

Re: [petsc-users] MUMPS error reporting in PETSc-3.7

2016-06-05 Thread Barry Smith
> On Jun 5, 2016, at 12:12 PM, David Knezevic > wrote: > > In PETSc-3.6, when I hit an error using the MUMPS solver, I would get an > error like this: > > [0]PETSC ERROR: - Error Message >

[petsc-users] MUMPS error reporting in PETSc-3.7

2016-06-05 Thread David Knezevic
In PETSc-3.6, when I hit an error using the MUMPS solver, I would get an error like this: [0]PETSC ERROR: - Error Message -- [0]PETSC ERROR: Error in external library [0]PETSC ERROR: Error reported by MUMPS in

Re: [petsc-users] MUMPS error and superLU error

2015-06-22 Thread venkatesh g
Hi I have restructured my matrix eigenvalue problem to see why B is singular as you suggested by changing the governing equations in different form. Now my matrix B is not singular. Both A and B are invertible in Ax=lambda Bx. Still I receive error in MUMPS as it uses large memory (attached is

Re: [petsc-users] MUMPS error and superLU error

2015-06-22 Thread Hong
Venkatesh, You may also test superlu_dist, which may use less memory. Hong On Mon, Jun 22, 2015 at 12:43 PM, Barry Smith bsm...@mcs.anl.gov wrote: There is nothing we can really do to help on the PETSc side. I do note from the output REDISTRIB: TOTAL DATA LOCAL/SENT = 328575589

Re: [petsc-users] MUMPS error and superLU error

2015-06-22 Thread Barry Smith
There is nothing we can really do to help on the PETSc side. I do note from the output REDISTRIB: TOTAL DATA LOCAL/SENT = 328575589 1437471711 GLOBAL TIME FOR MATRIX DISTRIBUTION =206.6792 ** Memory relaxation parameter ( ICNTL(14) ):35 ** Rank of

Re: [petsc-users] MUMPS error and superLU error

2015-05-31 Thread venkatesh g
Hi all, I tried to run my Generalized Eigenproblem in 120 x 24 = 2880 cores. The matrix size of A = 20GB and B = 5GB. It got killed after 7 Hrs of run time. Please see the mumps error log. Why must it fail ? I gave the command: aprun -n 240 -N 24 ./ex7 -f1 a110t -f2 b110t -st_type sinvert

[petsc-users] MUMPS error and superLU error

2015-05-31 Thread Dave May
It failed due to a lack of memory. OOM stands for out of memory. OOM killer terminated your job means you ran out of memory. On Sunday, 31 May 2015, venkatesh g venkateshg...@gmail.com javascript:_e(%7B%7D,'cvml','venkateshg...@gmail.com'); wrote: Hi all, I tried to run my Generalized

Re: [petsc-users] MUMPS error and superLU error

2015-05-31 Thread venkatesh g
Sorry the command was: aprun -n 2880 -N 24 ./ex7 -f1 a110t -f2 b110t -st_type sinvert -eps_nev 1 -st_ksp_type preonly -st_pc_type lu -st_pc_factor_mat_solver_package mumps -mat_mumps_cntl_1 1e-6 -mat_mumps_icntl_4 2 cheers, Venkatesh On Sun, May 31, 2015 at 6:07 PM, venkatesh g

Re: [petsc-users] MUMPS error and superLU error

2015-05-31 Thread Hong
venkatesh, As we discussed previously, even on smaller problems, both mumps and superlu_dist failed, although Mumps gave OOM error in numerical factorization. You acknowledged that B is singular, which may need additional reformulation for your eigenvalue problems. The option '-st_type sinvert'

Re: [petsc-users] MUMPS error and superLU error

2015-05-29 Thread venkatesh g
On Tue, May 26, 2015 at 9:02 PM, Hong hzh...@mcs.anl.gov wrote: 'A serial job in MATLAB for the same matrices takes 60GB. ' Can you run this case in serial? If so, try petsc, superlu or mumps to make sure the matrix is non-singular. B matrix is singular but I get my result in Petsc and Mumps

Re: [petsc-users] MUMPS error and superLU error

2015-05-29 Thread Hong
venkatesh: On Tue, May 26, 2015 at 9:02 PM, Hong hzh...@mcs.anl.gov wrote: 'A serial job in MATLAB for the same matrices takes 60GB. ' Can you run this case in serial? If so, try petsc, superlu or mumps to make sure the matrix is non-singular. B matrix is singular but I get my result in

Re: [petsc-users] MUMPS error and superLU error

2015-05-26 Thread Hong
'A serial job in MATLAB for the same matrices takes 60GB. ' Can you run this case in serial? If so, try petsc, superlu or mumps to make sure the matrix is non-singular. Both mumps and superlu_dist show crash in MatFactorNumeric(). Mumps gives error [16]PETSC ERROR: Error reported by MUMPS in

Re: [petsc-users] MUMPS error

2015-05-25 Thread Matthew Knepley
On Mon, May 25, 2015 at 2:54 AM, venkatesh g venkateshg...@gmail.com wrote: Ok this will load the matrices in parallel correct ? Yes Matt On Sun, May 24, 2015 at 7:36 PM, Matthew Knepley knep...@gmail.com wrote: On Sun, May 24, 2015 at 8:57 AM, venkatesh g venkateshg...@gmail.com

Re: [petsc-users] MUMPS error

2015-05-25 Thread venkatesh g
Ok this will load the matrices in parallel correct ? On Sun, May 24, 2015 at 7:36 PM, Matthew Knepley knep...@gmail.com wrote: On Sun, May 24, 2015 at 8:57 AM, venkatesh g venkateshg...@gmail.com wrote: I am using Matload option as in the ex7.c code given by the Slepc. ierr =

Re: [petsc-users] MUMPS error

2015-05-25 Thread venkatesh g
Ok now I have restructured the matrix. Now I will not use mumps and use only Petsc and Slepc. I am getting some error. I will open another thread. Venkatesh On Mon, May 25, 2015 at 2:36 PM, Matthew Knepley knep...@gmail.com wrote: On Mon, May 25, 2015 at 2:54 AM, venkatesh g

Re: [petsc-users] MUMPS error

2015-05-24 Thread venkatesh g
I am using Matload option as in the ex7.c code given by the Slepc. ierr = MatLoad(A,viewer);CHKERRQ(ierr); There is no problem here right ? or any additional option is required for very large matrices while running the eigensolver in parallel ? cheers, Venkatesh On Sat, May 23, 2015 at 5:43

Re: [petsc-users] MUMPS error

2015-05-24 Thread Matthew Knepley
On Sun, May 24, 2015 at 8:57 AM, venkatesh g venkateshg...@gmail.com wrote: I am using Matload option as in the ex7.c code given by the Slepc. ierr = MatLoad(A,viewer);CHKERRQ(ierr); There is no problem here right ? or any additional option is required for very large matrices while running

Re: [petsc-users] MUMPS error

2015-05-23 Thread venkatesh g
Hi, The same eigenproblem runs with 120 GB RAM in a serial machine in Matlab. In Cray I fired with 240*4 GB RAM in parallel. So it has to go in right ? And for small matrices it is having negative scaling i.e 24 core is running faster. I have attached the submission script. Pls see.. Kindly

Re: [petsc-users] MUMPS error

2015-05-23 Thread Matthew Knepley
On Sat, May 23, 2015 at 6:44 AM, venkatesh g venkateshg...@gmail.com wrote: Hi, The same eigenproblem runs with 120 GB RAM in a serial machine in Matlab. In Cray I fired with 240*4 GB RAM in parallel. So it has to go in right ? I do not know how MUMPS allocates memory, but the message is

Re: [petsc-users] MUMPS error

2015-05-23 Thread venkatesh g
Hi, Thanks. Per node it has 24 cores and each core has 4 GB RAM. And the job was submitted in 10 nodes. So, does it mean it requires 10G for one core ? or for 1 node ? cheers, Venkatesh On Sat, May 23, 2015 at 5:17 PM, Matthew Knepley knep...@gmail.com wrote: On Sat, May 23, 2015 at 6:44 AM,

Re: [petsc-users] MUMPS error

2015-05-23 Thread Matthew Knepley
On Sat, May 23, 2015 at 2:39 AM, venkatesh g venkateshg...@gmail.com wrote: Hi again, I have installed the Petsc and Slepc in Cray with intel compilers with Mumps. I am getting this error when I solve eigenvalue problem with large matrices: [201]PETSC ERROR: Error reported by MUMPS in

Re: [petsc-users] MUMPS error

2015-05-19 Thread venkatesh g
Hi, I have attached the log of the command which I gave in the master node: make streams NPMAX=32 I dont know why it says 'It appears you have only 1 node'. But other codes run in parallel with good scaling on 8 nodes. Kindly let me know. Venkatesh On Mon, May 18, 2015 at 11:21 PM, Barry

Re: [petsc-users] MUMPS error

2015-05-19 Thread Matthew Knepley
On Tue, May 19, 2015 at 1:04 AM, venkatesh g venkateshg...@gmail.com wrote: Hi, I have attached the log of the command which I gave in the master node: make streams NPMAX=32 I dont know why it says 'It appears you have only 1 node'. But other codes run in parallel with good scaling on 8

Re: [petsc-users] MUMPS error

2015-05-18 Thread Matthew Knepley
On Mon, May 18, 2015 at 8:29 AM, venkatesh g venkateshg...@gmail.com wrote: Hi I have attached the performance logs for 2 jobs on different processors. I had to increase the workspace icntl(14) when I submit on more cores since it is failing with small value of icntl(14). 1.

Re: [petsc-users] MUMPS error

2015-05-18 Thread venkatesh g
Hi, I have emailed the mumps-user list. Actually the cluster has 8 nodes with 16 cores, and other codes scale well. I wanted to ask if this job takes much time, then if I submit on more cores, I have to increase the icntl(14).. which would again take long time. So is there another way ? cheers,

Re: [petsc-users] MUMPS error

2015-05-18 Thread Barry Smith
Run the streams benchmark on this system and send the results. http://www.mcs.anl.gov/petsc/documentation/faq.html#computers On May 18, 2015, at 11:14 AM, venkatesh g venkateshg...@gmail.com wrote: Hi, I have emailed the mumps-user list. Actually the cluster has 8 nodes with 16

Re: [petsc-users] MUMPS error

2015-05-18 Thread Hong
venkatesh: The differences are in MatSolve() as Matt said. The questions are 1. why np=2 calls 7624 MatSolve(), while np=8 calls 9600 times? With slepc/sinvert/lu, would both runs use same algorithm? 2. with np=2, each MatSolve takes 0.027 sec, while np=8 takes 0.065 sec. Would MatSolve be

Re: [petsc-users] MUMPS error

2015-05-17 Thread venkatesh g
Hi, Thanks for the information. I now increased the workspace by adding '-mat_mumps_icntl_14 100' It works. However, the problem is, if I submit in 1 core I get the answer in 200 secs, but with 4 cores and '-mat_mumps_icntl_14 100' it takes 3500secs. My command line is: 'mpiexec -np 4 ./ex7 -f1

Re: [petsc-users] MUMPS error

2015-05-17 Thread Matthew Knepley
On Sun, May 17, 2015 at 1:38 AM, venkatesh g venkateshg...@gmail.com wrote: Hi, Thanks for the information. I now increased the workspace by adding '-mat_mumps_icntl_14 100' It works. However, the problem is, if I submit in 1 core I get the answer in 200 secs, but with 4 cores and

Re: [petsc-users] MUMPS error

2015-05-16 Thread Matthew Knepley
On Sat, May 16, 2015 at 7:08 AM, venkatesh g venkateshg...@gmail.com wrote: Hi, I am trying to solving AX=lambda BX eigenvalue problem. A and B are of sizes 3600x3600 I run with this command : 'mpiexec -np 4 ./ex7 -f1 a2 -f2 b2 -eps_nev 1 -st_type sinvert -eps_max_it 5000 -st_ksp_type

Re: [petsc-users] MUMPS error

2015-05-16 Thread David Knezevic
On Sat, May 16, 2015 at 8:08 AM, venkatesh g venkateshg...@gmail.com wrote: Hi, I am trying to solving AX=lambda BX eigenvalue problem. A and B are of sizes 3600x3600 I run with this command : 'mpiexec -np 4 ./ex7 -f1 a2 -f2 b2 -eps_nev 1 -st_type sinvert -eps_max_it 5000 -st_ksp_type

[petsc-users] MUMPS error

2015-05-16 Thread venkatesh g
Hi, I am trying to solving AX=lambda BX eigenvalue problem. A and B are of sizes 3600x3600 I run with this command : 'mpiexec -np 4 ./ex7 -f1 a2 -f2 b2 -eps_nev 1 -st_type sinvert -eps_max_it 5000 -st_ksp_type preonly -st_pc_type lu -st_pc_factor_mat_solver_package mumps' I get this error: (I

[petsc-users] MUMPS error

2015-04-21 Thread siddhesh godbole
Hello, i wanted to use mumps package for parallel direct linear solver in my spectral transformation call to ksp during EPSSolve. i ran this command as per section 3.4.1 in slepc guide. though this message tell me that i need to download mumps i still see the mumps package in my

Re: [petsc-users] MUMPS error

2015-04-21 Thread Jose E. Roman
Proceed as the error message says: add --download-mumps in the configure line, together with additional dependency packages, e.g.: $ ./configure --download-mumps --download-scalapack --download-metis --download-parmetis Jose El 21/04/2015, a las 10:21, siddhesh godbole escribió: Hello,

Re: [petsc-users] MUMPS error -13 info

2014-04-18 Thread Hong Zhang
Gong Ding, The reported bug is fixed in petsc-release https://bitbucket.org/petsc/petsc/commits/151787a63101f3b7c1ee9a4abd20e4c3fe8caf18?at=maint Thanks for your contribution! Hong On Sun, Apr 13, 2014 at 2:42 AM, Gong Ding gd...@ustc.edu wrote: Dear Sir, I see an error: Transient compute

Re: [petsc-users] MUMPS error -13 info

2014-04-13 Thread Hong Zhang
Thanks for reporting it. Do you agree with the fix below? If so, I'll patch petsc-release. Hong $ git diff ../../../../mat/impls/aij/mpi/mumps/mumps.c diff --git a/src/mat/impls/aij/mpi/mumps/mumps.c b/src/mat/impls/aij/mpi/mumps/mumps.c index f07788a..e4b2e44 100644 ---

[petsc-users] MUMPS error

2013-02-23 Thread Hui Zhang
Error reported by MUMPS in numerical factorization phase: Cannot allocate required memory 70397575 megabytes The required memory is too large. Is that normal? I'm solving a Helmholtz equation on 100^3 mesh partitioned to 4^3 processors.

[petsc-users] Mumps Error

2011-09-01 Thread Gabriele Jost
-bounces at mcs.anl.gov] on behalf of Max Rudolph [rudo...@berkeley.edu] Sent: Wednesday, August 31, 2011 3:34 PM To: petsc-users at mcs.anl.gov Subject: Re: [petsc-users] Mumps Error I'm not sure if you figured out a solution yet, but I think that you might want to run with -mat_mumps_icntl_14 100

[petsc-users] Mumps Error

2011-09-01 Thread Hong Zhang
of Max Rudolph [rudolph at berkeley.edu] Sent: Wednesday, August 31, 2011 3:34 PM To: petsc-users at mcs.anl.gov Subject: Re: [petsc-users] Mumps Error I'm not sure if you figured out a solution yet, but I think that you might want to run with -mat_mumps_icntl_14 100 Max Dear Users

[petsc-users] Mumps Error

2011-08-31 Thread Adam Byrd
Dear Users, I'm having trouble figuring out why the MUMPS solver is failing on a specific range of one of my parameters. When using the PETSc direct solver on a single processor I have no issues. There error is: [0]PETSC ERROR: - Error Message

[petsc-users] Mumps Error

2011-08-31 Thread Barry Smith
Adam, You need to use the MUMPS documentation to interpret what the error INFO(1)=-9, INFO(2)=13 means and then go from that. Barry On Aug 31, 2011, at 10:13 AM, Adam Byrd wrote: Dear Users, I'm having trouble figuring out why the MUMPS solver is failing on a specific range

[petsc-users] Mumps Error

2011-08-31 Thread Max Rudolph
I'm not sure if you figured out a solution yet, but I think that you might want to run with -mat_mumps_icntl_14 100 Max Dear Users, I'm having trouble figuring out why the MUMPS solver is failing on a specific range of one of my parameters. When using the PETSc direct solver on a