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