Blame tests/mpq/io.c

Packit 5c3484
/* Test conversion and I/O using mpq_out_str and mpq_inp_str.
Packit 5c3484
Packit 5c3484
Copyright 1993, 1994, 1996, 2000, 2001, 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 "config.h"
Packit 5c3484
Packit 5c3484
#include <stdio.h>
Packit 5c3484
#include <stdlib.h>
Packit 5c3484
#if HAVE_UNISTD_H
Packit 5c3484
#include <unistd.h>		/* for unlink */
Packit 5c3484
#endif
Packit 5c3484
Packit 5c3484
#include "gmp.h"
Packit 5c3484
#include "gmp-impl.h"
Packit 5c3484
#include "tests.h"
Packit 5c3484
Packit 5c3484
#define FILENAME  "io.tmp"
Packit 5c3484
Packit 5c3484
void
Packit 5c3484
debug_mp (mpq_t x, int base)
Packit 5c3484
{
Packit 5c3484
  mpq_out_str (stdout, base, x); fputc ('\n', stdout);
Packit 5c3484
}
Packit 5c3484
Packit 5c3484
int
Packit 5c3484
main (int argc, char **argv)
Packit 5c3484
{
Packit 5c3484
  mpq_t  op1, op2;
Packit 5c3484
  mp_size_t size;
Packit 5c3484
  int i;
Packit 5c3484
  int reps = 10000;
Packit 5c3484
  FILE *fp;
Packit 5c3484
  int base;
Packit 5c3484
  gmp_randstate_ptr rands;
Packit 5c3484
  mpz_t bs;
Packit 5c3484
  unsigned long bsi, size_range;
Packit 5c3484
  size_t nread;
Packit 5c3484
Packit 5c3484
  tests_start ();
Packit 5c3484
  rands = RANDS;
Packit 5c3484
Packit 5c3484
  mpz_init (bs);
Packit 5c3484
Packit 5c3484
  if (argc == 2)
Packit 5c3484
    reps = atoi (argv[1]);
Packit 5c3484
Packit 5c3484
  mpq_init (op1);
Packit 5c3484
  mpq_init (op2);
Packit 5c3484
Packit 5c3484
  fp = fopen (FILENAME, "w+");
Packit 5c3484
Packit 5c3484
  for (i = 0; i < reps; i++)
Packit 5c3484
    {
Packit 5c3484
      mpz_urandomb (bs, rands, 32);
Packit 5c3484
      size_range = mpz_get_ui (bs) % 10 + 2;
Packit 5c3484
Packit 5c3484
      mpz_urandomb (bs, rands, size_range);
Packit 5c3484
      size = mpz_get_ui (bs);
Packit 5c3484
      mpz_errandomb (mpq_numref(op1), rands, 512L);
Packit 5c3484
      mpz_errandomb_nonzero (mpq_denref(op1), rands, 512L);
Packit 5c3484
      mpq_canonicalize (op1);
Packit 5c3484
Packit 5c3484
      mpz_urandomb (bs, rands, 1);
Packit 5c3484
      bsi = mpz_get_ui (bs);
Packit 5c3484
      if ((bsi & 1) != 0)
Packit 5c3484
	mpq_neg (op1, op1);
Packit 5c3484
Packit 5c3484
      mpz_urandomb (bs, rands, 16);
Packit 5c3484
      bsi = mpz_get_ui (bs);
Packit 5c3484
      base = bsi % 36 + 1;
Packit 5c3484
      if (base == 1)
Packit 5c3484
	base = 0;
Packit 5c3484
Packit 5c3484
      rewind (fp);
Packit 5c3484
      if (mpq_out_str (fp, base, op1) == 0
Packit 5c3484
	  || putc (' ', fp) == EOF
Packit 5c3484
	  || fflush (fp) != 0)
Packit 5c3484
	{
Packit 5c3484
	  printf ("mpq_out_str write error\n");
Packit 5c3484
	  abort ();
Packit 5c3484
	}
Packit 5c3484
Packit 5c3484
      rewind (fp);
Packit 5c3484
      nread = mpq_inp_str (op2, fp, base);
Packit 5c3484
      if (nread == 0)
Packit 5c3484
	{
Packit 5c3484
	  if (ferror (fp))
Packit 5c3484
	    printf ("mpq_inp_str stream read error\n");
Packit 5c3484
	  else
Packit 5c3484
	    printf ("mpq_inp_str data conversion error\n");
Packit 5c3484
	  abort ();
Packit 5c3484
	}
Packit 5c3484
Packit 5c3484
      if (nread != ftell(fp))
Packit 5c3484
	{
Packit 5c3484
	  printf ("mpq_inp_str nread doesn't match ftell\n");
Packit 5c3484
	  printf ("  nread  %lu\n", (unsigned long) nread);
Packit 5c3484
	  printf ("  ftell  %ld\n", ftell(fp));
Packit 5c3484
	  abort ();
Packit 5c3484
	}
Packit 5c3484
Packit 5c3484
      if (mpq_cmp (op1, op2))
Packit 5c3484
	{
Packit 5c3484
	  printf ("ERROR\n");
Packit 5c3484
	  printf ("op1  = "); debug_mp (op1, -16);
Packit 5c3484
	  printf ("op2  = "); debug_mp (op2, -16);
Packit 5c3484
	  printf ("base = %d\n", base);
Packit 5c3484
	  abort ();
Packit 5c3484
	}
Packit 5c3484
    }
Packit 5c3484
Packit 5c3484
  fclose (fp);
Packit 5c3484
Packit 5c3484
  unlink (FILENAME);
Packit 5c3484
Packit 5c3484
  mpz_clear (bs);
Packit 5c3484
  mpq_clear (op1);
Packit 5c3484
  mpq_clear (op2);
Packit 5c3484
Packit 5c3484
  tests_end ();
Packit 5c3484
  exit (0);
Packit 5c3484
}