Calculation of the first restriction of a critical point.
is the smallest eigen-value for the matrix:
Where And should be the eigen-vector corresponding to the smallest eigen-value when
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(ArModel), | intent(in) | :: | model | |||
real(kind=pr), | intent(in) | :: | X(4) |
Vector of variables |
||
real(kind=pr), | intent(in) | :: | s |
Distance between the two fluids compositions to the main composition |
||
real(kind=pr), | intent(in) | :: | z0(:) |
Molar fractions of the first fluid |
||
real(kind=pr), | intent(in) | :: | zi(:) |
Molar fractions of the second fluid |
||
real(kind=pr), | intent(in) | :: | u(:) |
Eigen-vector that defines the direction between the two compositions |
||
real(kind=pr), | intent(out), | optional | :: | u_new(:) |
Eigen-vector corresponding to the smallest eigenvalue of the matrix |
|
real(kind=pr), | intent(out), | optional | :: | P |
Pressure of the system [bar] |
Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|
real(kind=pr), | private | :: | M(size(z0),size(z0)) | ||||
real(kind=pr), | private | :: | Pin | ||||
real(kind=pr), | private | :: | T | ||||
real(kind=pr), | private | :: | V | ||||
real(kind=pr), | private | :: | a | ||||
real(kind=pr), | private | :: | dlnf_dn(size(z0),size(z0)) | ||||
integer, | private | :: | i | ||||
integer, | private | :: | j | ||||
real(kind=pr), | private | :: | lambda(size(z0)) | ||||
real(kind=pr), | private | :: | n(size(z0)) | ||||
integer, | private | :: | nc | ||||
real(kind=pr), | private | :: | vectors(size(z0),size(z0)) | ||||
real(kind=pr), | private | :: | z(size(z0)) |