From 944ffdf44a232a9b1aa087755455f211443b605f Mon Sep 17 00:00:00 2001
From: Robert Izzard <r.izzard@surrey.ac.uk>
Date: Sat, 7 Sep 2019 14:18:39 +0100
Subject: [PATCH] add a failure when discs are run with a solver other than
forward euler
---
src/disc/evolve_disc.c | 8 +++++++-
1 file changed, 7 insertions(+), 1 deletion(-)
diff --git a/src/disc/evolve_disc.c b/src/disc/evolve_disc.c
index a491703dd..305879a9f 100644
--- a/src/disc/evolve_disc.c
+++ b/src/disc/evolve_disc.c
@@ -25,12 +25,18 @@ void evolve_disc(struct stardata_t * const stardata,
Discdebug(1,"EV t=%g dt=%g\n",
disc->lifetime/YEAR_LENGTH_IN_SECONDS,
dtsecs/YEAR_LENGTH_IN_SECONDS);
+
+ if(stardata->preferences->solver != SOLVER_FORWARD_EULER)
+ {
+ Exit_binary_c(BINARY_C_ALGORITHM_OUT_OF_RANGE,
+ "The discs algorithm has not been tested with solvers other than Forward-Euler. Please only use this solver.");
+ }
disc->first = disc->first==TRUE ? TRUE : Boolean_(Is_zero(disc->lifetime));
disc->dT = 0.0;
Clear_disc_feedback(disc);
-
+
if(Disc_is_disc(disc))
{
/*
--
GitLab