This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

PETSc error code is: 72.

0 votes

Hello everybody

I ' am trying to solve a problem depending on velocity, pressure and temperature but it does not work, here is the error PETSc error code is: 72.

Solving nonlinear variational problem.
Newton iteration 0: r (abs) = 1.483e+02 (tol = 1.000e-10) r (rel) = 1.000e+00 (tol = 1.000e-09)
[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.3.0, Patch 3, Wed Aug 29 11:26:24 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: Unknown Name on a darwin14. named Ndongos-MacBook-Air.local by ndongofall Fri May 22 11:26:20 2015
[0]PETSC ERROR: Libraries linked from /Users/johannr/fenics-1.5.0/local/lib
[0]PETSC ERROR: Configure run at Mon Jan 12 22:17:07 2015
[0]PETSC ERROR: Configure options --prefix=/Users/johannr/fenics-1.5.0/local COPTFLAGS=-O2 --with-debugging=0 --with-clanguage=cxx --with-c-support=1 --with-blas-lapack-dir=/usr --with-umfpack=1 --with-umfpack-include=/Users/johannr/fenics-1.5.0/local/include/suitesparse --with-umfpack-lib="[/Users/johannr/fenics-1.5.0/local/lib/libumfpack.a,/Users/johannr/fenics-1.5.0/local/lib/libamd.a]" --with-spooles=1 --with-spooles-include=/Users/johannr/fenics-1.5.0/local/include --with-spooles-lib=/Users/johannr/fenics-1.5.0/local/lib/libspooles.a --with-ptscotch=1 --with-ptscotch-dir=/Users/johannr/fenics-1.5.0/local --with-ml=1 --with-ml-include=/Users/johannr/fenics-1.5.0/local/include/trilinos --with-ml-lib=/Users/johannr/fenics-1.5.0/local/lib/libml.dylib --with-hdf5=1 --with-hdf5-dir=/Users/johannr/fenics-1.5.0/local --with-x=0 -with-x11=0 --with-fortran=0 --with-shared-libraries=1 PETSC_DIR=/Users/johannr/fenics-1.5.0/fenics-superbuild/build-fenics/CMakeExternals/src/PETSc PETSC_ARCH=darwin14.0.0-cxx-opt
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: VecNorm() line 169 in /Users/johannr/fenics-1.5.0/fenics-superbuild/build-fenics/CMakeExternals/src/PETSc/src/vec/vec/interface/rvector.c
Traceback (most recent call last):
File "steady.py", line 5, in
problem.steady_state([50])
File "/Users/ndongofall/Desktop/Stage/2D/ns_Cartesian.py", line 44, in steady_state
solver.solve()
RuntimeError:

*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at


*** fenics@fenicsproject.org


*** Remember to include the error message listed below and, if possible,
*** include a minimal running example to reproduce the error.


*** -------------------------------------------------------------------------
*** Error: Unable to successfully call PETSc function 'VecNorm'.
*** Reason: PETSc error code is: 72.
*** Where: This error was encountered inside /Users/johannr/fenics-1.5.0/fenics-superbuild/build-fenics/CMakeExternals/src/DOLFIN/dolfin/la/PETScVector.cpp.
*** Process: unknown


*** DOLFIN version: 1.5.0
*** Git changeset: f467b66dcfd821ec20e9f9070c7cef5a991dbc42
*** -------------------------------------------------

I really appreciate a help !

asked May 22, 2015 by fall FEniCS Novice (230 points)

1 Answer

+1 vote

Hi,

It seems that there is a division by 0 or something similar in the Newton solver implied by the way you implemented the variational problem. Could you provide some code sample in order to see where the error comes from?

answered May 22, 2015 by MathieuFV FEniCS User (1,490 points)

could you send me your e-mail address in order to receive the entire code, it's so long, i try to make a code which can determine a steady of a however 2D physical problem (like boussinesq, flow past a cylinder etc ...)( ndantyfall@icloud.com)

Maybe you could send it to me by private message? Not that I don't want to provide my email address but it definitely will go faster on the fenics q&a platform...

...