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