test_zcomplexabs.f - Man Page

program zabs
zabs tests the robustness and precision of the intrinsic ABS for double complex

zabs tests the robustness and precision of the intrinsic ABS for double complex

Author

Weslley S. Pereira, University of Colorado Denver, U.S.

 Real values for test:
 (1) x = 2**m, where m = MINEXPONENT-DIGITS, ..., MINEXPONENT-1. Stop on the first success.
     Mind that not all platforms might implement subnormal numbers.
 (2) x = 2**m, where m = MINEXPONENT, ..., 0. Stop on the first success.
 (3) x = OV, where OV is the overflow threshold. OV^2 overflows but the norm is OV.
 (4) x = 2**m, where m = MAXEXPONENT-1, ..., 1. Stop on the first success.

 Tests:
 (a) y = x + 0 * I, |y| = x
 (b) y = 0 + x * I, |y| = x
 (c) y = (3/4)*x + x * I, |y| = (5/4)*x whenever (3/4)*x and (5/4)*x can be exactly stored
 (d) y = (1/2)*x + (1/2)*x * I, |y| = (1/2)*x*sqrt(2) whenever (1/2)*x can be exactly stored

 Special cases:

 (i) Inf propagation
    (1) y = Inf + 0 * I, |y| is Inf.
    (2) y =-Inf + 0 * I, |y| is Inf.
    (3) y = 0 + Inf * I, |y| is Inf.
    (4) y = 0 - Inf * I, |y| is Inf.
    (5) y = Inf + Inf * I, |y| is Inf.

 (n) NaN propagation
    (1) y = NaN + 0 * I, |y| is NaN.
    (2) y = 0 + NaN * I, |y| is NaN.
    (3) y = NaN + NaN * I, |y| is NaN.

Definition at line 44 of file test_zcomplexabs.f.