lib/silcmath/mp_tfm.c

   1 /*
   2
   3   mp_tfm.c
   4
   5   Author: Pekka Riikonen <priikone@silcnet.org>
   6
   7   Copyright (C) 2005 - 2006 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 "silc.h"
  22 #include "mp_tfm.h"
  23
  24 void silc_mp_init(SilcMPInt *mp)
  25 {
  26   fp_init(mp);
  27 }
  28
  29 void silc_mp_uninit(SilcMPInt *mp)
  30 {
  31   fp_zero(mp);
  32 }
  33
  34 size_t silc_mp_size(SilcMPInt *mp)
  35 {
  36   return fp_unsigned_bin_size(mp);
  37 }
  38
  39 size_t silc_mp_sizeinbase(SilcMPInt *mp, int base)
  40 {
  41   int size = 0;
  42   fp_radix_size(mp, base, &size);
  43   if (size > 1)
  44     size--;
  45   return size;
  46 }
  47
  48 void silc_mp_set(SilcMPInt *dst, SilcMPInt *src)
  49 {
  50   fp_copy(src, dst);
  51 }
  52
  53 void silc_mp_set_ui(SilcMPInt *dst, SilcUInt32 ui)
  54 {
  55   fp_set(dst, ui);
  56 }
  57
  58 void silc_mp_set_si(SilcMPInt *dst, SilcInt32 si)
  59 {
  60   fp_set(dst, si);
  61 }
  62
  63 void silc_mp_set_str(SilcMPInt *dst, const char *str, int base)
  64 {
  65   fp_read_radix(dst, str, base);
  66 }
  67
  68 SilcUInt32 silc_mp_get_ui(SilcMPInt *mp)
  69 {
  70   fp_int *tmp = mp;
  71   return tmp->used > 0 ? tmp->dp[0] : 0;
  72 }
  73
  74 char *silc_mp_get_str(char *str, SilcMPInt *mp, int base)
  75 {
  76   if (fp_toradix(mp, str, base) != MP_OKAY)
  77     return NULL;
  78   return str;
  79 }
  80
  81 void silc_mp_add(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
  82 {
  83   fp_add(mp1, mp2, dst);
  84 }
  85
  86 void silc_mp_add_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
  87 {
  88   mp_add_d(mp1, (mp_digit)ui, dst);
  89 }
  90
  91 void silc_mp_sub(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
  92 {
  93   fp_sub(mp1, mp2, dst);
  94 }
  95
  96 void silc_mp_sub_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
  97 {
  98   fp_sub_d(mp1, (mp_digit)ui, dst);
  99 }
 100
 101 void silc_mp_mul(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 102 {
 103   fp_mul(mp1, mp2, dst);
 104 }
 105
 106 void silc_mp_mul_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 107 {
 108   fp_mul_d(mp1, (mp_digit)ui, dst);
 109 }
 110
 111 void silc_mp_mul_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
 112 {
 113   fp_mul_2d(mp1, exp, dst);
 114 }
 115
 116 void silc_mp_sqrt(SilcMPInt *dst, SilcMPInt *src)
 117 {
 118   fp_sqrt(src, dst);
 119 }
 120
 121 void silc_mp_div(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 122 {
 123   fp_div(mp1, mp2, dst, NULL);
 124 }
 125
 126 void silc_mp_div_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 127 {
 128   fp_div_d(mp1, (mp_digit)ui, dst, NULL);
 129 }
 130
 131 void silc_mp_div_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
 132                     SilcMPInt *mp2)
 133 {
 134   fp_div(mp1, mp2, q, r);
 135 }
 136
 137 void silc_mp_div_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
 138 {
 139   fp_div_2d(mp1, exp, dst, NULL);
 140 }
 141
 142 void silc_mp_div_2exp_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
 143                          SilcUInt32 exp)
 144 {
 145   fp_div_2d(mp1, exp, q, r);
 146 }
 147
 148 void silc_mp_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 149 {
 150   fp_mod(mp1, mp2, dst);
 151 }
 152
 153 void silc_mp_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 154 {
 155   fp_digit d;
 156   fp_mod_d(mp1, ui, &d);
 157   silc_mp_set_ui(dst, d);
 158 }
 159
 160 void silc_mp_mod_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 161 {
 162   fp_mod_2d(mp1, ui, dst);
 163 }
 164
 165 void silc_mp_pow(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp)
 166 {
 167   SILC_NOT_IMPLEMENTED("silc_mp_pow");
 168   assert(FALSE);
 169 }
 170
 171 void silc_mp_pow_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
 172 {
 173   SILC_NOT_IMPLEMENTED("silc_mp_pow_ui");
 174   assert(FALSE);
 175 }
 176
 177 void silc_mp_pow_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp,
 178                      SilcMPInt *mod)
 179 {
 180   fp_exptmod(mp1, exp, mod, dst);
 181 }
 182
 183 void silc_mp_pow_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp,
 184                         SilcMPInt *mod)
 185 {
 186   SilcMPInt tmp;
 187   silc_mp_init(&tmp);
 188   silc_mp_set_ui(&tmp, exp);
 189   silc_mp_pow_mod(dst, mp1, &tmp, mod);
 190   silc_mp_uninit(&tmp);
 191 }
 192
 193 void silc_mp_gcd(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 194 {
 195   fp_gcd(mp1, mp2, dst);
 196 }
 197
 198 void silc_mp_gcdext(SilcMPInt *g, SilcMPInt *s, SilcMPInt *t, SilcMPInt *mp1,
 199                     SilcMPInt *mp2)
 200 {
 201   SILC_NOT_IMPLEMENTED("silc_mp_gcdext");
 202   assert(FALSE);
 203 }
 204
 205 int silc_mp_cmp(SilcMPInt *mp1, SilcMPInt *mp2)
 206 {
 207   return fp_cmp(mp1, mp2);
 208 }
 209
 210 int silc_mp_cmp_si(SilcMPInt *mp1, SilcInt32 si)
 211 {
 212   return fp_cmp_d(mp1, si);
 213 }
 214
 215 int silc_mp_cmp_ui(SilcMPInt *mp1, SilcUInt32 ui)
 216 {
 217   return fp_cmp_d(mp1, ui);
 218 }
 219
 220 void silc_mp_abs(SilcMPInt *dst, SilcMPInt *src)
 221 {
 222   fp_abs(src, dst);
 223 }
 224
 225 void silc_mp_neg(SilcMPInt *dst, SilcMPInt *src)
 226 {
 227   fp_neg(src, dst);
 228 }
 229
 230 void silc_mp_and(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 231 {
 232   SILC_NOT_IMPLEMENTED("silc_mp_and");
 233   assert(FALSE);
 234 }
 235
 236 void silc_mp_or(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 237 {
 238   SILC_NOT_IMPLEMENTED("silc_mp_or");
 239   assert(FALSE);
 240 }
 241
 242 void silc_mp_xor(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 243 {
 244   SILC_NOT_IMPLEMENTED("silc_mp_xor");
 245   assert(FALSE);
 246 }