3 * $Id: mpbarrett.h,v 1.3 2001/06/16 12:58:47 mdw Exp $
5 * Barrett modular reduction
7 * (c) 1999 Straylight/Edgeware
10 /*----- Licensing notice --------------------------------------------------*
12 * This file is part of Catacomb.
14 * Catacomb is free software; you can redistribute it and/or modify
15 * it under the terms of the GNU Library General Public License as
16 * published by the Free Software Foundation; either version 2 of the
17 * License, or (at your option) any later version.
19 * Catacomb is distributed in the hope that it will be useful,
20 * but WITHOUT ANY WARRANTY; without even the implied warranty of
21 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
22 * GNU Library General Public License for more details.
24 * You should have received a copy of the GNU Library General Public
25 * License along with Catacomb; if not, write to the Free
26 * Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
30 /*----- Revision history --------------------------------------------------*
32 * $Log: mpbarrett.h,v $
33 * Revision 1.3 2001/06/16 12:58:47 mdw
34 * Added simultaneous exponentiation with Barrett reduction.
36 * Revision 1.2 2000/10/08 12:03:44 mdw
37 * (mpbarrett_reduce): Cope with negative numbers.
39 * Revision 1.1 1999/12/10 23:22:00 mdw
40 * Barrett reduction support: works with even moduli.
44 /*----- Notes on Barrett reduction ----------------------------------------*
46 * Barrett reduction is a technique for computing modular residues. Unlike
47 * Montgomery reduction, it doesn't have restrictions on the modulus (except
48 * that it be positive) and doesn't confuse matters by putting an extra
49 * factor all the way through your computation.
51 * It's useful for slightly less heavy-duty work than Montgomery reduction
52 * because the precomputation phase is rather simpler, involving a single
55 * Sometimes it's useful to exponentiate modulo an even number, so there's a
56 * modexp routine provided which uses Barrett reduction rather than
57 * Montgomery reduction. This is handy when you're working on indices in an
58 * even-order cyclic group or something.
61 #ifndef CATACOMB_MPBARRETT_H
62 #define CATACOMB_MPBARRETT_H
68 /*----- Header files ------------------------------------------------------*/
74 /*----- Data structures ---------------------------------------------------*/
76 typedef struct mpbarrett
{
82 /*----- Functions provided ------------------------------------------------*/
84 /* --- @mpbarrett_create@ --- *
86 * Arguments: @mpbarrett *mb@ = pointer to Barrett reduction context
87 * @mp *m@ = modulus to work to
92 * Use: Initializes a Barrett reduction context ready for use.
95 extern void mpbarrett_create(mpbarrett */
*mb*/
, mp */
*m*/
);
97 /* --- @mpbarrett_destroy@ --- *
99 * Arguments: @mpbarrett *mb@ = pointer to Barrett reduction context
103 * Use: Destroys a Barrett reduction context releasing any resources
107 extern void mpbarrett_destroy(mpbarrett */
*mb*/
);
109 /* --- @mpbarrett_reduce@ --- *
111 * Arguments: @mpbarrett *mb@ = pointer to Barrett reduction context
112 * @mp *d@ = destination for result
113 * @mp *m@ = number to reduce
115 * Returns: The residue of @m@ modulo the number in the reduction
118 * Use: Performs an efficient modular reduction.
121 extern mp
*mpbarrett_reduce(mpbarrett */
*mb*/
, mp */
*d*/
, mp */
*m*/
);
123 /* --- @mpbarrett_exp@ --- *
125 * Arguments: @mpbarrett *mb@ = pointer to Barrett reduction context
126 * @mp *d@ = fake destination
130 * Returns: Result, %$a^e \bmod m$%.
133 extern mp
*mpbarrett_exp(mpbarrett */
*mb*/
, mp */
*d*/
, mp */
*a*/
, mp */
*e*/
);
135 /* --- @mpbarrett_mexp@ --- *
137 * Arguments: @mpbarrett *mb@ = pointer to Barrett reduction context
138 * @mp *d@ = fake destination
139 * @mp_expfactor *f@ = pointer to array of factors
140 * @size_t n@ = number of factors supplied
142 * Returns: If the bases are %$g_0, g_1, \ldots, g_{n-1}$% and the
143 * exponents are %$e_0, e_1, \ldots, e_{n-1}$% then the result
146 * %$g_0^{e_0} g_1^{e_1} \ldots g_{n-1}^{e_{n-1}} \bmod m$%
149 extern mp
*mpbarrett_mexp(mpbarrett */
*mb*/
, mp */
*d*/
,
150 mp_expfactor */
*f*/
, size_t /*n*/);
152 /*----- That's all, folks -------------------------------------------------*/