1  /* Test file for mpfr_log1p.
       2  
       3  Copyright 2001-2023 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, 9,
     143               4, 128, 800, 40);
     144  
     145    tests_end_mpfr ();
     146    return 0;
     147  }