diff --git a/src/runner_doiact_functions_sinks.h b/src/runner_doiact_functions_sinks.h index fd293d886c2b6cfa250718bba28ebb73eee52669..e5cd24769b8537b508f40677f9ae65e957c52dd4 100644 --- a/src/runner_doiact_functions_sinks.h +++ b/src/runner_doiact_functions_sinks.h @@ -292,10 +292,26 @@ void DO_NONSYM_PAIR1_SINKS_NAIVE(struct runner *r, struct cell *restrict ci, #ifdef SWIFT_DEBUG_CHECKS /* Check that particles have been drifted to the current time */ - if (si->ti_drift != e->ti_current) - error("Particle si not drifted to current time"); - if (sj->ti_drift != e->ti_current) - error("Particle sj not drifted to current time"); + if (si->ti_drift != e->ti_current) { + scheduler_write_cell_dependencies_debug(&r->e->sched, e->verbose, e->step, ci); + scheduler_write_cell_dependencies_debug(&r->e->sched, e->verbose, e->step, cj); + error("Particle si not drifted to current time. si->id = %lld, sj->id = %lld | i: hydro super = %lld, grav super = %lld," + " top = %lld, c = %lld | j hydro super = %lld, grav super = %lld," + " top = %lld, c = %lld", + si->id, sj->id, ci->hydro.super->cellID, ci->grav.super->cellID, ci->top->cellID, + ci->cellID, cj->hydro.super->cellID, cj->grav.super->cellID, cj->top->cellID, + cj->cellID); + } + if (sj->ti_drift != e->ti_current) { + scheduler_write_cell_dependencies_debug(&r->e->sched, e->verbose, e->step, ci); + scheduler_write_cell_dependencies_debug(&r->e->sched, e->verbose, e->step, cj); + error("Particle sj not drifted to current time. si->id = %lld, sj->id = %lld | i: hydro super = %lld, grav super = %lld," + " top = %lld, c = %lld | j hydro super = %lld, grav super = %lld," + " top = %lld, c = %lld", + si->id, sj->id, ci->hydro.super->cellID, ci->grav.super->cellID, ci->top->cellID, + ci->cellID, cj->hydro.super->cellID, cj->grav.super->cellID, cj->top->cellID, + cj->cellID); + } #endif if (r2 < hig2 || r2 < hjg2) {