mpi.h
上传用户:lyxiangda
上传日期:2007-01-12
资源大小:3042k
文件大小:10k
- /*
- * mpi.h
- *
- * Arbitrary precision integer arithmetic library
- *
- * The contents of this file are subject to the Mozilla Public
- * License Version 1.1 (the "License"); you may not use this file
- * except in compliance with the License. You may obtain a copy of
- * the License at http://www.mozilla.org/MPL/
- *
- * Software distributed under the License is distributed on an "AS
- * IS" basis, WITHOUT WARRANTY OF ANY KIND, either express or
- * implied. See the License for the specific language governing
- * rights and limitations under the License.
- *
- * The Original Code is the MPI Arbitrary Precision Integer Arithmetic
- * library.
- *
- * The Initial Developer of the Original Code is Michael J. Fromberger.
- * Portions created by Michael J. Fromberger are
- * Copyright (C) 1998, 1999, 2000 Michael J. Fromberger.
- * All Rights Reserved.
- *
- * Contributor(s):
- * Netscape Communications Corporation
- *
- * Alternatively, the contents of this file may be used under the
- * terms of the GNU General Public License Version 2 or later (the
- * "GPL"), in which case the provisions of the GPL are applicable
- * instead of those above. If you wish to allow use of your
- * version of this file only under the terms of the GPL and not to
- * allow others to use your version of this file under the MPL,
- * indicate your decision by deleting the provisions above and
- * replace them with the notice and other provisions required by
- * the GPL. If you do not delete the provisions above, a recipient
- * may use your version of this file under either the MPL or the
- * GPL.
- *
- * $Id: mpi.h,v 1.14.2.1 2000/11/21 03:32:40 nelsonb%netscape.com Exp $
- */
- #ifndef _H_MPI_
- #define _H_MPI_
- #include "mpi-config.h"
- #if MP_DEBUG
- #undef MP_IOFUNC
- #define MP_IOFUNC 1
- #endif
- #if MP_IOFUNC
- #include <stdio.h>
- #include <ctype.h>
- #endif
- #include <limits.h>
- #ifdef macintosh
- #include <Types.h>
- #else
- #include <sys/types.h>
- #endif
- #define MP_NEG 1
- #define MP_ZPOS 0
- #define MP_OKAY 0 /* no error, all is well */
- #define MP_YES 0 /* yes (boolean result) */
- #define MP_NO -1 /* no (boolean result) */
- #define MP_MEM -2 /* out of memory */
- #define MP_RANGE -3 /* argument out of range */
- #define MP_BADARG -4 /* invalid parameter */
- #define MP_UNDEF -5 /* answer is undefined */
- #define MP_LAST_CODE MP_UNDEF
- typedef unsigned int mp_sign;
- typedef unsigned int mp_size;
- typedef int mp_err;
- #define MP_32BIT_MAX 4294967295U
- #if !defined(ULONG_MAX)
- #error "ULONG_MAX not defined"
- #elif !defined(UINT_MAX)
- #error "UINT_MAX not defined"
- #elif !defined(USHRT_MAX)
- #error "USHRT_MAX not defined"
- #endif
- #if !defined(MP_USE_UINT_DIGIT) && ULONG_MAX > MP_32BIT_MAX
- typedef unsigned long mp_digit;
- #define MP_DIGIT_MAX ULONG_MAX
- #define MP_HALF_DIGIT_MAX UINT_MAX
- #undef MP_NO_MP_WORD
- #define MP_NO_MP_WORD 1
- #else
- typedef unsigned int mp_digit;
- #define MP_DIGIT_MAX UINT_MAX
- #define MP_HALF_DIGIT_MAX USHRT_MAX
- #endif
- #ifndef MP_NO_MP_WORD
- #if defined(ULONG_LONG_MAX) /* GCC, HPUX */
- #define MP_ULONG_LONG_MAX ULONG_LONG_MAX
- #elif defined(ULLONG_MAX) /* Solaris */
- #define MP_ULONG_LONG_MAX ULLONG_MAX
- /* MP_ULONG_LONG_MAX was defined to be ULLONG_MAX */
- #elif defined(ULONGLONG_MAX) /* IRIX, AIX */
- #define MP_ULONG_LONG_MAX ULONGLONG_MAX
- #endif
- #if defined(MP_ULONG_LONG_MAX)
- #if (MP_ULONG_LONG_MAX > UINT_MAX) || defined(SOLARIS)
- #if MP_ULONG_LONG_MAX == ULONG_MAX || (defined(SOLARIS) && defined(NSS_USE_64))
- typedef unsigned long mp_word;
- typedef long mp_sword;
- #define MP_WORD_MAX ULONG_MAX
- #else
- typedef unsigned long long mp_word;
- typedef long long mp_sword;
- #define MP_WORD_MAX MP_ULONG_LONG_MAX
- #endif
- #else
- /* MP_ULONG_LONG_MAX <= UINT_MAX */
- #define MP_NO_MP_WORD 1
- #endif
- #else
- /* MP_ULONG_LONG_MAX not defined */
- #define MP_NO_MP_WORD 1
- #endif
- #endif /* !MP_NO_MP_WORD */
- #if !defined(MP_WORD_MAX) && defined(MP_DEFINE_SMALL_WORD)
- typedef unsigned int mp_word;
- typedef int mp_sword;
- #define MP_WORD_MAX UINT_MAX
- #endif
- #define MP_DIGIT_BIT (CHAR_BIT*sizeof(mp_digit))
- #define MP_WORD_BIT (CHAR_BIT*sizeof(mp_word))
- #define MP_RADIX (1+(mp_word)MP_DIGIT_MAX)
- #define MP_HALF_DIGIT_BIT (MP_DIGIT_BIT/2)
- #define MP_HALF_RADIX (1+(mp_digit)MP_HALF_DIGIT_MAX)
- /* MP_HALF_RADIX really ought to be called MP_SQRT_RADIX, but it's named
- ** MP_HALF_RADIX because it's the radix for MP_HALF_DIGITs, and it's
- ** consistent with the other _HALF_ names.
- */
- #if MP_DIGIT_MAX == USHRT_MAX
- #define MP_DIGIT_FMT "%04X" /* printf() format for 1 digit */
- #elif MP_DIGIT_MAX == UINT_MAX
- #define MP_DIGIT_FMT "%08X" /* printf() format for 1 digit */
- #elif MP_DIGIT_MAX == ULONG_MAX
- #define MP_DIGIT_FMT "%016lX" /* printf() format for 1 digit */
- #else
- #define MP_DIGIT_FMT "%016llX" /* printf() format for 1 digit */
- #endif
- /* Macros for accessing the mp_int internals */
- #define MP_SIGN(MP) ((MP)->sign)
- #define MP_USED(MP) ((MP)->used)
- #define MP_ALLOC(MP) ((MP)->alloc)
- #define MP_DIGITS(MP) ((MP)->dp)
- #define MP_DIGIT(MP,N) (MP)->dp[(N)]
- /* This defines the maximum I/O base (minimum is 2) */
- #define MP_MAX_RADIX 64
- typedef struct {
- mp_sign sign; /* sign of this quantity */
- mp_size alloc; /* how many digits allocated */
- mp_size used; /* how many digits used */
- mp_digit *dp; /* the digits themselves */
- } mp_int;
- /* Default precision */
- mp_size mp_get_prec(void);
- void mp_set_prec(mp_size prec);
- /* Memory management */
- mp_err mp_init(mp_int *mp);
- mp_err mp_init_size(mp_int *mp, mp_size prec);
- mp_err mp_init_copy(mp_int *mp, const mp_int *from);
- mp_err mp_copy(const mp_int *from, mp_int *to);
- void mp_exch(mp_int *mp1, mp_int *mp2);
- void mp_clear(mp_int *mp);
- void mp_zero(mp_int *mp);
- void mp_set(mp_int *mp, mp_digit d);
- mp_err mp_set_int(mp_int *mp, long z);
- #define mp_set_long(mp,z) mp_set_int(mp,z)
- mp_err mp_set_ulong(mp_int *mp, unsigned long z);
- /* Single digit arithmetic */
- mp_err mp_add_d(const mp_int *a, mp_digit d, mp_int *b);
- mp_err mp_sub_d(const mp_int *a, mp_digit d, mp_int *b);
- mp_err mp_mul_d(const mp_int *a, mp_digit d, mp_int *b);
- mp_err mp_mul_2(const mp_int *a, mp_int *c);
- mp_err mp_div_d(const mp_int *a, mp_digit d, mp_int *q, mp_digit *r);
- mp_err mp_div_2(const mp_int *a, mp_int *c);
- mp_err mp_expt_d(const mp_int *a, mp_digit d, mp_int *c);
- /* Sign manipulations */
- mp_err mp_abs(const mp_int *a, mp_int *b);
- mp_err mp_neg(const mp_int *a, mp_int *b);
- /* Full arithmetic */
- mp_err mp_add(const mp_int *a, const mp_int *b, mp_int *c);
- mp_err mp_sub(const mp_int *a, const mp_int *b, mp_int *c);
- mp_err mp_mul(const mp_int *a, const mp_int *b, mp_int *c);
- #if MP_SQUARE
- mp_err mp_sqr(const mp_int *a, mp_int *b);
- #else
- #define mp_sqr(a, b) mp_mul(a, a, b)
- #endif
- mp_err mp_div(const mp_int *a, const mp_int *b, mp_int *q, mp_int *r);
- mp_err mp_div_2d(const mp_int *a, mp_digit d, mp_int *q, mp_int *r);
- mp_err mp_expt(mp_int *a, mp_int *b, mp_int *c);
- mp_err mp_2expt(mp_int *a, mp_digit k);
- mp_err mp_sqrt(const mp_int *a, mp_int *b);
- /* Modular arithmetic */
- #if MP_MODARITH
- mp_err mp_mod(const mp_int *a, const mp_int *m, mp_int *c);
- mp_err mp_mod_d(const mp_int *a, mp_digit d, mp_digit *c);
- mp_err mp_addmod(const mp_int *a, const mp_int *b, const mp_int *m, mp_int *c);
- mp_err mp_submod(const mp_int *a, const mp_int *b, const mp_int *m, mp_int *c);
- mp_err mp_mulmod(const mp_int *a, const mp_int *b, const mp_int *m, mp_int *c);
- #if MP_SQUARE
- mp_err mp_sqrmod(const mp_int *a, const mp_int *m, mp_int *c);
- #else
- #define mp_sqrmod(a, m, c) mp_mulmod(a, a, m, c)
- #endif
- mp_err mp_exptmod(const mp_int *a, const mp_int *b, const mp_int *m, mp_int *c);
- mp_err mp_exptmod_d(const mp_int *a, mp_digit d, const mp_int *m, mp_int *c);
- #endif /* MP_MODARITH */
- /* Comparisons */
- int mp_cmp_z(const mp_int *a);
- int mp_cmp_d(const mp_int *a, mp_digit d);
- int mp_cmp(const mp_int *a, const mp_int *b);
- int mp_cmp_mag(mp_int *a, mp_int *b);
- int mp_cmp_int(const mp_int *a, long z);
- int mp_isodd(const mp_int *a);
- int mp_iseven(const mp_int *a);
- /* Number theoretic */
- #if MP_NUMTH
- mp_err mp_gcd(mp_int *a, mp_int *b, mp_int *c);
- mp_err mp_lcm(mp_int *a, mp_int *b, mp_int *c);
- mp_err mp_xgcd(const mp_int *a, const mp_int *b, mp_int *g, mp_int *x, mp_int *y);
- mp_err mp_invmod(const mp_int *a, const mp_int *m, mp_int *c);
- mp_err mp_invmod_xgcd(const mp_int *a, const mp_int *m, mp_int *c);
- #endif /* end MP_NUMTH */
- /* Input and output */
- #if MP_IOFUNC
- void mp_print(mp_int *mp, FILE *ofp);
- #endif /* end MP_IOFUNC */
- /* Base conversion */
- mp_err mp_read_raw(mp_int *mp, char *str, int len);
- int mp_raw_size(mp_int *mp);
- mp_err mp_toraw(mp_int *mp, char *str);
- mp_err mp_read_radix(mp_int *mp, const char *str, int radix);
- int mp_radix_size(mp_int *mp, int radix);
- mp_err mp_toradix(mp_int *mp, char *str, int radix);
- int mp_tovalue(char ch, int r);
- #define mp_tobinary(M, S) mp_toradix((M), (S), 2)
- #define mp_tooctal(M, S) mp_toradix((M), (S), 8)
- #define mp_todecimal(M, S) mp_toradix((M), (S), 10)
- #define mp_tohex(M, S) mp_toradix((M), (S), 16)
- /* Error strings */
- const char *mp_strerror(mp_err ec);
- /* Octet string conversion functions */
- mp_err mp_read_unsigned_octets(mp_int *mp, const unsigned char *str, mp_size len);
- int mp_unsigned_octet_size(const mp_int *mp);
- mp_err mp_to_unsigned_octets(const mp_int *mp, unsigned char *str, mp_size maxlen);
- mp_err mp_to_signed_octets(const mp_int *mp, unsigned char *str, mp_size maxlen);
- mp_err mp_to_fixlen_octets(const mp_int *mp, unsigned char *str, mp_size len);
- /* Miscellaneous */
- mp_size mp_trailing_zeros(const mp_int *mp);
- #define MP_CHECKOK(x) if (MP_OKAY > (res = (x))) goto CLEANUP
- #define MP_CHECKERR(x) if (MP_OKAY > (res = (x))) goto CLEANUP
- #if defined(MP_API_COMPATIBLE)
- #define NEG MP_NEG
- #define ZPOS MP_ZPOS
- #define DIGIT_MAX MP_DIGIT_MAX
- #define DIGIT_BIT MP_DIGIT_BIT
- #define DIGIT_FMT MP_DIGIT_FMT
- #define RADIX MP_RADIX
- #define MAX_RADIX MP_MAX_RADIX
- #define SIGN(MP) MP_SIGN(MP)
- #define USED(MP) MP_USED(MP)
- #define ALLOC(MP) MP_ALLOC(MP)
- #define DIGITS(MP) MP_DIGITS(MP)
- #define DIGIT(MP,N) MP_DIGIT(MP,N)
- #if MP_ARGCHK == 1
- #define ARGCHK(X,Y) {if(!(X)){return (Y);}}
- #elif MP_ARGCHK == 2
- #include <assert.h>
- #define ARGCHK(X,Y) assert(X)
- #else
- #define ARGCHK(X,Y) /* */
- #endif
- #endif /* defined MP_API_COMPATIBLE */
- #endif /* end _H_MPI_ */