MPFR v4.0.2
[mpfr.git] / tests / tlog1p.c
1 /* Test file for mpfr_log1p.
2
3 Copyright 2001-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 #ifdef CHECK_EXTERNAL
26 static int
27 test_log1p (mpfr_ptr a, mpfr_srcptr b, mpfr_rnd_t rnd_mode)
28 {
29   int res;
30   int ok = rnd_mode == MPFR_RNDN && mpfr_number_p (b) && mpfr_get_prec (a)>=53;
31   if (ok)
32     {
33       mpfr_print_raw (b);
34     }
35   res = mpfr_log1p (a, b, rnd_mode);
36   if (ok)
37     {
38       printf (" ");
39       mpfr_print_raw (a);
40       printf ("\n");
41     }
42   return res;
43 }
44 #else
45 #define test_log1p mpfr_log1p
46 #endif
47
48 #define TEST_FUNCTION test_log1p
49 #define TEST_RANDOM_EMAX 80
50 #include "tgeneric.c"
51
52 static void
53 special (void)
54 {
55   mpfr_t x;
56   int inex;
57
58   mpfr_init (x);
59
60   mpfr_set_nan (x);
61   mpfr_clear_flags ();
62   inex = test_log1p (x, x, MPFR_RNDN);
63   MPFR_ASSERTN (mpfr_nan_p (x) && inex == 0);
64   MPFR_ASSERTN (__gmpfr_flags == MPFR_FLAGS_NAN);
65
66   mpfr_set_inf (x, -1);
67   mpfr_clear_flags ();
68   inex = test_log1p (x, x, MPFR_RNDN);
69   MPFR_ASSERTN (mpfr_nan_p (x) && inex == 0);
70   MPFR_ASSERTN (__gmpfr_flags == MPFR_FLAGS_NAN);
71
72   mpfr_set_inf (x, 1);
73   mpfr_clear_flags ();
74   inex = test_log1p (x, x, MPFR_RNDN);
75   MPFR_ASSERTN (mpfr_inf_p (x) && mpfr_sgn (x) > 0 && inex == 0);
76   MPFR_ASSERTN (__gmpfr_flags == 0);
77
78   mpfr_set_ui (x, 0, MPFR_RNDN);
79   mpfr_clear_flags ();
80   inex = test_log1p (x, x, MPFR_RNDN);
81   MPFR_ASSERTN (mpfr_cmp_ui (x, 0) == 0 && MPFR_IS_POS (x) && inex == 0);
82   MPFR_ASSERTN (__gmpfr_flags == 0);
83   mpfr_neg (x, x, MPFR_RNDN);
84   mpfr_clear_flags ();
85   inex = test_log1p (x, x, MPFR_RNDN);
86   MPFR_ASSERTN (mpfr_cmp_ui (x, 0) == 0 && MPFR_IS_NEG (x) && inex == 0);
87   MPFR_ASSERTN (__gmpfr_flags == 0);
88
89   mpfr_set_si (x, -1, MPFR_RNDN);
90   mpfr_clear_flags ();
91   inex = test_log1p (x, x, MPFR_RNDN);
92   MPFR_ASSERTN (mpfr_inf_p (x) && mpfr_sgn (x) < 0 && inex == 0);
93   MPFR_ASSERTN (__gmpfr_flags == MPFR_FLAGS_DIVBY0);
94
95   mpfr_set_si (x, -2, MPFR_RNDN);
96   mpfr_clear_flags ();
97   inex = test_log1p (x, x, MPFR_RNDN);
98   MPFR_ASSERTN (mpfr_nan_p (x) && inex == 0);
99   MPFR_ASSERTN (__gmpfr_flags == MPFR_FLAGS_NAN);
100
101   mpfr_clear (x);
102 }
103
104 static void
105 other (void)
106 {
107   mpfr_t x, y;
108
109   /* Bug reported by Guillaume Melquiond on 2006-08-14. */
110   mpfr_init2 (x, 53);
111   mpfr_set_str (x, "-1.5e4f72873ed9a@-100", 16, MPFR_RNDN);
112   mpfr_init2 (y, 57);
113   mpfr_log1p (y, x, MPFR_RNDU);
114   if (mpfr_cmp (x, y) != 0)
115     {
116       printf ("Error in tlog1p for x = ");
117       mpfr_out_str (stdout, 16, 0, x, MPFR_RNDN);
118       printf (", rnd = MPFR_RNDU\nExpected ");
119       mpfr_out_str (stdout, 16, 15, x, MPFR_RNDN);
120       printf ("\nGot      ");
121       mpfr_out_str (stdout, 16, 15, y, MPFR_RNDN);
122       printf ("\n");
123       exit (1);
124     }
125
126   mpfr_clear (y);
127   mpfr_clear (x);
128   return;
129 }
130
131 int
132 main (int argc, char *argv[])
133 {
134   tests_start_mpfr ();
135
136   special ();
137   other ();
138
139   test_generic (MPFR_PREC_MIN, 100, 50);
140
141   data_check ("data/log1p", mpfr_log1p, "mpfr_log1p");
142   bad_cases (mpfr_log1p, mpfr_expm1, "mpfr_log1p", 256, -64, 40,
143              4, 128, 800, 40);
144
145   tests_end_mpfr ();
146   return 0;
147 }