[Gmp-commit] /var/hg/gmp: Add another mpn test.

[mercurial at gmplib.org]

details:   /var/hg/gmp/rev/7fb2a8774242
changeset: 14597:7fb2a8774242
user:      Torbjorn Granlund <tege at gmplib.org>
date:      Wed Feb 01 23:09:23 2012 +0100
description:
Add another mpn test.

diffstat:

 ChangeLog             |    5 ++
 tests/mpn/Makefile.am |   13 +++--
 tests/mpn/t-mul.c     |  101 ++++++++++++++++++++++++++++++++++++++++++++++++++
 3 files changed, 113 insertions(+), 6 deletions(-)

diffs (146 lines):

diff -r 9bcb33a7abf2 -r 7fb2a8774242 ChangeLog
--- a/ChangeLog	Wed Feb 01 18:40:56 2012 +0100
+++ b/ChangeLog	Wed Feb 01 23:09:23 2012 +0100
@@ -1,3 +1,8 @@
+2012-02-01  Torbjorn Granlund  <tege at gmplib.org>
+
+	* tests/mpn/t-mul.c: New file.
+	* tests/mpn/Makefile.am: Compile it.
+
 2012-02-01  Marc Glisse  <marc.glisse at inria.fr>
 
 	* gmpxx.h: Remove check for g++ older than 2.91.
diff -r 9bcb33a7abf2 -r 7fb2a8774242 tests/mpn/Makefile.am
--- a/tests/mpn/Makefile.am	Wed Feb 01 18:40:56 2012 +0100
+++ b/tests/mpn/Makefile.am	Wed Feb 01 23:09:23 2012 +0100
@@ -1,6 +1,7 @@
 ## Process this file with automake to generate Makefile.in
 
-# Copyright 2001, 2002, 2003, 2009, 2010 Free Software Foundation, Inc.
+# Copyright 2001, 2002, 2003, 2009, 2010, 2011, 2012 Free Software Foundation,
+# Inc.
 #
 # This file is part of the GNU MP Library.
 #
@@ -22,11 +23,11 @@
 LDADD = $(top_builddir)/tests/libtests.la $(top_builddir)/libgmp.la
 
 check_PROGRAMS = t-asmtype t-aors_1 t-divrem_1 t-mod_1 t-fat t-get_d	\
-  t-instrument t-iord_u t-mp_bases t-perfsqr t-scan		\
-  t-toom22 t-toom32 t-toom33 t-toom42 t-toom43 t-toom44 	\
-  t-toom52 t-toom53 t-toom62 t-toom63 t-toom6h t-toom8h		\
-  t-hgcd t-hgcd_appr t-matrix22 t-mullo t-mulmod_bnm1 t-sqrmod_bnm1		\
-  t-mulmid t-invert t-div t-bdiv
+  t-instrument t-iord_u t-mp_bases t-perfsqr t-scan			\
+  t-toom22 t-toom32 t-toom33 t-toom42 t-toom43 t-toom44			\
+  t-toom52 t-toom53 t-toom62 t-toom63 t-toom6h t-toom8h			\
+  t-mul t-mullo t-mulmod_bnm1 t-sqrmod_bnm1 t-mulmid			\
+  t-hgcd t-hgcd_appr t-matrix22 t-invert t-div t-bdiv
 
 EXTRA_DIST = toom-shared.h
 
diff -r 9bcb33a7abf2 -r 7fb2a8774242 tests/mpn/t-mul.c
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tests/mpn/t-mul.c	Wed Feb 01 23:09:23 2012 +0100
@@ -0,0 +1,101 @@
+/* Test mpn_mul function for all sizes up to a selected limit.
+
+Copyright 2011, 2012 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 <stdlib.h>
+#include <stdio.h>
+
+#include "gmp.h"
+#include "gmp-impl.h"
+#include "tests.h"
+
+static unsigned
+isqrt (unsigned t)
+{
+  unsigned s, b;
+
+  for (b = 0, s = t;  b++, s >>= 1; )
+    ;
+
+  s = 1 << (b >> 1);
+  if (b & 1)
+    s += s >> 1;
+
+  do
+    {
+      b = t / s;
+      s = (s + b) >> 1;
+    }
+  while (b < s);
+
+  return s;
+}
+
+int
+main (int argc, char **argv)
+{
+  mp_ptr ap, bp, rp, refp;
+  mp_size_t max_n, an, bn, rn;
+  gmp_randstate_ptr rands;
+  int reps;
+  TMP_DECL;
+  TMP_MARK;
+
+  reps = 1;
+
+  tests_start ();
+  TESTS_REPS (reps, argv, argc);
+
+  rands = RANDS;
+
+  /* Re-interpret reps argument as a size argument.  */
+  max_n = isqrt (reps * 25000);
+
+  ap = TMP_ALLOC_LIMBS (max_n + 1);
+  bp = TMP_ALLOC_LIMBS (max_n + 1);
+  rp = TMP_ALLOC_LIMBS (2 * max_n);
+  refp = TMP_ALLOC_LIMBS (2 * max_n);
+
+  for (an = 1; an <= max_n; an += 1)
+    {
+      for (bn = 1; bn <= an; bn += 1)
+	{
+	  mpn_random2 (ap, an + 1);
+	  mpn_random2 (bp, bn + 1);
+
+	  refmpn_mul (refp, ap, an, bp, bn);
+	  mpn_mul (rp, ap, an, bp, bn);
+
+	  rn = an + bn;
+	  if (mpn_cmp (refp, rp, rn))
+	    {
+	      printf ("ERROR, an = %d, bn = %d, rn = %d\n",
+		      (int) an, (int) bn, (int) rn);
+	      printf ("a: "); mpn_dump (ap, an);
+	      printf ("b: "); mpn_dump (bp, bn);
+	      printf ("r:   "); mpn_dump (rp, rn);
+	      printf ("ref: "); mpn_dump (refp, rn);
+	      abort();
+	    }
+	}
+    }
+  TMP_FREE;
+  tests_end ();
+  return 0;
+}