intermediate release
[gnupg.git] / mpi / mpiutil.c
1 /* mpiutilac  -  Utility functions for MPI
2  *      Copyright (C) 1998 Free Software Foundation, Inc.
3  *
4  * This file is part of GNUPG.
5  *
6  * GNUPG is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 2 of the License, or
9  * (at your option) any later version.
10  *
11  * GNUPG is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with this program; if not, write to the Free Software
18  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA
19  */
20
21 #include <config.h>
22 #include <stdio.h>
23 #include <stdlib.h>
24 #include <assert.h>
25
26 #include "mpi.h"
27 #include "mpi-internal.h"
28 #include "memory.h"
29 #include "util.h"
30
31
32 #ifdef M_DEBUG
33   #undef mpi_alloc
34   #undef mpi_alloc_secure
35   #undef mpi_free
36 #endif
37
38 /****************
39  * fixme: It was a bad idea to use the number of limbs to allocate
40  *        because on a alpha the limbs are large but we normally need
41  *        integers of n bits - So we should chnage this to bits (or bytes).
42  *
43  *        But mpi_alloc is used in a lot of places :-)
44  */
45 MPI
46 #ifdef M_DEBUG
47 mpi_debug_alloc( unsigned nlimbs, const char *info )
48 #else
49 mpi_alloc( unsigned nlimbs )
50 #endif
51 {
52     MPI a;
53
54     if( DBG_MEMORY )
55         log_debug("mpi_alloc(%u)\n", nlimbs*BITS_PER_MPI_LIMB );
56   #ifdef M_DEBUG
57     a = m_debug_alloc( sizeof *a, info );
58     a->d = nlimbs? mpi_debug_alloc_limb_space( nlimbs, 0, info ) : NULL;
59   #else
60     a = m_alloc( sizeof *a );
61     a->d = nlimbs? mpi_alloc_limb_space( nlimbs, 0 ) : NULL;
62   #endif
63     a->alloced = nlimbs;
64     a->nlimbs = 0;
65     a->sign = 0;
66     a->flags = 0;
67     a->nbits = 0;
68     return a;
69 }
70
71 void
72 mpi_m_check( MPI a )
73 {
74     m_check(a);
75     m_check(a->d);
76 }
77
78 MPI
79 #ifdef M_DEBUG
80 mpi_debug_alloc_secure( unsigned nlimbs, const char *info )
81 #else
82 mpi_alloc_secure( unsigned nlimbs )
83 #endif
84 {
85     MPI a;
86
87     if( DBG_MEMORY )
88         log_debug("mpi_alloc_secure(%u)\n", nlimbs*BITS_PER_MPI_LIMB );
89   #ifdef M_DEBUG
90     a = m_debug_alloc( sizeof *a, info );
91     a->d = nlimbs? mpi_debug_alloc_limb_space( nlimbs, 1, info ) : NULL;
92   #else
93     a = m_alloc( sizeof *a );
94     a->d = nlimbs? mpi_alloc_limb_space( nlimbs, 1 ) : NULL;
95   #endif
96     a->alloced = nlimbs;
97     a->flags = 1;
98     a->nlimbs = 0;
99     a->sign = 0;
100     a->nbits = 0;
101     return a;
102 }
103
104
105 mpi_ptr_t
106 #ifdef M_DEBUG
107 mpi_debug_alloc_limb_space( unsigned nlimbs, int secure, const char *info )
108 #else
109 mpi_alloc_limb_space( unsigned nlimbs, int secure )
110 #endif
111 {
112     size_t len = nlimbs * sizeof(mpi_limb_t);
113
114     if( DBG_MEMORY )
115         log_debug("mpi_alloc_limb_space(%u)\n", (unsigned)len*8 );
116   #ifdef M_DEBUG
117     return secure? m_debug_alloc_secure(len, info):m_debug_alloc( len, info );
118   #else
119     return secure? m_alloc_secure( len ):m_alloc( len );
120   #endif
121 }
122
123 void
124 #ifdef M_DEBUG
125 mpi_debug_free_limb_space( mpi_ptr_t a, const char *info )
126 #else
127 mpi_free_limb_space( mpi_ptr_t a )
128 #endif
129 {
130     if( !a )
131         return;
132     if( DBG_MEMORY )
133         log_debug("mpi_free_limb_space of size %lu\n", (ulong)m_size(a)*8 );
134     m_free(a);
135 }
136
137
138 void
139 mpi_assign_limb_space( MPI a, mpi_ptr_t ap, unsigned nlimbs )
140 {
141     mpi_free_limb_space(a->d);
142     a->d = ap;
143     a->alloced = nlimbs;
144 }
145
146
147
148 /****************
149  * Resize the array of A to NLIMBS. the additional space is cleared
150  * (set to 0) [done by m_realloc()]
151  */
152 void
153 #ifdef M_DEBUG
154 mpi_debug_resize( MPI a, unsigned nlimbs, const char *info )
155 #else
156 mpi_resize( MPI a, unsigned nlimbs )
157 #endif
158 {
159     if( nlimbs <= a->alloced )
160         return; /* no need to do it */
161     /* FIXME: add realloc_secure based on a->secure */
162   #ifdef M_DEBUG
163     if( a->d )
164         a->d = m_debug_realloc(a->d, nlimbs * sizeof(mpi_limb_t), info );
165     else
166         a->d = m_debug_alloc_clear( nlimbs * sizeof(mpi_limb_t), info );
167   #else
168     if( a->d )
169         a->d = m_realloc(a->d, nlimbs * sizeof(mpi_limb_t) );
170     else
171         a->d = m_alloc_clear( nlimbs * sizeof(mpi_limb_t) );
172   #endif
173     a->alloced = nlimbs;
174 }
175
176 void
177 mpi_clear( MPI a )
178 {
179     a->nlimbs = 0;
180     a->nbits = 0;
181     a->flags = 0;
182 }
183
184
185 void
186 #ifdef M_DEBUG
187 mpi_debug_free( MPI a, const char *info )
188 #else
189 mpi_free( MPI a )
190 #endif
191 {
192     if( !a )
193         return;
194     if( DBG_MEMORY )
195         log_debug("mpi_free\n" );
196   #ifdef M_DEBUG
197     mpi_debug_free_limb_space(a->d, info);
198   #else
199     mpi_free_limb_space(a->d);
200   #endif
201     if( a->flags & ~3 )
202         log_bug("invalid flag value in mpi\n");
203     m_free(a);
204 }
205
206
207 void
208 mpi_set_secure( MPI a )
209 {
210     mpi_ptr_t ap, bp;
211
212     if( (a->flags & 1) )
213         return;
214     a->flags |= 1;
215     ap = a->d;
216     if( !a->nlimbs ) {
217         assert(!ap);
218         return;
219     }
220   #ifdef M_DEBUG
221     bp = mpi_debug_alloc_limb_space( a->nlimbs, 1, "set_secure" );
222   #else
223     bp = mpi_alloc_limb_space( a->nlimbs, 1 );
224   #endif
225     MPN_COPY( bp, ap, a->nlimbs );
226     a->d = bp;
227   #ifdef M_DEBUG
228     mpi_debug_free_limb_space(ap, "set_secure");
229   #else
230     mpi_free_limb_space(ap);
231   #endif
232 }
233
234
235 /****************
236  * Note: This copy function should not interpret the MPI
237  *       but copy it transparently.
238  */
239 MPI
240 #ifdef M_DEBUG
241 mpi_debug_copy( MPI a, const char *info )
242 #else
243 mpi_copy( MPI a )
244 #endif
245 {
246     int i;
247     MPI b;
248
249     if( a ) {
250       #ifdef M_DEBUG
251         b = mpi_is_secure(a)? mpi_debug_alloc_secure( a->nlimbs, info )
252                             : mpi_debug_alloc( a->nlimbs, info );
253       #else
254         b = mpi_is_secure(a)? mpi_alloc_secure( a->nlimbs )
255                             : mpi_alloc( a->nlimbs );
256       #endif
257         b->nlimbs = a->nlimbs;
258         b->sign = a->sign;
259         b->flags  = a->flags;
260         b->nbits = a->nbits;
261         for(i=0; i < b->nlimbs; i++ )
262             b->d[i] = a->d[i];
263     }
264     else
265         b = NULL;
266     return b;
267 }
268
269
270 void
271 mpi_set( MPI w, MPI u)
272 {
273     mpi_ptr_t wp, up;
274     mpi_size_t usize = u->nlimbs;
275     int usign = u->sign;
276
277     RESIZE_IF_NEEDED(w, usize);
278     wp = w->d;
279     up = u->d;
280     MPN_COPY( wp, up, usize );
281     w->nlimbs = usize;
282     w->nbits = u->nbits;
283     w->flags = u->flags;
284     w->sign = usign;
285 }
286
287
288 void
289 mpi_set_ui( MPI w, unsigned long u)
290 {
291     RESIZE_IF_NEEDED(w, 1);
292     w->d[0] = u;
293     w->nlimbs = u? 1:0;
294     w->sign = 0;
295     w->nbits = 0;
296     w->flags = 0;
297 }
298
299
300 MPI
301 mpi_alloc_set_ui( unsigned long u)
302 {
303   #ifdef M_DEBUG
304     MPI w = mpi_debug_alloc(1,"alloc_set_ui");
305   #else
306     MPI w = mpi_alloc(1);
307   #endif
308     w->d[0] = u;
309     w->nlimbs = u? 1:0;
310     w->sign = 0;
311     return w;
312 }
313
314
315 void
316 mpi_swap( MPI a, MPI b)
317 {
318     struct mpi_struct tmp;
319
320     tmp = *a; *a = *b; *b = tmp;
321 }
322
323