central_diff Subroutine

private subroutine central_diff(f, x, h, dfdx)

Compute the first derivative using a 2-point central difference [-h,h].

Arguments

Type IntentOptional Attributes Name
procedure(func) :: f
complex(kind=wp), intent(in) :: x
real(kind=wp), intent(in) :: h
real(kind=wp), intent(out) :: dfdx

Called by

proc~~central_diff~~CalledByGraph proc~central_diff complex_step_module::central_diff proc~complex_step_test complex_step_module::complex_step_test proc~complex_step_test->proc~central_diff

Source Code

    subroutine central_diff(f,x,h,dfdx)

    implicit none

    procedure(func)        :: f
    complex(wp),intent(in) :: x
    real(wp),intent(in)    :: h
    real(wp),intent(out)   :: dfdx

    dfdx = (f(x+h) - f(x-h)) / (2.0_wp*h)

    end subroutine central_diff