X-Git-Url: https://oss.titaniummirror.com/gitweb/?a=blobdiff_plain;f=mpfr%2Ftests%2Ftstckintc.c;fp=mpfr%2Ftests%2Ftstckintc.c;h=83fd63f85efbd7a8fc36eaf2562833419810ee9c;hb=6fed43773c9b0ce596dca5686f37ac3fc0fa11c0;hp=0000000000000000000000000000000000000000;hpb=27b11d56b743098deb193d510b337ba22dc52e5c;p=msp430-gcc.git diff --git a/mpfr/tests/tstckintc.c b/mpfr/tests/tstckintc.c new file mode 100644 index 00000000..83fd63f8 --- /dev/null +++ b/mpfr/tests/tstckintc.c @@ -0,0 +1,183 @@ +/* Test file for mpfr_custom_* + +Copyright 2005, 2006, 2007, 2008, 2009 Free Software Foundation, Inc. +Contributed by the Arenaire and Cacao projects, INRIA. + +This file is part of the GNU MPFR Library. + +The GNU 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 GNU 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 GNU MPFR Library; see the file COPYING.LIB. If not, write to +the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, +MA 02110-1301, USA. */ + +#include +#include +#include +#include + +#include "mpfr-test.h" + +#define BUFFER_SIZE 250 +#define PREC_TESTED 200 + +long Buffer[BUFFER_SIZE]; +char *stack = (char *) Buffer; +mp_prec_t p = PREC_TESTED; + +#define ALIGNED(s) (((s) + sizeof (long) - 1) / sizeof (long) * sizeof (long)) + +static void * +new_st (size_t s) +{ + void *p = (void *) stack; + stack += ALIGNED (s); + if (MPFR_UNLIKELY (stack > (char *) &Buffer[BUFFER_SIZE])) + { + printf ("Stack overflow.\n"); + exit (1); + } + return p; +} + + /* Alloc a new mpfr_t on the main stack */ +static mpfr_ptr +new_mpfr (mp_prec_t p) +{ + mpfr_ptr x = (mpfr_ptr) new_st (sizeof (mpfr_t)); + void *mantissa = new_st (mpfr_custom_get_size (p)); + mpfr_custom_init (mantissa, p); + mpfr_custom_init_set (x, 0, 0, p, mantissa); + return x; +} + + /* Garbage the stack by keeping only x */ +static mpfr_ptr +return_mpfr (mpfr_ptr x, char *old_stack) +{ + void *mantissa = mpfr_custom_get_mantissa (x); + size_t size_mantissa = mpfr_custom_get_size (mpfr_get_prec (x)); + mpfr_ptr newx; + + memmove (old_stack, x, sizeof (mpfr_t)); + memmove (old_stack + ALIGNED (sizeof (mpfr_t)), mantissa, size_mantissa); + newx = (mpfr_ptr) old_stack; + mpfr_custom_move (newx, old_stack + ALIGNED (sizeof (mpfr_t))); + stack = old_stack + ALIGNED (sizeof (mpfr_t)) + ALIGNED (size_mantissa); + return newx; +} + +static void +test1 (void) +{ + mpfr_ptr x, y; + char *org; + + org = stack; + x = new_mpfr (p); + y = new_mpfr (p); + mpfr_set_ui (x, 42, GMP_RNDN); + mpfr_set_ui (y, 17, GMP_RNDN); + mpfr_add (y, x, y, GMP_RNDN); + y = return_mpfr (y, org); + if (y != x || mpfr_cmp_ui (y, 59) != 0) + { + printf ("Compact (1) failed!\n"); + exit (1); + } + stack = org; +} + +/* We build the MPFR variable each time it is needed */ +/* a[0] is the kind, a[1] is the exponent, &a[2] is the mantissa */ +static long * +dummy_new (void) +{ + long *r; + + r = (long *) new_st (ALIGNED (2 * sizeof (long)) + + ALIGNED (mpfr_custom_get_size (p))); + MPFR_ASSERTN (r != NULL); + (mpfr_custom_init) (&r[2], p); + r[0] = (int) MPFR_NAN_KIND; + r[1] = 0; + return r; +} + +static long * +dummy_set_si (long si) +{ + mpfr_t x; + long * r = dummy_new (); + (mpfr_custom_init_set) (x, 0, 0, p, &r[2]); + mpfr_set_si (x, si, GMP_RNDN); + r[0] = mpfr_custom_get_kind (x); + r[1] = mpfr_custom_get_exp (x); + return r; +} + +static long * +dummy_add (long *a, long *b) +{ + mpfr_t x, y, z; + long *r = dummy_new (); + mpfr_custom_init_set (x, 0, 0, p, &r[2]); + (mpfr_custom_init_set) (y, a[0], a[1], p, &a[2]); + mpfr_custom_init_set (z, b[0], b[1], p, &b[2]); + mpfr_add (x, y, z, GMP_RNDN); + r[0] = (mpfr_custom_get_kind) (x); + r[1] = (mpfr_custom_get_exp) (x); + return r; +} + +static long * +dummy_compact (long *r, char *org_stack) +{ + memmove (org_stack, r, + ALIGNED (2*sizeof (long)) + ALIGNED ((mpfr_custom_get_size) (p))); + return (long *) org_stack; +} + +static void +test2 (void) +{ + mpfr_t x; + char *org = stack; + long *a, *b, *c; + + a = dummy_set_si (42); + b = dummy_set_si (17); + c = dummy_add (a, b); + c = dummy_compact (c, org); + (mpfr_custom_init_set) (x, c[0], c[1], p, &c[2]); + if (c != a || mpfr_cmp_ui (x, 59) != 0) + { + printf ("Compact (2) failed! c=%p a=%p\n", c, a); + mpfr_dump (x); + exit (1); + } + stack = org; +} + +int +main (void) +{ + tests_start_mpfr (); + /* We test iff long = mp_limb_t */ + if (sizeof (long) == sizeof (mp_limb_t)) + { + test1 (); + test2 (); + } + tests_end_mpfr (); + return 0; +}