inline void diffuseUV(const double _diff)
{
const double a = _dt * _diff * NX * NY;
linearSolverUV(a, 1.0 + 4 * a);
}
void ciMsaFluidSolver::diffuseUV( float diff )
{
float a = _dt * diff * _NX * _NY;
linearSolverUV( a, 1.0f + 4 * a );
}
void FluidSolver::diffuseUV( float diff )
{
float a = deltaT * diff * _NX * _NY;
linearSolverUV( a, 1.0 + 4 * a );
}