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

Make the drift functions ignore the inhibited particles.

parent 18b6cce3
......@@ -2895,6 +2895,9 @@ void cell_drift_part(struct cell *c, const struct engine *e, int force) {
struct part *const p = &parts[k];
struct xpart *const xp = &xparts[k];
/* Ignore inhibited particles */
if (part_is_inhibited(p, e)) continue;
/* Drift... */
drift_part(p, xp, dt_drift, dt_kick_hydro, dt_kick_grav, dt_therm,
ti_old_part, ti_current);
......@@ -3037,6 +3040,9 @@ void cell_drift_gpart(struct cell *c, const struct engine *e, int force) {
/* Get a handle on the gpart. */
struct gpart *const gp = &gparts[k];
/* Ignore inhibited particles */
if (gpart_is_inhibited(gp, e)) continue;
/* Drift... */
drift_gpart(gp, dt_drift, ti_old_gpart, ti_current);
......@@ -3082,6 +3088,9 @@ void cell_drift_gpart(struct cell *c, const struct engine *e, int force) {
/* Get a handle on the spart. */
struct spart *const sp = &sparts[k];
/* Ignore inhibited particles */
if (spart_is_inhibited(sp, e)) continue;
/* Drift... */
drift_spart(sp, dt_drift, ti_old_gpart, ti_current);
......
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