skiehl / lcsim Goto Github PK

Parallelize the simulation of multiple light curves and/or the Emmanoulopoulos algorithm.

The multi_logn method implements a custom distribution. It is not an integral part of the LightCurveSimulator class and should not be part of it. Better move it to a separate module.

Currently, power at lower frequencies can be included by setting nlcs>1 in the sim_tk() or sim_emp() methods. With e.g. nlcs=10, a single long light curve is simulated ten times as long as the requested total time. This is then split into 10 independent light curves, that have the requested total time, but include power at frequencies ten times lower than the one corresponding to the total time of the light curves. With this setup it is not possible to produce more than 10 light curves with a single method call, if this is the low frequency range that is wished for. Either even lower frequencies need to be included by increasing nlcs, which increases the run time as the long light curve grows in length, or multiple method calls are required.
The low frequency coverage should be separated from the number of light curves.

Add verbosity option to the main methods to get some feedback about the progress.

Add a method to check the PDF of simulated light curves against a target PDF and accept/reject simulations based on a matching criterion.