pest_newton_FRdFR Function

private function pest_newton_FRdFR(sf, x) result(A_FRdFR)

Function for 1D newton root search for PEST

Type Bound

t_newton_Root1D_FdF_pest

Arguments

Type IntentOptional Attributes Name
class(t_newton_Root1D_FdF_pest), intent(in) :: sf
real(kind=wp), intent(in) :: x

Return Value real(kind=wp), (2)


Calls

proc~~pest_newton_frdfr~~CallsGraph proc~pest_newton_frdfr t_newton_Root1D_FdF_pest%pest_newton_FRdFR proc~fbase_evaldof_x t_fBase%fBase_evalDOF_x proc~pest_newton_frdfr->proc~fbase_evaldof_x proc~fbase_eval t_fBase%fBase_eval proc~fbase_evaldof_x->proc~fbase_eval proc~fbase_eval_xn t_fBase%fBase_eval_xn proc~fbase_eval->proc~fbase_eval_xn

Source Code

FUNCTION pest_newton_FRdFR(sf, x) RESULT(A_FRdFR)
  !uses current zeta where newton is called, and A from subroutine above
  CLASS(t_newton_Root1D_FdF_pest), INTENT(IN) :: sf
  REAL(wp), INTENT(IN) :: x ! theta_iter
  REAL(wp) :: A_FRdFR(2) !output function and derivative
  !---------------------------------------------------
  A_FRdFR(1)=x      + sf%LA_fbase_in%evalDOF_x((/x,sf%zeta/),         0,sf%LA_in) !theta_iter+lambda = thet* (rhs)
  A_FRdFR(2)=1.0_wp + sf%LA_fbase_in%evalDOF_x((/x,sf%zeta/),DERIV_THET,sf%LA_in) !1+dlambda/dtheta
END FUNCTION pest_newton_FRdFR