xref: /titanic_51/usr/src/common/mpi/mpprime.h (revision f9fbec18f5b458b560ecf45d3db8e8bd56bf6942)
1*f9fbec18Smcpowers /*
2*f9fbec18Smcpowers  *  mpprime.h
3*f9fbec18Smcpowers  *
4*f9fbec18Smcpowers  *  Utilities for finding and working with prime and pseudo-prime
5*f9fbec18Smcpowers  *  integers
6*f9fbec18Smcpowers  *
7*f9fbec18Smcpowers  * ***** BEGIN LICENSE BLOCK *****
8*f9fbec18Smcpowers  * Version: MPL 1.1/GPL 2.0/LGPL 2.1
9*f9fbec18Smcpowers  *
10*f9fbec18Smcpowers  * The contents of this file are subject to the Mozilla Public License Version
11*f9fbec18Smcpowers  * 1.1 (the "License"); you may not use this file except in compliance with
12*f9fbec18Smcpowers  * the License. You may obtain a copy of the License at
13*f9fbec18Smcpowers  * http://www.mozilla.org/MPL/
14*f9fbec18Smcpowers  *
15*f9fbec18Smcpowers  * Software distributed under the License is distributed on an "AS IS" basis,
16*f9fbec18Smcpowers  * WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
17*f9fbec18Smcpowers  * for the specific language governing rights and limitations under the
18*f9fbec18Smcpowers  * License.
19*f9fbec18Smcpowers  *
20*f9fbec18Smcpowers  * The Original Code is the MPI Arbitrary Precision Integer Arithmetic library.
21*f9fbec18Smcpowers  *
22*f9fbec18Smcpowers  * The Initial Developer of the Original Code is
23*f9fbec18Smcpowers  * Michael J. Fromberger.
24*f9fbec18Smcpowers  * Portions created by the Initial Developer are Copyright (C) 1997
25*f9fbec18Smcpowers  * the Initial Developer. All Rights Reserved.
26*f9fbec18Smcpowers  *
27*f9fbec18Smcpowers  * Contributor(s):
28*f9fbec18Smcpowers  *
29*f9fbec18Smcpowers  * Alternatively, the contents of this file may be used under the terms of
30*f9fbec18Smcpowers  * either the GNU General Public License Version 2 or later (the "GPL"), or
31*f9fbec18Smcpowers  * the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
32*f9fbec18Smcpowers  * in which case the provisions of the GPL or the LGPL are applicable instead
33*f9fbec18Smcpowers  * of those above. If you wish to allow use of your version of this file only
34*f9fbec18Smcpowers  * under the terms of either the GPL or the LGPL, and not to allow others to
35*f9fbec18Smcpowers  * use your version of this file under the terms of the MPL, indicate your
36*f9fbec18Smcpowers  * decision by deleting the provisions above and replace them with the notice
37*f9fbec18Smcpowers  * and other provisions required by the GPL or the LGPL. If you do not delete
38*f9fbec18Smcpowers  * the provisions above, a recipient may use your version of this file under
39*f9fbec18Smcpowers  * the terms of any one of the MPL, the GPL or the LGPL.
40*f9fbec18Smcpowers  *
41*f9fbec18Smcpowers  * ***** END LICENSE BLOCK ***** */
42*f9fbec18Smcpowers /*
43*f9fbec18Smcpowers  * Copyright 2007 Sun Microsystems, Inc.  All rights reserved.
44*f9fbec18Smcpowers  * Use is subject to license terms.
45*f9fbec18Smcpowers  *
46*f9fbec18Smcpowers  * Sun elects to use this software under the MPL license.
47*f9fbec18Smcpowers  */
48*f9fbec18Smcpowers 
49*f9fbec18Smcpowers #ifndef _MP_PRIME_H
50*f9fbec18Smcpowers #define _MP_PRIME_H
51*f9fbec18Smcpowers 
52*f9fbec18Smcpowers #pragma ident	"%Z%%M%	%I%	%E% SMI"
53*f9fbec18Smcpowers 
54*f9fbec18Smcpowers #include "mpi.h"
55*f9fbec18Smcpowers 
56*f9fbec18Smcpowers extern const int prime_tab_size;   /* number of primes available */
57*f9fbec18Smcpowers extern const mp_digit prime_tab[];
58*f9fbec18Smcpowers 
59*f9fbec18Smcpowers /* Tests for divisibility    */
60*f9fbec18Smcpowers mp_err  mpp_divis(mp_int *a, mp_int *b);
61*f9fbec18Smcpowers mp_err  mpp_divis_d(mp_int *a, mp_digit d);
62*f9fbec18Smcpowers 
63*f9fbec18Smcpowers /* Random selection          */
64*f9fbec18Smcpowers mp_err  mpp_random(mp_int *a);
65*f9fbec18Smcpowers mp_err  mpp_random_size(mp_int *a, mp_size prec);
66*f9fbec18Smcpowers 
67*f9fbec18Smcpowers /* Pseudo-primality testing  */
68*f9fbec18Smcpowers mp_err  mpp_divis_vector(mp_int *a, const mp_digit *vec, int size, int *which);
69*f9fbec18Smcpowers mp_err  mpp_divis_primes(mp_int *a, mp_digit *np);
70*f9fbec18Smcpowers mp_err  mpp_fermat(mp_int *a, mp_digit w);
71*f9fbec18Smcpowers mp_err mpp_fermat_list(mp_int *a, const mp_digit *primes, mp_size nPrimes);
72*f9fbec18Smcpowers mp_err  mpp_pprime(mp_int *a, int nt);
73*f9fbec18Smcpowers mp_err mpp_sieve(mp_int *trial, const mp_digit *primes, mp_size nPrimes,
74*f9fbec18Smcpowers 		 unsigned char *sieve, mp_size nSieve);
75*f9fbec18Smcpowers mp_err mpp_make_prime(mp_int *start, mp_size nBits, mp_size strong,
76*f9fbec18Smcpowers 		      unsigned long * nTries);
77*f9fbec18Smcpowers 
78*f9fbec18Smcpowers #endif /* _MP_PRIME_H */
79