CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   Pow in lib64tlslibmso6 SigFpe when running coodles with SpalartAllmaras (http://www.cfd-online.com/Forums/openfoam-bugs/62526-pow-lib64tlslibmso6-sigfpe-when-running-coodles-spalartallmaras.html)

lillberg December 6, 2007 02:03

Hi forum and Henry! I've ru
 
Hi forum and Henry!

I've run into this problem when running coodles in parallel using the Spalart-Allmaras DES model.

[3] #0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpe
nFOAM.so"
[3] #1 Foam::sigFpe::sigFpeHandler(int) in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.s
o"
[3] #2 ?? in "/lib64/tls/libc.so.6"
[3] #3 ?? in "/lib64/tls/libm.so.6"
[3] #4 ?? in "/lib64/tls/libm.so.6"
[3] #5 pow in "/lib64/tls/libm.so.6"
[3] #6 Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) in "/local/users/exterli/OpenFOAM/OpenFO
AM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
[3] #7 void Foam::pow<foam::fvpatchfield,>(Foam::GeometricFiel d<double,>&,
Foam::GeometricField<double,> const&, Foam::dimensioned<double> const&) in "/local/users
/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleLESmodels.so "
[3] #8 Foam::tmp<foam::geometricfield<double,> > Foam::pow<Foam::fvPatchField, Foam::vol
Mesh>(Foam::GeometricField<double,> const&, Foam::dimensioned<double> const&) in "/local/
users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleLESmode ls.so"
[3] #9 Foam::tmp<foam::geometricfield<double,> > Foam::pow<Foam::fvPatchField, Foam::vol
Mesh>(Foam::GeometricField<double,> const&, double const&) in "/local/users/exterli/OpenF
OAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibleLESmodels.so"
[3] #10 Foam::compressible::LESmodels::SpalartAllmaras::fw (Foam::GeometricField<double,>
const&) const in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibl eLESmodels.so"
[3] #11 Foam::compressible::LESmodels::SpalartAllmaras::co rrect(Foam::tmp<foam::geometri cfield<foam::tensor<double>, Foam:
:fvPatchField, Foam::volMesh> > const&) in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libcompressibl
eLESmodels.so"
[3] #12 Foam::compressible::LESmodel::correct() in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libco
mpressibleLESmodels.so"
[3] #13 main in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/c oodles"
[3] #14 __libc_start_main in "/lib64/tls/libc.so.6"
[3] #15 Foam::lduMatrix::lduMatrix(Foam::lduMatrix&, bool) in "/local/users/exterli/OpenFOAM/OpenFOAM-1.4.1/applications/b
in/linux64GccDPOpt/coodles"
srun: error: n107: task3: Floating point exception (core dumped)

Running the same mesh with the same code using any other LES model works just fine. Tried setting floatTransfer to 0 but it didn't change anything. Printed the max/min values of what goes into fw but couldn't see anything unusual. The simulation runs for some 10.000 timesteps and the solution look Ok.

The cluster has an older kernel and glibc and the problem does not occur on my laptop which is running openSuSE 10.3. Hence, this is likely a system specific problem. Anyhow, I would like to hear if someone has any idea how to overcome this. (except upgrading the cluster nodes, 400 cpu's and I'm not admin :-( )

Best Regards

//Eric

henry December 6, 2007 03:06

Hi Eric, Try with this htt
 
Hi Eric,

Try with this http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif SpalartAllmaras.C

Henry

lillberg December 7, 2007 01:50

Works like a charm, thanks a l
 
Works like a charm, thanks a lot Henry!

Just a question: Why bounding r by 10?

//Eric

henry December 7, 2007 07:11

r must be clipped otherwise th
 
r must be clipped otherwise the (r^6)^6 causes a floating point exception. Clipping r with 10 ensures that (r^6)^6 does not cause over-flow in either single or double precision.

Henry

lillberg December 7, 2007 08:17

Thanks again Henry! //Eric
 
Thanks again Henry!

//Eric


All times are GMT -4. The time now is 09:41.