Calculate and for the given elliptic modulus
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=real32), | intent(in) | :: | k |
elliptic modulus |
||
| real(kind=real32), | intent(out) | :: | pw02_k |
|
||
| real(kind=real32), | 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(real32), intent(in) :: k !! elliptic modulus \( k \) real(real32), intent(out) :: pw02_k !! \( { k }^{ 2 } \) real(real32), intent(out) :: comp_k !! \( { k }^{ \prime } \) pw02_k = k * k comp_k = sqrt(1.0_real32 - pw02_k) end subroutine evaluate_modulus