lib/silcmath/mp_gmp.c

   1 /*
   2
   3   mp_gmp.c
   4
   5   Author: Pekka Riikonen <priikone@silcnet.org>
   6
   7   Copyright (C) 2001 - 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 <gmp.h>
  23
  24 void silc_mp_init(SilcMPInt *mp)
  25 {
  26   mpz_init(mp);
  27 }
  28
  29 void silc_mp_uninit(SilcMPInt *mp)
  30 {
  31   mpz_clear(mp);
  32 }
  33
  34 size_t silc_mp_size(SilcMPInt *mp)
  35 {
  36   return mpz_size(mp);
  37 }
  38
  39 size_t silc_mp_sizeinbase(SilcMPInt *mp, int base)
  40 {
  41   return mpz_sizeinbase(mp, base);
  42 }
  43
  44 void silc_mp_set(SilcMPInt *dst, SilcMPInt *src)
  45 {
  46   mpz_set(dst, src);
  47 }
  48
  49 void silc_mp_set_ui(SilcMPInt *dst, SilcUInt32 ui)
  50 {
  51   mpz_set_ui(dst, ui);
  52 }
  53
  54 void silc_mp_set_si(SilcMPInt *dst, SilcInt32 si)
  55 {
  56   mpz_set_si(dst, si);
  57 }
  58
  59 void silc_mp_set_str(SilcMPInt *dst, const char *str, int base)
  60 {
  61   mpz_set_str(dst, str, base);
  62 }
  63
  64 SilcUInt32 silc_mp_get_ui(SilcMPInt *mp)
  65 {
  66   return (SilcUInt32)mpz_get_ui(mp);
  67 }
  68
  69 char *silc_mp_get_str(char *str, SilcMPInt *mp, int base)
  70 {
  71   return mpz_get_str(str, base, mp);
  72 }
  73
  74 void silc_mp_add(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
  75 {
  76   mpz_add(dst, mp1, mp2);
  77 }
  78
  79 void silc_mp_add_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
  80 {
  81   mpz_add_ui(dst, mp1, ui);
  82 }
  83
  84 void silc_mp_sub(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
  85 {
  86   mpz_sub(dst, mp1, mp2);
  87 }
  88
  89 void silc_mp_sub_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
  90 {
  91   mpz_sub_ui(dst, mp1, ui);
  92 }
  93
  94 void silc_mp_mul(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
  95 {
  96   mpz_mul(dst, mp1, mp2);
  97 }
  98
  99 void silc_mp_mul_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 100 {
 101   mpz_mul_ui(dst, mp1, ui);
 102 }
 103
 104 void silc_mp_mul_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
 105 {
 106   mpz_mul_2exp(dst, mp1, exp);
 107 }
 108
 109 void silc_mp_sqrt(SilcMPInt *dst, SilcMPInt *src)
 110 {
 111   mpz_sqrt(dst, src);
 112 }
 113
 114 void silc_mp_div(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 115 {
 116   mpz_div(dst, mp1, mp2);
 117 }
 118
 119 void silc_mp_div_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 120 {
 121   mpz_div_ui(dst, mp1, ui);
 122 }
 123
 124 void silc_mp_div_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
 125                     SilcMPInt *mp2)
 126 {
 127   if (q && r)
 128     mpz_fdiv_qr(q, r, mp1, mp2);
 129   if (q && !r)
 130     mpz_div(q, mp1, mp2);
 131   if (!q && r)
 132     mpz_mod(r, mp1, mp2);
 133 }
 134
 135 void silc_mp_div_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
 136 {
 137   mpz_fdiv_q_2exp(dst, mp1, exp);
 138 }
 139
 140 void silc_mp_div_2exp_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
 141                          SilcUInt32 exp)
 142 {
 143   if (q)
 144     mpz_fdiv_q_2exp(q, mp1, exp);
 145   if (r)
 146     mpz_fdiv_r_2exp(r, mp1, exp);
 147 }
 148
 149 void silc_mp_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 150 {
 151   mpz_mod(dst, mp1, mp2);
 152 }
 153
 154 void silc_mp_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 155 {
 156   mpz_mod_ui(dst, mp1, ui);
 157 }
 158
 159 void silc_mp_mod_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
 160 {
 161   mpz_mod_2exp(dst, mp1, ui);
 162 }
 163
 164 void silc_mp_pow(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp)
 165 {
 166   SilcUInt32 uiexp = mpz_get_ui(exp);
 167   mpz_pow_ui(dst, mp1, uiexp);
 168 }
 169
 170 void silc_mp_pow_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
 171 {
 172   mpz_pow_ui(dst, mp1, exp);
 173 }
 174
 175 void silc_mp_pow_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp,
 176                      SilcMPInt *mod)
 177 {
 178   mpz_powm(dst, mp1, exp, mod);
 179 }
 180
 181 void silc_mp_pow_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp,
 182                         SilcMPInt *mod)
 183 {
 184   mpz_powm_ui(dst, mp1, exp, mod);
 185 }
 186
 187 void silc_mp_gcd(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 188 {
 189   mpz_gcd(dst, mp1, mp2);
 190 }
 191
 192 void silc_mp_gcdext(SilcMPInt *g, SilcMPInt *s, SilcMPInt *t, SilcMPInt *mp1,
 193                     SilcMPInt *mp2)
 194 {
 195   mpz_gcdext(g, s, t, mp1, mp2);
 196 }
 197
 198 int silc_mp_cmp(SilcMPInt *mp1, SilcMPInt *mp2)
 199 {
 200   return mpz_cmp(mp1, mp2);
 201 }
 202
 203 int silc_mp_cmp_si(SilcMPInt *mp1, SilcInt32 si)
 204 {
 205   return mpz_cmp_si(mp1, si);
 206 }
 207
 208 int silc_mp_cmp_ui(SilcMPInt *mp1, SilcUInt32 ui)
 209 {
 210   return mpz_cmp_ui(mp1, ui);
 211 }
 212
 213 void silc_mp_abs(SilcMPInt *dst, SilcMPInt *src)
 214 {
 215   mpz_abs(dst, src);
 216 }
 217
 218 void silc_mp_neg(SilcMPInt *dst, SilcMPInt *src)
 219 {
 220   mpz_neg(dst, src);
 221 }
 222
 223 void silc_mp_and(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 224 {
 225   mpz_and(dst, mp1, mp2);
 226 }
 227
 228 void silc_mp_or(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 229 {
 230   mpz_ior(dst, mp1, mp2);
 231 }
 232
 233 void silc_mp_xor(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
 234 {
 235   mpz_xor(dst, mp1, mp2);
 236 }