5 Author: Pekka Riikonen <priikone@silcnet.org>
7 Copyright (C) 2005 Pekka Riikonen
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.
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.
21 #include "silcincludes.h"
24 void silc_mp_init(SilcMPInt *mp)
29 void silc_mp_uninit(SilcMPInt *mp)
34 size_t silc_mp_size(SilcMPInt *mp)
36 return mp_unsigned_bin_size(mp);
39 size_t silc_mp_sizeinbase(SilcMPInt *mp, int base)
42 mp_radix_size(mp, base, &size);
48 void silc_mp_set(SilcMPInt *dst, SilcMPInt *src)
50 (void)mp_copy(src, dst);
53 void silc_mp_set_ui(SilcMPInt *dst, SilcUInt32 ui)
55 (void)mp_set_int(dst, ui);
58 void silc_mp_set_si(SilcMPInt *dst, SilcInt32 si)
60 (void)mp_set_int(dst, si);
63 void silc_mp_set_str(SilcMPInt *dst, const char *str, int base)
65 (void)mp_read_radix(dst, str, base);
68 SilcUInt32 silc_mp_get_ui(SilcMPInt *mp)
70 return (SilcUInt32)mp_get_int(mp);
73 char *silc_mp_get_str(char *str, SilcMPInt *mp, int base)
75 if (mp_toradix(mp, str, base) != MP_OKAY)
80 void silc_mp_add(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
82 (void)mp_add(mp1, mp2, dst);
85 void silc_mp_add_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
87 mp_add_d(mp1, (mp_digit)ui, dst);
90 void silc_mp_sub(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
92 (void)mp_sub(mp1, mp2, dst);
95 void silc_mp_sub_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
97 (void)mp_sub_d(mp1, (mp_digit)ui, dst);
100 void silc_mp_mul(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
102 (void)mp_mul(mp1, mp2, dst);
105 void silc_mp_mul_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
107 (void)mp_mul_d(mp1, (mp_digit)ui, dst);
110 void silc_mp_mul_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
112 (void)mp_mul_2d(mp1, exp, dst);
115 void silc_mp_sqrt(SilcMPInt *dst, SilcMPInt *src)
117 (void)mp_sqrt(src, dst);
120 void silc_mp_div(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
122 (void)mp_div(mp1, mp2, dst, NULL);
125 void silc_mp_div_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
127 (void)mp_div_d(mp1, (mp_digit)ui, dst, NULL);
130 void silc_mp_div_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
133 (void)mp_div(mp1, mp2, q, r);
136 void silc_mp_div_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
138 (void)mp_div_2d(mp1, exp, dst, NULL);
141 void silc_mp_div_2exp_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
144 (void)mp_div_2d(mp1, exp, q, r);
147 void silc_mp_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
149 (void)mp_mod(mp1, mp2, dst);
152 void silc_mp_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
155 (void)mp_mod_d(mp1, ui, &d);
156 silc_mp_set_ui(dst, d);
159 void silc_mp_mod_2exp(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 ui)
161 (void)mp_mod_2d(mp1, ui, dst);
164 void silc_mp_pow(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp)
166 SILC_NOT_IMPLEMENTED("silc_mp_pow");
170 void silc_mp_pow_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp)
172 (void)mp_expt_d(mp1, (mp_digit)exp, dst);
175 void silc_mp_pow_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp,
178 (void)mp_exptmod(mp1, exp, mod, dst);
181 void silc_mp_pow_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, SilcUInt32 exp,
186 silc_mp_set_ui(&tmp, exp);
187 silc_mp_pow_mod(dst, mp1, &tmp, mod);
188 silc_mp_uninit(&tmp);
191 void silc_mp_gcd(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
193 (void)mp_gcd(mp1, mp2, dst);
196 void silc_mp_gcdext(SilcMPInt *g, SilcMPInt *s, SilcMPInt *t, SilcMPInt *mp1,
199 (void)mp_exteuclid(mp1, mp2, s, t, g);
202 int silc_mp_cmp(SilcMPInt *mp1, SilcMPInt *mp2)
204 return mp_cmp(mp1, mp2);
207 int silc_mp_cmp_si(SilcMPInt *mp1, SilcInt32 si)
209 return mp_cmp_d(mp1, si);
212 int silc_mp_cmp_ui(SilcMPInt *mp1, SilcUInt32 ui)
214 return mp_cmp_d(mp1, ui);
217 void silc_mp_abs(SilcMPInt *dst, SilcMPInt *src)
219 (void)mp_abs(src, dst);
222 void silc_mp_neg(SilcMPInt *dst, SilcMPInt *src)
224 (void)mp_neg(src, dst);
227 void silc_mp_and(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
229 (void)mp_and(mp1, mp2, dst);
232 void silc_mp_or(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
234 (void)mp_or(mp1, mp2, dst);
237 void silc_mp_xor(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
239 (void)mp_xor(mp1, mp2, dst);