From fcdafd94a6f044fdff3610ea0443ab2c1b38a35a Mon Sep 17 00:00:00 2001 From: Yohai Meiron Date: Fri, 14 Jan 2022 17:54:38 -0500 Subject: [PATCH] Fixed PN force adjustment before PN parameters set --- phigrape.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/phigrape.cpp b/phigrape.cpp index d3a734e..4dd4865 100644 --- a/phigrape.cpp +++ b/phigrape.cpp @@ -449,11 +449,11 @@ int main(int argc, char *argv[]) black_hole_physics.set_softening(eps, config.live_smbh_custom_eps); black_hole_physics.adjust_softening(pot[0], pot[1], a[0], a[1], adot[0], adot[1]); } - if (config.binary_smbh_pn) black_hole_physics.adjust_post_newtonian(dt_min, a[0], a[1], adot[0], adot[1]); } if (config.binary_smbh_pn) { black_hole_physics.set_post_newtonian(config.pn_c, config.pn_usage.data()); if (config.pn_usage[6]) black_hole_physics.set_spins(config.smbh1_spin.data(), config.smbh2_spin.data()); + black_hole_physics.adjust_post_newtonian(dt_min, a[0], a[1], adot[0], adot[1]); } std::vector pot_ext(N, 0.);