diff options
author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2001-12-06 17:22:14 +0000 |
---|---|---|
committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2001-12-06 17:22:14 +0000 |
commit | 0fb01bc484b0b984de11dadaa5f4a12fb965c670 (patch) | |
tree | 401d4fed2025ee9ebaedf7feae4107a4083405da /tests/tconst_euler.c | |
parent | a00e540c4bac91bcf81c14a1387a10cad4798837 (diff) | |
download | mpfr-0fb01bc484b0b984de11dadaa5f4a12fb965c670.tar.gz |
Files renamed (constants).
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@1631 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'tests/tconst_euler.c')
-rw-r--r-- | tests/tconst_euler.c | 90 |
1 files changed, 90 insertions, 0 deletions
diff --git a/tests/tconst_euler.c b/tests/tconst_euler.c new file mode 100644 index 000000000..ab9f74321 --- /dev/null +++ b/tests/tconst_euler.c @@ -0,0 +1,90 @@ +/* Test file for mpfr_const_euler. + +Copyright (C) 2001 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" + +int +main (int argc, char *argv[]) +{ + mpfr_t gamma, y, z, t; + unsigned int err, prec, yprec, p0 = 2, p1 = 200; + mp_rnd_t rnd; + + prec = (argc < 2) ? 53 : atoi(argv[1]); + + if (argc > 1) + { + mpfr_init2 (gamma, prec); + mpfr_const_euler (gamma, GMP_RNDN); + printf("gamma="); mpfr_out_str (stdout, 10, 0, gamma, GMP_RNDD); + putchar ('\n'); + mpfr_clear (gamma); + return 0; + } + + mpfr_init (y); + mpfr_init (z); + mpfr_init (t); + + for (prec = p0; prec <= p1; prec++) + { + mpfr_set_prec (z, prec); + mpfr_set_prec (t, prec); + yprec = prec + 10; + + for (rnd=0; rnd<4; rnd++) + { + mpfr_set_prec (y, yprec); + mpfr_const_euler (y, rnd); + err = (rnd == GMP_RNDN) ? yprec + 1 : yprec; + if (mpfr_can_round (y, err, rnd, rnd, prec)) + { + mpfr_set (t, y, rnd); + mpfr_const_euler (z, rnd); + if (mpfr_cmp (t, z)) + { + printf ("results differ for prec=%u rnd_mode=%s\n", prec, + mpfr_print_rnd_mode (rnd)); + printf (" got "); + mpfr_out_str (stdout, 2, prec, z, GMP_RNDN); + putchar ('\n'); + printf (" expected "); + mpfr_out_str (stdout, 2, prec, t, GMP_RNDN); + putchar ('\n'); + printf (" approximation was "); + mpfr_print_binary (y); + putchar ('\n'); + exit (1); + } + } + } + } + + mpfr_clear (y); + mpfr_clear (z); + mpfr_clear (t); + + return 0; +} |