error: Floating point exception
 

Hi everyone,

I am trying to implement level Set Method in OpenFoam, but I got this error in pEqn.H.

I have not changed anything in "pEqn.H" or related to that.


Any advice?


"

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at /opt/openfoam220/src/finiteVolume/lnInclude/fvMatrixSolve.C:81
#8
at ~/OpenFOAM/suhas-2.2.0/applications/solvers/levelSetFoam/pEqn.H:37
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
in "/home/suhas/OpenFOAM/suhas-2.2.0/platforms/linux64GccDPOpt/bin/levelSetFoam"
Floating point exception

"


Thanks,

Hi,

I can recommend to build a OF-Version with debugging enabled. That helped me with the same kind of error

http://openfoamwiki.net/index.php/Ho...UG_-_libraries

your problem might be a division by zero, or infinity or a multiplication of fields of different size. I also thought the problem was in the pEqn.H but i turned out it was a lot earlier..

Greetings
Nicklas

Hi Nicklas,


Thank you for the reply. I found out the error using nemiver. There was a mistake in defining Heaviside function which gave an error in pEqn.H. Now everything's working fine!!..:)



Thanks,