CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Floating point error

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   March 12, 2009, 07:17
Default I tried to run snappyHexMesh o
Senior Member
Prapanch Nair
Join Date: Mar 2009
Location: Bangalore, India
Posts: 105
Rep Power: 8
prapanj is on a distinguished road
I tried to run snappyHexMesh on a Geometry of a building. And I got the following as output with a floating point exception error.

Determining initial surface intersections

Edge intersection testing:
Number of edges : 19875
Number of edges to retest : 19875
#0 Foam::error::printStack(Foam:stream&) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#2 Uninterpreted: [0xb8060400]
#3 Foam::triangleFuncs::intersectAxesBundle(Foam::Vec tor<double> const&, Foam::Vector<double> const&, Foam::Vector<double> const&, int, Foam::Field<foam::vector<double> > const&, double, Foam::Vector<double>&) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#4 Foam::triangleFuncs::intersectBb(Foam::Vector<doub le> const&, Foam::Vector<double> const&, Foam::Vector<double> const&, Foam::treeBoundBox const&) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#5 Foam::treeDataTriSurface::overlaps(int, Foam::treeBoundBox const&) const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#6 Foam::indexedOctree<foam::treedatatrisurface>::div ide(Foam::List<int> const&, Foam::treeBoundBox const&, Foam::List<foam::list<int> >&) const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#7 Foam::indexedOctree<foam::treedatatrisurface>::div ide(Foam::treeBoundBox const&, Foam::DynamicList<foam::list<int>, 0u, 2u, 1u>&, int) const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#8 Foam::indexedOctree<foam::treedatatrisurface>::spl itNodes(int, Foam::DynamicList<foam::indexedoctree<foam::treeda tatrisurface>::node, 0u, 2u, 1u>&, Foam::DynamicList<foam::list<int>, 0u, 2u, 1u>&) const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#9 Foam::indexedOctree<foam::treedatatrisurface>::ind exedOctree(Foam::treeDataTriSu rface const&, Foam::treeBoundBox const&, int, double, double) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#10 Foam::triSurfaceMesh::tree() const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#11 Foam::triSurfaceMesh::findLineAny(Foam::Field<foam ::vector<double> > const&, Foam::Field<foam::vector<double> > const&, Foam::List<foam::pointindexhit<foam::vector<double > > >&) const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#12 Foam::refinementSurfaces::findHigherIntersection(F oam::Field<foam::vector<double > > const&, Foam::Field<foam::vector<double> > const&, Foam::List<int> const&, Foam::List<int>&, Foam::List<int>&) const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#13 Foam::meshRefinement::updateIntersections(Foam::Li st<int> const&) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#14 Foam::meshRefinement::meshRefinement(Foam::fvMesh& , double, Foam::refinementSurfaces const&, Foam::shellSurfaces const&) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/"
#15 main in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/snappyHexMesh"
#16 __libc_start_main in "/lib/tls/i686/cmov/"
#17 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/snappyHexMesh"
Floating point exception

I had set FOAM_SIGFPE=0, but got the same error. Could someone please explain how I may set the value for FOAM_SIGFPE? I tried a different geometry in the same block and snappy was running fine. I had checked the STL in other softwares and found it OK. Please suggest a solution . Thank you
prapanj is offline   Reply With Quote

Old   March 12, 2009, 08:57
Default Well, I figured out the proble
Senior Member
Prapanch Nair
Join Date: Mar 2009
Location: Bangalore, India
Posts: 105
Rep Power: 8
prapanj is on a distinguished road
Well, I figured out the problem myself. SnappyHexMesh has a bug. If the stl file has a face that is coplanar with any coordinate plane, xy, xz, yz, Snappy crashes. Thank you
prapanj is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Floating point error - help sk FLUENT 2 December 17, 2008 11:32
Floating point error Chico Mbanu FLUENT 0 July 14, 2004 13:56
floating point error hussain FLUENT 2 June 15, 2004 03:12
floating point error Mavinakere FLUENT 3 March 1, 2002 18:22
floating point error sivakuamr Main CFD Forum 12 August 31, 2001 01:37

All times are GMT -4. The time now is 12:23.