1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
|
# Data file for mpc_atanh.
#
# Copyright (C) 2009 Philippe Th\'eveny
#
# This file is part of the MPC Library.
#
# The MPC Library is free software; you can redistribute it and/or modify
# it under the terms of the GNU Lesser General Public License as published by
# the Free Software Foundation; either version 2.1 of the License, or (at your
# option) any later version.
#
# The MPC Library is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
# or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
# License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with the MPC Library; see the file COPYING.LIB. If not, write to
# the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
# MA 02111-1307, USA.
#
# The line format respects the parameter order in function prototype as
# follow:
#
# 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
#
# where op = op_re + i * op_im, rop = rop_re + i * rop_im,
# rop_re is ROP_RE rounded to nearest to the precision of PREC_ROP_RE
# rop_im is ROP_IM rounded to nearest to the precision of PREC_ROP_IM
# op_re is OP_RE rounded to nearest to the precision of PREC_OP_RE
# op_im is OP_IM rounded to nearest to the precision of PREC_OP_IM
# ROP_RE is checked against Re(atan op) rounded to the precision PREC_ROP_RE
# in the direction RND_RE
# ROP_IM is checked against Im(atan op) rounded to the precision PREC_ROP_IM
# in the direction RND_IM
# INEX_RE is the ternary value for the real part with the following notation:
# "?" ternary value not checked
# "+" if ROP_RE is greater than the exact mathematical result
# "0" if ROP_RE is exactly the mathematical result
# "-" if ROP_RE is less than the exact mathematical result
# (m.m. INEX_IM)
# rounding modes notation:
# "N" is rounding to nearest
# "Z" is rounding towards zero
# "U" is rounding towards plus infinity
# "D" is rounding towards minus infinity
# Use prefixes "0b" for values in base two, "0x" for values in base sixteen,
# no prefix for value in base ten.
# In all bases, "nan" is NaN, "inf" is infinity;
# The sign of the result is checked with "+inf", "-inf", "-0", or "+0".
# special values (following ISO C99 standard)
0 + 53 -0 53 -0x1921FB54442D18p-52 53 -inf 53 -inf N N
0 + 53 -0 53 -0x1921FB54442D18p-52 53 -inf 53 -1 N N
0 + 53 -0 53 -0x1921FB54442D18p-52 53 -inf 53 -0 N N
0 - 53 -0 53 +0x1921FB54442D18p-52 53 -inf 53 +0 N N
0 - 53 -0 53 +0x1921FB54442D18p-52 53 -inf 53 +1 N N
0 - 53 -0 53 +0x1921FB54442D18p-52 53 -inf 53 +inf N N
0 0 53 -0 53 nan 53 -inf 53 nan N N
0 + 53 -0 53 -0x1921FB54442D18p-52 53 -6 53 -inf N N
0 - 53 -0 53 +0x1921FB54442D18p-52 53 -6 53 +inf N N
0 0 53 nan 53 nan 53 -6 53 nan N N
0 + 53 -0 53 -0x1921FB54442D18p-52 53 -0 53 -inf N N
0 0 53 -0 53 -0 53 -0 53 -0 N N
0 0 53 -0 53 +0 53 -0 53 +0 N N
0 - 53 -0 53 +0x1921FB54442D18p-52 53 -0 53 +inf N N
0 0 53 -0 53 nan 53 -0 53 nan N N
0 + 53 +0 53 -0x1921FB54442D18p-52 53 +0 53 -inf N N
0 0 53 +0 53 -0 53 +0 53 -0 N N
0 0 53 +0 53 +0 53 +0 53 +0 N N
0 - 53 +0 53 +0x1921FB54442D18p-52 53 +0 53 +inf N N
0 0 53 +0 53 nan 53 +0 53 nan N N
0 + 53 +0 53 -0x1921FB54442D18p-52 53 +6 53 -inf N N
0 - 53 +0 53 +0x1921FB54442D18p-52 53 +6 53 +inf N N
0 0 53 nan 53 nan 53 +6 53 nan N N
0 + 53 +0 53 -0x1921FB54442D18p-52 53 +inf 53 -inf N N
0 + 53 +0 53 -0x1921FB54442D18p-52 53 +inf 53 -1 N N
0 + 53 +0 53 -0x1921FB54442D18p-52 53 +inf 53 -0 N N
0 - 53 +0 53 +0x1921FB54442D18p-52 53 +inf 53 +0 N N
0 - 53 +0 53 +0x1921FB54442D18p-52 53 +inf 53 +1 N N
0 - 53 +0 53 +0x1921FB54442D18p-52 53 +inf 53 +inf N N
0 0 53 +0 53 nan 53 +inf 53 nan N N
0 + 53 0 53 -0x1921FB54442D18p-52 53 nan 53 -inf N N
0 0 53 nan 53 nan 53 nan 53 -1 N N
0 0 53 nan 53 nan 53 nan 53 -0 N N
0 0 53 nan 53 nan 53 nan 53 +0 N N
0 0 53 nan 53 nan 53 nan 53 +1 N N
0 - 53 0 53 +0x1921FB54442D18p-52 53 nan 53 +inf N N
0 0 53 nan 53 nan 53 nan 53 nan N N
# pure real argument
- + 53 -0x1E27076E2AF2E6p-57 53 -0x1921FB54442D18p-52 53 -17 53 -0 N N
- - 53 -0x1E27076E2AF2E6p-57 53 +0x1921FB54442D18p-52 53 -17 53 +0 N N
+ + 53 +0x1E27076E2AF2E6p-57 53 -0x1921FB54442D18p-52 53 +17 53 -0 N N
+ - 53 +0x1E27076E2AF2E6p-57 53 +0x1921FB54442D18p-52 53 +17 53 +0 N N
+ 0 53 -0x1F2272AE325A57p-53 53 -0 53 -.75 53 -0 N N
+ 0 53 -0x1F2272AE325A57p-53 53 +0 53 -.75 53 +0 N N
- 0 53 +0x1F2272AE325A57p-53 53 -0 53 +.75 53 -0 N N
- 0 53 +0x1F2272AE325A57p-53 53 +0 53 +.75 53 +0 N N
# pure imaginary argument
0 - 53 -0 53 -0x167D8863BC99BDp-52 53 -0 53 -6 N N
0 - 53 +0 53 -0x167D8863BC99BDp-52 53 +0 53 -6 N N
0 + 53 -0 53 +0x167D8863BC99BDp-52 53 -0 53 +6 N N
0 + 53 +0 53 +0x167D8863BC99BDp-52 53 +0 53 +6 N N
0 + 53 -0 53 -0x1F5B75F92C80DDp-55 53 -0 53 -.25 N N
0 + 53 +0 53 -0x1F5B75F92C80DDp-55 53 +0 53 -.25 N N
0 - 53 -0 53 +0x1F5B75F92C80DDp-55 53 -0 53 +.25 N N
0 - 53 +0 53 +0x1F5B75F92C80DDp-55 53 +0 53 +.25 N N
# IEEE-754 double precision
|
