Error codes from external packages

[Matthew Knepley]

On Mon, Jun 2, 2008 at 9:22 AM, Lars Rindorf
<Lars.Rindorf at teknologisk.dk> wrote:
> 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.

Ah, I was checking solve(), not factor(). This is an oversight in the
code. I am fixing it in the
dev version. We will have a release fairly soon, but you can always
get dev for these kinds of
bugs fixes quickly.

   Matt

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



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