[wenxu]

Dear everyone.

Now I am calculate B-PDF (or B distribution), when it run for a while, it crashed. I do not know it's a OF problem or a C++ problem.

The Equation is as follows:

253 {
254 Info<<" Gamma: "<<Gamma(pdfBeta[i])<<endl;
255 Info<<"pdfAlpha[i]: "<<pdfAlpha[i]<<" pdfBeta[i]: "<<pdfBeta[i]<<endl;
256 PDF[j] = std::Pow(hZ_[j],(pdfAlpha[i]-1e0)) * std::Pow((1e0 - hZ_[j]),(pdfBeta[i]-1.0)) * min(Gam ma(pdfAlpha[i] + pdfBeta[i]),1e17)/(min(Gamma(pdfAlpha[i]),1e17)*min(Gamma(pdfBeta[i]),1e17));
257 }

While the error is:

pdfAlpha[i]: 0.781371 pdfBeta[i]: 1.52441
Gamma: 0.887262
pdfAlpha[i]: 0.781371 pdfBeta[i]: 1.52441
Gamma: 0.887262
pdfAlpha[i]: 0.781371 pdfBeta[i]: 1.52441
#0 Foam::error::PrintStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 powl in "/lib/i386-linux-gnu/libm.so.6"
#4
at /usr/include/c++/4.6/cmath:354
#5
at /usr/include/c++/4.6/cmath:378
#6

Could anyone give me some hints, any suggestion is appreciate!

regards,
wenxu