/*

Copyright 1999-2001, 2004, 2009, 2011 Free Software Foundation, Inc.

This file is part of the GNU MP Library test suite.

The GNU MP Library test suite is free software; you can redistribute it

and/or modify it under the terms of the GNU 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 test suite 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 General

Public License for more details.

You should have received a copy of the GNU General Public License along with

the GNU MP Library test suite.  If not, see https://www.gnu.org/licenses/.  */

#include <stdlib.h>

#include <string.h>

#include <stdio.h>

#include "gmp.h"

#include "gmp-impl.h"

#include "tests.h"

#ifdef OPERATION_copyi

#define func MPN_COPY_INCR

#define reffunc refmpn_copyi

#define funcname "MPN_COPY_INCR"

#endif

#ifdef OPERATION_copyd

#define func MPN_COPY_DECR

#define reffunc refmpn_copyd

#define funcname "MPN_COPY_DECR"

#endif

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

#include <time.h>

int

cputime ()

{

  if (CLOCKS_PER_SEC < 100000)

    return clock () * 1000 / CLOCKS_PER_SEC;

  return clock () / (CLOCKS_PER_SEC / 1000);

}

#else

#include <sys/types.h>

#include <sys/time.h>

#include <sys/resource.h>

int

cputime ()

{

  struct rusage rus;

  getrusage (0, &rus;;

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

}

#endif

static void print_posneg (mp_limb_t);

static void mpn_print (mp_ptr, mp_size_t);

#define LXW ((int) (2 * sizeof (mp_limb_t)))

#define M * 1000000

#ifndef CLOCK

#error "Don't know CLOCK of your machine"

#endif

#ifndef OPS

#define OPS (CLOCK/5)

#endif

#ifndef SIZE

#define SIZE 496

#endif

#ifndef TIMES

#define TIMES OPS/(SIZE+1)

#endif

int

main (int argc, char **argv)

{

  mp_ptr s1, dx, dy;

  int i;

  long t0, t;

  unsigned int test;

  mp_size_t size;

  unsigned int ntests;

  s1 = malloc (SIZE * sizeof (mp_limb_t));

  dx = malloc ((SIZE + 2) * sizeof (mp_limb_t));

  dy = malloc ((SIZE + 2) * sizeof (mp_limb_t));

  ntests = ~(unsigned) 0;

  if (argc == 2)

    ntests = strtol (argv[1], 0, 0);

  for (test = 1; test <= ntests; test++)

    {

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

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

	{

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

	  fflush (stdout);

	}

#endif

#ifdef RANDOM

      size = random () % SIZE + 1;

#else

      size = SIZE;

#endif

      dx[0] = 0x87654321;

      dy[0] = 0x87654321;

      dx[size+1] = 0x12345678;

      dy[size+1] = 0x12345678;

#if TIMES != 1

      mpn_random (s1, size);

      t0 = cputime();

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

	func (dx+1, s1, size);

      t = cputime() - t0;

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

	      t, ((double) t * CLOCK) / (TIMES * size * 1000.0));

#endif

#ifndef NOCHECK

      mpn_random2 (s1, size);

#ifdef PRINT

      mpn_print (s1, size);

#endif

      /* Put garbage in the destination.  */

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

	{

	  dx[i+1] = 0xdead;

	  dy[i+1] = 0xbeef;

	}

      reffunc (dx+1, s1, size);

      func (dy+1, s1, size);

#ifdef PRINT

      mpn_print (dx+1, size);

      mpn_print (dy+1, size);

#endif

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

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

	{

	  mp_size_t s, e;

	  for (s = 0;; s++)

	    if ((unsigned long long) (dx+1)[s] != (unsigned long long) (dy+1)[s])

	      break;

	  for (e = size - 1;; e--)

	    if ((unsigned long long) (dx+1)[e] != (unsigned long long) (dy+1)[e])

	      break;

#ifndef PRINT

	  for (i = s; i <= e; i++)

	    {

	      printf ("%6d: ", i);

	      printf ("%0*llX ", LXW, (unsigned long long) (dx+1)[i]);

	      printf ("%0*llX ", LXW, (unsigned long long) (dy+1)[i]);

	      print_posneg ((dy+1)[i] - (dx+1)[i]);

	      printf ("\n");

	    }

#endif

	  printf ("\n");

	  if (dy[0] != 0x87654321)

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

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

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

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

	  abort();

	}

#endif

    }

  exit (0);

}

static void

print_posneg (mp_limb_t d)

{

  char buf[LXW + 2];

  if (d == 0)

    printf (" %*X", LXW, 0);

  else if (-d < d)

    {

      sprintf (buf, "%llX", (unsigned long long) -d);

      printf ("%*s-%s", LXW - (int) strlen (buf), "", buf);

    }

  else

    {

      sprintf (buf, "%llX", (unsigned long long) d);

      printf ("%*s+%s", LXW - (int) strlen (buf), "", buf);

    }

}

static void

mpn_print (mp_ptr p, mp_size_t size)

{

  mp_size_t i;

  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] >> (GMP_LIMB_BITS/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 (" ");

#endif

    }

  puts ("");

}