Changeset View
Changeset View
Standalone View
Standalone View
intern/elbeem/intern/solver_main.cpp
| Context not available. | |||||
| #include "globals.h" | #include "globals.h" | ||||
| #include <stdlib.h> | #include <stdlib.h> | ||||
| #include <cmath> | |||||
| using std::isfinite; | |||||
| /*****************************************************************************/ | /*****************************************************************************/ | ||||
| /*! perform a single LBM step */ | /*! perform a single LBM step */ | ||||
| Context not available. | |||||
| #ifndef WIN32 | #ifndef WIN32 | ||||
| // good indicator for instabilities | // good indicator for instabilities | ||||
| if( (!finite(mMxvx)) || (!finite(mMxvy)) || (!finite(mMxvz)) ) { CAUSE_PANIC; } | if( (!isfinite(mMxvx)) || (!isfinite(mMxvy)) || (!isfinite(mMxvz)) ) { CAUSE_PANIC; } | ||||
| if( (!finite(mCurrentMass)) || (!finite(mCurrentVolume)) ) { CAUSE_PANIC; } | if( (!isfinite(mCurrentMass)) || (!isfinite(mCurrentVolume)) ) { CAUSE_PANIC; } | ||||
| #endif // WIN32 | #endif // WIN32 | ||||
| // output total step time | // output total step time | ||||
| Context not available. | |||||
| if(!this->mPanic){ FSGR_FORIJK_BOUNDS(mMaxRefine) { \ | if(!this->mPanic){ FSGR_FORIJK_BOUNDS(mMaxRefine) { \ | ||||
| if(RFLAG(mMaxRefine,i,j,k,mLevel[mMaxRefine].setCurr)&(CFFluid|CFInter)) { \ | if(RFLAG(mMaxRefine,i,j,k,mLevel[mMaxRefine].setCurr)&(CFFluid|CFInter)) { \ | ||||
| for(int l=0;l<dTotalNum;l++) { \ | for(int l=0;l<dTotalNum;l++) { \ | ||||
| if(!finite(QCELL(mMaxRefine,i,j,k,mLevel[mMaxRefine].setCurr,l))) { errMsg("NNOFIN"," "<<str<<" at "<<PRINT_IJK<<" l"<<l<<" "); }\ | if(!isfinite(QCELL(mMaxRefine,i,j,k,mLevel[mMaxRefine].setCurr,l))) { errMsg("NNOFIN"," "<<str<<" at "<<PRINT_IJK<<" l"<<l<<" "); }\ | ||||
| }/*for*/ \ | }/*for*/ \ | ||||
| }/*if*/ \ | }/*if*/ \ | ||||
| } } | } } | ||||
| Context not available. | |||||