MPFR v4.0.2
[mpfr.git] / tests / tli2.c
1 /* mpfr_tli2 -- test file for dilogarithm function
2
3 Copyright 2007-2019 Free Software Foundation, Inc.
4 Contributed by the AriC and Caramba projects, INRIA.
5
6 This file is part of the GNU MPFR Library.
7
8 The GNU MPFR Library is free software; you can redistribute it and/or modify
9 it under the terms of the GNU Lesser General Public License as published by
10 the Free Software Foundation; either version 3 of the License, or (at your
11 option) any later version.
12
13 The GNU MPFR Library is distributed in the hope that it will be useful, but
14 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
15 or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
16 License for more details.
17
18 You should have received a copy of the GNU Lesser General Public License
19 along with the GNU MPFR Library; see the file COPYING.LESSER.  If not, see
20 https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc.,
21 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA. */
22
23 #include "mpfr-test.h"
24
25 #define TEST_FUNCTION mpfr_li2
26 #include "tgeneric.c"
27
28 static void
29 special (void)
30 {
31   mpfr_t x, y;
32   mpfr_init (x);
33   mpfr_init (y);
34
35   mpfr_set_nan (x);
36   mpfr_li2 (y, x, MPFR_RNDN);
37   if (!mpfr_nan_p (y))
38     {
39       printf ("Error for li2(NaN)\n");
40       exit (1);
41     }
42
43   mpfr_set_inf (x, -1);
44   mpfr_li2 (y, x, MPFR_RNDN);
45   if (!MPFR_IS_INF (y) || MPFR_IS_POS (y))
46     {
47       printf ("Error for li2(-Inf)\n");
48       exit (1);
49     }
50
51   mpfr_set_inf (x, 1);
52   mpfr_li2 (y, x, MPFR_RNDN);
53   if (!MPFR_IS_INF (y) || MPFR_IS_POS (y))
54     {
55       printf ("Error for li2(+Inf)\n");
56       exit (1);
57     }
58
59   mpfr_set_ui (x, 0, MPFR_RNDN);
60   mpfr_li2 (y, x, MPFR_RNDN);
61   if (!MPFR_IS_ZERO (y) || MPFR_IS_NEG (y))
62     {
63       printf ("Error for li2(+0)\n");
64       exit (1);
65     }
66
67   mpfr_set_ui (x, 0, MPFR_RNDN);
68   mpfr_neg (x, x, MPFR_RNDN);
69   mpfr_li2 (y, x, MPFR_RNDN);
70   if (!MPFR_IS_ZERO (y) || MPFR_IS_POS (y))
71     {
72       printf ("Error for li2(-0)\n");
73       exit (1);
74     }
75
76   mpfr_clear (x);
77   mpfr_clear (y);
78 }
79
80 static void
81 normal (void)
82 {
83   int inexact;
84   mpfr_t x, y;
85   mpfr_init (x);
86   mpfr_init (y);
87
88   /* x1 = 2^-3 */
89   mpfr_set_str (x, "1p-3", 2, MPFR_RNDD);
90   mpfr_li2 (x, x, MPFR_RNDN);
91   if (mpfr_cmp_str (x, "0x1087a7a9e42141p-55", 16, MPFR_RNDN) != 0)
92     {
93       printf ("Error for li2(x1)\n");
94       exit (1);
95     }
96
97   /* check MPFR_FAST_COMPUTE_IF_SMALL_INPUT */
98   mpfr_set_prec (x, 2);
99   mpfr_set_prec (y, 20);
100   mpfr_set_ui_2exp (x, 1, -21, MPFR_RNDN);
101   mpfr_li2 (y, x, MPFR_RNDN);
102   MPFR_ASSERTN(mpfr_cmp (y, x) == 0);
103
104   mpfr_set_si_2exp (x, -1, -21, MPFR_RNDN);
105   mpfr_li2 (y, x, MPFR_RNDN);
106   MPFR_ASSERTN(mpfr_cmp (y, x) == 0);
107
108   /* worst case */
109   /* x2 = 0x7F18EA6537E00E983196CDDC6EFAC57Fp-129
110      Li2(x2) = 2^-2 + 2^-6 + 2^-120 */
111   mpfr_set_prec (x, 128);
112   mpfr_set_str (x, "7F18EA6537E00E983196CDDC6EFAC57Fp-129", 16, MPFR_RNDN);
113
114   /* round to nearest mode and 4 bits of precision,
115      it should be rounded to 2^-2 + 2^-5 and */
116   mpfr_set_prec (y, 4);
117   inexact = mpfr_li2 (y, x, MPFR_RNDN);
118   if (inexact != 1 || mpfr_cmp_str (y, "0.1001p-1", 2, MPFR_RNDN) != 0)
119     {
120       printf ("Error for li2(x2, RNDN)\n");
121       exit (1);
122     }
123
124   /* round toward zero mode and 5 bits of precision,
125      it should be rounded to 2^-2 + 2^-6 */
126   mpfr_set_prec (y, 5);
127   inexact = mpfr_li2 (y, x, MPFR_RNDZ);
128   if (inexact != -1 || mpfr_cmp_str (y, "0.10001p-1", 2, MPFR_RNDN) != 0)
129     {
130       printf ("Error for li2(x2, RNDZ)\n");
131       exit (1);
132     }
133
134   /* round away from zero mode and 5 bits of precision,
135      it should be rounded to 2^-2 + 2^-5 */
136   inexact = mpfr_li2 (y, x, MPFR_RNDU);
137   if (inexact != 1 || mpfr_cmp_str (y, "0.10010p-1", 2, MPFR_RNDN) != 0)
138     {
139       printf ("Error for li2(x2, RNDU)\n");
140       exit (1);
141     }
142
143   mpfr_clear (x);
144   mpfr_clear (y);
145 }
146
147 static void
148 bug20091013 (void)
149 {
150   mpfr_t x, y;
151   int inex;
152
153   mpfr_init2 (x, 17);
154   mpfr_init2 (y, 2);
155   mpfr_set_str_binary (x, "0.10000000000000000E-16");
156   inex = mpfr_li2 (y, x, MPFR_RNDN);
157   if (mpfr_cmp_ui_2exp (y, 1, -17) != 0)
158     {
159       printf ("Error in bug20091013()\n");
160       printf ("expected 2^(-17)\n");
161       printf ("got      ");
162       mpfr_dump (y);
163       exit (1);
164     }
165   if (inex >= 0)
166     {
167       printf ("Error in bug20091013()\n");
168       printf ("expected negative ternary value, got %d\n", inex);
169       exit (1);
170     }
171   mpfr_clear (x);
172   mpfr_clear (y);
173 }
174
175 int
176 main (int argc, char *argv[])
177 {
178   tests_start_mpfr ();
179
180   bug20091013 ();
181
182   special ();
183
184   normal ();
185
186   test_generic (MPFR_PREC_MIN, 100, 2);
187
188   data_check ("data/li2", mpfr_li2, "mpfr_li2");
189
190   tests_end_mpfr ();
191   return 0;
192 }