From 1f006bf77cf57ec0f52ab3a747e71a8f3b268907 Mon Sep 17 00:00:00 2001
From: Stefan Arridge <stefan.arridge@durham.ac.uk>
Date: Tue, 1 Nov 2016 17:36:03 +0000
Subject: [PATCH] In stats_collect_part_mapper in statistics.c, we don't
collect gpe if gpart does not exist
---
src/statistics.c | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/statistics.c b/src/statistics.c
index f1d8f61385..847bf24cb7 100644
--- a/src/statistics.c
+++ b/src/statistics.c
@@ -155,7 +155,8 @@ void stats_collect_part_mapper(void *map_data, int nr_parts, void *extra_data) {
/* Collect energies. */
stats.E_kin += 0.5f * m * (v[0] * v[0] + v[1] * v[1] + v[2] * v[2]);
stats.E_pot_self += 0.f;
- stats.E_pot_ext +=
+ if (gp != NULL)
+ stats.E_pot_ext +=
m * external_gravity_get_potential_energy(potential, phys_const, gp);
stats.E_int += m * hydro_get_internal_energy(p, dt);
stats.E_rad += cooling_get_radiated_energy(xp);
--
GitLab