lib/silcmath/mp_mpi.c

   1 /*
   2
   3   mp_mpi.c
   4
   5   Author: Pekka Riikonen <priikone@silcnet.org>
   6
   7   Copyright (C) 2001 Pekka Riikonen
   8
   9   This program is free software; you can redistribute it and/or modify
  10   it under the terms of the GNU General Public License as published by
  11   the Free Software Foundation; version 2 of the License.
  12
  13   This program is distributed in the hope that it will be useful,
  14   but WITHOUT ANY WARRANTY; without even the implied warranty of
  15   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16   GNU General Public License for more details.
  17
  18 */
  19 /* $Id$ */
  20
  21 #include "silcincludes.h"
  22 #include "mpi.h"
  23 #include "mplogic.h"
  24
  25 void silc_mp_init(SilcMPInt *mp)
  26 {
  27   (void)mp_init(mp);
  28 }
  29
  30 void silc_mp_uninit(SilcMPInt *mp)
  31 {
  32   (void)mp_clear(mp);
  33 }
  34
  35 size_t silc_mp_size(SilcMPInt *mp)
  36 {
  37   return mp_raw_size(mp);
  38 }
  39
  40 size_t silc_mp_sizeinbase(SilcMPInt *mp, int base)
  41 {
  42   size_t sib = mp_unsigned_octet_size(mp);
  43   sib = s_mp_outlen(sib * CHAR_BIT, base);
  44   if (sib > 2)
  45     sib -= 2;                   /* Looks like MPI returns extra zero
  46                                    bytes for C-strings. */
  47   return sib;
  48 }
  49
  50 void silc_mp_set(SilcMPInt *dst, SilcMPInt *src)
  51 {
  52   (void)mp_copy(src, dst);
  53 }
  54
  55 void silc_mp_set_ui(SilcMPInt *dst, SilcUInt32 ui)
  56 {
  57   mp_set(dst, ui);
  58 }
  59
  60 void silc_mp_set_si(SilcMPInt *dst, SilcInt32 si)
  61 {
  62   (void)mp_set_int(dst, si);
  63 }
  64
  65 void silc_mp_set_str(SilcMPInt *dst, const char *str, int base)
  66 {
  67   (void)mp_read_variable_radix(dst, str, base);
  68 }
  69
  70 SilcUInt32 silc_mp_get_ui(SilcMPInt *mp)
  71 {
  72   return (SilcUInt32)MP_DIGIT(mp, 0);
  73 }
  74
  75 char *silc_mp_get_str(char *str, SilcMPInt *mp, int base)
  76 {
  77   if (mp_toradix(mp, str, base) != MP_OKAY)
  78     return NULL;
  79   return str;
  80 }
  81
  82 void silc_mp_add(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
  83 {
  84   (void)mp_add(mp1, mp2, dst);
  85 }
  86
  87 void silc_mp_add_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
  88 {
  89   mp_add_d(mp1, ui, dst);
  90 }
  91
  92 void silc_mp_sub(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
  93 {
  94   (void)mp_sub(mp1, mp2, dst);
  95 }
  96
  97 void silc_mp_sub_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
  98 {
  99   (void)mp_sub_d(mp1, (mp_digit)ui, dst);
 100 }
 101
 102 void silc_mp_mul(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 103 {
 104   (void)mp_mul(mp1, mp2, dst);
 105 }
 106
 107 void silc_mp_mul_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 108 {
 109   (void)mp_mul_d(mp1, (mp_digit)ui, dst);
 110 }
 111
 112 void silc_mp_mul_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
 113 {
 114   SilcMPInt tmp;
 115   silc_mp_init(&tmp);
 116   (void)mp_2expt(&tmp, (mp_digit)exp);
 117   (void)mp_mul(mp1, &tmp, dst);
 118   silc_mp_uninit(&tmp);
 119 }
 120
 121 void silc_mp_sqrt(SilcMPInt *dst, SilcMPInt *src)
 122 {
 123   (void)mp_sqrt(src, dst);
 124 }
 125
 126 void silc_mp_div(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 127 {
 128   (void)mp_div(mp1, mp2, dst, NULL);
 129 }
 130
 131 void silc_mp_div_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 132 {
 133   (void)mp_div_d(mp1, (mp_digit)ui, dst, NULL);
 134 }
 135
 136 void silc_mp_div_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
 137                     SilcMPInt *mp2)
 138 {
 139   (void)mp_div(mp1, mp2, q, r);
 140 }
 141
 142 void silc_mp_div_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
 143 {
 144   SilcMPInt tmp;
 145   silc_mp_init(&tmp);
 146   (void)mp_2expt(&tmp, (mp_digit)exp);
 147   (void)mp_div(mp1, &tmp, dst, NULL);
 148   silc_mp_uninit(&tmp);
 149 }
 150
 151 void silc_mp_div_2exp_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
 152                          SilcUInt32 exp)
 153 {
 154   if (q) {
 155     (void)mp_2expt(q, (mp_digit)exp);
 156     (void)mp_div(mp1, q, q, r);
 157   }
 158 }
 159
 160 void silc_mp_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 161 {
 162   (void)mp_mod(mp1, mp2, dst);
 163 }
 164
 165 void silc_mp_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 166 {
 167   mp_digit uidst;
 168   (void)mp_mod_d(mp1, (mp_digit)ui, &uidst);
 169   mp_set(dst, uidst);
 170 }
 171
 172 void silc_mp_mod_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 173 {
 174   SilcMPInt tmp;
 175   silc_mp_init(&tmp);
 176   (void)mp_2expt(&tmp, (mp_digit)ui);
 177   (void)mp_mod(mp1, &tmp, dst);
 178   silc_mp_uninit(&tmp);
 179 }
 180
 181 void silc_mp_pow(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp)
 182 {
 183   (void)mp_expt(mp1, exp, dst);
 184 }
 185
 186 void silc_mp_pow_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
 187 {
 188   (void)mp_expt_d(mp1, (mp_digit)exp, dst);
 189 }
 190
 191 void silc_mp_pow_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp,
 192                      SilcMPInt *mod)
 193 {
 194   (void)mp_exptmod(mp1, exp, mod, dst);
 195 }
 196
 197 void silc_mp_pow_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp,
 198                         SilcMPInt *mod)
 199 {
 200   (void)mp_exptmod_d(mp1, (mp_digit)exp, mod, dst);
 201 }
 202
 203 void silc_mp_gcd(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 204 {
 205   (void)mp_gcd(mp1, mp2, dst);
 206 }
 207
 208 void silc_mp_gcdext(SilcMPInt *g, SilcMPInt *s, SilcMPInt *t, SilcMPInt *mp1,
 209                     SilcMPInt *mp2)
 210 {
 211   (void)mp_xgcd(mp1, mp2, g, s, t);
 212 }
 213
 214 int silc_mp_cmp(SilcMPInt *mp1, SilcMPInt *mp2)
 215 {
 216   return mp_cmp(mp1, mp2);
 217 }
 218
 219 int silc_mp_cmp_si(SilcMPInt *mp1, SilcInt32 si)
 220 {
 221   return mp_cmp_int(mp1, (long)si);
 222 }
 223
 224 int silc_mp_cmp_ui(SilcMPInt *mp1, SilcUInt32 ui)
 225 {
 226   return mp_cmp_d(mp1, ui);
 227 }
 228
 229 void silc_mp_abs(SilcMPInt *dst, SilcMPInt *src)
 230 {
 231   mp_abs(src, dst);
 232 }
 233
 234 void silc_mp_neg(SilcMPInt *dst, SilcMPInt *src)
 235 {
 236   mp_neg(src, dst);
 237 }
 238
 239 void silc_mp_and(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 240 {
 241   mpl_and(mp1, mp2, dst);
 242 }
 243
 244 void silc_mp_or(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 245 {
 246   mpl_or(mp1, mp2, dst);
 247 }
 248
 249 void silc_mp_xor(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 250 {
 251   mpl_xor(mp1, mp2, dst);
 252 }