delta_for_zero_omega Function

private function delta_for_zero_omega(omega_) result(delta_)

Checks if the input value equals 0; used in the calculation of off-diagonal elements of the centrifugal matrix; see Eq. 5 in "Coupling Matrix" section

Arguments

Type IntentOptional Attributes Name
integer(kind=int32), intent(in) :: omega_

input value which is to be compared with 0

Return Value real(kind=dp)

(output) 1 if omega_ = 0, 0 otherwise


Called by

proc~~delta_for_zero_omega~~CalledByGraph proc~delta_for_zero_omega delta_for_zero_omega proc~calculate_centrifugal_matrix calculate_centrifugal_matrix proc~calculate_centrifugal_matrix->proc~delta_for_zero_omega proc~initial_setup initial_setup proc~initial_setup->proc~calculate_centrifugal_matrix proc~numerov numerov proc~numerov->proc~initial_setup program~scattering SCATTERING program~scattering->proc~numerov

Contents

Source Code


Source Code

      function delta_for_zero_omega(omega_) result(delta_)
         !! Checks if the input value equals 0; used in the calculation
         !! of off-diagonal elements of the centrifugal matrix; see
         !! Eq. 5 in "Coupling Matrix" section
         !---------------------------------------------------------------------!
         integer(int32), intent(in) :: omega_
            !! input value which is to be compared with 0
         real(dp) :: delta_
            !! (output) 1 if omega_ = 0, 0 otherwise
         !---------------------------------------------------------------------!
         if (omega_ == 0) then
           delta_ = 1.0_dp
         else
           delta_ = 0.0_dp
         endif
         !---------------------------------------------------------------------!
      end function delta_for_zero_omega