Fix a bug: was only computing when there was a mutation AND we were at a...
Fix a bug: was only computing when there was a mutation AND we were at a periodic time step ( time()%period) Now computes for each required periodic timestep
Tested on my data only, not extensively (50 runs in 6 conditions, 2mio generations)