/* -*-c-*-
*
- * $Id: mpx.c,v 1.2 1999/11/13 01:50:59 mdw Exp $
+ * $Id: mpx.c,v 1.3 1999/11/13 01:57:31 mdw Exp $
*
* Low-level multiprecision arithmetic
*
/*----- Revision history --------------------------------------------------*
*
* $Log: mpx.c,v $
+ * Revision 1.3 1999/11/13 01:57:31 mdw
+ * Remove stray debugging code.
+ *
* Revision 1.2 1999/11/13 01:50:59 mdw
* Multiprecision routines finished and tested.
*
/* --- Work out the square at this point in the proceedings --- */
{
- mpw d = *dvv;
mpd x = (mpd)a * (mpd)a + *dvv;
*dvv++ = MPW(x);
c = MPW(x >> MPW_BITS);
avv++;
while (dvv < dvl && avv < avl) {
- mpw aa = *avv;
mpd x = (mpd)a * (mpd)*avv++;
mpd y = ((x << 1) & MPW_MAX) + c + *dvv;
c = (x >> (MPW_BITS - 1)) + (y >> MPW_BITS);