CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Floating point error (https://www.cfd-online.com/Forums/openfoam-solving/57801-floating-point-error.html)

prapanj March 12, 2009 06:17

I tried to run snappyHexMesh o
 
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) 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 March 12, 2009 07:57

Well, I figured out the proble
 
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


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