Commit 99219705 authored by David Deepwell's avatar David Deepwell
Browse files

Remove double execution of analysis function

A second calling of analysis was previously added to run on the
initial configuration, but this was repeated for each subsequent
time-step. This repetition has been fixed.
parent ecfdaf83
......@@ -70,9 +70,6 @@ namespace NSIntegrator {
double desttime = times[0] + the_timestep;
// Analyze the initial conditions
usercode->analysis(times[0],us[0],vs[0],ws[0],tracers_now,pressure);
/* Now, take timesteps until we reach our one-timestep destination*/
bool starting_step = false;
......@@ -87,6 +84,9 @@ namespace NSIntegrator {
total time integration error. So, if we're starting up,
we take a small timestep and progressively double it. */
if (times[-1] == times[-2]) {
// Analyze the initial conditions
usercode->analysis(times[0],us[0],vs[0],ws[0],tracers_now,pressure);
/* Startup, so take a tiny timestep */
my_step = the_timestep/64;
starting_step = true;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment