See ChangeLog: Mon Jan 24 22:24:38 CET 2000 Werner Koch
[gnupg.git] / mpi / mpi-add.c
index 047a2fa..3dd22eb 100644 (file)
@@ -1,14 +1,15 @@
 /* mpi-add.c  -  MPI functions
- *     Copyright (c) 1997 by Werner Koch (dd9jn)
+ *     Copyright (C) 1998 Free Software Foundation, Inc.
+ *     Copyright (C) 1994, 1996 Free Software Foundation, Inc.
  *
- * This file is part of G10.
+ * This file is part of GnuPG.
  *
- * G10 is free software; you can redistribute it and/or modify
+ * GnuPG is free software; you can redistribute it and/or modify
  * it under the terms of the GNU General Public License as published by
  * the Free Software Foundation; either version 2 of the License, or
  * (at your option) any later version.
  *
- * G10 is distributed in the hope that it will be useful,
+ * GnuPG is distributed in the hope that it will be useful,
  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  * GNU General Public License for more details.
  * You should have received a copy of the GNU General Public License
  * along with this program; if not, write to the Free Software
  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA
+ *
+ * Note: This code is heavily based on the GNU MP Library.
+ *      Actually it's the same code with only minor changes in the
+ *      way the data is stored; this is to support the abstraction
+ *      of an optional secure memory allocation which may be used
+ *      to avoid revealing of sensitive data due to paging etc.
+ *      The GNU MP Library itself is published under the LGPL;
+ *      however I decided to publish this code under the plain GPL.
  */
 
 #include <config.h>
@@ -68,7 +77,7 @@ mpi_add_ui(MPI w, MPI u, unsigned long v )
        else {
            mpihelp_sub_1(wp, up, usize, v);
            /* Size can decrease with at most one limb. */
-           wsize = (usize - (wp[usize-1]? 0:1));
+           wsize = usize - (wp[usize-1]==0);
            wsign = 1;
        }
     }
@@ -85,27 +94,30 @@ mpi_add(MPI w, MPI u, MPI v)
     mpi_size_t usize, vsize, wsize;
     int usign, vsign, wsign;
 
-    usize = u->nlimbs;
-    vsize = v->nlimbs;
-    usign = u->sign;
-    vsign = v->sign;
-
-    if( usize < vsize ) { /* Swap U and V. */
-       { MPI t; t = u; u = v; v = t; }
-       { mpi_size_t t = usize; usize = vsize; vsize = t; }
-       { int t = usign; usign = vsign; vsign = t; }
+    if( u->nlimbs < v->nlimbs ) { /* Swap U and V. */
+       usize = v->nlimbs;
+       usign = v->sign;
+       vsize = u->nlimbs;
+       vsign = u->sign;
+       wsize = usize + 1;
+       RESIZE_IF_NEEDED(w, wsize);
+       /* These must be after realloc (u or v may be the same as w).  */
+       up    = v->d;
+       vp    = u->d;
+    }
+    else {
+       usize = u->nlimbs;
+       usign = u->sign;
+       vsize = v->nlimbs;
+       vsign = v->sign;
+       wsize = usize + 1;
+       RESIZE_IF_NEEDED(w, wsize);
+       /* These must be after realloc (u or v may be the same as w).  */
+       up    = u->d;
+       vp    = v->d;
     }
-
-    /* If not space for w (and possible carry), increase space.  */
-    wsize = usize + 1;
-    if( w->alloced < wsize )
-       mpi_resize(w, wsize);
-    wsign = 0;
-
-    /* These must be after realloc (u or v may be the same as w).  */
-    up = u->d;
-    vp = v->d;
     wp = w->d;
+    wsign = 0;
 
     if( !vsize ) {  /* simple */
        MPN_COPY(wp, up, usize );
@@ -140,7 +152,7 @@ mpi_add(MPI w, MPI u, MPI v)
        wp[usize] = cy;
        wsize = usize + cy;
        if( usign )
-           wsize = 1;
+           wsign = 1;
     }
 
     w->nlimbs = wsize;
@@ -193,7 +205,7 @@ mpi_sub_ui(MPI w, MPI u, unsigned long v )
        else {
            mpihelp_sub_1(wp, up, usize, v);
            /* Size can decrease with at most one limb. */
-           wsize = (usize - (wp[usize-1]? 1:0));
+           wsize = usize - (wp[usize-1]==0);
        }
     }
 
@@ -208,7 +220,7 @@ mpi_sub(MPI w, MPI u, MPI v)
        MPI vv = mpi_copy(v);
        vv->sign = !vv->sign;
        mpi_add( w, u, vv );
-       m_free(vv);
+       mpi_free(vv);
     }
     else {
        /* fixme: this is not thread-save (we temp. modify v) */
@@ -219,3 +231,17 @@ mpi_sub(MPI w, MPI u, MPI v)
 }
 
 
+void
+mpi_addm( MPI w, MPI u, MPI v, MPI m)
+{
+    mpi_add(w, u, v);
+    mpi_fdiv_r( w, w, m );
+}
+
+void
+mpi_subm( MPI w, MPI u, MPI v, MPI m)
+{
+    mpi_sub(w, u, v);
+    mpi_fdiv_r( w, w, m );
+}
+