5 Author: Pekka Riikonen <priikone@poseidon.pspt.fi>
7 Copyright (C) 2001 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; either version 2 of the License, or
12 (at your option) any later version.
14 This program is distributed in the hope that it will be useful,
15 but WITHOUT ANY WARRANTY; without even the implied warranty of
16 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17 GNU General Public License for more details.
22 #include "silcincludes.h"
25 void silc_mp_init(SilcMPInt *mp)
30 void silc_mp_uninit(SilcMPInt *mp)
35 size_t silc_mp_size(SilcMPInt *mp)
40 size_t silc_mp_sizeinbase(SilcMPInt *mp, int base)
42 return mpz_sizeinbase(mp, base);
45 void silc_mp_set(SilcMPInt *dst, SilcMPInt *src)
50 void silc_mp_set_ui(SilcMPInt *dst, uint32 ui)
55 void silc_mp_set_si(SilcMPInt *dst, int32 si)
60 void silc_mp_set_str(SilcMPInt *dst, const char *str, int base)
62 mpz_set_str(dst, str, base);
65 uint32 silc_mp_get_ui(SilcMPInt *mp)
67 return (uint32)mpz_get_ui(mp);
70 char *silc_mp_get_str(char *str, SilcMPInt *mp, int base)
72 return mpz_get_str(str, base, mp);
75 void silc_mp_add(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
77 mpz_add(dst, mp1, mp2);
80 void silc_mp_add_ui(SilcMPInt *dst, SilcMPInt *mp1, uint32 ui)
82 mpz_add_ui(dst, mp1, ui);
85 void silc_mp_sub(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
87 mpz_sub(dst, mp1, mp2);
90 void silc_mp_sub_ui(SilcMPInt *dst, SilcMPInt *mp1, uint32 ui)
92 mpz_sub_ui(dst, mp1, ui);
95 void silc_mp_mul(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
97 mpz_mul(dst, mp1, mp2);
100 void silc_mp_mul_ui(SilcMPInt *dst, SilcMPInt *mp1, uint32 ui)
102 mpz_mul_ui(dst, mp1, ui);
105 void silc_mp_mul_2exp(SilcMPInt *dst, SilcMPInt *mp1, uint32 exp)
107 mpz_mul_2exp(dst, mp1, exp);
110 void silc_mp_sqrt(SilcMPInt *dst, SilcMPInt *src)
115 void silc_mp_div(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
117 mpz_div(dst, mp1, mp2);
120 void silc_mp_div_ui(SilcMPInt *dst, SilcMPInt *mp1, uint32 ui)
122 mpz_div_ui(dst, mp1, ui);
125 void silc_mp_div_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
129 mpz_fdiv_qr(q, r, mp1, mp2);
131 mpz_div(q, mp1, mp2);
133 mpz_mod(r, mp1, mp2);
136 void silc_mp_div_2exp(SilcMPInt *dst, SilcMPInt *mp1, uint32 exp)
138 mpz_fdiv_q_2exp(dst, mp1, exp);
141 void silc_mp_div_2exp_qr(SilcMPInt *q, SilcMPInt *r, SilcMPInt *mp1,
145 mpz_fdiv_q_2exp(q, mp1, exp);
147 mpz_fdiv_r_2exp(r, mp1, exp);
150 void silc_mp_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
152 mpz_mod(dst, mp1, mp2);
155 void silc_mp_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, uint32 ui)
157 mpz_mod_ui(dst, mp1, ui);
160 void silc_mp_mod_2exp(SilcMPInt *dst, SilcMPInt *mp1, uint32 ui)
162 mpz_mod_2exp(dst, mp1, ui);
165 void silc_mp_pow(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp)
167 uint32 uiexp = mpz_get_ui(exp);
168 mpz_pow_ui(dst, mp1, uiexp);
171 void silc_mp_pow_ui(SilcMPInt *dst, SilcMPInt *mp1, uint32 exp)
173 mpz_pow_ui(dst, mp1, exp);
176 void silc_mp_pow_mod(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *exp,
179 mpz_powm(dst, mp1, exp, mod);
182 void silc_mp_pow_mod_ui(SilcMPInt *dst, SilcMPInt *mp1, uint32 exp,
185 mpz_powm_ui(dst, mp1, exp, mod);
188 void silc_mp_gcd(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
190 mpz_gcd(dst, mp1, mp2);
193 void silc_mp_gcdext(SilcMPInt *g, SilcMPInt *s, SilcMPInt *t, SilcMPInt *mp1,
196 mpz_gcdext(g, s, t, mp1, mp2);
199 int silc_mp_cmp(SilcMPInt *mp1, SilcMPInt *mp2)
201 return mpz_cmp(mp1, mp2);
204 int silc_mp_cmp_si(SilcMPInt *mp1, int32 si)
206 return mpz_cmp_si(mp1, si);
209 int silc_mp_cmp_ui(SilcMPInt *mp1, uint32 ui)
211 return mpz_cmp_ui(mp1, ui);
214 void silc_mp_abs(SilcMPInt *dst, SilcMPInt *src)
219 void silc_mp_neg(SilcMPInt *dst, SilcMPInt *src)
224 void silc_mp_and(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
226 mpz_and(dst, mp1, mp2);
229 void silc_mp_or(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
231 mpz_ior(dst, mp1, mp2);
234 void silc_mp_xor(SilcMPInt *dst, SilcMPInt *mp1, SilcMPInt *mp2)
236 mpz_xor(dst, mp1, mp2);