summaryrefslogtreecommitdiff
path: root/tests/tset_q.c
blob: 7d6528f7bfd763bc74e67885e715e74b842e0df1 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
/* Test file for mpfr_set_q.

Copyright 2000, 2001, 2002, 2003 Free Software Foundation, Inc.

This file is part of the MPFR Library.

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

The MPFR Library 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 Lesser General Public
License for more details.

You should have received a copy of the GNU Lesser General Public License
along with the MPFR Library; see the file COPYING.LIB.  If not, write to
the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
MA 02111-1307, USA. */

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include "gmp.h"
#include "mpfr.h"
#include "mpfr-test.h"

static void
check (long int n, long int d, mp_rnd_t rnd, double y)
{
  mpq_t q;
  mpfr_t x, t;
  double z;
  int inexact, compare;

  mpfr_init2 (x, 53);
  mpfr_init2 (t, mpfr_get_prec (x) + mp_bits_per_limb);
  mpq_init (q);
  mpq_set_si (q, n, d);
  inexact = mpfr_set_q (x, q, rnd);
  z = mpfr_get_d1 (x);

  /* check values */
  if (y != z)
    {
      fprintf (stderr, "Error for q=%ld/%lu and rnd=%s\n", n, d, 
               mpfr_print_rnd_mode (rnd));
      fprintf (stderr, "correct result is %1.20e, mpfr_set_q gives %1.20e\n",
               y, z);
      exit (1);
    }

  /* check inexact flag */
  if (mpfr_mul_ui (t, x, (d < 0) ? (-d) : d, rnd))
    {
      fprintf (stderr, "t <- x * d should be exact\n");
      exit (1);
    }
  compare = mpfr_cmp_si (t, n);
  if (((inexact == 0) && (compare != 0)) ||
      ((inexact < 0) && (compare >= 0)) ||
      ((inexact > 0) && (compare <= 0)))
    {
      fprintf (stderr, "wrong inexact flag: expected %d, got %d\n", compare,
	       inexact);
      exit (1);
    }

  mpfr_clear (x);
  mpfr_clear (t);
  mpq_clear (q);
}

int
main (void)
{
  tests_start_mpfr ();

  check(-1647229822, 40619231, GMP_RNDZ, -4.055295438754120596e1);
  check(-148939696, 1673285490, GMP_RNDZ, -8.9010331404953485501e-2);
  check(-441322590, 273662545, GMP_RNDZ, -1.6126525096812205362);
  check(-1631156895, 1677687197, GMP_RNDU, -9.722652100563177191e-1);
  check(2141332571, 3117601, GMP_RNDZ, 6.8685267004982347316e2);
  check(75504803, 400207282, GMP_RNDU, 1.8866424074712365155e-1);
  check(643562308, 23100894, GMP_RNDD, 2.7858762002890447462e1);
  check(632549085, 1831935802, GMP_RNDN, 3.4528998467600230393e-1);
  check (1, 1, GMP_RNDN, 1.0);

  tests_end_mpfr ();
  return 0;
}