The infamous floating point exception
I tried to modify the sonicFoam solver and it went through wmake with no errors. When I run it, I get:
Code:
/*---------------------------------------------------------------------------*\ Any help will be greatly appreciated. :confused: |
Hi Costas
It is very hard to say, as it is a custom solver, however look at error #4, and your are explicitly told that it has something to do with the division operator. You are probably dividing by zero somewhere in the source. Bests, Niels |
Hi every1,
any concrete solution to the above problem? I have almost a similar problem but a little different. I am posting the failure below.. Create time Create fluid mesh for region fluid for time = 0 Create solid mesh for region aluminium for time = 0 *** Reading fluid mesh thermophysical properties for region fluid Adding to thermoFluid Selecting thermodynamics package hsRhoThermo<pureMixture<constTransport<specieTherm o<hConstThermo<perfectGas>>>>> #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 __restore_rt at sigaction.c:0 #3 Foam::hsRhoThermo<Foam::pureMixture<Foam::constTra nsport<Foam::specieThermo<Foam::hConstThermo<Foam: :perfectGas> > > > >::calculate() in "/opt/openfoam/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #4 Foam::hsRhoThermo<Foam::pureMixture<Foam::constTra nsport<Foam::specieThermo<Foam::hConstThermo<Foam: :perfectGas> > > > >::hsRhoThermo(Foam::fvMesh const&) in "/opt/openfoam/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #5 Foam::basicRhoThermo::addfvMeshConstructorToTable< Foam::hsRhoThermo<Foam::pureMixture<Foam::constTra nsport<Foam::specieThermo<Foam::hConstThermo<Foam: :perfectGas> > > > > >::New(Foam::fvMesh const&) in "/opt/openfoam/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #6 Foam::basicRhoThermo::New(Foam::fvMesh const&) in "/opt/openfoam/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so" #7 main in "/opt/openfoam/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" #8 __libc_start_main in "/lib64/libc.so.6" #9 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/opt/openfoam/OpenFOAM-2.0.1/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam" Fließkommafehler I tried the 'unset' fro $FOAM_SIGFPE but doesnt seem to work at all. Any ideas? IS there a problem with the libraries or the installation? |
Quote:
Nevertheless when I look at the stacktrace my guess is that we're looking at the "YouCan'tHavePressureZeroInACompressibleCase"-problem (so it is a setup problem) So: start the compilation of a Debug-version. While that is compiling you'll have time to look at your initial conditions ... |
Thank you for your response Mr.gschaider .... I am not trying to write a sovler myself, nut unfortunately I am receivng the same error as the one stated above, one of the reasons why I posted my situation here.
I use chtMultiRegionSimpleFoam as the std solver without modifications and also I use an incompressible case. I will definitely try with the debug version |
Quote:
Your problem may look similar but is not the same as the other: the stack trace is not the problem but a tool to help find the location of the actual problem and your stack trace and the one of the original poster show totally different locations in the important frames #3 and #4 |
Hi raghu
were you able to get rid of the floating point exception? what was the problem? Thanks ak |
Quote:
if you are refering to the problem I stated (please have a look at the stack trace once again as suggested by gschaider), then the solution to that is to input the perfect boundary conditions for temperature (if constant wall heat flux doesnt work, try constant wall temperature) in 0 folder. And please make doubly sure how the thermophysical properties have to be used (the syntax too). And yes, my problem is solved by doing the above. If all the above does not work, try to reduce the under-relaxation factor for temp. |
All times are GMT -4. The time now is 20:53. |