Skip to content
Snippets Groups Projects
Commit d079d64c authored by Loic Hausammann's avatar Loic Hausammann
Browse files

test if particle are active in runner_doiact_stars

parent fb480b90
No related branches found
No related tags found
3 merge requests!787Eagle stellar evolution matthieu,!781Eagle stellar evolution,!778Fix default models
...@@ -100,6 +100,9 @@ void DOSELF1_STARS(struct runner *r, struct cell *c, int timer) { ...@@ -100,6 +100,9 @@ void DOSELF1_STARS(struct runner *r, struct cell *c, int timer) {
/* Get a hold of the ith spart in ci. */ /* Get a hold of the ith spart in ci. */
struct spart *restrict si = &sparts[sid]; struct spart *restrict si = &sparts[sid];
if (!spart_is_active(si, e) || spart_is_inhibited(si, e))
continue;
const float hi = si->h; const float hi = si->h;
const float hig2 = hi * hi * kernel_gamma2; const float hig2 = hi * hi * kernel_gamma2;
const float six[3] = {(float)(si->x[0] - c->loc[0]), const float six[3] = {(float)(si->x[0] - c->loc[0]),
...@@ -180,6 +183,8 @@ void DO_NONSYM_PAIR1_STARS(struct runner *r, struct cell *restrict ci, ...@@ -180,6 +183,8 @@ void DO_NONSYM_PAIR1_STARS(struct runner *r, struct cell *restrict ci,
/* Get a hold of the ith spart in ci. */ /* Get a hold of the ith spart in ci. */
struct spart *restrict si = &sparts_i[sid]; struct spart *restrict si = &sparts_i[sid];
if (!spart_is_active(si, e) || spart_is_inhibited(si, e))
continue;
const float hi = si->h; const float hi = si->h;
const float hig2 = hi * hi * kernel_gamma2; const float hig2 = hi * hi * kernel_gamma2;
const float six[3] = {(float)(si->x[0] - (cj->loc[0] + shift[0])), const float six[3] = {(float)(si->x[0] - (cj->loc[0] + shift[0])),
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment