| Index: gcc/gmp/tests/mpz/t-fib_ui.c
|
| diff --git a/gcc/gmp/tests/mpz/t-fib_ui.c b/gcc/gmp/tests/mpz/t-fib_ui.c
|
| deleted file mode 100644
|
| index 1e84a288e3068d094e41a10f38561a8df1120ef7..0000000000000000000000000000000000000000
|
| --- a/gcc/gmp/tests/mpz/t-fib_ui.c
|
| +++ /dev/null
|
| @@ -1,156 +0,0 @@
|
| -/* Test mpz_fib_ui and mpz_fib2_ui.
|
| -
|
| -Copyright 2000, 2001, 2002 Free Software Foundation, Inc.
|
| -
|
| -This file is part of the GNU MP Library.
|
| -
|
| -The GNU MP 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 3 of the License, or (at your
|
| -option) any later version.
|
| -
|
| -The GNU MP 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 GNU MP Library. If not, see http://www.gnu.org/licenses/. */
|
| -
|
| -#include <stdio.h>
|
| -#include <stdlib.h>
|
| -#include "gmp.h"
|
| -#include "gmp-impl.h"
|
| -#include "tests.h"
|
| -
|
| -
|
| -/* Usage: t-fib_ui [x|num]
|
| -
|
| - Run with no arguments, tests goes up to the initial value of "limit"
|
| - below. With a number argument tests are carried up that far, or with a
|
| - literal "x" tests are continued without limit (this being only meant for
|
| - development purposes).
|
| -
|
| - The size tests performed are designed to partially replicate what will be
|
| - going on in mpz_fib_ui. There's plenty of ASSERTs there, but of course
|
| - they're not normally enabled.
|
| -
|
| - Misfeatures:
|
| -
|
| - The tests on MPN_FIB2_SIZE are a bit useless, since that macro includes a
|
| - +2 for the internal purposes of mpn_fib2_ui. It's probably better to
|
| - give mpn_fib2_ui a run with assertion checking enabled. */
|
| -
|
| -
|
| -#define MPZ_FIB_SIZE_FLOAT(n) \
|
| - ((mp_size_t) ((n) * 0.6942419 / GMP_NUMB_BITS + 1))
|
| -
|
| -
|
| -void
|
| -check_fib_table (void)
|
| -{
|
| - int i;
|
| - mp_limb_t want;
|
| -
|
| - ASSERT_ALWAYS (FIB_TABLE(-1) == 1);
|
| - ASSERT_ALWAYS (FIB_TABLE(0) == 0);
|
| -
|
| - for (i = 1; i <= FIB_TABLE_LIMIT; i++)
|
| - {
|
| - want = FIB_TABLE(i-1) + FIB_TABLE(i-2);
|
| - if (FIB_TABLE(i) != want)
|
| - {
|
| - printf ("FIB_TABLE(%d) wrong\n", i);
|
| - gmp_printf (" got %#Nx\n", &FIB_TABLE(i), 1);
|
| - gmp_printf (" want %#Nx\n", &want, 1);
|
| - abort ();
|
| - }
|
| - }
|
| -}
|
| -
|
| -
|
| -int
|
| -main (int argc, char *argv[])
|
| -{
|
| - unsigned long n;
|
| - unsigned long limit = 100 * BITS_PER_MP_LIMB;
|
| - mpz_t want_fn, want_fn1, got_fn, got_fn1;
|
| -
|
| - tests_start ();
|
| - mp_trace_base = -16;
|
| - if (argc > 1 && argv[1][0] == 'x')
|
| - limit = ULONG_MAX;
|
| - else if (argc > 1)
|
| - limit = atoi (argv[1]);
|
| -
|
| - check_fib_table ();
|
| -
|
| - /* start at n==0 */
|
| - mpz_init_set_ui (want_fn1, 1); /* F[-1] */
|
| - mpz_init_set_ui (want_fn, 0); /* F[0] */
|
| - mpz_init (got_fn);
|
| - mpz_init (got_fn1);
|
| -
|
| - for (n = 0; n < limit; n++)
|
| - {
|
| - /* check our float formula seems right */
|
| - if (MPZ_FIB_SIZE_FLOAT (n) < SIZ(want_fn))
|
| - {
|
| - printf ("MPZ_FIB_SIZE_FLOAT wrong at n=%lu\n", n);
|
| - printf (" MPZ_FIB_SIZE_FLOAT %ld\n", MPZ_FIB_SIZE_FLOAT (n));
|
| - printf (" SIZ(want_fn) %d\n", SIZ(want_fn));
|
| - abort ();
|
| - }
|
| -
|
| - /* check MPN_FIB2_SIZE seems right, compared to actual size and
|
| - compared to our float formula */
|
| - if (MPN_FIB2_SIZE (n) < MPZ_FIB_SIZE_FLOAT (n))
|
| - {
|
| - printf ("MPN_FIB2_SIZE wrong at n=%lu\n", n);
|
| - printf (" MPN_FIB2_SIZE %ld\n", MPN_FIB2_SIZE (n));
|
| - printf (" MPZ_FIB_SIZE_FLOAT %ld\n", MPZ_FIB_SIZE_FLOAT (n));
|
| - abort ();
|
| - }
|
| - if (MPN_FIB2_SIZE (n) < SIZ(want_fn))
|
| - {
|
| - printf ("MPN_FIB2_SIZE wrong at n=%lu\n", n);
|
| - printf (" MPN_FIB2_SIZE %ld\n", MPN_FIB2_SIZE (n));
|
| - printf (" SIZ(want_fn) %d\n", SIZ(want_fn));
|
| - abort ();
|
| - }
|
| -
|
| - mpz_fib2_ui (got_fn, got_fn1, n);
|
| - MPZ_CHECK_FORMAT (got_fn);
|
| - MPZ_CHECK_FORMAT (got_fn1);
|
| - if (mpz_cmp (got_fn, want_fn) != 0 || mpz_cmp (got_fn1, want_fn1) != 0)
|
| - {
|
| - printf ("mpz_fib2_ui(%lu) wrong\n", n);
|
| - mpz_trace ("want fn ", want_fn);
|
| - mpz_trace ("got fn ", got_fn);
|
| - mpz_trace ("want fn1", want_fn1);
|
| - mpz_trace ("got fn1", got_fn1);
|
| - abort ();
|
| - }
|
| -
|
| - mpz_fib_ui (got_fn, n);
|
| - MPZ_CHECK_FORMAT (got_fn);
|
| - if (mpz_cmp (got_fn, want_fn) != 0)
|
| - {
|
| - printf ("mpz_fib_ui(%lu) wrong\n", n);
|
| - mpz_trace ("want fn", want_fn);
|
| - mpz_trace ("got fn", got_fn);
|
| - abort ();
|
| - }
|
| -
|
| - mpz_add (want_fn1, want_fn1, want_fn); /* F[n+1] = F[n] + F[n-1] */
|
| - mpz_swap (want_fn1, want_fn);
|
| - }
|
| -
|
| - mpz_clear (want_fn);
|
| - mpz_clear (want_fn1);
|
| - mpz_clear (got_fn);
|
| - mpz_clear (got_fn1);
|
| -
|
| - tests_end ();
|
| - exit (0);
|
| -}
|
|
|