Gather up another utility.
[u/mdw/catacomb] / mpbarrett-exp.h
CommitLineData
f7ea7527 1/* -*-c-*-
2 *
b817bfc6 3 * $Id: mpbarrett-exp.h,v 1.3 2004/04/08 01:36:15 mdw Exp $
f7ea7527 4 *
5 * Exponentiation operations for Barrett reduction
6 *
7 * (c) 2001 Straylight/Edgeware
8 */
9
10/*----- Licensing notice --------------------------------------------------*
11 *
12 * This file is part of Catacomb.
13 *
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.
18 *
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.
23 *
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,
27 * MA 02111-1307, USA.
28 */
29
f7ea7527 30#ifndef CATACOMB_MPBARRETT_EXP_H
31#define CATACOMB_MPBARRETT_EXP_H
32
33#ifdef __cplusplus
34 extern "C" {
35#endif
36
37/*----- Exponentation definitions -----------------------------------------*/
38
39#define EXP_TYPE mp *
40
41#define EXP_COPY(d, x) d = MP_COPY(x)
42#define EXP_DROP(x) MP_DROP(x)
43
44#define EXP_MUL(a, x) do { \
45 mp *t = mp_mul(spare, a, x); \
46 spare = a; \
47 a = mpbarrett_reduce(mb, t, t); \
48} while (0)
49
50#define EXP_SQR(a) do { \
51 mp *t = mp_sqr(spare, a); \
52 spare = a; \
53 a = mpbarrett_reduce(mb, t, t); \
54} while (0)
55
34350998 56#define EXP_FIX(x)
57
f7ea7527 58#define EXP_SETMUL(d, x, y) do { \
59 d = mp_mul(MP_NEW, x, y); \
60 d = mpbarrett_reduce(mb, d, d); \
61} while (0)
62
63#define EXP_SETSQR(d, x) do { \
64 d = mp_sqr(MP_NEW, x); \
65 d = mpbarrett_reduce(mb, d, d); \
66} while (0)
67
68#include "exp.h"
69
70/*----- That's all, folks -------------------------------------------------*/
71
72#ifdef __cplusplus
73 }
74#endif
75
76#endif