/home66/gary/public_html/cloudy/c13_branch/source/grains_qheat.cpp File Reference

#include "cddefines.h"
#include "physconst.h"
#include "rfield.h"
#include "phycon.h"
#include "dense.h"
#include "hmi.h"
#include "thermal.h"
#include "trace.h"
#include "thirdparty.h"
#include "iterations.h"
#include "grainvar.h"
#include "grains.h"
Include dependency graph for grains_qheat.cpp:

Go to the source code of this file.

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 *, long, size_t, bool *)
STATIC double log_integral (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)
void GrainMakeDiffuse (void)
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 (void)

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.03
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]

Enumeration Type Documentation

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 

Definition at line 27 of file grains_qheat.cpp.


Function Documentation

STATIC double DebyeDeriv ( double  x,
long  n 
)

Definition at line 2532 of file grains_qheat.cpp.

References ASSERT, DEBUG_ENTRY, gauss_init(), gauss_legendre(), MAX2, POW2, POW3, powi(), and TotalInsanity().

Referenced by uderiv().

Here is the call graph for this function:

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 
)
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 
)
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  
)

Referenced by qheat().

void GrainMakeDiffuse ( void   ) 
void InitEnthalpy ( void   ) 

initialize interpolation arrays for grain enthalpy

Definition at line 2469 of file grains_qheat.cpp.

References GrainVar::bin, DEBUG_ENTRY, GrainVar::dsttmp, GRAIN_TMIN, gv, log_integral(), NDEMS, spline(), tlim, and uderiv().

Referenced by GrainsInit().

Here is the call graph for this function:

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

Definition at line 2443 of file grains_qheat.cpp.

References ASSERT, GrainVar::bin, cdEXIT, DEBUG_ENTRY, GrainVar::dsttmp, EXIT_FAILURE, gv, ioQQQ, NDEMS, and splint_safe().

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

Here is the call graph for this function:

STATIC double log_integral ( double  xx1,
double  yy1,
double  xx2,
double  yy2 
)

Definition at line 1785 of file grains_qheat.cpp.

References ASSERT, and DEBUG_ENTRY.

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

double no_atoms ( size_t  nd  )  [inline]

Definition at line 17 of file grains_qheat.cpp.

References ATOMIC_MASS_UNIT, GrainVar::bin, and gv.

Referenced by uderiv().

void qheat ( vector< double > &  qtemp,
vector< double > &  qprob,
long int *  qnbin,
size_t  nd 
)
STATIC void qheat_init ( size_t  nd,
vector< double > &  phiTilde,
double *  check 
)
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 
)
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 
)
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,
long  index,
size_t  nd,
bool *  lgBoundFail 
)

Definition at line 1612 of file grains_qheat.cpp.

References t_rfield::anu, ASSERT, GrainVar::bin, DEBUG_ENTRY, GrainVar::dsttmp, EN1RYD, gv, inv_ufunct(), log_integral(), MAX2, NDEMS, NQGRID, PROB_CUTOFF_LO, QHEAT_TOL, rfield, and splint_safe().

Referenced by GetProbDistr_LowLimit().

Here is the call graph for this function:

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

Definition at line 2417 of file grains_qheat.cpp.

References ASSERT, GrainVar::bin, cdEXIT, DEBUG_ENTRY, GrainVar::dsttmp, EXIT_FAILURE, gv, ioQQQ, NDEMS, and splint_safe().

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

Here is the call graph for this function:


Variable Documentation

const double cval[4] [static]
const double DEF_FAC = 3. [static]

Definition at line 69 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit(), and qheat().

const double DEN_SIL = 3.30 [static]

Definition at line 118 of file grains_qheat.cpp.

const double FWHM_RATIO = 0.1 [static]

Definition at line 80 of file grains_qheat.cpp.

Referenced by qheat().

const double FWHM_RATIO2 = 0.007 [static]

Definition at line 83 of file grains_qheat.cpp.

Referenced by qheat().

const long LOOP_MAX = 20L [static]

Definition at line 66 of file grains_qheat.cpp.

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

const double MAX_EVENTS = 150. [static]

Definition at line 62 of file grains_qheat.cpp.

Referenced by qheat().

const long MAX_LOOP = 2*NQGRID [static]

Definition at line 86 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit().

const double MW_SIL = 24.6051 [static]

Definition at line 121 of file grains_qheat.cpp.

const long NQMIN = 10L [static]
const long NQTEST = 500L [static]

Definition at line 76 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit().

const long NSTARTUP = 5L [static]

Definition at line 58 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit().

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

Definition at line 125 of file grains_qheat.cpp.

Referenced by uderiv().

const double PROB_CUTOFF_HI = 1.e-20 [static]

Definition at line 54 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), and GetProbDistr_LowLimit().

const double PROB_CUTOFF_LO = 1.e-15 [static]
const double PROB_TOL = 0.02 [static]

Definition at line 72 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), and GetProbDistr_LowLimit().

const double QHEAT_TOL = 5.e-3 [static]

Definition at line 89 of file grains_qheat.cpp.

Referenced by GetProbDistr_LowLimit(), and TryDoubleStep().

const double QT_RATIO = 1.03 [static]

Definition at line 103 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), and RebinQHeatResults().

const double RELAX = 15. [static]

Definition at line 96 of file grains_qheat.cpp.

Referenced by qheat().

const double SAFETY = 1.e+8 [static]
const double STRICT = 1. [static]

Definition at line 95 of file grains_qheat.cpp.

Referenced by GetProbDistr_HighLimit(), and qheat().

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

Definition at line 124 of file grains_qheat.cpp.

Referenced by InitEnthalpy(), and uderiv().

const long WIDE_FAIL_MAX = 3 [static]

Definition at line 92 of file grains_qheat.cpp.

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

 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines

Generated on 15 Nov 2012 for cloudy by  doxygen 1.6.1