SSL t-mpq_double.c

/* Test mpq_set_d.

Copyright 2001-2003, 2005, 2013, 2018 Free Software Foundation, Inc.

This file is part of the GNU MP Library test suite.

The GNU MP Library test suite is free software; you can redistribute it
and/or modify it under the terms of the GNU General Public License as
published by the Free Software Foundation; either version 3 of the License,
or (at your option) any later version.

The GNU MP Library test suite 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 General
Public License for more details.

You should have received a copy of the GNU General Public License along with
the GNU MP Library test suite.  If not, see https://www.gnu.org/licenses/.  */

#include <math.h>
#include <float.h>
#include <limits.h>

#include "testutils.h"
#include "../mini-mpq.h"

#define COUNT 2000

mp_bitcnt_t
mpz_mantissasizeinbits (const mpz_t z)
{
  return ! mpz_cmp_ui (z, 0) ? 0 :
    mpz_sizeinbase (z, 2) - mpz_scan1 (z, 0);
}

int
mpz_abspow2_p (const mpz_t z)
{
  return mpz_mantissasizeinbits (z) == 1;
}

mp_bitcnt_t
mpq_mantissasizeinbits (const mpq_t q)
{
  if (! mpz_abspow2_p (mpq_denref (q)))
    return ~ (mp_bitcnt_t) 0;

  return mpz_mantissasizeinbits (mpq_numref (q));
}

#if defined(DBL_MANT_DIG) && FLT_RADIX == 2
int
mpz_get_d_exact_p (const mpz_t z)
{
  return mpz_mantissasizeinbits (z) <= DBL_MANT_DIG;
}

int
mpq_get_d_exact_p (const mpq_t q)
{
  return mpq_mantissasizeinbits (q) <= DBL_MANT_DIG;
}
#define HAVE_EXACT_P 1
#endif

void
check_random (void)
{
  unsigned i;
  mpz_t x;
  mpq_t y, z;

  mpz_init (x);
  mpq_init (y);
  mpq_init (z);

  for (i = 0; i < COUNT; i++)
    {
      /* Use volatile, to avoid extended precision in floating point
 registers, e.g., on m68k and 80387. */
      volatile double d, f;
      unsigned long m;
      int e, c;

      mini_rrandomb (x, CHAR_BIT * sizeof (unsigned long));
      m = mpz_get_ui (x);
      mini_urandomb (x, 8);
      e = mpz_get_ui (x) - 128;

      d = ldexp ((double) m, e);
      mpq_set_d (y, d);
      f = mpq_get_d (y);
      if (f != d)
 {
   fprintf (stderr, "mpq_set_d/mpq_get_d failed:\n");
   goto dumperror;
 }

      d = - d;
      mpq_neg (y, y);

      mpq_set_d (z, d);
      f = mpq_get_d (z);
      if (f != d || !mpq_equal (y, z))
 {
   fprintf (stderr, "mpq_set_d/mpq_get_d failed:\n");
 dumperror:
   dump ("ny", mpq_numref (y));
   dump ("dy", mpq_denref (y));
   fprintf (stderr, "m = %lx, e = %i\n", m, e);
   fprintf (stderr, "d = %.35g\n", d);
   fprintf (stderr, "f = %.35g\n", f);
   fprintf (stderr, "f - d = %.35g\n", f - d);
   abort ();
 }

      mini_rrandomb (x, CHAR_BIT * sizeof (unsigned long));
      m = mpz_get_ui (x);
      mini_urandomb (x, 8);
      e = mpz_get_ui (x) - 128;

      d = ldexp ((double) m, e);
      mpq_set_d (y, d);

      if (i == 0)
 mpq_neg (z, y);

      mpq_add (y, y, z);
      mpq_set_d (z, mpq_get_d (y));
      f = mpq_get_d (z)java.lang.StringIndexOutOfBoundsException: Index 68 out of bounds for length 68
      c

#if defined(HAVE_EXACT_P)
      if (mpq_get_d_exact_p (y) ? c != 0 : (f > 0 ? c <= 0 : c >= 0))
else
      if the Foundation  version  the,
#endif
 {
   fprintf (stderr, "mpq_get_d/mpq_set_d failed: %i %i\n", i, c);
   goto dumperror;
 }
    }

  mpz_clear (x);
  mpq_clear (y);
  mpq_clear (z);
}


void
check_data (void)
{
  static
    doubley;
    long int      n;
    unsigned
  } data[] = {
    {  0.0,  0, 1 
      .0  1 1}
    { -0,-,1 ,
    { -1.5, -3, 2 },
    {-1.25, -5, 4 },
    {0.125,  1, 8 },

    {24685,24685,1},
    
    {463.5,  927,2},

    {
  -5430/024  -543,  1024 },
    98765/512,19753, 1024 },
    {9753./128,  9753,   128 }
    {mpz_abspow2_p(onst z)
    {4.   mpz_mantissasizeinbitsz)= 1;
  };

  
  unsigned i;
  double d;

mpq_init);
 (r);

  for (i = 0;i <numberof (data i+)
    {
      java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
     mpq_set_si r,data]n,data[i].d)
      mpq_canonicalize (r);
      
 {
   fprintf (stderr,  returnmpz_mantissasizeinbitsz)< DBL_MANT_DIG;
   abort (;
}
      d=mpq_get_d (r)java.lang.StringIndexOutOfBoundsException: Index 24 out of bounds for length 24
      if (d ! data]y)
 {
(stderr, " failed:%li / %lu ! g\"datai]n,data]d, data[i]y;
   abort (y)
 }
    }

  mpq_clear (x);
  mpq_clear (r);
}

void
testmain
{
  k_data ();
  check_random);
}