/* -*-c-*-
*
- * $Id: mpx.c,v 1.9 2000/06/26 07:52:50 mdw Exp $
+ * $Id: mpx.c,v 1.10 2000/10/08 12:06:12 mdw Exp $
*
* Low-level multiprecision arithmetic
*
/*----- Revision history --------------------------------------------------*
*
* $Log: mpx.c,v $
+ * Revision 1.10 2000/10/08 12:06:12 mdw
+ * Provide @mpx_ueq@ for rapidly testing equality of two integers.
+ *
* Revision 1.9 2000/06/26 07:52:50 mdw
* Portability fix for the bug fix.
*
MPX_UADDN(dv, dvl, 1);
}
+/* --- @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.
+ */
+
+int mpx_ueq(const mpw *av, const mpw *avl, const mpw *bv, const mpw *bvl)
+{
+ MPX_SHRINK(av, avl);
+ MPX_SHRINK(bv, bvl);
+ if (avl - av != bvl - bv)
+ return (0);
+ while (av < avl) {
+ if (*av++ != *bv++)
+ return (0);
+ }
+ return (1);
+}
+
/* --- @mpx_ucmp@ --- *
*
* Arguments: @const mpw *av, *avl@ = first argument vector base and limit
}
return (0);
}
-
+
/* --- @mpx_uadd@ --- *
*
* Arguments: @mpw *dv, *dvl@ = destination vector base and limit
ALLOC(d, dl, al - a + (n + MPW_BITS - 1) / MPW_BITS);
mpx_lsl(d, dl, a, al, n);
- if (MPX_UCMP(d, dl, !=, c, cl)) {
+ if (!mpx_ueq(d, dl, c, cl)) {
fprintf(stderr, "\n*** lsl(%i) failed\n", n);
dumpmp(" a", a, al);
dumpmp("expected", c, cl);
ALLOC(d, dl, al - a + (n + MPW_BITS - 1) / MPW_BITS + 1);
mpx_lsr(d, dl, a, al, n);
- if (MPX_UCMP(d, dl, !=, c, cl)) {
+ if (!mpx_ueq(d, dl, c, cl)) {
fprintf(stderr, "\n*** lsr(%i) failed\n", n);
dumpmp(" a", a, al);
dumpmp("expected", c, cl);
ALLOC(d, dl, MAX(al - a, bl - b) + 1);
mpx_uadd(d, dl, a, al, b, bl);
- if (MPX_UCMP(d, dl, !=, c, cl)) {
+ if (!mpx_ueq(d, dl, c, cl)) {
fprintf(stderr, "\n*** uadd failed\n");
dumpmp(" a", a, al);
dumpmp(" b", b, bl);
ALLOC(d, dl, al - a);
mpx_usub(d, dl, a, al, b, bl);
- if (MPX_UCMP(d, dl, !=, c, cl)) {
+ if (!mpx_ueq(d, dl, c, cl)) {
fprintf(stderr, "\n*** usub failed\n");
dumpmp(" a", a, al);
dumpmp(" b", b, bl);
ALLOC(d, dl, (al - a) + (bl - b));
mpx_umul(d, dl, a, al, b, bl);
- if (MPX_UCMP(d, dl, !=, c, cl)) {
+ if (!mpx_ueq(d, dl, c, cl)) {
fprintf(stderr, "\n*** umul failed\n");
dumpmp(" a", a, al);
dumpmp(" b", b, bl);
ALLOC(d, dl, 2 * (al - a));
mpx_usqr(d, dl, a, al);
- if (MPX_UCMP(d, dl, !=, c, cl)) {
+ if (!mpx_ueq(d, dl, c, cl)) {
fprintf(stderr, "\n*** usqr failed\n");
dumpmp(" a", a, al);
dumpmp("expected", c, cl);
ALLOC(s, sl, (bl - b) + 1);
mpx_udiv(qq, qql, a, al, b, bl, s, sl);
- if (MPX_UCMP(qq, qql, !=, q, ql) ||
- MPX_UCMP(a, al, !=, r, rl)) {
+ if (!mpx_ueq(qq, qql, q, ql) ||
+ !mpx_ueq(a, al, r, rl)) {
fprintf(stderr, "\n*** udiv failed\n");
dumpmp(" divisor", b, bl);
dumpmp("expect r", r, rl);
return (0);
}
-
#endif
/*----- That's all, folks -------------------------------------------------*/
/* -*-c-*-
*
- * $Id: mpx.h,v 1.9 1999/12/22 15:49:07 mdw Exp $
+ * $Id: mpx.h,v 1.10 2000/10/08 12:06:12 mdw Exp $
*
* Low level multiprecision arithmetic
*
/*----- Revision history --------------------------------------------------*
*
* $Log: mpx.h,v $
+ * 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.
*
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