Skip to content

Conversation

@claudiodsf
Copy link
Member

This PR introduces support for user-specified attenuation models, allowing either a fixed value or a frequency-dependent Q model.

By defining Q_model in the configuration file, the t_star parameter is no longer inverted.

See discussion in #73 for more details and context.

Remaining tasks:

  • Support for grid search (GS) and importance sampling algorithms (IS)
  • Documentation

…tor during inversion.

It only works if `inv_algorithm` is set to `TNC`, `LM`, `BH`.

Spectral figures are broken when using `Q_model`:
- modelled spectra are not plotted without t_star values.
- stacked spectra are not plotted since the summary spectrum cannot
  be computed without t_star values.

TODO:
- Add support for `Q_model` with `GS` and `IS` algorithms.
- Fix figures when using `Q_model`.
Also, better handling of logarithmic averages.
This is used to correctly compute the spectral integral for radiated energy,
when a frequency-dependent t_star model is used.

It will probably also be used later for plotting functions.

t_star_model is removed from Spectrum.stats before saving to HDF5,
as HDF5 does not support saving functions.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants