f2ba9e96867935dc624ece52573c174612f72825Dina K Nimeh/* BEGIN CSTYLED */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/*
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * mpi.h
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Arbitrary precision integer arithmetic library
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * ***** BEGIN LICENSE BLOCK *****
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Version: MPL 1.1/GPL 2.0/LGPL 2.1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * The contents of this file are subject to the Mozilla Public License Version
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * 1.1 (the "License"); you may not use this file except in compliance with
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * the License. You may obtain a copy of the License at
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * http://www.mozilla.org/MPL/
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Software distributed under the License is distributed on an "AS IS" basis,
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * for the specific language governing rights and limitations under the
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * License.
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * The Original Code is the MPI Arbitrary Precision Integer Arithmetic library.
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * The Initial Developer of the Original Code is
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Michael J. Fromberger.
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Portions created by the Initial Developer are Copyright (C) 1998
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * the Initial Developer. All Rights Reserved.
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Contributor(s):
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Netscape Communications Corporation
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Alternatively, the contents of this file may be used under the terms of
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * either the GNU General Public License Version 2 or later (the "GPL"), or
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * in which case the provisions of the GPL or the LGPL are applicable instead
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * of those above. If you wish to allow use of your version of this file only
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * under the terms of either the GPL or the LGPL, and not to allow others to
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * use your version of this file under the terms of the MPL, indicate your
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * decision by deleting the provisions above and replace them with the notice
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * and other provisions required by the GPL or the LGPL. If you do not delete
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * the provisions above, a recipient may use your version of this file under
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * the terms of any one of the MPL, the GPL or the LGPL.
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * ***** END LICENSE BLOCK ***** */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/*
f2ba9e96867935dc624ece52573c174612f72825Dina K Nimeh * Copyright (c) 2007, 2010, Oracle and/or its affiliates. All rights reserved.
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers *
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Sun elects to use this software under the MPL license.
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef _MPI_H
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define _MPI_H
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* $Id: mpi.h,v 1.22 2004/04/27 23:04:36 gerv%gerv.net Exp $ */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include "mpi-config.h"
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include <sys/param.h>
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifdef _KERNEL
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include <sys/debug.h>
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include <sys/systm.h>
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define assert ASSERT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define labs(a) (a >= 0 ? a : -a)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define UCHAR_MAX 255
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define memset(s, c, n) bzero(s, n)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define memcpy(a,b,c) bcopy((caddr_t)b, (caddr_t)a, c)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/*
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers * Generic #define's to cover missing things in the kernel
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef isdigit
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define isdigit(x) ((x) >= '0' && (x) <= '9')
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef isupper
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define isupper(x) (((unsigned)(x) >= 'A') && ((unsigned)(x) <= 'Z'))
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef islower
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define islower(x) (((unsigned)(x) >= 'a') && ((unsigned)(x) <= 'z'))
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef isalpha
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define isalpha(x) (isupper(x) || islower(x))
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef toupper
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define toupper(x) (islower(x) ? (x) - 'a' + 'A' : (x))
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef tolower
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define tolower(x) (isupper(x) ? (x) + 'a' - 'A' : (x))
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef isspace
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define isspace(x) (((x) == ' ') || ((x) == '\r') || ((x) == '\n') || \
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers ((x) == '\t') || ((x) == '\b'))
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif /* _KERNEL */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if MP_DEBUG
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_IOFUNC
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_IOFUNC 1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if MP_IOFUNC
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include <stdio.h>
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include <ctype.h>
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef _KERNEL
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include <limits.h>
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if defined(BSDI)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef ULLONG_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if defined( macintosh )
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include <Types.h>
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#elif defined( _WIN32_WCE)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* #include <sys/types.h> What do we need here ?? */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#else
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include <sys/types.h>
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_NEG 1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_ZPOS 0
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_OKAY 0 /* no error, all is well */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_YES 0 /* yes (boolean result) */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_NO -1 /* no (boolean result) */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_MEM -2 /* out of memory */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_RANGE -3 /* argument out of range */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_BADARG -4 /* invalid parameter */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_UNDEF -5 /* answer is undefined */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_LAST_CODE MP_UNDEF
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef unsigned int mp_sign;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef unsigned int mp_size;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef int mp_err;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef int mp_flag;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_32BIT_MAX 4294967295U
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if !defined(ULONG_MAX)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#error "ULONG_MAX not defined"
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#elif !defined(UINT_MAX)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#error "UINT_MAX not defined"
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#elif !defined(USHRT_MAX)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#error "USHRT_MAX not defined"
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if defined(ULONG_LONG_MAX) /* GCC, HPUX */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_ULONG_LONG_MAX ULONG_LONG_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#elif defined(ULLONG_MAX) /* Solaris */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_ULONG_LONG_MAX ULLONG_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* MP_ULONG_LONG_MAX was defined to be ULLONG_MAX */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#elif defined(ULONGLONG_MAX) /* IRIX, AIX */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_ULONG_LONG_MAX ULONGLONG_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* We only use unsigned long for mp_digit iff long is more than 32 bits. */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if !defined(MP_USE_UINT_DIGIT) && ULONG_MAX > MP_32BIT_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef unsigned long mp_digit;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_DIGIT_MAX ULONG_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_DIGIT_FMT "%016lX" /* printf() format for 1 digit */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_HALF_DIGIT_MAX UINT_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_NO_MP_WORD
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_NO_MP_WORD 1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_USE_LONG_DIGIT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_USE_LONG_DIGIT 1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_USE_LONG_LONG_DIGIT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#elif !defined(MP_USE_UINT_DIGIT) && defined(MP_ULONG_LONG_MAX)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef unsigned long long mp_digit;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_DIGIT_MAX MP_ULONG_LONG_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_DIGIT_FMT "%016llX" /* printf() format for 1 digit */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_HALF_DIGIT_MAX UINT_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_NO_MP_WORD
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_NO_MP_WORD 1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_USE_LONG_LONG_DIGIT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_USE_LONG_LONG_DIGIT 1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_USE_LONG_DIGIT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#else
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef unsigned int mp_digit;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_DIGIT_MAX UINT_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_DIGIT_FMT "%08X" /* printf() format for 1 digit */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_HALF_DIGIT_MAX USHRT_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_USE_UINT_DIGIT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_USE_UINT_DIGIT 1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_USE_LONG_LONG_DIGIT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#undef MP_USE_LONG_DIGIT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if !defined(MP_NO_MP_WORD)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if defined(MP_USE_UINT_DIGIT) && \
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers (defined(MP_ULONG_LONG_MAX) || (ULONG_MAX > UINT_MAX))
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if (ULONG_MAX > UINT_MAX)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef unsigned long mp_word;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef long mp_sword;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_WORD_MAX ULONG_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#else
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef unsigned long long mp_word;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef long long mp_sword;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_WORD_MAX MP_ULONG_LONG_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#else
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_NO_MP_WORD 1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif /* !defined(MP_NO_MP_WORD) */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if !defined(MP_WORD_MAX) && defined(MP_DEFINE_SMALL_WORD)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef unsigned int mp_word;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef int mp_sword;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_WORD_MAX UINT_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifndef CHAR_BIT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define CHAR_BIT 8
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_DIGIT_BIT (CHAR_BIT*sizeof(mp_digit))
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_WORD_BIT (CHAR_BIT*sizeof(mp_word))
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_RADIX (1+(mp_word)MP_DIGIT_MAX)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_HALF_DIGIT_BIT (MP_DIGIT_BIT/2)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_HALF_RADIX (1+(mp_digit)MP_HALF_DIGIT_MAX)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* MP_HALF_RADIX really ought to be called MP_SQRT_RADIX, but it's named
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers** MP_HALF_RADIX because it's the radix for MP_HALF_DIGITs, and it's
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers** consistent with the other _HALF_ names.
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers*/
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Macros for accessing the mp_int internals */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_FLAG(MP) ((MP)->flag)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_SIGN(MP) ((MP)->sign)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_USED(MP) ((MP)->used)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_ALLOC(MP) ((MP)->alloc)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_DIGITS(MP) ((MP)->dp)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_DIGIT(MP,N) (MP)->dp[(N)]
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* This defines the maximum I/O base (minimum is 2) */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_MAX_RADIX 64
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowerstypedef struct {
f2ba9e96867935dc624ece52573c174612f72825Dina K Nimeh mp_flag flag; /* KM_SLEEP/KM_NOSLEEP */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers mp_sign sign; /* sign of this quantity */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers mp_size alloc; /* how many digits allocated */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers mp_size used; /* how many digits used */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers mp_digit *dp; /* the digits themselves */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers} mp_int;
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Default precision */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_size mp_get_prec(void);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersvoid mp_set_prec(mp_size prec);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Memory management */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_init(mp_int *mp, int kmflag);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_init_size(mp_int *mp, mp_size prec, int kmflag);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_init_copy(mp_int *mp, const mp_int *from);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_copy(const mp_int *from, mp_int *to);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersvoid mp_exch(mp_int *mp1, mp_int *mp2);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersvoid mp_clear(mp_int *mp);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersvoid mp_zero(mp_int *mp);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersvoid mp_set(mp_int *mp, mp_digit d);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_set_int(mp_int *mp, long z);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define mp_set_long(mp,z) mp_set_int(mp,z)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_set_ulong(mp_int *mp, unsigned long z);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Single digit arithmetic */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_add_d(const mp_int *a, mp_digit d, mp_int *b);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_sub_d(const mp_int *a, mp_digit d, mp_int *b);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_mul_d(const mp_int *a, mp_digit d, mp_int *b);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_mul_2(const mp_int *a, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_div_d(const mp_int *a, mp_digit d, mp_int *q, mp_digit *r);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_div_2(const mp_int *a, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_expt_d(const mp_int *a, mp_digit d, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Sign manipulations */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_abs(const mp_int *a, mp_int *b);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_neg(const mp_int *a, mp_int *b);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Full arithmetic */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_add(const mp_int *a, const mp_int *b, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_sub(const mp_int *a, const mp_int *b, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_mul(const mp_int *a, const mp_int *b, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if MP_SQUARE
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_sqr(const mp_int *a, mp_int *b);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#else
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define mp_sqr(a, b) mp_mul(a, a, b)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_div(const mp_int *a, const mp_int *b, mp_int *q, mp_int *r);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_div_2d(const mp_int *a, mp_digit d, mp_int *q, mp_int *r);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_expt(mp_int *a, mp_int *b, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_2expt(mp_int *a, mp_digit k);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_sqrt(const mp_int *a, mp_int *b);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Modular arithmetic */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if MP_MODARITH
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_mod(const mp_int *a, const mp_int *m, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_mod_d(const mp_int *a, mp_digit d, mp_digit *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_addmod(const mp_int *a, const mp_int *b, const mp_int *m, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_submod(const mp_int *a, const mp_int *b, const mp_int *m, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_mulmod(const mp_int *a, const mp_int *b, const mp_int *m, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if MP_SQUARE
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_sqrmod(const mp_int *a, const mp_int *m, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#else
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define mp_sqrmod(a, m, c) mp_mulmod(a, a, m, c)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_exptmod(const mp_int *a, const mp_int *b, const mp_int *m, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_exptmod_d(const mp_int *a, mp_digit d, const mp_int *m, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif /* MP_MODARITH */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Comparisons */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_cmp_z(const mp_int *a);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_cmp_d(const mp_int *a, mp_digit d);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_cmp(const mp_int *a, const mp_int *b);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_cmp_mag(mp_int *a, mp_int *b);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_cmp_int(const mp_int *a, long z, int kmflag);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_isodd(const mp_int *a);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_iseven(const mp_int *a);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Number theoretic */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if MP_NUMTH
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_gcd(mp_int *a, mp_int *b, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_lcm(mp_int *a, mp_int *b, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_xgcd(const mp_int *a, const mp_int *b, mp_int *g, mp_int *x, mp_int *y);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_invmod(const mp_int *a, const mp_int *m, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_invmod_xgcd(const mp_int *a, const mp_int *m, mp_int *c);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif /* end MP_NUMTH */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Input and output */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if MP_IOFUNC
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersvoid mp_print(mp_int *mp, FILE *ofp);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif /* end MP_IOFUNC */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Base conversion */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_read_raw(mp_int *mp, char *str, int len);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_raw_size(mp_int *mp);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_toraw(mp_int *mp, char *str);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_read_radix(mp_int *mp, const char *str, int radix);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_read_variable_radix(mp_int *a, const char * str, int default_radix);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_radix_size(mp_int *mp, int radix);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_toradix(mp_int *mp, char *str, int radix);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_tovalue(char ch, int r);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define mp_tobinary(M, S) mp_toradix((M), (S), 2)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define mp_tooctal(M, S) mp_toradix((M), (S), 8)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define mp_todecimal(M, S) mp_toradix((M), (S), 10)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define mp_tohex(M, S) mp_toradix((M), (S), 16)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Error strings */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersconst char *mp_strerror(mp_err ec);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Octet string conversion functions */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_read_unsigned_octets(mp_int *mp, const unsigned char *str, mp_size len);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersint mp_unsigned_octet_size(const mp_int *mp);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_to_unsigned_octets(const mp_int *mp, unsigned char *str, mp_size maxlen);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_to_signed_octets(const mp_int *mp, unsigned char *str, mp_size maxlen);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_err mp_to_fixlen_octets(const mp_int *mp, unsigned char *str, mp_size len);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers/* Miscellaneous */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowersmp_size mp_trailing_zeros(const mp_int *mp);
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_CHECKOK(x) if (MP_OKAY > (res = (x))) goto CLEANUP
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MP_CHECKERR(x) if (MP_OKAY > (res = (x))) goto CLEANUP
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if defined(MP_API_COMPATIBLE)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define NEG MP_NEG
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define ZPOS MP_ZPOS
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define DIGIT_MAX MP_DIGIT_MAX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define DIGIT_BIT MP_DIGIT_BIT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define DIGIT_FMT MP_DIGIT_FMT
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define RADIX MP_RADIX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define MAX_RADIX MP_MAX_RADIX
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define FLAG(MP) MP_FLAG(MP)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define SIGN(MP) MP_SIGN(MP)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define USED(MP) MP_USED(MP)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define ALLOC(MP) MP_ALLOC(MP)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define DIGITS(MP) MP_DIGITS(MP)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define DIGIT(MP,N) MP_DIGIT(MP,N)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#if MP_ARGCHK == 1
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define ARGCHK(X,Y) {if(!(X)){return (Y);}}
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#elif MP_ARGCHK == 2
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#ifdef _KERNEL
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define ARGCHK(X,Y) ASSERT(X)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#else
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#include <assert.h>
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define ARGCHK(X,Y) assert(X)
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#else
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#define ARGCHK(X,Y) /* */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif /* defined MP_API_COMPATIBLE */
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers
f9fbec18f5b458b560ecf45d3db8e8bd56bf6942mcpowers#endif /* _MPI_H */
f2ba9e96867935dc624ece52573c174612f72825Dina K Nimeh/* END CSTYLED */