A true signum: -|a| if x<0, +|a| if x>0, 0 if x==0 (MOM6 signum).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | a | |||
| real(kind=wp), | intent(in) | :: | x |
pure function redi_signum(a, x) result(s) !$acc routine seq real(wp), intent(in) :: a, x real(wp) :: s s = sign(a, x) if (x == 0.0_wp) s = 0.0_wp end function redi_signum