solve_critical_point Subroutine

private subroutine solve_critical_point(model, nc, i_comp, V_init, T_init, Vc, Tc, Pc, converged)

Arguments

Type IntentOptional Attributes Name
class(ArModel), intent(in) :: model
integer, intent(in) :: nc
integer, intent(in) :: i_comp
real(kind=pr), intent(in) :: V_init
real(kind=pr), intent(in) :: T_init
real(kind=pr), intent(out) :: Vc
real(kind=pr), intent(out) :: Tc
real(kind=pr), intent(out) :: Pc
logical, intent(out) :: converged

Calls

proc~~solve_critical_point~~CallsGraph proc~solve_critical_point solve_critical_point proc~jacobian_critical_eqs jacobian_critical_eqs proc~solve_critical_point->proc~jacobian_critical_eqs proc~residual_critical_eqs residual_critical_eqs proc~solve_critical_point->proc~residual_critical_eqs proc~jacobian_critical_eqs->proc~residual_critical_eqs proc~pressure ArModel%pressure proc~residual_critical_eqs->proc~pressure residual_helmholtz residual_helmholtz proc~pressure->residual_helmholtz

Called by

proc~~solve_critical_point~~CalledByGraph proc~solve_critical_point solve_critical_point proc~find_critical_points_all_components find_critical_points_all_components proc~find_critical_points_all_components->proc~solve_critical_point proc~init_pcsaft init_pcsaft proc~init_pcsaft->proc~find_critical_points_all_components

Variables

Type Visibility Attributes Name Initial
real(kind=pr), private :: J(2,2)
real(kind=pr), private :: J_inv(2,2)
real(kind=pr), private :: P
real(kind=pr), private :: T
real(kind=pr), private :: V
real(kind=pr), private :: dX(2)
real(kind=pr), private :: det
real(kind=pr), private :: f1
real(kind=pr), private :: f2
integer, private :: iter
real(kind=pr), private :: norm_f