Create a templated experiment

[richelbilderbeek]

To ensure that we all do our experiment in the same way, a template function should be written, ensuring we do things in the same way.

This must be done by the team as a whole. I suggest I write a first draft, then let the team modify.

[richelbilderbeek]

I'll put this in the jkr namespace.

[richelbilderbeek]

Plan:

• [x] Monday: create simplest version
• [ ] Tuesday: teach students how to adopt this. Decide on seperated generations y/n. Implement this
• [ ] Wednesday: let students adopt this. Decide on offspring creation. Implement this
• [ ] Thursday: let students adopt this. Decide on viability. Implement this
• [ ] Friday: let students adopt this
• [ ] Monday: Decide on results. Implement this
• [ ] Tuesday: let students adopt this

[richelbilderbeek]

I suggest from

template <class parameters, class simulation, class results>
void do_experiment(const parameters& p)
{
  simulation s = create_simulation(p);
  run(s);
  save_ltt_plot(get_results(s), get_ltt_plot_filename(p));
}

to

template <class parameters, class simulation, class results>
void do_experiment(const parameters& p)
{
  simulation s = create_simulation(p);
  const int n_generations = get_n_generations(p);
  for (int t=0; t!=n_generations; ++t)
  {
    const auto next_population = create_next_population(s);
    set_population(s, next_population);
  }
  save_ltt_plot(get_results(s), get_ltt_plot_filename(p));
}

Sure, the get/set can be done in one line, but, thanks to auto, the rigidity of the function would not change from that.

Agree on this?

• [x] @JorikdeBoer
• [x] @Lumphie
• [x] @richelbilderbeek

[richelbilderbeek]

Suggestion for next iteration:

• enforce create_next_population only reads from simulation
• extract the RNG engine

template <class parameters, class simulation, class results>
void do_experiment(const parameters& p)
{
  simulation s = create_simulation(p);
  const int n_generations = get_n_generations(p);
  std::mt19937 rng_engine(get_rng_seed(p));
  for (int t=0; t!=n_generations; ++t)
  {
    const auto next_population
      = create_next_population(
        const_cast<const simulation&>(s), //Enforce only reading the simulation
        rng_engine
      );
    set_population(s, next_population);
  }
  save_ltt_plot(get_results(s), get_ltt_plot_filename(p));
}

Agree on this?

• [x] @JorikdeBoer
• [ ] @Lumphie

[richelbilderbeek]

This goes just fine, closing this Issue for now.