Commit de688b82 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Simplified drift

parent ba4e72b1
......@@ -100,8 +100,9 @@ void printgParticle(struct gpart *gparts, long long int id, size_t N) {
void printParticle_single(struct part *p, struct xpart *xp) {
printf("## Particle: id=%lld", p->id);
printf("## Particle: id=%lld ", p->id);
hydro_debug_particle(p, xp);
printf("\n");
}
#ifdef HAVE_METIS
......
......@@ -82,7 +82,7 @@ __attribute__((always_inline))
p->h *= expf(w1);
/* Predict density */
const float w2 = -3.0f * p->h_dt * h_inv * dt;
const float w2 = -3.0f * w1;
if (fabsf(w2) < 0.2f)
p->rho *= approx_expf(w2); /* 4th order expansion of exp(w) */
else
......
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