diff options
| author | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2008-12-09 19:29:03 +0000 |
|---|---|---|
| committer | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2008-12-09 19:29:03 +0000 |
| commit | 1100e2290732e9fd296aafdd0fecd5e0b6b2d411 (patch) | |
| tree | 44f50a7a47acf840a62d3cfdbccf6a6a0c9be2d5 /tests/sqrt.dat | |
| parent | fb12b07ad5c9b4c181c54455f438eaf448c3ca03 (diff) | |
| download | mpc-1100e2290732e9fd296aafdd0fecd5e0b6b2d411.tar.gz | |
[sqrt.c] introduced ternary return value for special, purely real
and purely imaginary arguments
[sqrt.dat] added ternary return value for such arguments
[read_data.c] output '?' instead of 'x' for unchecked return value
git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@402 211d60ee-9f03-0410-a15a-8952a2c7a4e4
Diffstat (limited to 'tests/sqrt.dat')
| -rw-r--r-- | tests/sqrt.dat | 76 |
1 files changed, 38 insertions, 38 deletions
diff --git a/tests/sqrt.dat b/tests/sqrt.dat index c1ba750..066819d 100644 --- a/tests/sqrt.dat +++ b/tests/sqrt.dat @@ -1,6 +1,6 @@ # Data file for mpc_sqrt. # -# Copyright (C) 2008 Philippe Th\'eveny +# Copyright (C) 2008 Philippe Th\'eveny, Andreas Enge # # This file is part of the MPC Library. # @@ -73,50 +73,50 @@ 0 0 53 +inf 53 +inf 53 nan 53 +inf N N 0 0 53 nan 53 nan 53 nan 53 nan N N -# pure real argument +# purely real argument # sqrt(x +i*0) = sqrt(x) +i*0, when x>0 # sqrt(x -i*0) = sqrt(x) -i*0, when x>0 # sqrt(x +i*0) = +0 +i*sqrt(-x) +i*0, when x<0 # sqrt(x -i*0) = +0 -i*sqrt(-x) +i*0, when x<0 -? ? 53 0x16a09e667f3bcdp-52 53 +0 53 2 17 +0 N N -? ? 53 +0 53 0x16a09e667f3bcdp-52 54 -2 16 +0 Z N -? ? 53 0x16a09e667f3bcdp-52 53 -0 55 2 15 -0 U N -? ? 53 +0 53 -0x16a09e667f3bcdp-52 56 -2 14 -0 D N -? ? 53 0x5a827999fcef30p-54 53 +0 57 2 13 +0 Z Z -? ? 53 +0 53 0x5a827999fcef30p-54 58 -2 12 +0 U Z -? ? 53 0x5a827999fcef30p-54 53 -0 59 2 11 -0 D Z -? ? 53 +0 53 -0x5a827999fcef30p-54 60 -2 10 -0 N Z -? ? 53 0x16a09e667f3bcdp-52 53 +0 61 2 9 +0 U U -? ? 53 +0 53 0x16a09e667f3bcdp-52 62 -2 8 +0 D U -? ? 53 0x16a09e667f3bcdp-52 53 -0 63 2 7 -0 N U -? ? 53 +0 53 -0x5a827999fcef30p-54 64 -2 6 -0 Z U -? ? 53 0x5a827999fcef30p-54 53 +0 65 2 5 +0 D D -? ? 53 +0 53 0x5a827999fcef30p-54 66 -2 4 +0 N D -? ? 53 0x5a827999fcef30p-54 53 -0 67 2 3 -0 Z D -? ? 53 +0 53 -0x16a09e667f3bcdp-52 68 -2 2 -0 U D ++ 0 53 0x16a09e667f3bcdp-52 53 +0 53 2 17 +0 N N +0 + 53 +0 53 0x16a09e667f3bcdp-52 54 -2 16 +0 Z N ++ 0 53 0x16a09e667f3bcdp-52 53 -0 55 2 15 -0 U N +0 - 53 +0 53 -0x16a09e667f3bcdp-52 56 -2 14 -0 D N +- 0 53 0x5a827999fcef30p-54 53 +0 57 2 13 +0 Z Z +0 - 53 +0 53 0x5a827999fcef30p-54 58 -2 12 +0 U Z +- 0 53 0x5a827999fcef30p-54 53 -0 59 2 11 -0 D Z +0 + 53 +0 53 -0x5a827999fcef30p-54 60 -2 10 -0 N Z ++ 0 53 0x16a09e667f3bcdp-52 53 +0 61 2 9 +0 U U +0 + 53 +0 53 0x16a09e667f3bcdp-52 62 -2 8 +0 D U ++ 0 53 0x16a09e667f3bcdp-52 53 -0 63 2 7 -0 N U +0 + 53 +0 53 -0x5a827999fcef30p-54 64 -2 6 -0 Z U +- 0 53 0x5a827999fcef30p-54 53 +0 65 2 5 +0 D D +0 - 53 +0 53 0x5a827999fcef30p-54 66 -2 4 +0 N D +- 0 53 0x5a827999fcef30p-54 53 -0 67 2 3 -0 Z D +0 - 53 +0 53 -0x16a09e667f3bcdp-52 68 -2 2 -0 U D -# pure imaginary argument -# sqrt(+/-0 +i*y) = sqrt(y/2) * (1 +i), when y >0 +# purely imaginary argument +# sqrt(+/-0 +i*y) = sqrt(y/2) * (1 +i), when y >0 # sqrt(+/-0 +i*y) = sqrt(-y/2) * (1 -i), when y < 0 -? ? 53 0x16a09e667f3bcdp-52 53 0x16a09e667f3bcdp-52 53 +0 53 4 N N -? ? 53 0x5a827999fcef30p-54 53 0x16a09e667f3bcdp-52 51 -0 54 4 Z N -? ? 53 0x16a09e667f3bcdp-52 53 -0x16a09e667f3bcdp-52 49 +0 55 -4 U N -? ? 53 0x5a827999fcef30p-54 53 -0x16a09e667f3bcdp-52 47 -0 56 -4 D N -? ? 53 0x5a827999fcef30p-54 53 0x5a827999fcef30p-54 45 +0 57 4 Z Z -? ? 53 0x16a09e667f3bcdp-52 53 0x5a827999fcef30p-54 43 -0 58 4 U Z -? ? 53 0x5a827999fcef30p-54 53 -0x5a827999fcef30p-54 41 +0 59 -4 D Z -? ? 53 0x16a09e667f3bcdp-52 53 -0x5a827999fcef30p-54 39 -0 60 -4 N Z -? ? 53 0x16a09e667f3bcdp-52 53 0x16a09e667f3bcdp-52 37 +0 61 4 U U -? ? 53 0x5a827999fcef30p-54 53 0x16a09e667f3bcdp-52 35 -0 62 4 D U -? ? 53 0x16a09e667f3bcdp-52 53 -0x5a827999fcef30p-54 33 +0 63 -4 N U -? ? 53 0x5a827999fcef30p-54 53 -0x5a827999fcef30p-54 31 -0 64 -4 Z U -? ? 53 0x5a827999fcef30p-54 53 0x5a827999fcef30p-54 29 +0 65 4 D D -? ? 53 0x16a09e667f3bcdp-52 53 0x5a827999fcef30p-54 27 -0 66 4 N D -? ? 53 0x5a827999fcef30p-54 53 -0x16a09e667f3bcdp-52 25 +0 67 -4 Z D -? ? 53 0x16a09e667f3bcdp-52 53 -0x16a09e667f3bcdp-52 23 -0 68 -4 U D ++ + 53 0x16a09e667f3bcdp-52 53 0x16a09e667f3bcdp-52 53 +0 53 4 N N +- + 53 0x5a827999fcef30p-54 53 0x16a09e667f3bcdp-52 51 -0 54 4 Z N ++ - 53 0x16a09e667f3bcdp-52 53 -0x16a09e667f3bcdp-52 49 +0 55 -4 U N +- - 53 0x5a827999fcef30p-54 53 -0x16a09e667f3bcdp-52 47 -0 56 -4 D N +- - 53 0x5a827999fcef30p-54 53 0x5a827999fcef30p-54 45 +0 57 4 Z Z ++ - 53 0x16a09e667f3bcdp-52 53 0x5a827999fcef30p-54 43 -0 58 4 U Z +- + 53 0x5a827999fcef30p-54 53 -0x5a827999fcef30p-54 41 +0 59 -4 D Z ++ + 53 0x16a09e667f3bcdp-52 53 -0x5a827999fcef30p-54 39 -0 60 -4 N Z ++ + 53 0x16a09e667f3bcdp-52 53 0x16a09e667f3bcdp-52 37 +0 61 4 U U +- + 53 0x5a827999fcef30p-54 53 0x16a09e667f3bcdp-52 35 -0 62 4 D U ++ + 53 0x16a09e667f3bcdp-52 53 -0x5a827999fcef30p-54 33 +0 63 -4 N U +- + 53 0x5a827999fcef30p-54 53 -0x5a827999fcef30p-54 31 -0 64 -4 Z U +- - 53 0x5a827999fcef30p-54 53 0x5a827999fcef30p-54 29 +0 65 4 D D ++ - 53 0x16a09e667f3bcdp-52 53 0x5a827999fcef30p-54 27 -0 66 4 N D +- - 53 0x5a827999fcef30p-54 53 -0x16a09e667f3bcdp-52 25 +0 67 -4 Z D ++ - 53 0x16a09e667f3bcdp-52 53 -0x16a09e667f3bcdp-52 23 -0 68 -4 U D # bugs fixed in r160 2008-07-15 ? ? 19 0b11101001001001001100p+39 19 -0b1010110101100111011p-236 19 0b1.101010001010100000p+117 19 -0b1.001110111101100001p-158 N Z -? ? 2 0b11p+100 2 -0b11p+100 2 -0 2 -0b11p+203 N Z -? ? 2 0 2 -0b10p+117 2 -0b11p+235 2 -0 N Z +- + 2 0b11p+100 2 -0b11p+100 2 -0 2 -0b11p+203 N Z +0 + 2 0 2 -0b10p+117 2 -0b11p+235 2 -0 N Z |