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 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 1437471711 >

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-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 th

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'

[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 > wrote: > Hi all, > > I tried to run my Generalized Eigenproblem in 120 x 24 = 2880 cores. > The matrix size of A = 20GB and B =

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 wrote: > Hi all,

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 -eps_n

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

2015-05-29 Thread venkatesh g
Thanks. So I what other direct solver I can use for this singular B ? QR is a good option but what to do ? I will read the manual on dumping out null space. Venkatesh On Fri, May 29, 2015 at 7:42 PM, Hong wrote: > venkatesh: > >> On Tue, May 26, 2015 at 9:02 PM, Hong wrote: >> >>> 'A serial j

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

2015-05-29 Thread Hong
venkatesh: > On Tue, May 26, 2015 at 9:02 PM, Hong 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 an

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

2015-05-29 Thread venkatesh g
On Tue, May 26, 2015 at 9:02 PM, Hong 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 for small mat

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 num