Skip to content
Snippets Groups Projects
Commit b9760d6d authored by Jacob Kegerreis's avatar Jacob Kegerreis
Browse files

Undo same-position error check

parent b9509b00
Branches
Tags
1 merge request!545Add support for equations of state related to planetary physics
......@@ -707,16 +707,6 @@ static INLINE void runner_doself_grav_pp_full(struct runner *r,
swift_align_information(float, ci_cache->m, SWIFT_CACHE_ALIGNMENT);
swift_assume_size(gcount_padded, VEC_SIZE);
// int index = -gparts[pid].id_or_neg_offset; ///###
// if (e->s->parts[index].id == 7708) { ///###
// message("a");
// printParticle_single(&e->s->parts[index], &e->s->xparts[index]);
// }
// if (e->s->parts[index].id == 7727) { ///###
// message("b");
// printParticle_single(&e->s->parts[index], &e->s->xparts[index]);
// }
/* Loop over every other particle in the cell. */
for (int pjd = 0; pjd < gcount_padded; pjd++) {
......@@ -735,19 +725,6 @@ static INLINE void runner_doself_grav_pp_full(struct runner *r,
const float dz = z_i - z_j;
const float r2 = dx * dx + dy * dy + dz * dz;
if (r2 == 0) { ///###
message("Particles at identical positions!");
int index;
index = -gparts[pid].id_or_neg_offset;
printParticle_single(&e->s->parts[index], &e->s->xparts[index]);
index = -gparts[pjd].id_or_neg_offset;
printParticle_single(&e->s->parts[index], &e->s->xparts[index]);
error("");
}
#ifdef SWIFT_DEBUG_CHECKS
if (r2 == 0.f && h2_i == 0.)
error("Interacting particles with 0 distance and 0 softening.");
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment