lib/silcmath/mp_tma.c

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