4.3BSD-UWisc/man/cat3/mp.3x

MP(3X)              UNIX Programmer's Manual               MP(3X)



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

SYNOPSIS
     #include <mp.h>
     #include <stdio.h>

     typedef struct mint { int len; short *val; } MINT;

     madd(a, b, c)
     msub(a, b, c)
     mult(a, b, c)
     mdiv(a, b, q, r)
     pow(a, b, m, c)
     gcd(a, b, c)
     invert(a, b, c)
     rpow(a, n, c)
     msqrt(a, b, r)
     mcmp(a, b)
     move(a, b)
     min(a)
     omin(a)
     fmin(a, f)
     m_in(a, n, f)
     mout(a)
     omout(a)
     fmout(a, f)
     m_out(a, n, f)
     MINT *a, *b, *c, *m, *q, *r;
     FILE *f;
     int n;

     sdiv(a, n, q, r)
     MINT *a, *q;
     short n;
     short *r;

     MINT *itom(n)

DESCRIPTION
     These routines perform arithmetic on integers of arbitrary
     length.  The integers are stored using the defined type
     _M_I_N_T.  Pointers to a _M_I_N_T can be initialized using the func-
     tion _i_t_o_m which sets the initial value to _n.  After that,
     space is managed automatically by the routines.

     _m_a_d_d, _m_s_u_b and _m_u_l_t assign to _c the sum, difference and pro-
     duct, respectively, of _a and _b.  _m_d_i_v assigns to _q and _r the
     quotient and remainder obtained from dividing _a by _b.  _s_d_i_v
     is like _m_d_i_v except that the divisor is a short integer _n



Printed 12/27/86          June 4, 1986                          1






MP(3X)              UNIX Programmer's Manual               MP(3X)



     and the remainder is placed in a short whose address is
     given as _r.  _m_s_q_r_t produces the integer square root of _a in
     _b and places the remainder in _r.  _r_p_o_w calculates in _c the
     value of _a raised to the (``regular'' integral) power _n,
     while _p_o_w calculates this with a full multiple precision
     exponent _b and the result is reduced modulo _m.  _g_c_d returns
     the greatest common denominator of _a and _b in _c, and _i_n_v_e_r_t
     computes _c such that _a*_c mod _b = 1, for _a and _b relatively
     prime.  _m_c_m_p returns a negative, zero or positive integer
     value when _a is less than, equal to or greater than _b,
     respectively.  _m_o_v_e copies _a to _b.  _m_i_n and _m_o_u_t do decimal
     input and output while _o_m_i_n and _o_m_o_u_t do octal input and
     output.  More generally, _f_m_i_n and _f_m_o_u_t do decimal input and
     output using file _f, and _m__i_n and _m__o_u_t do I/O with arbi-
     trary radix _n.  On input, records should have the form of
     strings of digits terminated by a newline; output records
     have a similar form.

     Programs which use the multiple-precision arithmetic library
     must be loaded using the loader flag -_l_m_p.

FILES
     /usr/include/mp.h   include file
     /usr/lib/libmp.a    object code library

SEE ALSO
     dc(1), bc(1)

DIAGNOSTICS
     Illegal operations and running out of memory produce mes-
     sages and core images.

BUGS
     Bases for input and output should be <= 10.

     _d_c(1) and _b_c(1) don't use this library.

     The input and output routines are a crock.

     _p_o_w is also the name of a standard math library routine.















Printed 12/27/86          June 4, 1986                          2