SV: Error codes from external packages

[Lars Rindorf]

Hi again

I solved the problem! I changed the petsc source code, to show the error code. It is a memory allocation error in malloc inside umfpack. 

Thanks anyway.
Lars 

-----Oprindelig meddelelse-----
Fra: Lars Rindorf 
Sendt: 2. juni 2008 16:22
Til: 'petsc-users at mcs.anl.gov'
Emne: SV: Error codes from external packages

Hi Matthew

I have included the -mat_umfpack_prl parameter, but it does not make any difference. I have checked the spelling with petsc manual. When umfpack crashes it (umfpack_di_numeric) returns an error code. I want to access that code. 

Here is what petsc returns (with/without -mat_umfpack_prl option):
[lhr at localhost notch_patch]$ getdp patch_notch.pro -msh patch_notch_29000.msh -pre res_static -cal  -ksp_type preonly -pc_type lu -mat_type umfpack -mat_umfpack_prl P r e - P r o c e s s i n g . . .
E n d   P r e - P r o c e s s i n g
P r o c e s s i n g . . .
Operation : Generate[A]
Info      : Setting System {A,b} to zero
Resources : cpu 6.41402 s
Operation : Solve[A]
PETSc     : N: 136188
[0]PETSC ERROR: --------------------- Error Message ------------------------------------
[0]PETSC ERROR: Error in external library!
[0]PETSC ERROR: umfpack_di_numeric failed!
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 13, Thu May 15 17:29:26 CDT 2008 HG revision: 4466c6289a0922df26e20626fd4a0b4dd03c8124
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: patch_notch.pro on a linux-gnu named localhost.localdomain by lhr Mon Jun  2 16:03:50 2008 [0]PETSC ERROR: Libraries linked from /home/lhr/Desktop/getdp/petsc/petsc-2.3.3-p13/lib/linux-gnu-c-debug
[0]PETSC ERROR: Configure run at Mon Jun  2 15:58:42 2008 [0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=g77 --with-blas-lapack-dir=/opt/intel/mkl/10.0.1.014/lib/em64t/ --download-mpich=ifneeded --download-umfpack=ifneeded --with-shared=0 [0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: MatLUFactorNumeric_UMFPACK() line 129 in src/mat/impls/aij/seq/umfpack/umfpack.c
[0]PETSC ERROR: MatLUFactorNumeric() line 2227 in src/mat/interface/matrix.c [0]PETSC ERROR: PCSetUp_LU() line 280 in src/ksp/pc/impls/factor/lu/lu.c [0]PETSC ERROR: PCSetUp() line 787 in src/ksp/pc/interface/precon.c [0]PETSC ERROR: KSPSetUp() line 234 in src/ksp/ksp/interface/itfunc.c [0]PETSC ERROR: KSPSolve() line 347 in src/ksp/ksp/interface/itfunc.c [0]PETSC ERROR: User provided function() line 1472 in unknowndirectory/LinAlg_PETSC.c
[unset]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 76) - process 0

KR, Lars


-----Oprindelig meddelelse-----
Fra: owner-petsc-users at mcs.anl.gov [mailto:owner-petsc-users at mcs.anl.gov] På vegne af Matthew Knepley
Sendt: 2. juni 2008 15:07
Til: petsc-users at mcs.anl.gov
Emne: Re: Error codes from external packages

On Mon, Jun 2, 2008 at 7:17 AM, Lars Rindorf <Lars.Rindorf at teknologisk.dk> wrote:
> Dear all
>
> If I want to know the error codes of an external package that crashes, 
> then what can I do?

We call

  umfpack_zl_report_status()

with the UMFPACK status code in the event of a failure. Printing for this is controlled by the -mat_umfpack_prl option I believe.

   Matt

> The problem arises with umfpack when the size of the matrix is more 
> than 118000x118000, corresponding to 2.4 Gb memory consumption. It 
> simply returns "umfpack_di_numeric" (factorization of a real matrix) failed.
> Has anybody else experienced this?
>
> KR, Lars
>
>



--
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener