1 # Data file for mpc_sqrt.
3 # Copyright (C) 2008 Philippe Th\'eveny, Andreas Enge
5 # This file is part of the MPC Library.
7 # The MPC Library is free software; you can redistribute it and/or modify
8 # it under the terms of the GNU Lesser General Public License as published by
9 # the Free Software Foundation; either version 2.1 of the License, or (at your
10 # option) any later version.
12 # The MPC Library is distributed in the hope that it will be useful, but
13 # WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
14 # or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
15 # License for more details.
17 # You should have received a copy of the GNU Lesser General Public License
18 # along with the MPC Library; see the file COPYING.LIB. If not, write to
19 # the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
22 # The line format respects the parameter order in function prototype as
25 # INEX_RE INEX_IM PREC_ROP_RE ROP_RE PREC_ROP_IM ROP_IM PREC_OP_RE OP_RE PREC_OP_IM OP_IM RND_RE RND_IM
27 # see sin.dat for precisions
29 # special values (following ISO C99 standard)
30 0 0 53 +inf 53 -inf 53 -inf 53 -inf N N
31 0 0 53 +0 53 -inf 53 -inf 53 -1 N N
32 0 0 53 +0 53 -inf 53 -inf 53 -0 N N
33 0 0 53 +0 53 +inf 53 -inf 53 +0 N N
34 0 0 53 +0 53 +inf 53 -inf 53 +1 N N
35 0 0 53 +inf 53 +inf 53 -inf 53 +inf N N
36 0 0 53 nan 53 inf 53 -inf 53 nan N N
37 0 0 53 +inf 53 -inf 53 -1 53 -inf N N
38 0 0 53 +0 53 -1 53 -1 53 -0 N N
39 0 0 53 +0 53 +1 53 -1 53 +0 N N
40 0 0 53 +inf 53 +inf 53 -1 53 +inf N N
41 0 0 53 nan 53 nan 53 -1 53 nan N N
42 0 0 53 +inf 53 -inf 53 -0 53 -inf N N
43 0 0 53 +1 53 -1 53 -0 53 -2 N N
44 0 0 53 +0 53 -0 53 -0 53 -0 N N
45 0 0 53 +0 53 +0 53 -0 53 +0 N N
46 0 0 53 +1 53 +1 53 -0 53 +2 N N
47 0 0 53 +inf 53 +inf 53 -0 53 +inf N N
48 0 0 53 nan 53 nan 53 -0 53 nan N N
49 0 0 53 +inf 53 -inf 53 +0 53 -inf N N
50 0 0 53 +1 53 -1 53 +0 53 -2 N N
51 0 0 53 +0 53 -0 53 +0 53 -0 N N
52 0 0 53 +0 53 +0 53 +0 53 +0 N N
53 0 0 53 +1 53 +1 53 +0 53 +2 N N
54 0 0 53 +inf 53 +inf 53 +0 53 +inf N N
55 0 0 53 nan 53 nan 53 +0 53 nan N N
56 0 0 53 +inf 53 -inf 53 +1 53 -inf N N
57 0 0 53 +1 53 -0 53 +1 53 -0 N N
58 0 0 53 +1 53 +0 53 +1 53 +0 N N
59 0 0 53 +inf 53 +inf 53 +1 53 +inf N N
60 0 0 53 nan 53 nan 53 +1 53 nan N N
61 0 0 53 +inf 53 -inf 53 +inf 53 -inf N N
62 0 0 53 +inf 53 -0 53 +inf 53 -1 N N
63 0 0 53 +inf 53 -0 53 +inf 53 -0 N N
64 0 0 53 +inf 53 +0 53 +inf 53 +0 N N
65 0 0 53 +inf 53 +0 53 +inf 53 +1 N N
66 0 0 53 +inf 53 +inf 53 +inf 53 +inf N N
67 0 0 53 +inf 53 nan 53 +inf 53 nan N N
68 0 0 53 +inf 53 -inf 53 nan 53 -inf N N
69 0 0 53 nan 53 nan 53 nan 53 -1 N N
70 0 0 53 nan 53 nan 53 nan 53 -0 N N
71 0 0 53 nan 53 nan 53 nan 53 +0 N N
72 0 0 53 nan 53 nan 53 nan 53 +1 N N
73 0 0 53 +inf 53 +inf 53 nan 53 +inf N N
74 0 0 53 nan 53 nan 53 nan 53 nan N N
76 # purely real argument
77 # sqrt(x +i*0) = sqrt(x) +i*0, when x>0
78 # sqrt(x -i*0) = sqrt(x) -i*0, when x>0
79 # sqrt(x +i*0) = +0 +i*sqrt(-x) +i*0, when x<0
80 # sqrt(x -i*0) = +0 -i*sqrt(-x) +i*0, when x<0
81 + 0 53 0x16a09e667f3bcdp-52 53 +0 53 2 17 +0 N N
82 0 + 53 +0 53 0x16a09e667f3bcdp-52 54 -2 16 +0 Z N
83 + 0 53 0x16a09e667f3bcdp-52 53 -0 55 2 15 -0 U N
84 0 - 53 +0 53 -0x16a09e667f3bcdp-52 56 -2 14 -0 D N
85 - 0 53 0x5a827999fcef30p-54 53 +0 57 2 13 +0 Z Z
86 0 - 53 +0 53 0x5a827999fcef30p-54 58 -2 12 +0 U Z
87 - 0 53 0x5a827999fcef30p-54 53 -0 59 2 11 -0 D Z
88 0 + 53 +0 53 -0x5a827999fcef30p-54 60 -2 10 -0 N Z
89 + 0 53 0x16a09e667f3bcdp-52 53 +0 61 2 9 +0 U U
90 0 + 53 +0 53 0x16a09e667f3bcdp-52 62 -2 8 +0 D U
91 + 0 53 0x16a09e667f3bcdp-52 53 -0 63 2 7 -0 N U
92 0 + 53 +0 53 -0x5a827999fcef30p-54 64 -2 6 -0 Z U
93 - 0 53 0x5a827999fcef30p-54 53 +0 65 2 5 +0 D D
94 0 - 53 +0 53 0x5a827999fcef30p-54 66 -2 4 +0 N D
95 - 0 53 0x5a827999fcef30p-54 53 -0 67 2 3 -0 Z D
96 0 - 53 +0 53 -0x16a09e667f3bcdp-52 68 -2 2 -0 U D
98 # purely imaginary argument
99 # sqrt(+/-0 +i*y) = sqrt(y/2) * (1 +i), when y >0
100 # sqrt(+/-0 +i*y) = sqrt(-y/2) * (1 -i), when y < 0
101 + + 53 0x16a09e667f3bcdp-52 53 0x16a09e667f3bcdp-52 53 +0 53 4 N N
102 - + 53 0x5a827999fcef30p-54 53 0x16a09e667f3bcdp-52 51 -0 54 4 Z N
103 + - 53 0x16a09e667f3bcdp-52 53 -0x16a09e667f3bcdp-52 49 +0 55 -4 U N
104 - - 53 0x5a827999fcef30p-54 53 -0x16a09e667f3bcdp-52 47 -0 56 -4 D N
105 - - 53 0x5a827999fcef30p-54 53 0x5a827999fcef30p-54 45 +0 57 4 Z Z
106 + - 53 0x16a09e667f3bcdp-52 53 0x5a827999fcef30p-54 43 -0 58 4 U Z
107 - + 53 0x5a827999fcef30p-54 53 -0x5a827999fcef30p-54 41 +0 59 -4 D Z
108 + + 53 0x16a09e667f3bcdp-52 53 -0x5a827999fcef30p-54 39 -0 60 -4 N Z
109 + + 53 0x16a09e667f3bcdp-52 53 0x16a09e667f3bcdp-52 37 +0 61 4 U U
110 - + 53 0x5a827999fcef30p-54 53 0x16a09e667f3bcdp-52 35 -0 62 4 D U
111 + + 53 0x16a09e667f3bcdp-52 53 -0x5a827999fcef30p-54 33 +0 63 -4 N U
112 - + 53 0x5a827999fcef30p-54 53 -0x5a827999fcef30p-54 31 -0 64 -4 Z U
113 - - 53 0x5a827999fcef30p-54 53 0x5a827999fcef30p-54 29 +0 65 4 D D
114 + - 53 0x16a09e667f3bcdp-52 53 0x5a827999fcef30p-54 27 -0 66 4 N D
115 - - 53 0x5a827999fcef30p-54 53 -0x16a09e667f3bcdp-52 25 +0 67 -4 Z D
116 + - 53 0x16a09e667f3bcdp-52 53 -0x16a09e667f3bcdp-52 23 -0 68 -4 U D
118 # bugs fixed in r160 2008-07-15
119 - + 19 0b11101001001001001100p+39 19 -0b1010110101100111011p-236 19 0b1.101010001010100000p+117 19 -0b1.001110111101100001p-158 N Z
120 - + 2 0b11p+100 2 -0b11p+100 2 -0 2 -0b11p+203 N Z
121 0 + 2 0 2 -0b10p+117 2 -0b11p+235 2 -0 N Z
projects / platform / upstream / mpc.git / blob