19 double PumpDamp, PumpTau;
21 double operator() (
double x)
25 double opfun_v =
sexp(PumpTau*v)*v;
96 double t1 = 3.02*pow(t.
Emis().
damp()*1e3,-0.064 );
106 value = core + wings;
112 value =
MIN2(1., value );
119 else if(
rt.nLineContShield == 0 )
142 return (0.98925439 + 0.084594094*t)/(1. + t*(0.64794212 + t*0.44743976));
178 yinc1 = opfun.
sum( 0.,
BREAK, func );
179 yinc2 = opfun.
sum(
BREAK, 100., func );
181 a0 = 0.886227*(1. + func.PumpDamp);
182 conpmp_v = (yinc1 + yinc2)/
a0;
sys_float sexp(sys_float x)
NORETURN void TotalInsanity(void)
sys_float SDIV(sys_float x)
#define DEBUG_ENTRY(funcname)
double sum(double min, double max, Integrand func)
EmissionList::reference Emis() const
double esc_PRD_1side(double tau, double a)
#define LINE_CONT_SHIELD_FERLAND
#define LINE_CONT_SHIELD_PESC
double esca0k2(double taume)
#define LINE_CONT_SHIELD_FEDERMAN
double esc_CRDwing_1side(double tau, double a)
STATIC double FITTED(double t)
double RT_continuum_shield_fcn(const TransitionProxy &t)
STATIC double conpmp(const TransitionProxy &t)
static const double BREAK
void VoigtH(realnum a, const realnum v[], realnum y[], int n)