third_party/gmp/tests/mpf/t-set.c

/* Test mpf_set, mpf_init_set.

Copyright 2004, 2012 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 <stdio.h>
#include <stdlib.h>
#include "gmp-impl.h"
#include "tests.h"

void
check_reuse (void)
{
  /* Try mpf_set(f,f) when f is bigger than prec.  In the past this had
     resulted in an MPN_COPY with invalid operand overlap. */
  mpf_t  f;
  mp_size_t      limbs = 20;
  unsigned long  bits = limbs * GMP_NUMB_BITS;
  mpf_init2 (f, bits);
  refmpf_fill (f, limbs, GMP_NUMB_MAX);
  mpf_set_prec_raw (f, bits / 2);
  mpf_set (f, f);
  MPF_CHECK_FORMAT (f);
  mpf_set_prec_raw (f, bits);
  mpf_clear (f);
}

void
check_random (long reps)
{
  unsigned long test;
  gmp_randstate_ptr rands;
  mpf_t a, b;
  mpz_t z;
  int precbits;

#define PRECBITS 10

  rands = RANDS;

  mpz_init (z);
  mpf_init2 (a, 1 << PRECBITS);

  for (test = 0; test < reps; test++)
    {
      mpz_urandomb (z, rands, PRECBITS + 1);
      precbits = mpz_get_ui (z) + 1;
      mpz_urandomb (z, rands, precbits);
      mpz_setbit (z, precbits  - 1);	/* make sure msb is set */
      mpf_set_z (a, z);
      if (precbits & 1)
	mpf_neg (a, a);
      mpz_urandomb (z, rands, PRECBITS);
      mpf_div_2exp (a, a, mpz_get_ui (z) + 1);
      mpz_urandomb (z, rands, PRECBITS);
      precbits -= mpz_get_ui (z);
      if (precbits <= 0)
	precbits = 1 - precbits;
      mpf_set_default_prec (precbits);

      mpf_init_set (b, a);
      MPF_CHECK_FORMAT (b);
      if (!mpf_eq (a, b, precbits))
	{
	  printf ("mpf_init_set wrong.\n");
	  abort();
	}

      mpf_set_ui (b, 0);
      mpf_set (b, a);
      MPF_CHECK_FORMAT (b);
      if (!mpf_eq (a, b, precbits))
	{
	  printf ("mpf_set wrong.\n");
	  abort();
	}

      mpf_clear (b);
    }

  mpf_clear (a);
  mpz_clear (z);
}

int
main (int argc, char *argv[])
{
  long reps = 10000;

  tests_start ();
  TESTS_REPS (reps, argv, argc);

  check_reuse ();
  check_random (reps);

  tests_end ();
  exit (0);
}