mp(3)							     Library Functions Manual							     mp(3)

NAME

       madd,  msub,  mult,  mdiv, gcd, invert, rpow, msqrt, mcmp, min, omin, fmin, m_in, mout, omout, fmout, m_out, sdiv, itom - Performs multiple
       precision integer arithmetic

LIBRARY

       Object Code Library (libmp.a)

SYNOPSIS

       #include <mp.h> #include <stdio.h>
       typedef struct mint { int len; short *val; } MINT; madd(      MINT *a,	   MINT *b,	 MINT *c );
       msub(	  MINT *a,	MINT *b,      MINT *c );
       mult(	  MINT *a,	MINT *b,      MINT *c );
       mdiv(	  MINT *a,	MINT *b,      MINT *q,	    MINT *r );
       pow(	 MINT *a,      MINT *b,      MINT *m,	   MINT *c );
       gcd(	 MINT *a,      MINT *b,      MINT *c );
       invert(	    MINT *a,	  MINT *b,	MINT *c );
       rpow(	  MINT *a,	int n,	    MINT *c );
       msqrt(	   MINT *a,	 MINT *b,      MINT *r );
       mcmp(	  MINT *a,	MINT *b );
       move(	  MINT *a,	MINT *b );
       min(	 MINT *a );
       omin(	  MINT *a );
       fmin(	  MINT *a,	FILE *f );
       m_in(	  MINT *a,	int n,	    FILE *f );
       mout(	  MINT *a );
       omout(	   MINT *a );

       fmout(	   MINT *a,	 FILE *f );
       m_out(	   MINT *a,	 int n,      FILE *f );
       sdiv(	  MINT *a,	short n,      MINT *q,	    short *r );
       *itom(	   short n );

DESCRIPTION

       These functions perform arithmetic on integers of arbitrary length.  The integers are stored using the defined type MINT.   Pointers  to  a
       MINT  can  be  initialized using the itom() function, which sets the initial value to n.  After that, space is managed automatically by the
       routines.

       The madd(), msub() , and mult() functions assign to c the sum, difference, and product, respectively, of a and b.

       The mdiv() function assigns to q and r the quotient and remainder obtained from dividing a by b.  The sdiv() function is  like  the  mdiv()
       function except that the divisor is a short integer n and the remainder is placed in a short integer whose address is given as r.

       The msqrt() function produces the integer square root of a in b and places the remainder in r.

       The  rpow()  function calculates in c the value of a raised to the (``regular'' integral) power n, while the pow() function calculates this
       with a full multiple precision exponent b and the result is reduced modulo m.

       The gcd() function returns the greatest common denominator of a and b in c, and the invert() function computes c such that a*c mod b  =	1,
       for a and b relatively prime.

       The mcmp() function returns a negative, zero, or positive integer value when a is less than, equal to, or greater than b, respectively.

       The move() function copies a to b.

       The  min()  and mout() functions do decimal input and output while the omin() and omout() functions do octal input and output.  More gener-
       ally, the fmin() and fmout() functions do decimal input and output using file f, and m_in() and m_out do input and  output  with  arbitrary
       radix n.

       On input, records should have the form of strings of digits terminated by a newline; output records have a similar form.

NOTES

       Programs which use the multiple-precision arithmetic library must be compiled with -lmp.  delim off

																	     mp(3)