xref: /netbsd-src/external/lgpl3/mpc/dist/tests/pow_ui.dat (revision 90a8ff2142ed565a73c3c0859f0b1e7d216aeb7b)
1# Data file for mpc_pow_ui.
2#
3# Copyright (C) 2010, 2011, 2012 INRIA
4#
5# This file is part of GNU MPC.
6#
7# GNU MPC is free software; you can redistribute it and/or modify it under
8# the terms of the GNU Lesser General Public License as published by the
9# Free Software Foundation; either version 3 of the License, or (at your
10#o ption) any later version.
11#
12# GNU MPC is distributed in the hope that it will be useful, but WITHOUT ANY
13# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
14# FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
15# more details.
16#
17# You should have received a copy of the GNU Lesser General Public License
18# along with this program. If not, see http://www.gnu.org/licenses/ .
19#
20# The line format respects the parameter order in function prototype as
21# follow:
22#
23# INEX_RE INEX_IM PREC_ROP_RE  ROP_RE  PREC_ROP_IM  ROP_IM  PREC_OP1_RE  OP1_RE  PREC_OP1_IM  OP1_IM  OP2  RND_RE  RND_IM
24#
25# For further details, see add_fr.dat.
26
27# special cases, copied from pow.dat
280 0 53   +1 53    0    53  nan 53   +0  +0 N N
290 0 53  nan 53  nan    53  nan 53   +0  +1 N N
300 0 53  inf 53  nan    53 +inf 53   +0  +1 N N
310 0 53 +inf 53  nan    53 +inf 53   +1  +1 N N
320 0 53 +inf 53  nan    53 +inf 53   -1  +1 N N
330 0 53 +inf 53  nan    53 -inf 53   +0  +1 N N
340 0 53 +inf 53  nan    53 -inf 53   +1  +1 N N
350 0 53 +inf 53  nan    53 -inf 53   -1  +1 N N
36
370 0 53   +1 53   +0    53   +0 53   +0  +0 N N
380 0 53   +0 53   +0    53   +0 53   +0  +1 N N
39
400 0 53   +1 53   +0    53   +0 53   +1  +0 N N
410 0 53   +1 53   -0    53   +0 53   +1  +0 N D
420 0 53   +1 53   +0    53   -0 53   +1  +0 N N
430 0 53   +1 53   +0    53   -1 53   +0  +0 N N
440 0 53   +1 53   -0    53   -1 53   -0  +0 N N
450 0 53   +1 53   -0    53   -0 53   -1  +0 N N
460 0 53   +1 53   -0    53   +0 53   -1  +0 N N
47
480 0 53   +1 53   +0    53 +inf 53   +2  +0 N N
490 0 53   +1 53   +0    53 +inf 53   -0  +0 N N
500 0 53   +1 53   +0    53   +2 53 +inf  +0 N N
510 0 53   +1 53   +0    53   +2 53   +0  +0 N N
520 0 53   +1 53   +0    53   +0 53   +2  +0 N N
530 0 53   +1 53   +0    53   +0 53 +inf  +0 N N
540 0 53   +1 53   +0    53   -0 53   +2  +0 N N
550 0 53   +1 53   +0    53   -0 53 +inf  +0 N N
560 0 53   +1 53   +0    53   -5 53 +inf  +0 N N
570 0 53   +1 53   +0    53   -2 53   +0  +0 N N
580 0 53   +1 53   +0    53 -inf 53   +0  +0 N N
590 0 53   +1 53   +0    53 -inf 53   +3  +0 N N
600 0 53   +1 53   -0    53 +0.5 53 -0.5  +0 N N
610 0 53   +1 53   -0    53 +0.5 53   +0  +0 N N
620 0 53   +1 53   -0    53 +0.5 53   -0  +0 N N
630 0 53   +1 53   -0    53 -0.5 53   -0  +0 N N
640 0 53   +1 53   -0    53   +0 53 -0.5  +0 N N
650 0 53   +1 53   -0    53   -0 53 -0.5  +0 N N
660 0 53   +9 53   +0    53   +3 53   +0  +2 N N
670 0 53   +1 53   +0    53   +1 53   +0  +4 N N
680 0 53   +1 53   -0    53   +1 53   -0  +4 N N
690 0 53 0.25 53   -0    53 +0.5 53   -0  +2 N N
70
710 0 53    1 53    0    53   +2 53   -1  +0 N N
720 0 53    1 53    0    53   -2 53   -1  +0 N N
730 0 53    1 53    0    53   -2 53   -0  +0 N N
740 0 53    1 53    0    53 +0.5 53 +0.5  +0 N N
750 0 53    1 53    0    53 -0.5 53 +0.5  +0 N N
760 0 53    1 53    0    53 -0.5 53   +0  +0 N N
770 0 53    1 53    0    53   +0 53 +0.5  +0 N N
780 0 53    1 53    0    53   -0 53 +0.5  +0 N N
790 0 53    1 53    0    53   -0 53   -4  +0 N N
800 0 53    1 53    0    53   +0 53   -4  +0 N N
810 0 53    1 53    0    53   -1 53   -0  +0 N N
820 0 53    1 53    0    53   -1 53   +0  +0 N N
83
840 0 53    4 53    0    53   +2 53   -0  +2 N N
850 0 53    1 53    0    53   +1 53   +0  +2 N N
860 0 53    1 53    0    53   +1 53   -0  +2 N N
87
88# overflow: (a+a*I)^n is real for n divisible by 4
89? 0 53 +inf 53    0    53 1e100000000  53 1e100000000   1000000000 N N
90# underflow
91? ? 53    0 53    0    53 1e-100000000 53 1e-100000000  1000000000 N N
92# cannot round after one loop
93? ? 420  -0x1.c3fb41a71665f9a144927e70cbc2dc899e9e30880c0b5aa924ad8a538b4cd06e503f38bdbb7cfcfded29f7504fe0c91ecd4230984@-187 420 -0xc.82a09ac98133eb05b2643c98eb1c8e1a1609e75f682b14098176abd6c8b4b3c6c72dadaf8929f9bd87f8c78d03361bacb9fb13140@-292  420 0x1.cf13ce58adc4e639fd1c3063ffc9291433647999951bc04ba6797ec4de0335336ad0a28df18573d3b6322ebab662c08eadaed4a8e@-8 420 0x3.cf71d602ca6f754ebd6af522154f3ee1c46da0a52deb1f60016fca4b1e0b4b447b752169e837bb1866aa3734850cd158a7e3ca33c@-9  24 N N
94# exact result
950 0  5 -9 5 46  5 3 5 2  3 D U
96
97# special cases with exponents 1 and 2
980 0 53  3 53  4  53 3 53 4 1 N N
990 0 53 -7 53 24  53 3 53 4 2 N N
100
101# coverage test
102+ - 2 1 2 1 100 0x8ac78f4e6aec091921cbdc891p-99 100 0x4a5f2972ea662e1cc0caebb81p-100 3 N N
103