GitLab

Update the nebular models from Akio Inoue. Now log U is sampled in steps of 0.1 dex down to -4. Also deviations from the 10000K case B assumption are taken into account.

Limit the number of threads to 1 per process in order to prevent MKL from creating numerous threads leading to a strong oversubscription of the CPU/RAM as pcigale already does its own parallelisation.

Implement the computation of the mass-weighted age for the bc03 and m2005 module. This is done in the corresponding convolve() functions. The information is returned in a new variable called info_all, which also contains all the information when the entire SFH is included (rather than separating old and young populations).

Rather than calling convolve() twice, call it just once with the separation age as an argument et let it do the various different convolutions. While this is functionally equivalent, this allows to pass the full SFH to convolve(), which in turn is necessary to compute additional quantities. This patch is a preparatory work for the the implementation of mass-weighted ages.

Do not return the break strengths in the bc03 module as they are not computed corrected by a convolvution.

There is no need to return the wavelength grid after convolution in the data bc03 module, it is readily available in the bc03 sed module.

Homogenise the bc03 and m2005 data modules so they use the same variable names as far as possible. Comments have also been corrected and made consistent.

Version 0.10.0

In order to release a new version correcting some bugs, we remove the
work in progress correction for upper limits.

Correct the function adjusting the data and implements a new way to indicate
upper limits.

FIXME: This breaks other parts of the code that must be corrected because now
upper limits are indicated with a negative flux.

The test fails so that we can correct the function.

Rename (and correct a bug into) the SFH quenching module to sfh_quenching_trunk
and add a sfh_quenching_smooth version where the SFR is linearly decreasing
after the quenching.

The WFCAM Ks filter does not exist.
The filters J, and Y were wrong.
Add H, K, and Z filters.

Use a more human-friendly format for the ASCII output of savefluxes to facilitate visual inspection as a FITS output is now available for machine reading.

When on a grid, make sure that the models are always computed in the same order to facilitate comparisons between different runs.

Remove the output_file and output_format parameters from the savefluxes module. Now results are systematically saved as ASCII and FITS tables.

Clarify that if the parameters are given through a file, then the pdf_analysis module should not be used.

Correct the description of the radio module. It stated erroneously that it handled thermal emission whereas this is done by the nebular module.

As we loop over the bands to compute the χ², it is smarter that the bands are the slower varying index and the models the faster varying index. That way it is more likely that models at the correct redshift are prefetched into the cache.

Use a multiplication rather than a division to convert luminosities into fluxes. This should yield a small speedup.

When giving the list of parameters through a file, cigale did not compute properly what parameter changed between to successive models.

With the new sanity check of the input parameters, cigale did not handle the fact that the redshift could be given in the parameters file. Now this is handled properly.

A significant fraction of the total run time is spent computing integrals (e.g. fluxes in passbands). We can make the integration faster by rewriting the trapezoidal rule in terms of np.dot(). This allows to offload the computation to optimised libraries. The end result is that the integration is twice as fast, with a gain of ~10-15% on the total run time.

0.9.0 release