/* Test file for in-place operations.
Copyright 2000, 2001, 2002 Free Software Foundation.
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 "gmp.h"
#include "mpfr.h"
#include "mpfr-impl.h"
typedef void (*fct_t)();
fct_t testfunc;
void test3 _PROTO ((char *, mp_prec_t, mp_rnd_t));
void test4 _PROTO ((char *, mp_prec_t, mp_rnd_t));
void test3a _PROTO ((char *, mp_prec_t, mp_rnd_t));
void test2ui _PROTO ((char *, mp_prec_t, mp_rnd_t));
void testui2 _PROTO ((char *, mp_prec_t, mp_rnd_t));
void test2 _PROTO ((char *, mp_prec_t, mp_rnd_t));
void test2a _PROTO ((char *, mp_prec_t));
int mpfr_compare _PROTO ((mpfr_t, mpfr_t));
/* same than mpfr_cmp, but returns 0 for both NaN's */
int
mpfr_compare (mpfr_t a, mpfr_t b)
{
return (MPFR_IS_NAN(a)) ? !MPFR_IS_NAN(b) :
(MPFR_IS_NAN(b) || mpfr_cmp(a, b));
}
void
test3 (char *foo, mp_prec_t prec, mp_rnd_t rnd)
{
mpfr_t ref1, ref2, ref3;
mpfr_t res1;
int i;
#ifdef DEBUG
printf("checking %s\n", foo);
#endif
mpfr_init2 (ref1, prec);
mpfr_init2 (ref2, prec);
mpfr_init2 (ref3, prec);
mpfr_init2 (res1, prec);
/* for each variable, consider each of the following 6 possibilities:
NaN, +Infinity, -Infinity, +0, -0 or a random number */
for (i=0; i<36; i++) {
if (i%6==0) mpfr_set_nan (ref2);
if (i%6==1) mpfr_set_inf (ref2, 1);
if (i%6==2) mpfr_set_inf (ref2, -1);
if (i%6==3) mpfr_set_d (ref2, 0.0, GMP_RNDN);
if (i%6==4) mpfr_set_d (ref2, -0.0, GMP_RNDN);
if (i%6==5) mpfr_random (ref2);
if (i/6==0) mpfr_set_nan (ref3);
if (i/6==1) mpfr_set_inf (ref3, 1);
if (i/6==2) mpfr_set_inf (ref3, -1);
if (i/6==3) mpfr_set_d (ref3, 0.0, GMP_RNDN);
if (i/6==4) mpfr_set_d (ref3, -0.0, GMP_RNDN);
if (i/6==5) mpfr_random (ref3);
/* reference call: foo(a, b, c) */
testfunc (ref1, ref2, ref3, rnd);
/* foo(a, a, c) */
mpfr_set (res1, ref2, rnd); /* exact operation */
testfunc (res1, res1, ref3, rnd);
if (mpfr_compare (res1, ref1)) {
fprintf (stderr, "Error for %s(a, a, c) for a=%e, c=%e\n", foo,
mpfr_get_d1 (ref2), mpfr_get_d1 (ref3));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref1),
mpfr_get_d1 (res1));
exit (1);
}
/* foo(a, b, a) */
mpfr_set (res1, ref3, rnd);
testfunc (res1, ref2, res1, rnd);
if (mpfr_compare (res1, ref1)) {
fprintf (stderr, "Error for %s(a, b, a) for b=%e, a=%e\n", foo,
mpfr_get_d1 (ref2), mpfr_get_d1 (ref3));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref1),
mpfr_get_d1 (res1));
exit (1);
}
/* foo(a, a, a) */
mpfr_set (ref3, ref2, rnd);
testfunc (ref1, ref2, ref3, rnd);
mpfr_set (res1, ref2, rnd);
testfunc (res1, res1, res1, rnd);
if (mpfr_compare (res1, ref1)) {
fprintf (stderr, "Error for %s(a, a, a) for a=%e\n", foo,
mpfr_get_d1 (ref2));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref1),
mpfr_get_d1 (res1));
exit (1);
}
}
mpfr_clear (ref1);
mpfr_clear (ref2);
mpfr_clear (ref3);
mpfr_clear (res1);
}
void
test4 (char *foo, mp_prec_t prec, mp_rnd_t rnd)
{
mpfr_t ref, op1, op2, op3;
mpfr_t res;
int i, j, k;
#ifdef DEBUG
printf("checking %s\n", foo);
#endif
mpfr_init2 (ref, prec);
mpfr_init2 (op1, prec);
mpfr_init2 (op2, prec);
mpfr_init2 (op3, prec);
mpfr_init2 (res, prec);
/* for each variable, consider each of the following 6 possibilities:
NaN, +Infinity, -Infinity, +0, -0 or a random number */
for (i=0; i<6; i++)
{
MPFR_CLEAR_FLAGS(op1);
if (i==0) mpfr_set_nan (op1);
if (i==1) mpfr_set_inf (op1, 1);
if (i==2) mpfr_set_inf (op1, -1);
if (i==3) mpfr_set_d (op1, 0.0, GMP_RNDN);
if (i==4) mpfr_set_d (op1, -0.0, GMP_RNDN);
if (i==5) mpfr_random (op1);
for (j=0; j<6; j++)
{
MPFR_CLEAR_FLAGS(op2);
if (j==0) mpfr_set_nan (op2);
if (j==1) mpfr_set_inf (op2, 1);
if (j==2) mpfr_set_inf (op2, -1);
if (j==3) mpfr_set_d (op2, 0.0, GMP_RNDN);
if (j==4) mpfr_set_d (op2, -0.0, GMP_RNDN);
if (j==5) mpfr_random (op2);
for (k=0; k<6; k++)
{
MPFR_CLEAR_FLAGS(op3);
if (k==0) mpfr_set_nan (op3);
if (k==1) mpfr_set_inf (op3, 1);
if (k==2) mpfr_set_inf (op3, -1);
if (k==3) mpfr_set_d (op3, 0.0, GMP_RNDN);
if (k==4) mpfr_set_d (op3, -0.0, GMP_RNDN);
if (k==5) mpfr_random (op3);
/* reference call: foo(s, a, b, c) */
testfunc (ref, op1, op2, op3, rnd);
/* foo(a, a, b, c) */
mpfr_set (res, op1, rnd); /* exact operation */
testfunc (res, res, op2, op3, rnd);
if (mpfr_compare (res, ref))
{
fprintf (stderr,
"Error for %s(a, a, b, c) for a=%e, b=%e, c=%e\n",
foo,
mpfr_get_d1 (op1), mpfr_get_d1 (op2), mpfr_get_d1 (op3));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref),
mpfr_get_d1 (res));
exit (1);
}
/* foo(b, a, b, c) */
mpfr_set (res, op2, rnd);
testfunc (res, op1, res, op3, rnd);
if (mpfr_compare (res, ref))
{
fprintf (stderr,
"Error for %s(b, a, b, c) for a=%e, b=%e, c=%e\n",
foo,
mpfr_get_d1 (op1), mpfr_get_d1 (op2), mpfr_get_d1 (op3));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref),
mpfr_get_d1 (res));
exit (1);
}
/* foo(c, a, b, c) */
mpfr_set (res, op3, rnd);
testfunc (res, op1, op2, res, rnd);
if (mpfr_compare (res, ref))
{
fprintf (stderr,
"Error for %s(c, a, b, c) for a=%e, b=%e, c=%e\n",
foo,
mpfr_get_d1 (op1), mpfr_get_d1 (op2), mpfr_get_d1 (op3));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref),
mpfr_get_d1 (res));
exit (1);
}
/* foo(a, a, a,c) */
testfunc (ref, op1, op1, op3, rnd);
mpfr_set (res, op1, rnd);
testfunc (res, res, res, op3, rnd);
if (mpfr_compare (res, ref))
{
fprintf (stderr,
"Error for %s(a, a, a, c) for a=%e, a=%e, c=%e\n",
foo,
mpfr_get_d1 (op1), mpfr_get_d1 (op2), mpfr_get_d1 (op3));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref),
mpfr_get_d1 (res));
exit (1);
}
/* foo(a, a, b,a) */
testfunc (ref, op1, op2, op1, rnd);
mpfr_set (res, op1, rnd);
testfunc (res, res, op2, res, rnd);
if (mpfr_compare (res, ref))
{
fprintf (stderr,
"Error for %s(a, a, b, a) for a=%e, a=%e, c=%e\n",
foo,
mpfr_get_d1 (op1), mpfr_get_d1 (op2), mpfr_get_d1 (op3));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref),
mpfr_get_d1 (res));
exit (1);
}
/* foo(b, a, b, b) */
testfunc (ref, op1, op2, op2, rnd);
mpfr_set (res, op2, rnd);
testfunc (res, op1, res, res, rnd);
if (mpfr_compare (res, ref))
{
fprintf (stderr,
"Error for %s(b, a, b, b) for a=%e, a=%e, c=%e\n",
foo,
mpfr_get_d1 (op1), mpfr_get_d1 (op2), mpfr_get_d1 (op3));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref),
mpfr_get_d1 (res));
exit (1);
}
/* foo (a, a, a, a) */
testfunc (ref, op1, op1, op1 ,rnd);
mpfr_set (res, op1, rnd);
testfunc (res, res, res, res, rnd);
if (mpfr_compare (res, ref))
{
fprintf (stderr,
"Error for %s(a, a, a, a) for a=%e\n", foo,
mpfr_get_d1 (op1));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref),
mpfr_get_d1 (res));
exit (1);
}
}
}
}
mpfr_clear (ref);
mpfr_clear (op1);
mpfr_clear (op2);
mpfr_clear (op3);
mpfr_clear (res);
}
void
test2ui (char *foo, mp_prec_t prec, mp_rnd_t rnd)
{
mpfr_t ref1, ref2;
unsigned int ref3;
mp_limb_t c[1];
mpfr_t res1;
int i;
#ifdef DEBUG
printf("checking %s\n", foo);
#endif
mpfr_init2 (ref1, prec);
mpfr_init2 (ref2, prec);
mpfr_init2 (res1, prec);
/* ref2 can be NaN, +Inf, -Inf, +0, -0 or any number
ref3 can be 0 or any number */
for (i=0; i<12; i++)
{
if (i%6==0) mpfr_set_nan (ref2);
if (i%6==1) mpfr_set_inf (ref2, 1);
if (i%6==2) mpfr_set_inf (ref2, -1);
if (i%6==3) mpfr_set_d (ref2, 0.0, GMP_RNDN);
if (i%6==4) mpfr_set_d (ref2, -0.0, GMP_RNDN);
if (i%6==5) mpfr_random (ref2);
if (i/6==0)
ref3=0;
else
{
mpn_random (c, 1);
ref3 = (unsigned int) c[0];
}
/* reference call: foo(a, b, c) */
testfunc (ref1, ref2, ref3, rnd);
/* foo(a, a, c) */
mpfr_set (res1, ref2, rnd); /* exact operation */
testfunc (res1, res1, ref3, rnd);
if (mpfr_compare (res1, ref1))
{
fprintf (stderr, "Error for %s(a, a, c) for a=%e c=%u\n", foo,
mpfr_get_d1 (ref2), ref3);
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref1),
mpfr_get_d1 (res1));
exit (1);
}
}
mpfr_clear (ref1);
mpfr_clear (ref2);
mpfr_clear (res1);
}
void
testui2 (char *foo, mp_prec_t prec, mp_rnd_t rnd)
{
mpfr_t ref1, ref3;
unsigned int ref2;
mp_limb_t c[1];
mpfr_t res1;
int i;
#ifdef DEBUG
printf("checking %s\n", foo);
#endif
mpfr_init2 (ref1, prec);
mpfr_init2 (ref3, prec);
mpfr_init2 (res1, prec);
mpfr_random (ref3);
mpn_random (c, 1);
ref2 = (unsigned int) c[0];
for (i=0; i<12; i++) {
if (i%6==0) mpfr_set_nan (ref3);
if (i%6==1) mpfr_set_inf (ref3, 1);
if (i%6==2) mpfr_set_inf (ref3, -1);
if (i%6==3) mpfr_set_d (ref3, 0.0, GMP_RNDN);
if (i%6==4) mpfr_set_d (ref3, -0.0, GMP_RNDN);
if (i%6==5) mpfr_random (ref3);
if (i/6==0) ref2=0;
else {
mpn_random (c, 1);
ref2 = (unsigned int) c[0];
}
/* reference call: foo(a, b, c) */
testfunc (ref1, ref2, ref3, rnd);
/* foo(a, b, a) */
mpfr_set (res1, ref3, rnd); /* exact operation */
testfunc (res1, ref2, res1, rnd);
if (mpfr_compare (res1, ref1)) {
fprintf (stderr, "Error for %s(a, b, a) for b=%u a=%e\n", foo,
ref2, mpfr_get_d1 (ref3));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref1),
mpfr_get_d1 (res1));
exit (1);
}
}
mpfr_clear (ref1);
mpfr_clear (ref3);
mpfr_clear (res1);
}
/* foo(mpfr_ptr, mpfr_srcptr, mp_rndt) */
void
test2 (char *foo, mp_prec_t prec, mp_rnd_t rnd)
{
mpfr_t ref1, ref2;
mpfr_t res1;
int i;
#ifdef DEBUG
printf("checking %s\n", foo);
#endif
mpfr_init2 (ref1, prec);
mpfr_init2 (ref2, prec);
mpfr_init2 (res1, prec);
mpfr_random (ref2);
for (i=0; i<6; i++) {
if (i==0) mpfr_set_nan (ref2);
if (i==1) mpfr_set_inf (ref2, 1);
if (i==2) mpfr_set_inf (ref2, -1);
if (i==3) mpfr_set_d (ref2, 0.0, GMP_RNDN);
if (i==4) mpfr_set_d (ref2, -0.0, GMP_RNDN);
if (i==5) mpfr_random (ref2);
/* reference call: foo(a, b) */
testfunc (ref1, ref2, rnd);
/* foo(a, a) */
mpfr_set (res1, ref2, rnd); /* exact operation */
testfunc (res1, res1, rnd);
if (mpfr_compare (res1, ref1)) {
fprintf (stderr, "Error for %s(a, a) for a=%e\n", foo, mpfr_get_d1 (ref2));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref1),
mpfr_get_d1 (res1));
exit (1);
}
}
mpfr_clear (ref1);
mpfr_clear (ref2);
mpfr_clear (res1);
}
/* foo(mpfr_ptr, mpfr_srcptr) */
void
test2a (char *foo, mp_prec_t prec)
{
mpfr_t ref1, ref2;
mpfr_t res1;
int i;
#ifdef DEBUG
printf("checking %s\n", foo);
#endif
mpfr_init2 (ref1, prec);
mpfr_init2 (ref2, prec);
mpfr_init2 (res1, prec);
mpfr_random (ref2);
for (i=0; i<6; i++) {
if (i==0) mpfr_set_nan (ref2);
if (i==1) mpfr_set_inf (ref2, 1);
if (i==2) mpfr_set_inf (ref2, -1);
if (i==3) mpfr_set_d (ref2, 0.0, GMP_RNDN);
if (i==4) mpfr_set_d (ref2, -0.0, GMP_RNDN);
if (i==5) mpfr_random (ref2);
/* reference call: foo(a, b) */
testfunc (ref1, ref2);
/* foo(a, a) */
mpfr_set (res1, ref2, GMP_RNDN); /* exact operation */
testfunc (res1, res1);
if (mpfr_compare (res1, ref1)) {
fprintf (stderr, "Error for %s(a, a) for a=%e\n", foo, mpfr_get_d1 (ref2));
fprintf (stderr, "expected %e, got %e\n", mpfr_get_d1 (ref1),
mpfr_get_d1 (res1));
exit (1);
}
}
mpfr_clear (ref1);
mpfr_clear (ref2);
mpfr_clear (res1);
}
/* one operand, two results */
void
test3a (char *foo, mp_prec_t prec, mp_rnd_t rnd)
{
mpfr_t ref1, ref2, ref3;
mpfr_t res1, res2;
int i;
#ifdef DEBUG
printf("checking %s\n", foo);
#endif
mpfr_init2 (ref1, prec);
mpfr_init2 (ref2, prec);
mpfr_init2 (ref3, prec);
mpfr_init2 (res1, prec);
mpfr_init2 (res2, prec);
mpfr_random (ref3);
for (i=0; i<6; i++) {
if (i==0) mpfr_set_nan (ref3);
if (i==1) mpfr_set_inf (ref3, 1);
if (i==2) mpfr_set_inf (ref3, -1);
if (i==3) mpfr_set_d (ref3, 0.0, GMP_RNDN);
if (i==4) mpfr_set_d (ref3, -0.0, GMP_RNDN);
if (i==5) mpfr_random (ref3);
/* reference call: foo(a, b, c) */
testfunc (ref1, ref2, ref3, rnd);
/* foo(a, b, a) */
mpfr_set (res1, ref3, rnd); /* exact operation */
testfunc (res1, res2, res1, rnd);
if (mpfr_compare (res1, ref1) || mpfr_compare (res2, ref2)) {
fprintf (stderr, "Error for %s(a, b, a) for a=%e\n", foo, mpfr_get_d1 (ref3));
fprintf (stderr, "expected (%e,%e), got (%e,%e)\n", mpfr_get_d1 (ref1),
mpfr_get_d1 (ref2), mpfr_get_d1 (res1), mpfr_get_d1 (res2));
exit (1);
}
/* foo(a, b, b) */
mpfr_set (res2, ref3, rnd); /* exact operation */
testfunc (res1, res2, res2, rnd);
if (mpfr_compare (res1, ref1) || mpfr_compare (res2, ref2)) {
fprintf (stderr, "Error for %s(a, b, b) for b=%e\n", foo, mpfr_get_d1 (ref3));
fprintf (stderr, "expected (%e,%e), got (%e,%e)\n", mpfr_get_d1 (ref1),
mpfr_get_d1 (ref2), mpfr_get_d1 (res1), mpfr_get_d1 (res2));
exit (1);
}
}
mpfr_clear (ref1);
mpfr_clear (ref2);
mpfr_clear (ref3);
mpfr_clear (res1);
mpfr_clear (res2);
}
int
main (void)
{
testfunc = (fct_t) mpfr_add; test3 ("mpfr_add", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_add_ui; test2ui ("mpfr_add_ui", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_agm; test3 ("mpfr_agm", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_ceil; test2 ("mpfr_ceil", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_div; test3 ("mpfr_div", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_div_2exp; test2ui ("mpfr_div_2exp", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_div_ui; test2ui ("mpfr_div_ui", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_exp; test2 ("mpfr_exp", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_floor; test2 ("mpfr_floor", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_log; test2 ("mpfr_log", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_mul; test3 ("mpfr_mul", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_mul_2exp; test2ui ("mpfr_mul_2exp", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_mul_ui; test2ui ("mpfr_mul_ui", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_neg; test2 ("mpfr_neg", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_pow_ui; test2ui ("mpfr_pow_ui", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_reldiff; test3 ("mpfr_reldiff", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_sub; test3 ("mpfr_sub", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_sub_ui; test2ui ("mpfr_sub_ui", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_sqrt; test2 ("mpfr_sqrt", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_ui_div; testui2 ("mpfr_ui_div", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_ui_sub; testui2 ("mpfr_ui_sub", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_trunc; test2 ("mpfr_trunc", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_asin; test2 ("mpfr_asin", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_acos; test2 ("mpfr_acos", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_atan; test2 ("mpfr_atan", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_sinh; test2 ("mpfr_sinh", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_cosh; test2 ("mpfr_cosh", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_tanh; test2 ("mpfr_tanh", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_asinh; test2 ("mpfr_asinh", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_acosh; test2 ("mpfr_acosh", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_atanh; test2 ("mpfr_atanh", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_exp2; test2 ("mpfr_exp2", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_cos; test2 ("mpfr_cos", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_sin; test2 ("mpfr_sin", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_tan; test2 ("mpfr_tan", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_log10; test2 ("mpfr_log10", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_log2; test2 ("mpfr_log2", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_pow; test3 ("mpfr_pow", 53, GMP_RNDN);
testfunc = (fct_t) mpfr_fma; test4 ("mpfr_fma", 53, GMP_RNDN);
return 0;
}
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to reuse.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / mpfr / tests