#include <conv.h>
the variables that deal with the convergence of the model
Definition at line 77 of file conv.h.
a convergence criteria
Definition at line 205 of file conv.h.
Referenced by ConvIterCheck(), ParseIterations(), and zero().
remember the average electron density error
Definition at line 127 of file conv.h.
Referenced by ConvInitSolution(), PrtFinal(), radius_increment(), and zero().
Definition at line 131 of file conv.h.
Referenced by ConvInitSolution(), PrtFinal(), radius_increment(), and zero().
Definition at line 135 of file conv.h.
Referenced by ConvInitSolution(), PrtFinal(), radius_increment(), and zero().
double t_conv::BadConvIoniz[2] |
when the lgConvIoniz flag is set false, the old and new numbers, the reason for the lack of convergence, should be set to following
Definition at line 98 of file conv.h.
Referenced by CO_drive(), CO_solve(), ConvBase(), ConvFail(), ConvIoniz(), GrainChargeTemp(), and lgIonizConverg().
remember the biggest electron density error as test of convergence quality
Definition at line 169 of file conv.h.
Referenced by ConvInitSolution(), PrtFinal(), radius_increment(), and zero().
remember the biggest and average heating-cooling error
Definition at line 130 of file conv.h.
Referenced by ConvInitSolution(), PrtFinal(), radius_increment(), and zero().
remember the biggest and average pressure error
Definition at line 134 of file conv.h.
Referenced by ConvInitSolution(), PrtFinal(), radius_increment(), and zero().
this says why the ionization did not converge, reasons can be a large change in the level of ionization, or in the heating
Definition at line 84 of file conv.h.
Referenced by CO_drive(), CO_solve(), ConvBase(), ConvEdenIoniz(), ConvFail(), ConvIoniz(), get_total_abundance_ions(), GrainChargeTemp(), H2_LevelPops(), lgConvEden(), and zero().
this gives the reason the model was declared not converged when 'iter to convergence' command is given
Definition at line 88 of file conv.h.
Referenced by ConvIterCheck(), PrtComment(), and zero().
char t_conv::chSolverEden[20] |
which electron density solver to use? set with set eden solver command, simple and new
Definition at line 191 of file conv.h.
Referenced by ConvEdenIoniz(), and zero().
char t_conv::chSolverTemp[20] |
which temperature density solver to use? set with set eden solver command, default and brent
Definition at line 195 of file conv.h.
Referenced by ConvTempEdenIoniz(), and zero().
double t_conv::dCmHdT |
numerical estimate of d(cooling-heating)/dT
Definition at line 222 of file conv.h.
Referenced by ConvTempEdenIoniz(), IterStart(), lgConvTemp(), and zero().
double t_conv::EdenErrorAllowed |
this is relative error in the electron density we want set in zero to 0.01 reset with set eden error command
Definition at line 210 of file conv.h.
Referenced by advection_set_default(), ConvBase(), GrainDrive(), GrainsInit(), H2_LevelPops(), hmole(), IonHydro(), lgConvEden(), ParseSet(), and zero().
failmx is largest relative error in heating cooling match
Definition at line 160 of file conv.h.
Referenced by ConvFail(), IterStart(), and zero().
allowed error in heating - cooling balance, set with TOLERANCE command default set in zerologic
Definition at line 219 of file conv.h.
Referenced by advection_set_default(), ConvBase(), ConvFail(), ConvInitSolution(), CoolEvaluate(), CoolHeatError(), GrainsInit(), H2_LevelPops(), lgConvTemp(), lgCoolHeatCheckConverge(), ParseSet(), and zero().
vector<double> t_conv::hist_pres_correct |
Definition at line 229 of file conv.h.
Referenced by PressureChange(), and SaveDo().
vector<double> t_conv::hist_pres_current |
Definition at line 229 of file conv.h.
Referenced by PressureChange(), and SaveDo().
vector<double> t_conv::hist_pres_density |
these are used to retain the density pressure history in current zone, can be output with save pressure history
Definition at line 229 of file conv.h.
Referenced by PressureChange(), and SaveDo().
long int t_conv::hist_pres_nzone |
Definition at line 230 of file conv.h.
Referenced by ConvInitSolution(), and PressureChange().
vector<double> t_conv::hist_temp_cool |
Definition at line 234 of file conv.h.
Referenced by CoolHeatError(), and SaveDo().
vector<double> t_conv::hist_temp_heat |
Definition at line 234 of file conv.h.
Referenced by CoolHeatError(), and SaveDo().
long int t_conv::hist_temp_nzone |
Definition at line 235 of file conv.h.
Referenced by ConvInitSolution(), and CoolHeatError().
vector<double> t_conv::hist_temp_temp |
these are used to retain the temp/heat/cooling history in current zone, can be output with save temperature history
Definition at line 234 of file conv.h.
Referenced by CoolHeatError(), and SaveDo().
long int t_conv::ifailz[12] |
zones where converge problems occurred
Definition at line 187 of file conv.h.
Referenced by ConvFail().
bool t_conv::lgAutoIt |
says "interate to convergence" command given
Definition at line 202 of file conv.h.
Referenced by ConvIterCheck(), ParseCommands(), ParseIterations(), PrtComment(), and zero().
bool t_conv::lgBadStop |
flag saying that calculation stopped for bad reason mostly set in lgEndfun
Definition at line 199 of file conv.h.
Referenced by iter_end_check(), PrtComment(), and PrtFinal().
bool t_conv::lgConvEden |
true when the electron density has converged
Definition at line 151 of file conv.h.
Referenced by ConvEdenIoniz(), ConvPresTempEdenIoniz(), CoolSave(), lgConvEden(), transition::outline(), SaveHeat(), and zero().
bool t_conv::lgConvIoniz |
this flag is used in ConvPresTempEdenIoniz to check that ionization has converged
Definition at line 91 of file conv.h.
Referenced by CO_drive(), CO_solve(), ConvBase(), ConvFail(), ConvIoniz(), ConvPresTempEdenIoniz(), CoolSave(), EdenError(), get_total_abundance_ions(), GrainChargeTemp(), transition::outline(), RT_OTS_PrtRate(), and SaveHeat().
bool t_conv::lgConvPops |
this flag is used in ConvPresTempEdenIoniz to check that populations have converged
Definition at line 94 of file conv.h.
Referenced by CO_drive(), CO_solve(), ConvFail(), ConvPresTempEdenIoniz(), H2_LevelPops(), hmole(), and transition::outline().
bool t_conv::lgConvPres |
true if pressure is converged, false if not
Definition at line 148 of file conv.h.
Referenced by ConvPresTempEdenIoniz(), CoolSave(), lgConvPres(), transition::outline(), PressureChange(), SaveDo(), SaveHeat(), and zero().
bool t_conv::lgConvTemp |
true if temperature is converged, false if not
Definition at line 145 of file conv.h.
Referenced by ConvFail(), ConvPresTempEdenIoniz(), ConvTempEdenIoniz(), CoolSave(), lgConvTemp(), transition::outline(), SaveDo(), SaveHeat(), and zero().
first sweep through solvers in this zone. also true is search phase
Definition at line 104 of file conv.h.
Referenced by ConvBase(), ConvInitSolution(), ConvPresTempEdenIoniz(), Fe2_cooling(), RT_line_all(), RT_line_escape(), RT_line_pumping(), RT_tau_inc(), and RT_tau_init().
flag set in ConvBase, saying whether ionization stage is trimmed down
Definition at line 138 of file conv.h.
Referenced by ConvBase(), ion_trim(), and RT_line_all().
last sweep through solvers in this zone, so update fine opacities
Definition at line 106 of file conv.h.
Referenced by ConvInitSolution(), ConvPresTempEdenIoniz(), Fe2_cooling(), RT_line_all(), RT_line_fine_opacity(), RT_line_one(), and RT_tau_inc().
bool t_conv::lgMap |
lgMap is option to map failures
Definition at line 184 of file conv.h.
Referenced by ConvFail(), ParseFail(), and zero().
bool t_conv::lgOscilOTS |
this is true if ots rates are oscillating, and this is why ionization is not converged
Definition at line 142 of file conv.h.
Referenced by ConvBase(), ConvInitSolution(), ConvIoniz(), and RT_tau_init().
bool t_conv::lgSearch |
conv.lgSearch is true during initial temp-ion search phase, false after first zone established
Definition at line 124 of file conv.h.
Referenced by atom_levelN(), CO_drive(), CO_solve(), ConvBase(), ConvEdenIoniz(), ConvInitSolution(), ConvIoniz(), CoolEvaluate(), dBase_solve(), Fe2_cooling(), get_total_abundance_ions(), GrainCharge(), GrainChargeTemp(), GrainDrive(), H21_cm_pops(), H2_Cooling(), H2_LevelPops(), hmole(), hmole_step(), ion_trim(), iso_level(), iter_end_check(), map_do(), MyAssert(), OpacityAddTotal(), PresTotCurrent(), and store_new_densities().
long int t_conv::LimFail |
LimFail is limit to number of te failures, set with "failures" cmnd
Definition at line 181 of file conv.h.
Referenced by ConvFail(), grid_do(), ParseFail(), and zero().
long int t_conv::limPres2Ioniz |
a limit to the above, in case where one zone takes forever to not converge, usually very large, set with SET PRESIONIZ command
Definition at line 110 of file conv.h.
Referenced by ConvBase(), ParseSet(), and zero().
long int t_conv::nGrainFail |
number of grain ionization balance failures
Definition at line 178 of file conv.h.
Referenced by ConvFail(), IterStart(), and zero().
long int t_conv::nIonFail |
nIonFail is number of ionization failures
Definition at line 172 of file conv.h.
Referenced by cdNwcns(), ConvFail(), IterStart(), and zero().
long int t_conv::nNeFail |
nNeFail is number of electron density failures
Definition at line 166 of file conv.h.
Referenced by cdNwcns(), ConvFail(), IterStart(), and zero().
long int t_conv::nPopFail |
nIonFail is number of level population failures
Definition at line 175 of file conv.h.
Referenced by ConvFail(), IterStart(), and zero().
long int t_conv::nPreFail |
nPreFail is number of pressure failures
Definition at line 163 of file conv.h.
Referenced by cdNwcns(), ConvFail(), IterStart(), and zero().
long int t_conv::nPres2Ioniz |
this will count the number of ionizations in one call from ConvPresTempEdenIoniz
Definition at line 101 of file conv.h.
Referenced by ConvBase(), ConvInitSolution(), ConvIoniz(), ConvPresTempEdenIoniz(), CoolEvaluate(), hmole_reactions(), ion_photo(), lgIonizConverg(), OpacityAddTotal(), PrtZone(), RT_line_all(), SaveDo(), and zero().
long int t_conv::nTeFail |
nTeFail number of temperature failures
Definition at line 157 of file conv.h.
Referenced by cdNwcns(), ConvFail(), IterStart(), ShowMe(), and zero().
long int t_conv::nTotalFailures |
total number of all falures, used to trigger abort
Definition at line 154 of file conv.h.
Referenced by ConvFail(), IterStart(), and zero().
long int t_conv::nTotalIoniz |
counts the number of calls to conv base after iteration starts can be used to determine very first pass through an iteration reset to zero at start of each iteration in a simulation
Definition at line 115 of file conv.h.
Referenced by ChargTranEval(), CO_solve(), ConvBase(), ConvInitSolution(), dBase_solve(), fill_array(), find_solution(), get_total_abundance_ions(), hmole(), ion_trim(), IonHelium(), iso_collide(), iso_continuum_lower(), iso_drive(), iso_ionize_recombine(), iso_level(), iso_radiative_recomb(), iso_update_rates(), iter_end_check(), lgCheckMonitors(), lgHomogeneousSource(), map_do(), PresTotCurrent(), PrtFinal(), RT_DestProb(), vib_evap(), and zero().
long int t_conv::nTotalIoniz_start |
the same counter but set to zero after the initial solution is converged, so this is a measure of the number of calls within true zones
Definition at line 120 of file conv.h.
Referenced by ConvInitSolution(), lgCheckMonitors(), and PrtFinal().
this is relative error in the pressure, initialized to 0.02 in reset with set pressure error command
Definition at line 215 of file conv.h.
Referenced by advection_set_default(), lgConvPres(), ParseSet(), PresTotCurrent(), and zero().
double t_conv::sigma_dCmHdT |
1-sigma uncertainty in numerical estimate of d(cooling-heating)/dT
Definition at line 225 of file conv.h.
Referenced by ConvTempEdenIoniz(), and lgConvTemp().