OpenXM_contrib/gmp/mpn/tests/rshift.c - diff

Return to rshift.c CVS log

Up to [local] / OpenXM_contrib / gmp / mpn / tests

Diff for /OpenXM_contrib/gmp/mpn/tests/Attic/rshift.c between version 1.1.1.1 and 1.1.1.2

version 1.1.1.1, 2000/01/10 15:35:26

version 1.1.1.2, 2000/09/09 14:12:46

Line 1

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 2.1 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; 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 "gmp.h"

#include "gmp-impl.h"

#include "longlong.h"

#ifndef USG

#if defined (USG) || defined (__SVR4) || defined (_UNICOS) || defined (__hpux)

#include <sys/time.h>

#include <time.h>

#include <sys/resource.h>

unsigned long

int

cputime ()

{

struct rusage rus;

if (CLOCKS_PER_SEC < 100000)

return clock () * 1000 / CLOCKS_PER_SEC;

getrusage (0, &rus);

return clock () / (CLOCKS_PER_SEC / 1000);

return rus.ru_utime.tv_sec * 1000 + rus.ru_utime.tv_usec / 1000;

}

#else

#include <time.h>

#include <sys/types.h>

#include <sys/time.h>

#include <sys/resource.h>

#ifndef CLOCKS_PER_SEC

int

#define CLOCKS_PER_SEC 1000000

#endif

#if CLOCKS_PER_SEC >= 10000

#define CLOCK_TO_MILLISEC(cl) ((cl) / (CLOCKS_PER_SEC / 1000))

#else

#define CLOCK_TO_MILLISEC(cl) ((cl) * 1000 / CLOCKS_PER_SEC)

#endif

unsigned long

cputime ()

{

return CLOCK_TO_MILLISEC (clock ());

struct rusage rus;

getrusage (0, &rus);

return rus.ru_utime.tv_sec * 1000 + rus.ru_utime.tv_usec / 1000;

}

#endif

Line 41 cputime ()

Line 54 cputime ()

#if defined (__m88k__)

#define CLOCK 20 M

#elif defined (__i386__)

#define CLOCK (16.666667 M)

#define CLOCK (16666667)

#elif defined (__m68k__)

#define CLOCK (20 M)

#elif defined (_IBMR2)

Line 62 cputime ()

Line 75 cputime ()

#endif

#ifndef OPS

#define OPS 10000000

#define OPS (CLOCK/5)

#endif

#ifndef SIZE

#define SIZE 496

Line 144 main (argc, argv)

Line 157 main (argc, argv)

for (test = 0; ; test++)

{

#if TIMES == 1 && ! defined (PRINT)

if (test % (SIZE > 10000 ? 1 : 10000 / SIZE) == 0)

{

printf ("\r%d", test);

fflush (stdout);

}

#endif

#if TIMES == 1

cnt = random () % (BITS_PER_MP_LIMB - 1) + 1;

#endif

#ifdef RANDOM

size = (random () % SIZE + 1);

size = random () % SIZE + 1;

#else

size = SIZE;

#endif

Line 156 main (argc, argv)

Line 181 main (argc, argv)

dx[0] = 0x87654321;

dy[0] = 0x87654321;

#ifdef PRINT

#if TIMES != 1

mpn_print (s1, size);

mpn_random (s1, size);

#endif

#ifndef NOCHECK

t0 = cputime();

for (i = 0; i < TIMES; i++)

cyx = refmpn_rshift (dx+1, s1, size, cnt);

refmpn_rshift (dx+1, s1, size, cnt);

t = cputime() - t0;

#if TIMES != 1

printf ("refmpn_rshift: %5ldms (%.2f cycles/limb)\n",

t, ((double) t * CLOCK) / (OPS * 1000.0));

((double) t * CLOCK) / (OPS * 1000.0));

#endif

#ifdef PRINT

printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyx); mpn_print (dx+1, size);

#endif

t0 = cputime();

for (i = 0; i < TIMES; i++)

cyy = mpn_rshift (dx+1, s1, size, cnt);

mpn_rshift (dx+1, s1, size, cnt);

t = cputime() - t0;

#if TIMES != 1

printf ("mpn_rshift: %5ldms (%.2f cycles/limb)\n",

t, ((double) t * CLOCK) / (OPS * 1000.0));

((double) t * CLOCK) / (OPS * 1000.0));

#endif

#ifndef NOCHECK

mpn_random2 (s1, size);

#ifdef PRINT

printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyy); mpn_print (dx+1, size);

printf ("%-*d ", (int) (2 * sizeof(mp_limb_t)), cnt); mpn_print (s1, size);

#endif

#ifndef NOCHECK

/* Put garbage in the destination. */

for (i = 1; i <= size; i++)

for (i = 0; i < size; i++)

{

dx[i] = 0x7654321;

dx[i+1] = 0xdead;

dy[i] = 0x1234567;

dy[i+1] = 0xbeef;

}

cyx = refmpn_rshift (dx+1, s1, size, cnt);

cyy = mpn_rshift (dy+1, s1, size, cnt);

#ifdef PRINT

printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyx);

mpn_print (dx+1, size);

printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyy);

mpn_print (dy+1, size);

#endif

if (cyx != cyy || mpn_cmp (dx, dy, size+2) != 0

|| dx[size+1] != 0x12345678 || dx[0] != 0x87654321)

|| dx[0] != 0x87654321 || dx[size+1] != 0x12345678)

{

#ifndef PRINT

printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyx);

Line 205 main (argc, argv)

Line 232 main (argc, argv)

printf ("%*lX ", (int) (2 * sizeof(mp_limb_t)), cyy);

mpn_print (dy+1, size);

#endif

printf ("TEST NUMBER %d\n", test);

if (dy[size+1] != 0x12345678)

printf ("clobbered at high end\n");

if (dy[0] != 0x87654321)

printf ("clobbered at low end\n");

abort();

}

#endif

Line 217 mpn_print (mp_ptr p, mp_size_t size)

Line 249 mpn_print (mp_ptr p, mp_size_t size)

for (i = size - 1; i >= 0; i--)

{

#ifdef _LONG_LONG_LIMB

printf ("%0*lX%0*lX", (int) (sizeof(mp_limb_t)),

(unsigned long) (p[i] >> (BITS_PER_MP_LIMB/2)),

(int) (sizeof(mp_limb_t)), (unsigned long) (p[i]));

#else

printf ("%0*lX", (int) (2 * sizeof(mp_limb_t)), p[i]);

#endif

#ifdef SPACE

if (i != 0)

printf (" ");

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>