diff --git a/src/debug.c b/src/debug.c index 8675355ca3a97eaa29d79a7bc0e064bb577168d3..f63ed557d908607e6cd53edf8467675a938245d7 100644 --- a/src/debug.c +++ b/src/debug.c @@ -393,7 +393,7 @@ void dumpCells(const char *prefix, int active, int mpiactive, int pactive, fclose(file); } -#ifdef HAVE_METIS +#if defined(WITH_MPI) && defined(HAVE_METIS) /** * @brief Dump the METIS graph in standard format, simple format and weights diff --git a/src/debug.h b/src/debug.h index 2cff37fa6f6f03185dc769bb770fe3a98a424ce1..5a646948204cc809c66430bb01c9ee90e21d720a 100644 --- a/src/debug.h +++ b/src/debug.h @@ -39,7 +39,7 @@ int checkCellhdxmax(const struct cell *c, int *depth); void dumpCells(const char *prefix, int active, int mpiactive, int pactive, struct space *s, int rank, int step); -#ifdef HAVE_METIS +#if defined(WITH_MPI) && defined(HAVE_METIS) #include "metis.h" void dumpMETISGraph(const char *prefix, idx_t nvtxs, idx_t ncon, idx_t *xadj, idx_t *adjncy, idx_t *vwgt, idx_t *vsize, idx_t *adjwgt); diff --git a/src/runner.c b/src/runner.c index 27ce272163bfe0335b334211cc30299d3b33699f..c08c425a7bed76d0cdc77593079af0d184b7d265 100644 --- a/src/runner.c +++ b/src/runner.c @@ -1001,6 +1001,21 @@ void runner_do_kick1(struct runner *r, struct cell *c, int timer) { /* do the kick */ kick_part(p, xp, ti_begin, ti_begin + ti_step / 2, timeBase); + + /* Update the accelerations to be used in the drift for hydro */ + if (p->gpart != NULL) { + + const float a_tot[3] = {p->a_hydro[0] + p->gpart->a_grav[0], + p->a_hydro[1] + p->gpart->a_grav[1], + p->a_hydro[2] + p->gpart->a_grav[2]}; + + p->a_hydro[0] = a_tot[0]; + p->a_hydro[1] = a_tot[1]; + p->a_hydro[2] = a_tot[2]; + p->gpart->a_grav[0] = a_tot[0]; + p->gpart->a_grav[1] = a_tot[1]; + p->gpart->a_grav[2] = a_tot[2]; + } } }