/* -*-c-*-
*
- * $Id: mpx.h,v 1.8 1999/12/11 10:57:43 mdw Exp $
+ * $Id: mpx.h,v 1.11 2000/10/08 15:48:35 mdw Exp $
*
* Low level multiprecision arithmetic
*
/*----- Revision history --------------------------------------------------*
*
* $Log: mpx.h,v $
+ * Revision 1.11 2000/10/08 15:48:35 mdw
+ * Rename Karatsuba constants now that we have @gfx_kmul@ too.
+ *
+ * Revision 1.10 2000/10/08 12:06:12 mdw
+ * Provide @mpx_ueq@ for rapidly testing equality of two integers.
+ *
+ * Revision 1.9 1999/12/22 15:49:07 mdw
+ * New function for division by a small integer.
+ *
* Revision 1.8 1999/12/11 10:57:43 mdw
* Karatsuba squaring algorithm.
*
extern void mpx_2c(mpw */*dv*/, mpw */*dvl*/,
const mpw */*v*/, const mpw */*vl*/);
+/* --- @mpx_ueq@ --- *
+ *
+ * Arguments: @const mpw *av, *avl@ = first argument vector base and limit
+ * @const mpw *bv, *bvl@ = second argument vector base and limit
+ *
+ * Returns: Nonzero if the two vectors are equal.
+ *
+ * Use: Performs an unsigned integer test for equality.
+ */
+
+extern int mpx_ueq(const mpw */*av*/, const mpw */*avl*/,
+ const mpw */*bv*/, const mpw */*bvl*/);
+
/* --- @mpx_ucmp@ --- *
*
* Arguments: @const mpw *av, *avl@ = first argument vector base and limit
const mpw */*dv*/, const mpw */*dvl*/,
mpw */*sv*/, mpw */*svl*/);
+/* --- @mpx_udivn@ --- *
+ *
+ * Arguments: @mpw *qv, *qvl@ = storage for the quotient (may overlap
+ * dividend)
+ * @const mpw *rv, *rvl@ = dividend
+ * @mpw d@ = single-precision divisor
+ *
+ * Returns: Remainder after divison.
+ *
+ * Use: Performs a single-precision division operation.
+ */
+
+extern mpw mpx_udivn(mpw */*qv*/, mpw */*qvl*/,
+ const mpw */*rv*/, const mpw */*rvl*/, mpw /*d*/);
+
/*----- Karatsuba multiplication algorithms -------------------------------*/
-/* --- @KARATSUBA_CUTOFF@ --- *
+/* --- @MPK_THRESH@ --- *
*
* This is the limiting length for using Karatsuba algorithms. It's best to
* use the simpler classical multiplication method on numbers smaller than
* this.
*/
-#define KARATSUBA_CUTOFF 16
+#define MPK_THRESH 16
-/* --- @KARATSUBA_SLOP@ --- *
+/* --- @MPK_SLOP@ --- *
*
* The extra number of words required as scratch space by the Karatsuba
* routines. This is a (generous) guess, since the actual amount of space
* required is proportional to the recursion depth.
*/
-#define KARATSUBA_SLOP 32
+#define MPK_SLOP 64
/* --- @mpx_kmul@ --- *
*
* The destination and scratch buffers must be twice as large as
* the larger argument. The scratch space must be twice as
* large as the larger argument, plus the magic number
- * @KARATSUBA_SLOP@.
+ * @MPK_SLOP@.
*/
extern void mpx_kmul(mpw */*dv*/, mpw */*dvl*/,
*
* The destination must be twice as large as the argument. The
* scratch space must be twice as large as the argument, plus
- * the magic number @KARATSUBA_SLOP@.
+ * the magic number @MPK_SLOP@.
*/
extern void mpx_ksqr(mpw */*dv*/, mpw */*dvl*/,