solve_generalized_point Subroutine

public subroutine solve_generalized_point(model, z, np, kinds_x, kind_w, X, ns1, S1, ns2, S2, max_iters, F, df, iters)

Uses

  • proc~~solve_generalized_point~~UsesGraph proc~solve_generalized_point solve_generalized_point module~yaeos__math yaeos__math proc~solve_generalized_point->module~yaeos__math module~yaeos__constants yaeos__constants module~yaeos__math->module~yaeos__constants module~yaeos__math_continuation yaeos__math_continuation module~yaeos__math->module~yaeos__math_continuation module~yaeos__math_linalg yaeos__math_linalg module~yaeos__math->module~yaeos__math_linalg iso_fortran_env iso_fortran_env module~yaeos__constants->iso_fortran_env module~yaeos__math_continuation->module~yaeos__constants module~yaeos__math_continuation->module~yaeos__math_linalg module~yaeos__math_linalg->module~yaeos__constants

Function to solve the multiphase flash problem.

Arguments

Type IntentOptional Attributes Name
class(ArModel), intent(in) :: model

Model to use.
real(kind=pr), intent(in) :: z(:)

Mixture global composition.
integer, intent(in) :: np

Number of x phases.
character(len=14), intent(in) :: kinds_x(np)

Kind of the x phases.
character(len=14), intent(in) :: kind_w

Kind of the w phase.
real(kind=pr), intent(inout) :: X(:)

Vector of variables.
integer, intent(in) :: ns1

Number of first specification.
real(kind=pr), intent(in) :: S1

First specification value.
integer, intent(in) :: ns2

Number of second specification.
real(kind=pr), intent(in) :: S2

Second specification value.
integer, intent(in) :: max_iters

Maximum number of iterations.
real(kind=pr), intent(out) :: F(size(X))

Vector of functions valuated.
real(kind=pr), intent(out), dimension(size(X), size(X)) :: df

Jacobian matrix.
integer, intent(out) :: iters

Number of iterations performed.

Calls

proc~~solve_generalized_point~~CallsGraph proc~solve_generalized_point solve_generalized_point proc~pt_f_np~3 pt_F_NP proc~solve_generalized_point->proc~pt_f_np~3 proc~solve_system solve_system proc~solve_generalized_point->proc~solve_system proc~lnphi_pt ArModel%lnphi_pt proc~pt_f_np~3->proc~lnphi_pt none~dgesv dgesv proc~solve_system->none~dgesv proc~lnphi_vt ArModel%lnphi_vt proc~lnphi_pt->proc~lnphi_vt proc~volume~3 ArModel%volume proc~lnphi_pt->proc~volume~3 residual_helmholtz residual_helmholtz proc~lnphi_vt->residual_helmholtz get_v0 get_v0 proc~volume~3->get_v0 interface~newton newton proc~volume~3->interface~newton proc~newton_1d newton_1d interface~newton->proc~newton_1d

Called by

proc~~solve_generalized_point~~CalledByGraph proc~solve_generalized_point solve_generalized_point proc~create_generalized_isoz_line create_generalized_isoz_line proc~create_generalized_isoz_line->proc~solve_generalized_point

Variables

Type Visibility Attributes Name Initial
real(kind=pr), public, dimension(size(X)) :: dX

Newton step vector.
integer, public :: i
integer, public :: iBetas(np+1)

Index of the betas in the vector X
integer, public :: nc

Number of components