macroatom.h

#ifndef MACROATOM_H
#define MACROATOM_H

#include "globals.h"
#include "packet.h"

void macroatom_open_file(int my_rank);
void macroatom_close_file();

void do_macroatom(Packet &pkt, const MacroAtomState &pktmastate);

#pragma omp declare simd
[[nodiscard]] auto rad_deexcitation_ratecoeff(int nonemptymgi, int element, int ion, int lower, double epsilon_trans,
                                              float A_ul, double upperstatweight, double nnlevelupper, double t_current)
    -> double;

#pragma omp declare simd
[[nodiscard]] auto rad_excitation_ratecoeff(int nonemptymgi, int element, int ion, int lower, int uptransindex,
                                            double epsilon_trans, double nnlevel_lower, int lineindex, double t_current)
    -> double;

#pragma omp declare simd
[[nodiscard]] auto rad_recombination_ratecoeff(float T_e, float nne, int element, int upperion, int upperionlevel,
                                               int lowerionlevel, int nonemptymgi) -> double;
#pragma omp declare simd
[[nodiscard]] auto stim_recombination_ratecoeff(float nne, int element, int upperion, int upper, int lower,
                                                int nonemptymgi) -> double;

#pragma omp declare simd
[[nodiscard]] auto col_recombination_ratecoeff(float T_e, float nne, int element, int upperion, int upper, int lower,
                                               double epsilon_trans) -> double;
#pragma omp declare simd
[[nodiscard]] auto col_ionization_ratecoeff(float T_e, float nne, int element, int ion, int lower, int phixstargetindex,
                                            double epsilon_trans) -> double;
#pragma omp declare simd
[[nodiscard]] auto col_deexcitation_ratecoeff(float T_e, float nne, double epsilon_trans, int element, int ion,
                                              int upper, const LevelTransition &downtransition) -> double;
#pragma omp declare simd
[[nodiscard]] auto col_excitation_ratecoeff(float T_e, float nne, int element, int ion, int lower, int uptransindex,
                                            double epsilon_trans, double lowerstatweight) -> double;

#endif  // MACROATOM_H