From d23ad7f7f2c93dc84acdb38697f31ca2272b024a Mon Sep 17 00:00:00 2001
From: wbhart <goodwillhart@googlemail.com>
Date: Wed, 1 Apr 2009 20:09:59 +0000
Subject: [PATCH] Went back to my optimised mul.c as it was faster on Windows.

---
 mpz/mul.c | 28 ++++++++++++++++++++--------
 1 file changed, 20 insertions(+), 8 deletions(-)

diff --git a/mpz/mul.c b/mpz/mul.c
index 78f1a5bd..6608c172 100644
--- a/mpz/mul.c
+++ b/mpz/mul.c
@@ -50,13 +50,7 @@ mult (mpz_srcptr u, mpz_srcptr v, mpz_ptr w)
   usize = ABS (usize);
   vsize = ABS (vsize);
 
-  if (usize < vsize)
-    {
-      MPZ_SRCPTR_SWAP (u, v);
-      MP_SIZE_T_SWAP (usize, vsize);
-    }
-
-  if (vsize == 0)
+  if ((usize == 0) || (vsize == 0))
     {
       SIZ(w) = 0;
       return;
@@ -92,6 +86,25 @@ mult (mpz_srcptr u, mpz_srcptr v, mpz_ptr w)
     }
 #endif
 
+  wsize = usize + vsize;
+  
+  if ((wsize <= 32) && (w != u) && (w != v))
+  {
+	  MPZ_REALLOC (w, wsize);
+     wp = PTR(w);
+     if (usize > vsize) mpn_mul_basecase(wp, PTR(u), usize, PTR(v), vsize);
+	  else mpn_mul_basecase(wp, PTR(v), vsize, PTR(u), usize);
+	  wsize -= (wp[wsize - 1] == 0);
+	  SIZ(w) = (sign_product >= 0 ? wsize : -wsize);
+     return;
+  }
+
+  if (usize < vsize)
+    {
+      MPZ_SRCPTR_SWAP (u, v);
+      MP_SIZE_T_SWAP (usize, vsize);
+    }
+
   TMP_MARK;
   free_me = NULL;
   up = u->_mp_d;
@@ -99,7 +112,6 @@ mult (mpz_srcptr u, mpz_srcptr v, mpz_ptr w)
   wp = w->_mp_d;
 
   /* Ensure W has space enough to store the result.  */
-  wsize = usize + vsize;
   if (w->_mp_alloc < wsize)
     {
       if (wp == up || wp == vp)