/* -*-c-*-
*
- * $Id: mp-sqrt.c,v 1.5 2004/04/08 01:36:15 mdw Exp $
+ * $Id$
*
* Compute integer square roots
*
* (c) 2000 Straylight/Edgeware
*/
-/*----- Licensing notice --------------------------------------------------*
+/*----- Licensing notice --------------------------------------------------*
*
* This file is part of Catacomb.
*
* it under the terms of the GNU Library General Public License as
* published by the Free Software Foundation; either version 2 of the
* License, or (at your option) any later version.
- *
+ *
* Catacomb is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Library General Public License for more details.
- *
+ *
* You should have received a copy of the GNU Library General Public
* License along with Catacomb; if not, write to the Free
* Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
/* --- Sanity preservation --- */
- assert(((void)"imaginary root in mp_sqrt", !(a->f & MP_NEG)));
+ assert(!MP_NEGP(a));
/* --- Deal with trivial cases --- */
MP_SHRINK(a);
- if (a->v == a->vl) {
+ if (MP_ZEROP(a)) {
mp_drop(d);
return (MP_ZERO);
}
for (;;) {
q = mp_sqr(q, d);
q = mp_sub(q, q, a);
- if (q->v == q->vl)
+ if (MP_ZEROP(q))
break;
- if (q->f & MP_NEG) {
+ if (MP_NEGP(q)) {
r = mp_lsl(r, d, 1);
r->f |= MP_NEG;
- if (MP_CMP(q, <=, r))
+ if (MP_CMP(q, >=, r))
break;
}
mp_div(&r, &q, q, d);
if (!MP_EQ(q, qq)) {
ok = 0;
fputs("\n*** sqrt failed", stderr);
- fputs("\n*** a = ", stderr); mp_writefile(a, stderr, 10);
+ fputs("\n*** a = ", stderr); mp_writefile(a, stderr, 10);
fputs("\n*** result = ", stderr); mp_writefile(q, stderr, 10);
fputs("\n*** expect = ", stderr); mp_writefile(qq, stderr, 10);
fputc('\n', stderr);