bug_end_time.c

Time condition
Iteration condition

This small example demonstrates the inability to change the date of the final event of the simulation during execution. This happens with both event end (t = end_time) (“time” condition) and event end (i = end_iteration) (“iteration” condition) types of criteria.

In both cases, the simulation ends when reaching the first definition of the final date (end_time or end_iteration), regardless if we increased or decreased it.

#include "run.h"

Time condition

#if 1

Let end_time be the initial date of the final event, which ends the run.

double end_time = 10.;


int main ()
{  
  run();
}


event init (i = 0)
{
  printf ("Before init, end_time is %g\n", end_time);

For some reason, let’s decide to change the final time. This can be done for instance if a convergence criterion hasn’t been reached and the run needs to last a little longer.

  end_time += 1.;
    
  printf ("After init, end_time is %g\n", end_time);
}


event end (t = end_time) {
  printf ("End event reached at t = %g, whereas end_time is %g\n", t, end_time);
}

The results we get in out are as follows :
Before init, end_time is 10
After init, end_time is 11
End event reached at t = 10, whereas end_time is 11
# Quadtree, 1 steps, 4.6e-05 CPU, 3.1e-05 real, 1.32e+08 points.step/s, 0 var

Iteration condition

#else

Now let’s try the same but this time the run ends when the number of iterations i reaches end_iteration.

int end_iteration = 10;

int main ()
{  
  run();
}

event init (i = 0) {
  printf ("Before init, end_iteration is %i\n", end_iteration);
  end_iteration -= 1;
  printf ("After init, end_iteration is %i\n", end_iteration);
}

event end (i = end_iteration) {
  printf ("End event reached at i = %i, whereas end_iteration is %i\n", i, end_iteration);
}

#endif

We get similar results :
Before init, end_iteration is 10
After init, end_iteration is 9
End event reached at i = 10, whereas end_iteration is 9
# Quadtree, 10 steps, 5.3e-05 CPU, 3.8e-05 real, 1.08e+09 points.step/s, 0 var