Blame mini-gmp/tests/t-div_2exp.c

Packit 5c3484
/*
Packit 5c3484
Packit 5c3484
Copyright 2012, Free Software Foundation, Inc.
Packit 5c3484
Packit 5c3484
This file is part of the GNU MP Library test suite.
Packit 5c3484
Packit 5c3484
The GNU MP Library test suite is free software; you can redistribute it
Packit 5c3484
and/or modify it under the terms of the GNU General Public License as
Packit 5c3484
published by the Free Software Foundation; either version 3 of the License,
Packit 5c3484
or (at your option) any later version.
Packit 5c3484
Packit 5c3484
The GNU MP Library test suite is distributed in the hope that it will be
Packit 5c3484
useful, but WITHOUT ANY WARRANTY; without even the implied warranty of
Packit 5c3484
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General
Packit 5c3484
Public License for more details.
Packit 5c3484
Packit 5c3484
You should have received a copy of the GNU General Public License along with
Packit 5c3484
the GNU MP Library test suite.  If not, see https://www.gnu.org/licenses/.  */
Packit 5c3484
Packit 5c3484
#include <assert.h>
Packit 5c3484
#include <stdlib.h>
Packit 5c3484
#include <stdio.h>
Packit 5c3484
Packit 5c3484
#include "testutils.h"
Packit 5c3484
Packit 5c3484
#define MAXBITS 400
Packit 5c3484
#define COUNT 10000
Packit 5c3484
Packit 5c3484
typedef void div_func (mpz_t, const mpz_t, mp_bitcnt_t);
Packit 5c3484
Packit 5c3484
void
Packit 5c3484
testmain (int argc, char **argv)
Packit 5c3484
{
Packit 5c3484
  unsigned i;
Packit 5c3484
  mpz_t a, res, ref;
Packit 5c3484
  mp_bitcnt_t b;
Packit 5c3484
Packit 5c3484
  mpz_init (a);
Packit 5c3484
  mpz_init (res);
Packit 5c3484
  mpz_init (ref);
Packit 5c3484
Packit 5c3484
  for (i = 0; i < COUNT; i++)
Packit 5c3484
    {
Packit 5c3484
      unsigned j;
Packit 5c3484
      for (j = 0; j < 6; j++)
Packit 5c3484
	{
Packit 5c3484
	  static const enum hex_random_op ops[6] =
Packit 5c3484
	    {
Packit 5c3484
	      OP_CDIV_Q_2, OP_CDIV_R_2,
Packit 5c3484
	      OP_FDIV_Q_2, OP_FDIV_R_2,
Packit 5c3484
	      OP_TDIV_Q_2, OP_TDIV_R_2
Packit 5c3484
	    };
Packit 5c3484
	  static const char *name[6] =
Packit 5c3484
	    {
Packit 5c3484
	      "cdiv_q", "cdiv_r",
Packit 5c3484
	      "fdiv_q", "fdiv_r",
Packit 5c3484
	      "tdiv_q", "tdiv_r"
Packit 5c3484
	    };
Packit 5c3484
	  static div_func * const div [6] =
Packit 5c3484
	    {
Packit 5c3484
	      mpz_cdiv_q_2exp, mpz_cdiv_r_2exp,
Packit 5c3484
	      mpz_fdiv_q_2exp, mpz_fdiv_r_2exp,
Packit 5c3484
	      mpz_tdiv_q_2exp, mpz_tdiv_r_2exp
Packit 5c3484
	    };
Packit 5c3484
Packit 5c3484
	  mini_random_bit_op (ops[j], MAXBITS, a, &b, ref);
Packit 5c3484
	  div[j] (res, a, b);
Packit 5c3484
	  if (mpz_cmp (ref, res))
Packit 5c3484
	    {
Packit 5c3484
	      fprintf (stderr, "mpz_%s_2exp failed:\n", name[j]);
Packit 5c3484
	      dump ("a", a);
Packit 5c3484
	      fprintf (stderr, "b: %lu\n", b);
Packit 5c3484
	      dump ("r", res);
Packit 5c3484
	      dump ("ref", ref);
Packit 5c3484
	      abort ();
Packit 5c3484
	    }
Packit 5c3484
	}
Packit 5c3484
    }
Packit 5c3484
  mpz_clear (a);
Packit 5c3484
  mpz_clear (res);
Packit 5c3484
  mpz_clear (ref);
Packit 5c3484
}