rho2_derivative Function

public pure function rho2_derivative(rho, deriv) result(rho2_prime)

evaluate the n-th derivative of (rho^2) with respect to rho ~sqrt(magnetic flux).

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: rho

rho position rho ~sqrt(magnetic flux)
integer, intent(in) :: deriv

derivative in rho

Return Value real(kind=wp)

n-th derivative of rho2 with respect to rho ~sqrt(magnetic flux).


Calls

proc~~rho2_derivative~~CallsGraph proc~rho2_derivative rho2_derivative proc~poly_derivative_prefactor poly_derivative_prefactor proc~rho2_derivative->proc~poly_derivative_prefactor

Called by

proc~~rho2_derivative~~CalledByGraph proc~rho2_derivative rho2_derivative proc~rprofile_drho2 c_rProfile%rProfile_drho2 proc~rprofile_drho2->proc~rho2_derivative proc~rprofile_drho3 c_rProfile%rProfile_drho3 proc~rprofile_drho3->proc~rho2_derivative proc~rprofile_drho4 c_rProfile%rProfile_drho4 proc~rprofile_drho4->proc~rho2_derivative proc~rprofile_eval_at_rho c_rProfile%rProfile_eval_at_rho proc~rprofile_eval_at_rho->proc~rho2_derivative proc~rprofile_eval_at_rho->proc~rprofile_drho2 proc~rprofile_eval_at_rho->proc~rprofile_drho3 proc~rprofile_eval_at_rho->proc~rprofile_drho4 proc~evaluate_profile evaluate_profile proc~evaluate_profile->proc~rprofile_eval_at_rho proc~init_boozer init_boozer proc~init_boozer->proc~rprofile_eval_at_rho proc~init_la_from_solution Init_LA_from_Solution proc~init_la_from_solution->proc~rprofile_eval_at_rho proc~initprofilesgp InitProfilesGP proc~initprofilesgp->proc~rprofile_eval_at_rho proc~initvmec InitVMEC proc~initvmec->proc~rprofile_eval_at_rho proc~vmec1d_visu VMEC1D_visu proc~vmec1d_visu->proc~rprofile_eval_at_rho proc~writestatetoascii WriteStateToASCII proc~writestatetoascii->proc~rprofile_eval_at_rho interface~writestate WriteState interface~writestate->proc~writestatetoascii proc~analyze Analyze proc~analyze->proc~vmec1d_visu proc~initsolutionmhd3d t_functional_mhd3d%InitSolutionMHD3D proc~initsolutionmhd3d->proc~init_la_from_solution proc~initsolutionmhd3d->proc~initprofilesgp proc~initsolutionmhd3d->interface~writestate program~gvec_post GVEC_POST program~gvec_post->proc~initprofilesgp proc~initsolution~2 InitSolution proc~initsolution~2->proc~initsolutionmhd3d proc~minimizemhd3d_descent t_minimizer_mhd3d%MinimizeMHD3D_descent proc~minimizemhd3d_descent->interface~writestate proc~rungvec rungvec proc~rungvec->proc~initsolutionmhd3d proc~minimizemhd3d t_functional_mhd3d%MinimizeMHD3D proc~rungvec->proc~minimizemhd3d proc~minimizemhd3d->proc~minimizemhd3d_descent proc~start_rungvec start_rungvec proc~start_rungvec->proc~rungvec program~gvec GVEC program~gvec->proc~rungvec proc~minimize minimize proc~minimize->proc~minimizemhd3d

Source Code

  PURE FUNCTION rho2_derivative(rho,deriv) RESULT(rho2_prime)
    ! INPUT VARIABLES -------------------------!
    REAL(wp), INTENT(IN) :: rho                !! rho position rho ~sqrt(magnetic flux)
    INTEGER, INTENT(IN)  :: deriv  !! derivative in rho
    ! OUTPUT VARIABLES -------------------------!
    REAL(wp) :: rho2_prime     !!n-th derivative of rho2 with respect to rho ~sqrt(magnetic flux).
    ! CODE --------------------------------------------------------------------------------------------------------------------------!
    IF (deriv>2) THEN
      rho2_prime = 0.0_wp
  ELSE
      rho2_prime = poly_derivative_prefactor(2,deriv)*rho**(2-deriv)
  END IF
  END FUNCTION rho2_derivative