* Arguments: @gfreduce *r@ = structure to fill in
* @mp *x@ = an integer
*
- * Returns: ---
+ * Returns: Zero if successful; nonzero on failure.
*
* Use: Initializes a context structure for reduction.
*/
-void mpreduce_create(mpreduce *r, mp *p)
+int mpreduce_create(mpreduce *r, mp *p)
{
mpscan sc;
enum { Z = 0, Z1 = 2, X = 4, X0 = 6 };
/* --- Fill in the easy stuff --- */
- assert(MP_POSP(p));
+ if (!MP_POSP(p))
+ return (-1);
d = mp_bits(p);
r->lim = d/MPW_BITS;
r->s = d%MPW_BITS;
}
}
if ((DA(&iv)[DA_LEN(&iv) - 1].op & ~1u) == MPRI_SUB) {
- fprintf(stderr,
- "mpreduce can't handle moduli of the form 2^m + x\n");
- abort();
+ mp_drop(r->p);
+ DA_DESTROY(&iv);
+ return (-1);
}
#undef INSTR
r->iv[i + r->in].argy = b;
}
}
- DA_DESTROY(&iv);
+ DA_DESTROY(&iv);
+ return (0);
}
/* --- @mpreduce_destroy@ --- *