new release
[gnupg.git] / mpi / mpih-sub.c
1 /* mpihelp-sub.c  -  MPI helper functions
2  *      Copyright (C) 1998 Free Software Foundation, Inc.
3  *      Copyright (C) 1994, 1996 Free Software Foundation, Inc.
4  *
5  * This file is part of GNUPG.
6  *
7  * GNUPG is free software; you can redistribute it and/or modify
8  * it under the terms of the GNU General Public License as published by
9  * the Free Software Foundation; either version 2 of the License, or
10  * (at your option) any later version.
11  *
12  * GNUPG is distributed in the hope that it will be useful,
13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
15  * GNU General Public License for more details.
16  *
17  * You should have received a copy of the GNU General Public License
18  * along with this program; if not, write to the Free Software
19  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA
20  *
21  * Note: This code is heavily based on the GNU MP Library.
22  *       Actually it's the same code with only minor changes in the
23  *       way the data is stored; this is to support the abstraction
24  *       of an optional secure memory allocation which may be used
25  *       to avoid revealing of sensitive data due to paging etc.
26  *       The GNU MP Library itself is published under the LGPL;
27  *       however I decided to publish this code under the plain GPL.
28  */
29
30 #include <config.h>
31 #include <stdio.h>
32 #include <stdlib.h>
33
34 #include "mpi-internal.h"
35
36 #ifndef __GNUC__
37 mpi_limb_t
38 mpihelp_sub_1(mpi_ptr_t res_ptr,  mpi_ptr_t s1_ptr,
39               mpi_size_t s1_size, mpi_limb_t s2_limb )
40 {
41     mpi_limb_t x;
42
43     x = *s1_ptr++;
44     s2_limb = x - s2_limb;
45     *res_ptr++ = s2_limb;
46     if( s2_limb > x ) {
47         while( --s1_size ) {
48             x = *s1_ptr++;
49             *res_ptr++ = x - 1;
50             if( x )
51                 goto leave;
52         }
53         return 1;
54     }
55
56   leave:
57     if( res_ptr != s1_ptr ) {
58         mpi_size_t i;
59         for( i=0; i < s1_size-1; i++ )
60             res_ptr[i] = s1_ptr[i];
61     }
62     return 0;
63 }
64
65
66 mpi_limb_t
67 mpihelp_sub( mpi_ptr_t res_ptr, mpi_ptr_t s1_ptr, mpi_size_t s1_size,
68                                 mpi_ptr_t s2_ptr, mpi_size_t s2_size)
69 {
70     mpi_limb_t cy = 0;
71
72     if( s2_size )
73         cy = mpihelp_sub_n(res_ptr, s1_ptr, s2_ptr, s2_size);
74
75     if( s1_size - s2_size )
76         cy = mpihelp_sub_1(res_ptr + s2_size, s1_ptr + s2_size,
77                                       s1_size - s2_size, cy);
78     return cy;
79 }
80 #endif
81