Modules
module	m_pressure_relaxation
	Pressure relaxation for the six-equation multi-component model via Newton–Raphson equilibration and volume-fraction correction.

Functions/Subroutines
subroutine, public	m_pressure_relaxation::s_pressure_relaxation_procedure (q_cons_vf)
	The main pressure relaxation procedure.
impure subroutine, public	m_pressure_relaxation::s_report_pressure_relaxation
	One line at the end of a run if the equilibration ever stopped on its iteration cap. Silence means every cell reached the tolerance.
logical function	m_pressure_relaxation::s_needs_pressure_relaxation (q_cons_vf, j, k, l)
	Check if pressure relaxation is needed for this cell.
subroutine	m_pressure_relaxation::s_correct_volume_fractions (q_cons_vf, j, k, l)
	Correct volume fractions to physical bounds.
subroutine	m_pressure_relaxation::s_equilibrate_pressure (q_cons_vf, j, k, l, hit_cap, resid, unusable)
	Main pressure equilibration using Newton-Raphson.
subroutine	m_pressure_relaxation::s_correct_internal_energies (q_cons_vf, j, k, l, rho, gamma, pi_inf, qv_mix)
	Correct internal energies using equilibrated pressure.

Variables
integer	m_pressure_relaxation::n_hit_cap = 0
	Cell updates in which the Newton below stopped on its iteration cap rather than on the tolerance. Counted because the answer then depends on that cap, and nothing else in the code says so.
real(wp)	m_pressure_relaxation::worst_residual = 0._wp
	Largest \|f\| left behind when the cap was reached.

Detailed Description

Modules