Commit 18986850 authored by Oliver Zier's avatar Oliver Zier
Browse files

Remove old TODOs

parent 76236103
...@@ -54,7 +54,6 @@ void sph_particle_data::set_velocity_gradients(void) ...@@ -54,7 +54,6 @@ void sph_particle_data::set_velocity_gradients(void)
} }
else else
{ {
// TODO check
DivVel = (dvel[0][0] + dvel[1][1]) / Density; DivVel = (dvel[0][0] + dvel[1][1]) / Density;
CurlVel = fabs((dvel[1][0] - dvel[0][1]) / Density); CurlVel = fabs((dvel[1][0] - dvel[0][1]) / Density);
} }
...@@ -142,7 +141,7 @@ void sph_particle_data::set_viscosity_coefficient(double dt) ...@@ -142,7 +141,7 @@ void sph_particle_data::set_viscosity_coefficient(double dt)
double shockIndicator = -dDivVel_dt > 0 ? -dDivVel_dt : 0; double shockIndicator = -dDivVel_dt > 0 ? -dDivVel_dt : 0;
double hsml = Hsml * All.cf_atime; double hsml = Hsml * All.cf_atime;
double Hsml2 = hsml * hsml; double Hsml2 = hsml * hsml;
double alpha_tar = (Hsml2 * shockIndicator) / (Hsml2 * shockIndicator + Csnd * Csnd) * All.ArtBulkViscConst; // TODO check max double alpha_tar = (Hsml2 * shockIndicator) / (Hsml2 * shockIndicator + Csnd * Csnd) * All.ArtBulkViscConst;
double DivVel2 = DivVel * DivVel; double DivVel2 = DivVel * DivVel;
double CurlVel2 = CurlVel * CurlVel; double CurlVel2 = CurlVel * CurlVel;
......
...@@ -567,20 +567,6 @@ void sim::hydro_force(int step_indicator) ...@@ -567,20 +567,6 @@ void sim::hydro_force(int step_indicator)
* routine) * routine)
*/ */
} }
#ifdef TIMEDEP_ART_VISC
/* double csnd_over_h = Sp.SphP[i].Csnd / Sp.SphP[i].Hsml;
double f = fabs(Sp.SphP[i].DivVel) / (fabs(Sp.SphP[i].DivVel) + Sp.SphP[i].CurlVel + 0.0001 * csnd_over_h /
NgbTree.fac_mu);
Sp.SphP[i].Dtalpha =
-(Sp.SphP[i].alpha - All.AlphaMin) * All.DecayTime * 0.5 * Sp.SphP[i].MaxSignalVel / (Sp.SphP[i].Hsml *
NgbTree.fac_mu) + f * All.ViscSource * std::max<double>(0.0, -Sp.SphP[i].DivVel); if(All.ComovingIntegrationOn)
{
Sp.SphP[i].Dtalpha /= (All.cf_hubble_a * All.Time * All.Time);
}*/
// TODO check
#endif
} }
if(step_indicator == SECOND_HALF_STEP) if(step_indicator == SECOND_HALF_STEP)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment