gauss_freq Function

private function gauss_freq(x)

Uses

  • proc~~gauss_freq~~UsesGraph proc~gauss_freq gauss_freq module~w90_constants w90_constants proc~gauss_freq->module~w90_constants

gauss_freq(x) = (1+erf(x/sqrt(2)))/2 = erfc(-x/sqrt(2))/2 - See comments in erf

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: x

Return Value real(kind=dp)


Calls

proc~~gauss_freq~~CallsGraph proc~gauss_freq gauss_freq proc~qe_erfc qe_erfc proc~gauss_freq->proc~qe_erfc proc~qe_erf qe_erf proc~qe_erfc->proc~qe_erf proc~qe_erf->proc~qe_erfc

Called by

proc~~gauss_freq~~CalledByGraph proc~gauss_freq gauss_freq proc~utility_wgauss utility_wgauss proc~utility_wgauss->proc~gauss_freq

Source Code

  function gauss_freq(x)
    !---------------------------------------------------------------------
    !
    !! gauss_freq(x) = (1+erf(x/sqrt(2)))/2 = erfc(-x/sqrt(2))/2
    !! - See comments in erf
    !
    use w90_constants, only: dp

    implicit none
    real(kind=dp), intent(in) :: x
    real(kind=dp)            :: gauss_freq
    real(kind=dp), parameter :: c = 0.7071067811865475_dp
    !        ( c= sqrt(1/2) )
    !
    gauss_freq = 0.5_dp*qe_erfc(-x*c)
    !
    return
  end function gauss_freq