X-Git-Url: https://git.distorted.org.uk/u/mdw/catacomb/blobdiff_plain/c76161cce326aee7a0835c095cd99b694db3d3d4..45c0fd363937c6e9b05da04a9167e9912c05ca0c:/mp-jacobi.c diff --git a/mp-jacobi.c b/mp-jacobi.c index 0e187cf..d0a67b6 100644 --- a/mp-jacobi.c +++ b/mp-jacobi.c @@ -1,13 +1,13 @@ /* -*-c-*- * - * $Id: mp-jacobi.c,v 1.4 2000/12/06 20:31:33 mdw Exp $ + * $Id$ * * Compute Jacobi symbol * * (c) 1999 Straylight/Edgeware */ -/*----- Licensing notice --------------------------------------------------* +/*----- Licensing notice --------------------------------------------------* * * This file is part of Catacomb. * @@ -15,35 +15,18 @@ * 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, * MA 02111-1307, USA. */ -/*----- Revision history --------------------------------------------------* - * - * $Log: mp-jacobi.c,v $ - * Revision 1.4 2000/12/06 20:31:33 mdw - * Add assertion to prevent crapness. - * - * Revision 1.3 2000/07/20 17:14:34 mdw - * Simplify by using @mp_odd@. - * - * Revision 1.2 1999/12/10 23:19:02 mdw - * Improve error-checking. - * - * Revision 1.1 1999/11/22 20:50:37 mdw - * Add support for computing Jacobi symbols. - * - */ - /*----- Header files ------------------------------------------------------*/ #include "mp.h" @@ -67,7 +50,7 @@ int mp_jacobi(mp *a, mp *n) { int s = 1; - assert(((void)"n must be odd in mp_jacobi", MP_LEN(n) && (n->v[0] & 1))); + assert(MP_ODDP(n)); /* --- Take copies of the arguments --- */ @@ -83,7 +66,7 @@ int mp_jacobi(mp *a, mp *n) /* --- Some simple special cases --- */ MP_SHRINK(a); - if (MP_LEN(a) == 0) { + if (MP_ZEROP(a)) { s = 0; goto done; }