#include "cddefines.h"
#include "rfield.h"
#include "phycon.h"
#include "dense.h"
#include "hmi.h"
#include "thermal.h"
#include "trace.h"
#include "iterations.h"
#include "vectorize.h"
#include "grainvar.h"
#include "grains.h"

Include dependency graph for grains_qheat.cpp:

Enumerations
enum	QH_Code { QH_OK , QH_ANALYTIC , QH_ANALYTIC_RELAX , QH_DELTA , QH_NEGRATE_FAIL , QH_LOOP_FAIL , QH_ARRAY_FAIL , QH_THIGH_FAIL , QH_RETRY , QH_CONV_FAIL , QH_BOUND_FAIL , QH_DELTA_FAIL , QH_NO_REBIN , QH_LOW_FAIL , QH_HIGH_FAIL , QH_STEP_FAIL , QH_FATAL , QH_WIDE_FAIL , QH_NBIN_FAIL , QH_REBIN_FAIL }

Functions
double	no_atoms (size_t nd)

STATIC void	qheat_init (size_t, vector< double > &, double *)

STATIC void	GetProbDistr_LowLimit (size_t, double, double, double, const vector< double > &, const vector< double > &, vector< double > &, vector< double > &, vector< double > &, long , double , long , QH_Code )

STATIC double	TryDoubleStep (vector< double > &, vector< double > &, vector< double > &, vector< double > &, vector< double > &, const vector< double > &, const vector< double > &, double, double , double, long, size_t, bool )

STATIC double	log_integral (double, double, double, double, double, double, double, double)

STATIC void	ScanProbDistr (const vector< double > &, const vector< double > &, long, double, long, double, long , long , long , long , QH_Code *)

STATIC long	RebinQHeatResults (size_t, long, long, vector< double > &, vector< double > &, vector< double > &, vector< double > &, vector< double > &, vector< double > &, vector< double > &, QH_Code *)

STATIC void	GetProbDistr_HighLimit (long, double, double, double, vector< double > &, vector< double > &, vector< double > &, double , long , double , QH_Code )

STATIC double	uderiv (double, size_t)

STATIC double	ufunct (double, size_t, bool *)

STATIC double	inv_ufunct (double, size_t, bool *)

STATIC double	DebyeDeriv (double, long)

realnum	fast_expm1 (realnum x)

STATIC long	GrainMakeDiffuseSingle (double Tgrain, double fracpop, avx_ptr< realnum > &flux, long nflux)

void	GrainMakeDiffuse ()

void	qheat (vector< double > &qtemp, vector< double > &qprob, long int *qnbin, size_t nd)

STATIC void	GetProbDistr_LowLimit (size_t nd, double rel_tol, double Umax, double fwhm, const vector< double > &Phi, const vector< double > &PhiDrv, vector< double > &qtemp, vector< double > &qprob, vector< double > &dPdlnT, long int qnbin, double new_tmin, long nWideFail, QH_Code ErrorCode)

void	InitEnthalpy ()

Variables
static const long	NQMIN = 10L

static const double	PROB_CUTOFF_LO = 1.e-15

static const double	PROB_CUTOFF_HI = 1.e-20

static const double	SAFETY = 1.e+8

static const long	NSTARTUP = 5L

static const double	MAX_EVENTS = 150.

static const long	LOOP_MAX = 20L

static const double	DEF_FAC = 3.

static const double	PROB_TOL = 0.02

static const long	NQTEST = 500L

static const double	FWHM_RATIO = 0.1

static const double	FWHM_RATIO2 = 0.007

static const long	MAX_LOOP = 2*NQGRID

static const double	QHEAT_TOL = 5.e-3

static const long	WIDE_FAIL_MAX = 3

static const double	STRICT = 1.

static const double	RELAX = 15.

static const double	QT_RATIO = 1.074

static const double	DEN_SIL = 3.30

static const double	MW_SIL = 24.6051

static const double	tlim [5] ={0.,50.,150.,500.,DBL_MAX}

static const double	ppower [4] ={2.00,1.30,0.68,0.00}

static const double	cval [4]

static const realnum	LIM1 = 2.e-4f

static const realnum	LIM2 = sqrtf(6.e-4f)

static const realnum	LIM3 = cbrtf(24.e-4f)

Enumeration Type Documentation

◆ QH_Code

enum QH_Code

Enumerator
QH_OK
QH_ANALYTIC
QH_ANALYTIC_RELAX
QH_DELTA
QH_NEGRATE_FAIL
QH_LOOP_FAIL
QH_ARRAY_FAIL
QH_THIGH_FAIL
QH_RETRY
QH_CONV_FAIL
QH_BOUND_FAIL
QH_DELTA_FAIL
QH_NO_REBIN
QH_LOW_FAIL
QH_HIGH_FAIL
QH_STEP_FAIL
QH_FATAL
QH_WIDE_FAIL
QH_NBIN_FAIL
QH_REBIN_FAIL

Function Documentation

◆ DebyeDeriv()

STATIC double DebyeDeriv	(	double	x,
		long	n )

References ASSERT, DEBUG_ENTRY, gauss_init(), gauss_legendre(), MAX2, pow2, pow3, powi(), STATIC, and TotalInsanity().

Referenced by uderiv().

Here is the call graph for this function:

◆ fast_expm1()

realnum fast_expm1 ( realnum x )

inline

References LIM1, LIM2, and LIM3.

Referenced by GrainMakeDiffuseSingle().

◆ GetProbDistr_HighLimit()

STATIC void GetProbDistr_HighLimit	(	long	nd,
		double	TolFac,
		double	Umax,
		double	fwhm,
		vector< double > &	qtemp,
		vector< double > &	qprob,
		vector< double > &	dPdlnT,
		double *	tol,
		long *	qnbin,
		double *	new_tmin,
		QH_Code *	ErrorCode )

References DEBUG_ENTRY, fprintf(), gv, inv_ufunct(), ioQQQ, log_integral(), MAX2, MIN2, NDEMS, NQGRID, NQMIN, nzone, pow(), pow2, PROB_CUTOFF_HI, PROB_CUTOFF_LO, PROB_TOL, QH_ANALYTIC, QH_ANALYTIC_RELAX, QH_BOUND_FAIL, QH_CONV_FAIL, QH_THIGH_FAIL, QH_WIDE_FAIL, QT_RATIO, SAFETY, splint_safe(), STATIC, STRICT, trace, uderiv(), ufunct(), and vlog().

Referenced by qheat().

Here is the call graph for this function:

◆ GetProbDistr_LowLimit() [1/2]

STATIC void GetProbDistr_LowLimit	(	size_t	nd,
		double	rel_tol,
		double	Umax,
		double	fwhm,
		const vector< double > &	Phi,
		const vector< double > &	PhiDrv,
		vector< double > &	qtemp,
		vector< double > &	qprob,
		vector< double > &	dPdlnT,
		long int *	qnbin,
		double *	new_tmin,
		long *	nWideFail,
		QH_Code *	ErrorCode )

References ASSERT, DEBUG_ENTRY, DEF_FAC, fprintf(), GRAIN_TMIN, gv, ioQQQ, max(), MAX2, MAX3, MAX_LOOP, MIN2, NDEMS, NQGRID, NQMIN, NQTEST, NSTARTUP, nzone, pow3, PROB_CUTOFF_HI, PROB_CUTOFF_LO, PROB_TOL, QH_ARRAY_FAIL, QH_BOUND_FAIL, QH_CONV_FAIL, QH_DELTA_FAIL, QH_HIGH_FAIL, QH_LOOP_FAIL, QH_LOW_FAIL, QH_NBIN_FAIL, QH_NO_REBIN, QH_STEP_FAIL, QH_THIGH_FAIL, QH_WIDE_FAIL, QHEAT_TOL, RebinQHeatResults(), SAFETY, ScanProbDistr(), spldrv_safe(), splint_safe(), STATIC, trace, TryDoubleStep(), uderiv(), ufunct(), and WIDE_FAIL_MAX.

Here is the call graph for this function:

◆ GetProbDistr_LowLimit() [2/2]

STATIC void GetProbDistr_LowLimit	(	size_t	,
		double	,
		double	,
		double	,
		const vector< double > &	,
		const vector< double > &	,
		vector< double > &	,
		vector< double > &	,
		vector< double > &	,
		long *	,
		double *	,
		long *	,
		QH_Code *	)

References STATIC.

Referenced by qheat().

◆ GrainMakeDiffuse()

void GrainMakeDiffuse ( )

main routine for generating the grain diffuse emission

References ASSERT, CONSERV_TOL, avx_ptr< T, lgBC >::data(), DEBUG_ENTRY, dense, get_ptr(), GrainMakeDiffuseSingle(), gv, hmi, ipHYDROGEN, max(), MAX2, NQGRID, pow2, qheat(), rfield, STRG_CAR, STRG_SIL, thermal, and TotalInsanity().

Referenced by RT_diffuse().

Here is the call graph for this function:

◆ GrainMakeDiffuseSingle()

STATIC long GrainMakeDiffuseSingle	(	double	Tgrain,
		double	fracpop,
		avx_ptr< realnum > &	flux,
		long	nflux )

References fast_expm1(), avx_ptr< T, lgBC >::ptr0(), rfield, STATIC, and vexpm1().

Referenced by GrainMakeDiffuse().

Here is the call graph for this function:

◆ InitEnthalpy()

void InitEnthalpy ( )

initialize interpolation arrays for grain enthalpy

References DEBUG_ENTRY, GRAIN_TMIN, gv, log_integral(), NDEMS, spline(), tlim, uderiv(), and vlog().

Referenced by GrainsInit().

Here is the call graph for this function:

◆ inv_ufunct()

STATIC double inv_ufunct	(	double	enthalpy,
		size_t	nd,
		bool *	lgBoundErr )

References ASSERT, cdEXIT, DEBUG_ENTRY, EXIT_FAILURE, fprintf(), gv, ioQQQ, NDEMS, splint_safe(), and STATIC.

Referenced by GetProbDistr_HighLimit(), qheat(), and TryDoubleStep().

Here is the call graph for this function:

◆ log_integral()

STATIC double log_integral	(	double	xx1,
		double	yy1,
		double	xx2,
		double	yy2,
		double	log_xx1,
		double	log_yy1,
		double	log_xx2,
		double	log_yy2 )

References DEBUG_ENTRY, and STATIC.

Referenced by GetProbDistr_HighLimit(), InitEnthalpy(), RebinQHeatResults(), and TryDoubleStep().

◆ no_atoms()

double no_atoms ( size_t nd )

inline

References gv.

Referenced by mie_read_form(), and uderiv().

◆ qheat()

void qheat	(	vector< double > &	qtemp,
		vector< double > &	qprob,
		long int *	qnbin,
		size_t	nd )

References ASSERT, CONSERV_TOL, DEBUG_ENTRY, DEF_FAC, fprintf(), FWHM_RATIO, FWHM_RATIO2, GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), GRAIN_TMIN, gv, inv_ufunct(), ioQQQ, iterations, LOOP_MAX, MAX2, MAX_EVENTS, MIN2, NDEMS, NQGRID, NULL, nzone, open_data(), pow2, PROB_CUTOFF_LO, QH_ANALYTIC, QH_CONV_FAIL, QH_DELTA, QH_DELTA_FAIL, QH_FATAL, QH_HIGH_FAIL, QH_LOW_FAIL, QH_NEGRATE_FAIL, QH_NO_REBIN, QH_OK, QH_RETRY, QH_WIDE_FAIL, qheat_init(), RELAX, rfield, spldrv_safe(), STRICT, TorF(), trace, uderiv(), ufunct(), and WIDE_FAIL_MAX.

Referenced by GrainMakeDiffuse().

Here is the call graph for this function:

◆ qheat_init()

STATIC void qheat_init	(	size_t	nd,
		vector< double > &	phiTilde,
		double *	check )

References ASSERT, CONSERV_TOL, DEBUG_ENTRY, dprintf(), ChargeBin::DustZ, ChargeBin::ehat, ChargeBin::fac1, ChargeBin::FracPop, gv, ChargeBin::HeatingRate2, ioQQQ, ChargeBin::ipThresInf, ChargeBin::ipThresInfVal, LIM2, LIM3, max(), MAX2, min(), MIN2, ChargeBin::nfill, phycon, ChargeBin::PotSurf, ChargeBin::PotSurfInc, powpq(), rfield, sign(), STATIC, ChargeBin::ThresInfInc, ChargeBin::ThresSurf, ChargeBin::ThresSurfVal, ChargeBin::yhat, and ChargeBin::yhat_primary.

Referenced by qheat().

Here is the call graph for this function:

◆ RebinQHeatResults()

STATIC long RebinQHeatResults	(	size_t	nd,
		long	nstart,
		long	nend,
		vector< double > &	p,
		vector< double > &	qtemp,
		vector< double > &	qprob,
		vector< double > &	dPdlnT,
		vector< double > &	u1,
		vector< double > &	delu,
		vector< double > &	Lambda,
		QH_Code *	ErrorCode )

References ASSERT, CONSERV_TOL, DEBUG_ENTRY, fprintf(), frac(), gv, ioQQQ, log_integral(), MAX2, MIN2, NDEMS, NQGRID, NQMIN, pow(), PROB_CUTOFF_LO, QH_CONV_FAIL, QH_REBIN_FAIL, QT_RATIO, SAFETY, splint_safe(), STATIC, trace, uderiv(), ufunct(), and vlog().

Referenced by GetProbDistr_LowLimit().

Here is the call graph for this function:

◆ ScanProbDistr()

STATIC void ScanProbDistr	(	const vector< double > &	u1,
		const vector< double > &	dPdlnT,
		long	nbin,
		double	maxVal,
		long	nmax,
		double	qtmin1,
		long *	nstart,
		long *	nstart2,
		long *	nend,
		long *	nWideFail,
		QH_Code *	ErrorCode )

References ASSERT, DEBUG_ENTRY, fprintf(), GRAIN_TMIN, ioQQQ, MAX2, NQMIN, QH_HIGH_FAIL, QH_LOW_FAIL, QH_NBIN_FAIL, QH_NO_REBIN, QH_WIDE_FAIL, STATIC, trace, and WIDE_FAIL_MAX.

Referenced by GetProbDistr_LowLimit().

Here is the call graph for this function:

◆ TryDoubleStep()

STATIC double TryDoubleStep	(	vector< double > &	u1,
		vector< double > &	delu,
		vector< double > &	p,
		vector< double > &	qtemp,
		vector< double > &	Lambda,
		const vector< double > &	Phi,
		const vector< double > &	PhiDrv,
		double	step,
		double *	cooling,
		double	p_max,
		long	index,
		size_t	nd,
		bool *	lgBoundFail )

References ASSERT, DEBUG_ENTRY, gv, inv_ufunct(), log_integral(), MAX2, NDEMS, NQGRID, PROB_CUTOFF_LO, QHEAT_TOL, rfield, splint_safe(), STATIC, and vlog().

Referenced by GetProbDistr_LowLimit().

Here is the call graph for this function:

◆ uderiv()

STATIC double uderiv	(	double	temp,
		size_t	nd )

References ASSERT, cdEXIT, cval, DEBUG_ENTRY, DebyeDeriv(), ENTH_CAR, ENTH_CAR2, ENTH_PAH, ENTH_PAH2, ENTH_SIC, ENTH_SIL, ENTH_SIL2, EXIT_FAILURE, exp10(), fprintf(), gv, ioQQQ, MAX2, MIN2, no_atoms(), pow(), pow2, ppower, STATIC, and tlim.

Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), InitEnthalpy(), qheat(), and RebinQHeatResults().

Here is the call graph for this function:

◆ ufunct()

STATIC double ufunct	(	double	temp,
		size_t	nd,
		bool *	lgBoundErr )

References ASSERT, cdEXIT, DEBUG_ENTRY, EXIT_FAILURE, fprintf(), gv, ioQQQ, NDEMS, splint_safe(), and STATIC.

Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), qheat(), and RebinQHeatResults().

Here is the call graph for this function:

Variable Documentation

◆ cval

const double cval[4]

static

Initial value:

={
        1.40e3/DEN_SIL*ATOMIC_MASS_UNIT*MW_SIL/EN1RYD,
        2.20e4/DEN_SIL*ATOMIC_MASS_UNIT*MW_SIL/EN1RYD,
        4.80e5/DEN_SIL*ATOMIC_MASS_UNIT*MW_SIL/EN1RYD,
        3.41e7/DEN_SIL*ATOMIC_MASS_UNIT*MW_SIL/EN1RYD}

Referenced by uderiv().

◆ DEF_FAC

const double DEF_FAC = 3.

static

Referenced by GetProbDistr_LowLimit(), and qheat().

◆ DEN_SIL

const double DEN_SIL = 3.30

static

◆ FWHM_RATIO

const double FWHM_RATIO = 0.1

static

Referenced by qheat().

◆ FWHM_RATIO2

const double FWHM_RATIO2 = 0.007

static

Referenced by qheat().

◆ LIM1

const realnum LIM1 = 2.e-4f

static

Referenced by fast_expm1().

◆ LIM2

const realnum LIM2 = sqrtf(6.e-4f)

static

Referenced by fast_expm1(), and qheat_init().

◆ LIM3

const realnum LIM3 = cbrtf(24.e-4f)

static

Referenced by fast_expm1(), and qheat_init().

◆ LOOP_MAX

const long LOOP_MAX = 20L

static

Referenced by cnewton(), lgCoolNetConverge(), and qheat().

◆ MAX_EVENTS

const double MAX_EVENTS = 150.

static

Referenced by qheat().

◆ MAX_LOOP

const long MAX_LOOP = 2*NQGRID

static

Referenced by GetProbDistr_LowLimit().

◆ MW_SIL

const double MW_SIL = 24.6051

static

◆ NQMIN

const long NQMIN = 10L

static

Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), RebinQHeatResults(), and ScanProbDistr().

◆ NQTEST

const long NQTEST = 500L

static

Referenced by GetProbDistr_LowLimit().

◆ NSTARTUP

const long NSTARTUP = 5L

static

Referenced by GetProbDistr_LowLimit().

◆ ppower

const double ppower[4] ={2.00,1.30,0.68,0.00}

static

Referenced by uderiv().

◆ PROB_CUTOFF_HI

const double PROB_CUTOFF_HI = 1.e-20

static

Referenced by GetProbDistr_HighLimit(), and GetProbDistr_LowLimit().

◆ PROB_CUTOFF_LO

const double PROB_CUTOFF_LO = 1.e-15

static

Referenced by GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), qheat(), RebinQHeatResults(), and TryDoubleStep().

◆ PROB_TOL

const double PROB_TOL = 0.02

static

Referenced by GetProbDistr_HighLimit(), and GetProbDistr_LowLimit().

◆ QHEAT_TOL

const double QHEAT_TOL = 5.e-3

static

Referenced by GetProbDistr_LowLimit(), and TryDoubleStep().

◆ QT_RATIO

const double QT_RATIO = 1.074

static

Referenced by GetProbDistr_HighLimit(), and RebinQHeatResults().

◆ RELAX

const double RELAX = 15.

static

Referenced by qheat().

◆ SAFETY

const double SAFETY = 1.e+8

static

Referenced by gauss_legendre(), GetProbDistr_HighLimit(), GetProbDistr_LowLimit(), radius_next(), RebinQHeatResults(), spldrv_safe(), and splint_safe().

◆ STRICT

const double STRICT = 1.

static

Referenced by GetProbDistr_HighLimit(), and qheat().

◆ tlim

const double tlim[5] ={0.,50.,150.,500.,DBL_MAX}

static

Referenced by InitEnthalpy(), and uderiv().

◆ WIDE_FAIL_MAX

const long WIDE_FAIL_MAX = 3

static

Referenced by GetProbDistr_LowLimit(), qheat(), and ScanProbDistr().

Enumerations

Functions

Variables

Enumeration Type Documentation

◆ QH_Code

Function Documentation

◆ DebyeDeriv()

◆ fast_expm1()

◆ GetProbDistr_HighLimit()

◆ GetProbDistr_LowLimit() [1/2]

◆ GetProbDistr_LowLimit() [2/2]

◆ GrainMakeDiffuse()

◆ GrainMakeDiffuseSingle()

◆ InitEnthalpy()

◆ inv_ufunct()

◆ log_integral()

◆ no_atoms()

◆ qheat()

◆ qheat_init()

◆ RebinQHeatResults()

◆ ScanProbDistr()

◆ TryDoubleStep()

◆ uderiv()

◆ ufunct()

Variable Documentation

◆ cval

◆ DEF_FAC

◆ DEN_SIL

◆ FWHM_RATIO

◆ FWHM_RATIO2

◆ LIM1

◆ LIM2

◆ LIM3

◆ LOOP_MAX

◆ MAX_EVENTS

◆ MAX_LOOP

◆ MW_SIL

◆ NQMIN

◆ NQTEST

◆ NSTARTUP

◆ ppower

◆ PROB_CUTOFF_HI

◆ PROB_CUTOFF_LO

◆ PROB_TOL

◆ QHEAT_TOL

◆ QT_RATIO

◆ RELAX

◆ SAFETY

◆ STRICT

◆ tlim

◆ WIDE_FAIL_MAX