Fix bug dividing small things by large ones.
[u/mdw/catacomb] / mpbarrett-exp.h
CommitLineData
f7ea7527 1/* -*-c-*-
2 *
3 * $Id: mpbarrett-exp.h,v 1.1 2001/06/16 12:58:12 mdw Exp $
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
30/*----- Revision history --------------------------------------------------*
31 *
32 * $Log: mpbarrett-exp.h,v $
33 * Revision 1.1 2001/06/16 12:58:12 mdw
34 * Parameters for generic exponentiation.
35 *
36 */
37
38#ifndef CATACOMB_MPBARRETT_EXP_H
39#define CATACOMB_MPBARRETT_EXP_H
40
41#ifdef __cplusplus
42 extern "C" {
43#endif
44
45/*----- Exponentation definitions -----------------------------------------*/
46
47#define EXP_TYPE mp *
48
49#define EXP_COPY(d, x) d = MP_COPY(x)
50#define EXP_DROP(x) MP_DROP(x)
51
52#define EXP_MUL(a, x) do { \
53 mp *t = mp_mul(spare, a, x); \
54 spare = a; \
55 a = mpbarrett_reduce(mb, t, t); \
56} while (0)
57
58#define EXP_SQR(a) do { \
59 mp *t = mp_sqr(spare, a); \
60 spare = a; \
61 a = mpbarrett_reduce(mb, t, t); \
62} while (0)
63
64#define EXP_SETMUL(d, x, y) do { \
65 d = mp_mul(MP_NEW, x, y); \
66 d = mpbarrett_reduce(mb, d, d); \
67} while (0)
68
69#define EXP_SETSQR(d, x) do { \
70 d = mp_sqr(MP_NEW, x); \
71 d = mpbarrett_reduce(mb, d, d); \
72} while (0)
73
74#include "exp.h"
75
76/*----- That's all, folks -------------------------------------------------*/
77
78#ifdef __cplusplus
79 }
80#endif
81
82#endif