Calculate and for the given elliptic modulus
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=real64), | intent(in) | :: | k | elliptic modulus | ||
| real(kind=real64), | intent(out) | :: | pw02_k | 
 | ||
| real(kind=real64), | intent(out) | :: | comp_k | 
 | 
elemental subroutine evaluate_modulus(k, pw02_k, comp_k) !! Calculate \( { k }^{ 2 } \) !! and \( { k }^{ \prime } := \sqrt{ 1 - { k }^{ 2 } } \) !! for the given elliptic modulus \( k \) real(real64), intent(in) :: k !! elliptic modulus \( k \) real(real64), intent(out) :: pw02_k !! \( { k }^{ 2 } \) real(real64), intent(out) :: comp_k !! \( { k }^{ \prime } \) pw02_k = k * k comp_k = sqrt(1.0_real64 - pw02_k) end subroutine evaluate_modulus