Annotation of OpenXM/src/kan96xx/gmp-2.0.2-ssh-2/mpz/jacobi.c, Revision 1.1
1.1 ! takayama 1: /* mpz_jacobi (op1, op2).
! 2: Contributed by Bennet Yee (bsy) at Carnegie-Mellon University
! 3:
! 4: Copyright (C) 1991, 1996 Free Software Foundation, Inc.
! 5:
! 6: This file is part of the GNU MP Library.
! 7:
! 8: The GNU MP Library is free software; you can redistribute it and/or modify
! 9: it under the terms of the GNU Library General Public License as published by
! 10: the Free Software Foundation; either version 2 of the License, or (at your
! 11: option) any later version.
! 12:
! 13: The GNU MP Library is distributed in the hope that it will be useful, but
! 14: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
! 15: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public
! 16: License for more details.
! 17:
! 18: You should have received a copy of the GNU Library General Public License
! 19: along with the GNU MP Library; see the file COPYING.LIB. If not, write to
! 20: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
! 21: MA 02111-1307, USA. */
! 22:
! 23: #include "gmp.h"
! 24:
! 25: /* Precondition: both p and q are positive */
! 26:
! 27: int
! 28: #if __STDC__
! 29: mpz_jacobi (mpz_srcptr pi, mpz_srcptr qi)
! 30: #else
! 31: mpz_jacobi (pi, qi)
! 32: mpz_srcptr pi, qi;
! 33: #endif
! 34: {
! 35: #if GCDCHECK
! 36: int retval;
! 37: mpz_t gcdval;
! 38:
! 39: mpz_init (gcdval);
! 40: mpz_gcd (gcdval, pi, qi);
! 41: if (!mpz_cmp_ui (gcdval, 1L))
! 42: {
! 43: /* J(ab,cb) = J(ab,c)J(ab,b) = J(ab,c)J(0,b) = J(ab,c)*0 */
! 44: retval = 0;
! 45: }
! 46: else
! 47: retval = mpz_legendre (pi, qi);
! 48: mpz_clear (gcdval);
! 49: return retval;
! 50: #else
! 51: return mpz_legendre (pi, qi);
! 52: #endif
! 53: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>