Raise exception if infinite_sites is called during a simulation. (#820)

Closes #769
molpopgen · Oct 15, 2021 · 4011bed · 4011bed
1 parent 95a3601
commit 4011bed
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 0 deletions.
diff --git a/doc/misc/changelog.md b/doc/misc/changelog.md
@@ -33,6 +33,9 @@ Behavior changes
   This change makes the timings consistent with all other events and also makes
   certain operations easier/feasible.
   {pr}`809`
+* Calling {func}`fwdpy11.infinite_sites` during a simulation now raises `RuntimeError`.
+  {pr}`820`
+  {issue}`769`
 
 New features
 

diff --git a/fwdpy11/src/ts/infinite_sites.cc b/fwdpy11/src/ts/infinite_sites.cc
@@ -39,6 +39,11 @@ init_infinite_sites(py::module& m)
         "infinite_sites",
         [](const fwdpy11::GSLrng_t& rng, fwdpy11::Population& pop,
            const double mu) -> unsigned {
+            if (pop.is_simulating)
+                {
+                    throw std::runtime_error(
+                        "infinite_sites cannot be called during a simulation");
+                }
             if (mu <= 0.0)
                 {
                     return 0u;