homotopy Subroutine

public subroutine homotopy(sub, x, tol, max_its, its)

Uses

  • proc~~homotopy~~UsesGraph proc~homotopy homotopy module~yaeos__math yaeos__math proc~homotopy->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

Arguments

Type IntentOptional Attributes Name
procedure(to_solve) :: sub
real(kind=pr), intent(inout) :: x(:)
real(kind=pr), intent(in) :: tol
integer, intent(in) :: max_its
integer, intent(out) :: its

Calls

proc~~homotopy~~CallsGraph proc~homotopy homotopy none~wrap wrap proc~homotopy->none~wrap proc~newton newton proc~homotopy->proc~newton proc~solve_system solve_system proc~newton->proc~solve_system none~dgesv dgesv proc~solve_system->none~dgesv

Variables

Type Visibility Attributes Name Initial
real(kind=pr), public :: F(size(X))
real(kind=pr), public :: G(size(X))
real(kind=pr), public :: H(size(X))
real(kind=pr), public :: X0(size(X))
real(kind=pr), public :: dF(size(X),size(X))
real(kind=pr), public :: dG(size(X),size(X))
real(kind=pr), public :: dH(size(X),size(X))
real(kind=pr), public :: dX(size(x))
real(kind=pr), public :: dt
integer, public :: i
integer, public :: info
integer, public :: nt
real(kind=pr), public :: t

Subroutines

subroutine fun(n, x, fvec, iflag)

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: n
real(kind=pr), intent(in) :: x(n)
real(kind=pr), intent(out) :: fvec(n)
integer, intent(inout) :: iflag

subroutine wrap(X, H, dH)

Arguments

Type IntentOptional Attributes Name
real(kind=pr), intent(in) :: X(:)
real(kind=pr), intent(out) :: H(:)
real(kind=pr), intent(out) :: dH(:,:)