Commit 53bc63dc authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Better name for temporary variable

parent 3348cd78
......@@ -664,7 +664,7 @@ void runner_dodrift(struct runner *r, struct cell *c, int timer) {
const float ti_current = r->e->ti_current;
struct part *restrict p, *restrict parts = c->parts;
struct xpart *restrict xp, *restrict xparts = c->xparts;
float dx_max = 0.f, h_max = 0.f;
float dx_max = 0.f, dx2_max = 0.f, h_max = 0.f;
float w;
TIMER_TIC
......@@ -710,17 +710,17 @@ void runner_dodrift(struct runner *r, struct cell *c, int timer) {
hydro_predict_extra(p, xp, ti_old, ti_current, timeBase);
/* Compute (square of) motion since last cell construction */
const float dx = (p->x[0] - xp->x_old[0]) * (p->x[0] - xp->x_old[0]) +
(p->x[1] - xp->x_old[1]) * (p->x[1] - xp->x_old[1]) +
(p->x[2] - xp->x_old[2]) * (p->x[2] - xp->x_old[2]);
dx_max = fmaxf(dx_max, dx);
const float dx2 = (p->x[0] - xp->x_old[0]) * (p->x[0] - xp->x_old[0]) +
(p->x[1] - xp->x_old[1]) * (p->x[1] - xp->x_old[1]) +
(p->x[2] - xp->x_old[2]) * (p->x[2] - xp->x_old[2]);
dx2_max = fmaxf(dx2_max, dx2);
/* Maximal smoothing length */
h_max = fmaxf(p->h, h_max);
}
/* Now, get the maximal particle motion from its square */
dx_max = sqrtf(dx_max);
dx_max = sqrtf(dx2_max);
}
/* Otherwise, aggregate data from children. */
......
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