/* -*-c-*-
*
- * $Id: mpbarrett.c,v 1.1 1999/12/10 23:21:59 mdw Exp $
+ * $Id: mpbarrett.c,v 1.2 1999/12/11 01:50:56 mdw Exp $
*
* Barrett modular reduction
*
/*----- Revision history --------------------------------------------------*
*
* $Log: mpbarrett.c,v $
+ * Revision 1.2 1999/12/11 01:50:56 mdw
+ * Improve initialization slightly.
+ *
* Revision 1.1 1999/12/10 23:21:59 mdw
* Barrett reduction support: works with even moduli.
*
mp_shrink(m);
mb->k = MP_LEN(m);
mb->m = MP_COPY(m);
- b = mp_lsl(MP_NEW, MP_ONE, 2 * MPW_BITS * mb->k);
+ b = mp_create(2 * mb->k + 1);
+ MPX_ZERO(b->v, b->vl - 1);
+ b->vl[-1] = 1;
mp_div(&b, 0, b, m);
mb->mu = b;
}