Differential D4649 Diff 14611 intern/elbeem/intern/solver_adap.cpp

Changeset View

Standalone View

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.