[petsc-users] Use of superlu_dis

Jed Brown jedbrown at mcs.anl.gov
Fri Jan 4 09:03:28 CST 2013


0. ALWAYS send the ENTIRE error message.

1. Use a --with-debugging=1 (the default) build, which will likely check
earlier for NaN

2. Run in a debugger with -fp_trap.


On Fri, Jan 4, 2013 at 8:54 AM, w_ang_temp <w_ang_temp at 163.com> wrote:

> Hello,
>     I have a problem about superlu_dis.
>    The result of Ax=b is NaN. The error information is:
>
>  ****** Warning from MC64A/AD. INFO(1) =  2
> Linear solve converged due to CONVERGED_ITS iterations 1
> [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [0]PETSC ERROR: Floating point exception!
> [0]PETSC ERROR: Infinite or not-a-number generated in norm!
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 7, Thu Mar 15 09:30:51
> CDT 2012
> [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: ./ex4f on a arch-linu named node4 by wang Fri Jan  4
> 20:48:16 2013
> [0]PETSC ERROR: Libraries linked from
> /public/soft/ddc/soft/petsc/petsc-3.2-p7/arch-linux2-c-opt/lib
> [0]PETS C ERROR: Configure run at Tue Jan  1 22:23:45 2013
> [0]PETSC ERROR: Configure options
> --with-mpi-dir=/public/soft/ddc/soft/mpich2/ --download-f-blas-lapack=1
> --download-hypre=1 --with-x=1 --with-debugging=0 --download-superlu_dist
> --download-parmetis --download-mumps --download-scalapack
> --download-blacs=/public/soft/ddc/soft/petsc/blacs-dev.tar.gz
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: VecNorm() line 167 in src/vec/vec/interface/rvector.c
>
>
>     Does it mean that the matrix is singular? Or the size of the matrix
> (in the project, it is 250000X250000) is too large?
>     Thanks.
>                                                           Jim
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130104/5830f108/attachment.html>


More information about the petsc-users mailing list