DPCHIP Sign-Testing Routine
Returns:
The object is to do this without multiplying ARG1*ARG2, to avoid possible over/underflow problems.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in) | :: | arg1 | |||
real(kind=wp), | intent(in) | :: | arg2 |
pure function dpchst (arg1, arg2) result(s) real(wp),intent(in) :: arg1 real(wp),intent(in) :: arg2 real(wp) :: s if ((arg1==zero) .or. (arg2==zero)) then s = zero else s = sign(one,arg1) * sign(one,arg2) end if end function dpchst