pressure_change.cpp File Reference

#include "cddefines.h"
#include "pressure_change.h"
#include "colden.h"
#include "conv.h"
#include "cosmology.h"
#include "dark_matter.h"
#include "dense.h"
#include "dynamics.h"
#include "geometry.h"
#include "phycon.h"
#include "pressure.h"
#include "radius.h"
#include "struc.h"
#include "thermal.h"
#include "trace.h"
#include "wind.h"

Include dependency graph for pressure_change.cpp:

Enumerations
enum	{   CPRE, SUBSONIC, SUPERSONIC, STRONGD,   ORIGINAL, SHOCK, ANTISHOCK, ANTISHOCK2 }

Functions
double	zoneDensity ()
STATIC double	stepDensity (const PresMode &presmode, solverState &st)
STATIC void	logPressureState ()
STATIC bool	lgTestPressureConvergence (double new_density)
STATIC double	limitedDensityScaling (double new_density, double dP_chng_factor)
void	PressureChange (double dP_chng_factor, const PresMode &presmode, solverState &st, bool &lgStable)
double	pressureZone (const PresMode &presmode)

Enumeration Type Documentation

anonymous enum

Enumerator
CPRE
SUBSONIC
SUPERSONIC
STRONGD
ORIGINAL
SHOCK
ANTISHOCK
ANTISHOCK2

Function Documentation

STATIC bool lgTestPressureConvergence

(

double

new_density

)

References conv, t_conv::PressureErrorAllowed, and scalingDensity().

Referenced by PressureChange().

Here is the call graph for this function:

STATIC double limitedDensityScaling	(	double	new_density,
		double	dP_chng_factor
	)

References conv, MAX2, t_conv::MaxFractionalDensityStepPerIteration, MIN2, and scalingDensity().

Referenced by PressureChange().

Here is the call graph for this function:

STATIC void logPressureState

(

)

References conv, dense, t_dense::gas_phase, t_conv::hist_pres_current, t_conv::hist_pres_density, t_conv::hist_pres_error, ipHYDROGEN, pressure, t_pressure::PresTotlCurr, and t_pressure::PresTotlError.

Referenced by PressureChange().

void PressureChange	(	double	dP_chng_factor,
		const PresMode &	presmode,
		solverState &	st,
		bool &	lgStable
	)

References conv, ConvTempEdenIoniz(), DEBUG_ENTRY, dense, t_radius::depth, t_geometry::FillFac, fixit, fprintf(), t_dense::gas_phase, geometry, ioQQQ, ipHYDROGEN, t_conv::lgConvPres, lgTestPressureConvergence(), t_trace::lgTrace, limitedDensityScaling(), logPressureState(), nzone, phycon, t_pressure::pinzon, t_pressure::pres_radiation_lines_curr, t_pressure::PresGasCurr, t_pressure::PresInteg, t_pressure::PresRamCurr, pressure, PresTotCurrent(), t_pressure::PresTotlCurr, t_pressure::PresTotlError, t_pressure::PresTotlInit, radius, ScaleAllDensities(), scalingDensity(), stepDensity(), SUBSONIC, SUPERSONIC, t_phycon::te, TempChange(), TorF(), trace, wind, Wind::windv, t_dense::xMassDensity, and PresMode::zone.

Referenced by ConvPresTempEdenIoniz().

Here is the call graph for this function:

double pressureZone

(

const PresMode &

presmode

)

References CPRE, dark, DEBUG_ENTRY, t_dynamics::DivergePresInteg, dynamics, fixit, PresMode::global, t_pressure::gravity_symmetry, t_pressure::IntegRhoGravity, t_pressure::lgContRadPresOn, t_dark_matter::lgNFW_Set, pow(), t_pressure::PresInteg, t_pressure::PresPowerlaw, pressure, t_pressure::PresTotlInit, radius, t_radius::Radius, and t_radius::rinner.

Referenced by ConvPresTempEdenIoniz().

Here is the call graph for this function:

STATIC double stepDensity	(	const PresMode &	presmode,
		solverState &	st
	)

returns updated best guess for equilibrium scalingDensity

References t_dynamics::chPresMode, CPRE, DEBUG_ENTRY, t_radius::depth, t_dynamics::DivergePresInteg, solverState::dp, DynaFlux(), dynamics, solverState::erp, fprintf(), PresMode::global, ioQQQ, solverState::lastzone, t_pressure::lgContRadPresOn, t_trace::lgTrace, t_dynamics::lgTracePrint, t_trace::nTrConvg, nzone, ORIGINAL, t_pressure::PresInteg, solverState::press, pressure, t_pressure::PresTotlCurr, t_pressure::PresTotlError, t_pressure::PresTotlInit, radius, scalingDensity(), SMALLFLOAT, SUBSONIC, SUPERSONIC, trace, wind, Wind::windv, and PresMode::zone.

Referenced by PressureChange().

Here is the call graph for this function:

double zoneDensity

(

)

References Wind::AccelGravity, Wind::AccelTotalOutward, cdEXIT, t_dense::cfirst, t_dense::chDenseLaw, colden, t_colden::colden, cosmology, t_dense::csecnd, DEBUG_ENTRY, t_dense::den0, dense, dense_fabden(), dense_parametric_wind(), t_dense::DensityPower, t_radius::depth, t_dense::DLW, t_radius::drad, Wind::emdot, EXIT_FAILURE, fixit, t_dense::flcPhase, t_dense::flong, fprintf(), t_dense::gas_phase, GetDensity(), t_radius::glbden, t_radius::glbdst, t_radius::glbpow, t_radius::glbrad, ioQQQ, ipCOL_HTOT, ipHYDROGEN, Wind::lgBallistic(), t_dense::lgDenFlucRadius, Wind::lgDisk, t_cosmology::lgDo, Wind::lgStatic(), t_trace::lgTrace, Wind::lgVelPos, t_trace::lgWind, nzone, pow(), POW2, pressure, t_pressure::PresTotlError, t_radius::r1r0sq, radius, t_radius::Radius, t_cosmology::redshift_current, t_radius::rinner, t_dense::rscale, scalingDensity(), ShowMe(), struc, DepthTable::tabval(), TorF(), TotalInsanity(), trace, wind, Wind::windv, and t_struc::windv.

Referenced by ConvPresTempEdenIoniz().

Here is the call graph for this function: