Changeset View
Changeset View
Standalone View
Standalone View
intern/elbeem/intern/solver_adap.cpp
| Context not available. | |||||
| #include "solver_relax.h" | #include "solver_relax.h" | ||||
| #include "particletracer.h" | #include "particletracer.h" | ||||
| #include <cmath> | |||||
| using std::isfinite; | |||||
| /*****************************************************************************/ | /*****************************************************************************/ | ||||
| //! coarse step functions | //! coarse step functions | ||||
| Context not available. | |||||
| uz += (dfDvecZ[l]*m); | uz += (dfDvecZ[l]*m); | ||||
| } | } | ||||
| #ifndef WIN32 | #ifndef WIN32 | ||||
| if (!finite(rho)) { | if (!isfinite(rho)) { | ||||
| errMsg("adaptTimestep","Brute force non-finite rho at"<<PRINT_IJK); // DEBUG! | errMsg("adaptTimestep","Brute force non-finite rho at"<<PRINT_IJK); // DEBUG! | ||||
| rho = 1.0; | rho = 1.0; | ||||
| ux = uy = uz = 0.0; | ux = uy = uz = 0.0; | ||||
| Context not available. | |||||