The standardisation of error bars must be done on the final estimates, not on the estimates for a single block.

Use a much more readable formula to compute the uncertainties. We use the fact that Var(X) = E(Var(X)) + Var(E(X)). At the same time, broadcast the weights array to the same shape as the other arrays in order to simplify the syntax for the computation.

Do not correct intensive physical properties for the difference in redshift between the observation and the grid.

Read the list of sed modules from pcigale.ini. This is more flexible than inferring this list from the input parameters file as some modules can be missing (for instance when using the default values for all the parameters of a given module).

Compute the best fit at the at the exact distance of the analysed object and not at the distance of the rounded redshift. This eliminates the need for a correction factor.

Give the redshift of the model rather than the list of redshifts to correct for the redshift difference between the model and the observation.

Add a results manager to store the results of the analysis (bayesian and best fit estimates, best chi², etc.). Several such managers can be merged to combine appropriately the analysis of the observations in different blocks. This is the base class. The implementation will come in subsequent patches.

More correct block size computation to ensure that all the models are computed and that each block contains all the redshifts for a given model.

Save the model id rather than the index within the block as this would give an incorrect id when using more than one block.

Switch savefluxes to using a models manager and make all the necessary adjustments with the workers.

Add a models manager that stores the fluxes, the physical properties, and all the necessary information to compute the models. Actual use will come in forthcoming patches.

Add the possibility to store the list of blocks of models to compute to the parameters manager. If the number of blocks is not given, we assume it must be computed in 1 block.

Implement an observations manager. For now it is mainly used as a simple container for the observations table but upcoming rearchitecturing will make full use of it.

Remove the TOLERANCE constant and replace it with 0. There is no good reason to have a tolerance and its value is quite arbitrary anyway.

Do not crash when trying to compute equivalent widths when the nebular component emission is missing.

Rename ParametersHandler to ParametersManager to be consistent with new managers that will come in the near future.

The output directory is a bit of a mess, with it being defined in different places. As we never make use of the possibility of changing it, let's simply hardcode it to out/. It makes filenames slightly easier to read too.

Move the backup_dir() function to the base class of the analysis modules as it makes sense for analysers to prepare the output directories.

Simplify the cache architecture as we always use the in-memory one. The removes one level of indirection and makes the model generation faster by 1%.

Implement the possibility of not caching a module in order to save memory. A typical use case is when having many different stellar populations at different metallicities. It can take a substantial amount of memory, especially when using many parallel processes.

Correct a bug that was preventing the module restframe_parameters from being listed among the available modules.

Tiny optimisation replacing a division with a subtraction. Subtractions are quite a bit faster than divisions.