CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Floating point error

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 12, 2009, 07:17
Default I tried to run snappyHexMesh o
  #1
Senior Member
 
Prapanch Nair
Join Date: Mar 2009
Location: Bangalore, India
Posts: 105
Rep Power: 17
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/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#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/libmeshTools.so"
#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/libmeshTools.so"
#5 Foam::treeDataTriSurface::overlaps(int, Foam::treeBoundBox const&) const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libmeshTools.so"
#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/libmeshTools.so"
#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/libmeshTools.so"
#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/libmeshTools.so"
#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/libmeshTools.so"
#10 Foam::triSurfaceMesh::tree() const in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libmeshTools.so"
#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/libmeshTools.so"
#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/libautoMesh.so"
#13 Foam::meshRefinement::updateIntersections(Foam::Li st<int> const&) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libautoMesh.so"
#14 Foam::meshRefinement::meshRefinement(Foam::fvMesh& , double, Foam::refinementSurfaces const&, Foam::shellSurfaces const&) in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libautoMesh.so"
#15 main in "/home/rwdi/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/snappyHexMesh"
#16 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#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
  #2
Senior Member
 
Prapanch Nair
Join Date: Mar 2009
Location: Bangalore, India
Posts: 105
Rep Power: 17
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

Reply


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 Off
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 14:56
floating point error hussain FLUENT 2 June 15, 2004 04:12
floating point error Mavinakere FLUENT 3 March 1, 2002 18:22
floating point error sivakuamr Main CFD Forum 12 August 31, 2001 02:37


All times are GMT -4. The time now is 04:10.