libtommath for Win-64

[jack]

made a win64 dll of LibTomMath https://www.libtom.net/

Code: QB64: [Select]

1. $Console:Only
2. _Dest _Console
3.  
4. Const MP_64BIT = 1
5. Type private_mp_word
6.     As String * 128 mpword
7. End Type
8. Const MP_DIGIT_BIT = 60
9.  
10. Type mp_int
11.     used As Long
12.     alloc As Long
13.     sign As Long
14.     dp As _Offset
15. End Type
16.  
17. Declare Dynamic Library "libtommath"
18.     Function mp_init& (a As mp_int) ' as mp_err
19.     Function mp_init_size& (a As mp_int, Byval size As Long) ' as mp_err
20.     Function mp_init_i32& (a As mp_int, Byval b As Long) ' as mp_err
21.     Function mp_init_l& (a As mp_int, Byval b As Long) ' ' as mp_err
22.     Function mp_init_u32& (a As mp_int, Byval b As _Unsigned Long) ' as mp_err
23.     Function mp_init_ul& (a As mp_int, Byval b As _Unsigned Long) ' as mp_err
24.     Function mp_init_i64& (a As mp_int, Byval b As _Integer64) ' as mp_err
25.     Function mp_init_ll& (a As mp_int, Byval b As _Integer64) ' as mp_err
26.     Function mp_init_u64& (a As mp_int, Byval b As _Unsigned _Integer64) ' as mp_err
27.     Function mp_init_ull& (a As mp_int, Byval b As _Unsigned _Integer64) ' as mp_err
28.     Function mp_init_set& (a As mp_int, Byval b As _Unsigned _Integer64) ' as mp_err
29.     Function mp_init_set_int& (a As mp_int, Byval b As _Unsigned Long) ' as mp_err
30.     Function mp_init_copy (a As mp_int, b As mp_int) ' as mp_err
31.     Sub mp_clear (a As mp_int)
32.     Sub mp_exch (a As mp_int, b As mp_int)
33.     Function mp_shrink& (a As mp_int) ' as mp_err
34.     Function mp_grow& (a As mp_int, Byval size As Long) ' as mp_err
35.     Function mp_iseven& (a As mp_int) ' as long
36.     Function mp_isodd& (a As mp_int) ' as long
37.     Sub mp_zero (a As mp_int)
38.     Function mp_get_double# (a As mp_int) ' as double
39.     Function mp_set_double& (a As mp_int, Byval b As Double) ' as mp_err
40.     Function mp_get_i32& (a As mp_int) ' as long
41.     Function mp_get_l& (a As mp_int) ' as long
42.     Function mp_get_int~& (a As mp_int) ' as ulong
43.     Function mp_get_long~& (a As mp_int) ' as ulong
44.     Function mp_get_i64&& (a As mp_int) ' as longint
45.     Function mp_get_ll&& (a As mp_int) ' as longint
46.     Function mp_get_long_long~&& (a As mp_int) ' as ulongint
47.     Sub mp_set_i32 (a As mp_int, Byval b As Long)
48.     Sub mp_set_l (a As mp_int, Byval b As Long)
49.     Function mp_set_long& (a As mp_int, Byval b As _Unsigned Long) ' as mp_err
50.     Sub mp_set_u32 (a As mp_int, Byval b As _Unsigned Long)
51.     Sub mp_set_ul (a As mp_int, Byval b As _Unsigned Long)
52.     Function mp_set_int& (a As mp_int, Byval b As _Unsigned Long) ' as mp_err
53.     Sub mp_set_i64 (a As mp_int, Byval b As _Integer64)
54.     Function mp_set_long_long& (a As mp_int, Byval b As _Unsigned _Integer64) ' as mp_err
55.     Sub mp_set_ll (a As mp_int, Byval b As _Integer64)
56.     Sub mp_set_u64 (a As mp_int, Byval b As _Unsigned _Integer64)
57.     Sub mp_set_ull (a As mp_int, Byval b As _Unsigned _Integer64)
58.     Sub mp_set (a As mp_int, Byval b As _Unsigned _Integer64)
59.     Function mp_get_mag_u32~& (a As mp_int) ' as ulong
60.     Function mp_get_mag_ul~& (a As mp_int) ' as ulong
61.     Function mp_get_mag_u64~&& (a As mp_int) ' as ulongint
62.     Function mp_get_mag_ull~&& (a As mp_int) ' as ulongint
63.     Function mp_copy (a As mp_int, b As mp_int) ' as mp_err
64.     Sub mp_clamp (a As mp_int)
65.     Function mp_export& (rop As _Offset, countp As _Unsigned Long, Byval order As Long, Byval size As _Unsigned Long, Byval endian As Long, Byval nails As _Unsigned Long, op As mp_int) ' as mp_err
66.     Function mp_import& (rop As mp_int, Byval count As _Unsigned Long, Byval order As Long, Byval size As _Unsigned Long, Byval endian As Long, Byval nails As _Unsigned Long, Byval op As _Offset) ' as mp_err
67.     Function mp_unpack& (rop As mp_int, Byval count As _Unsigned Long, Byval order As Long, Byval size As _Unsigned Long, Byval endian As Long, Byval nails As _Unsigned Long, op As _Offset) ' as mp_err
68.     Function mp_pack_count~& (a As mp_int, Byval nails As _Unsigned Long, Byval size As _Unsigned Long) ' as uinteger
69.     Function mp_pack& (rop As _Offset, Byval maxcount As _Unsigned Long, written As _Unsigned Long, Byval order As Long, Byval size As _Unsigned Long, Byval endian As Long, Byval nails As _Unsigned Long, op As mp_int) ' as mp_err
70.     Sub mp_rshd (a As mp_int, Byval b As Long)
71.     Function mp_lshd& (a As mp_int, Byval b As Long) ' as mp_err
72.     Function mp_div_2d& (a As mp_int, Byval b As Long, c As mp_int, d As mp_int) ' as mp_err
73.     Function mp_div_2& (a As mp_int, b As mp_int) ' as mp_err
74.     Function mp_div_3& (a As mp_int, c As mp_int, d As _Unsigned _Integer64) ' as mp_err
75.     Function mp_mul_2d& (a As mp_int, Byval b As Long, c As mp_int) ' as mp_err
76.     Function mp_mul_2& (a As mp_int, b As mp_int) ' as mp_err
77.     Function mp_mod_2d& (a As mp_int, Byval b As Long, c As mp_int) ' as mp_err
78.     Function mp_2expt& (a As mp_int, Byval b As Long) ' as mp_err
79.     Function mp_cnt_lsb& (a As mp_int) ' as long
80.     Function mp_rand& (a As mp_int, Byval digits As Long) ' as mp_err
81.     Function mp_rand_digit& (r As _Unsigned _Integer64) ' as mp_err
82.     Function mp_get_bit& (a As mp_int, Byval b As Long) ' as long
83.     Function mp_tc_xor& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
84.     Function mp_xor& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
85.     Function mp_tc_or& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
86.     Function mp_or& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
87.     Function mp_tc_and& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
88.     Function mp_and& (a As mp_int, b As mp_int, c As mp_int) 'as mp_err
89.     Function mp_complement& (a As mp_int, b As mp_int) ' as mp_err
90.     Function mp_tc_div_2d& (a As mp_int, Byval b As Long, c As mp_int) ' as mp_err
91.     Function mp_signed_rsh& (a As mp_int, Byval b As Long, c As mp_int) ' as mp_err
92.     Function mp_neg (a As mp_int, b As mp_int) ' as mp_err
93.     Function mp_abs& (a As mp_int, b As mp_int) ' as mp_err
94.     Function mp_cmp& (a As mp_int, b As mp_int) 'as mp_ord
95.     Function mp_cmp_mag& (a As mp_int, b As mp_int) ' as mp_ord
96.     Function mp_add& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
97.     Function mp_sub& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
98.     Function mp_mul& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
99.     Function mp_sqr& (a As mp_int, b As mp_int) ' as mp_err
100.     Function mp_div& (a As mp_int, b As mp_int, c As mp_int, d As mp_int) ' as mp_err
101.     Function mp_mod& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
102.     Function mp_incr& (a As mp_int) ' as mp_err
103.     Function mp_decr& (a As mp_int) ' as mp_err
104.     Function mp_cmp_d& (a As mp_int, Byval b As _Unsigned _Integer64) ' as mp_ord
105.     Function mp_add_d& (a As mp_int, Byval b As _Unsigned _Integer64, c As mp_int) ' as mp_err
106.     Function mp_sub_d& (a As mp_int, Byval b As _Unsigned _Integer64, c As mp_int) ' as mp_err
107.     Function mp_mul_d& (a As mp_int, Byval b As _Unsigned _Integer64, c As mp_int) ' as mp_err
108.     Function mp_div_d& (a As mp_int, Byval b As _Unsigned _Integer64, c As mp_int, d As _Unsigned _Integer64) ' as mp_err
109.     Function mp_mod_d& (a As mp_int, Byval b As _Unsigned _Integer64, c As _Unsigned _Integer64) ' as mp_err
110.     Function mp_addmod& (a As mp_int, b As mp_int, c As mp_int, d As mp_int) ' as mp_err
111.     Function mp_submod& (a As mp_int, b As mp_int, c As mp_int, d As mp_int) ' as mp_err
112.     Function mp_mulmod& (a As mp_int, b As mp_int, c As mp_int, d As mp_int) ' as mp_err
113.     Function mp_sqrmod& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
114.     Function mp_invmod& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
115.     Function mp_gcd& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
116.     Function mp_exteuclid& (a As mp_int, b As mp_int, U1 As mp_int, U2 As mp_int, U3 As mp_int) ' as mp_err
117.     Function mp_lcm& (a As mp_int, b As mp_int, c As mp_int) ' as mp_err
118.     Function mp_root_u32& (a As mp_int, Byval b As _Unsigned Long, c As mp_int) ' as mp_err
119.     Function mp_n_root& (a As mp_int, Byval b As _Unsigned _Integer64, c As mp_int) ' as mp_err
120.     Function mp_n_root_ex& (a As mp_int, Byval b As _Unsigned _Integer64, c As mp_int, Byval fast As Long) ' as mp_err
121.     Function mp_sqrt& (arg As mp_int, ret As mp_int) ' as mp_err
122.     Function mp_sqrtmod_prime& (n As mp_int, prime As mp_int, ret As mp_int) ' as mp_err
123.     Function mp_is_square& (arg As mp_int, ret As Long) ' as mp_err
124.     Function mp_jacobi& (a As mp_int, n As mp_int, c As Long) ' as mp_err
125.     Function mp_kronecker& (a As mp_int, p As mp_int, c As Long) ' as mp_err
126.     Function mp_reduce_setup& (a As mp_int, b As mp_int) ' as mp_err
127.     Function mp_reduce& (x As mp_int, m As mp_int, mu As mp_int) ' as mp_err
128.     Function mp_montgomery_setup& (n As mp_int, rho As _Unsigned _Integer64) ' as mp_err
129.     Function mp_montgomery_calc_normalization& (a As mp_int, b As mp_int) ' as mp_err
130.     Function mp_montgomery_reduce& (x As mp_int, n As mp_int, Byval rho As _Unsigned _Integer64) ' as mp_err
131.     Function mp_dr_is_modulus& (a As mp_int) ' as mp_bool
132.     Sub mp_dr_setup (a As mp_int, d As _Unsigned _Integer64)
133.     Function mp_dr_reduce& (x As mp_int, n As mp_int, Byval k As _Unsigned _Integer64) ' as mp_err
134.     Function mp_reduce_is_2k& (a As mp_int) ' as mp_bool
135.     Function mp_reduce_2k_setup& (a As mp_int, d As _Unsigned _Integer64) ' as mp_err
136.     Function mp_reduce_2k& (a As mp_int, n As mp_int, Byval d As _Unsigned _Integer64) ' as mp_err
137.     Function mp_reduce_is_2k_l& (a As mp_int) ' as mp_bool
138.     Function mp_reduce_2k_setup_l& (a As mp_int, d As mp_int) ' as mp_err
139.     Function mp_reduce_2k_l& (a As mp_int, n As mp_int, d As mp_int) ' as mp_err
140.     Function mp_exptmod& (G As mp_int, X As mp_int, P As mp_int, Y As mp_int) ' as mp_err
141.     Function mp_prime_is_divisible& (a As mp_int, result As Long) ' as mp_err
142.     Function mp_prime_fermat& (a As mp_int, b As mp_int, result As Long) ' as mp_err
143.     Function mp_prime_miller_rabin& (a As mp_int, b As mp_int, result As Long) ' as mp_err
144.     Function mp_prime_rabin_miller_trials& (ByVal size As Long) ' as long
145.     Function mp_prime_strong_lucas_selfridge& (a As mp_int, result As Long) ' as mp_err
146.     Function mp_prime_frobenius_underwood& (N As mp_int, result As Long) ' as mp_err
147.     Function mp_prime_is_prime& (a As mp_int, Byval t As Long, result As Long) ' as mp_err
148.     Function mp_prime_next_prime& (a As mp_int, Byval t As Long, Byval bbs_style As Long) ' as mp_err
149.     Function mp_prime_rand& (a As mp_int, Byval t As Long, Byval size As Long, Byval flags As Long) ' as mp_err
150.     Function mp_log_u32& (a As mp_int, Byval base As _Unsigned Long, c As _Unsigned Long) ' as mp_err
151.     Function mp_expt_u32& (a As mp_int, Byval b As _Unsigned Long, c As mp_int) ' as mp_err
152.     Function mp_expt_d& (a As mp_int, Byval b As _Unsigned _Integer64, c As mp_int) ' as mp_err
153.     Function mp_expt_d_ex& (a As mp_int, Byval b As _Unsigned _Integer64, c As mp_int, Byval fast As Long) ' as mp_err
154.     Function mp_count_bits& (a As mp_int) ' as long
155.     Function mp_unsigned_bin_size& (a As mp_int) ' as long
156.     Function mp_read_unsigned_bin& (a As mp_int, Byval b As _Offset, Byval c As Long) ' as mp_err
157.     Function mp_to_unsigned_bin& (a As mp_int, Byval b As _Offset) ' as mp_err
158.     Function mp_to_unsigned_bin_n& (a As mp_int, Byval b As _Offset, outlen As _Unsigned Long) ' as mp_err
159.     Function mp_signed_bin_size& (a As mp_int) ' as long
160.     Function mp_read_signed_bin& (a As mp_int, Byval b As _Offset, Byval c As Long) ' as mp_err
161.     Function mp_to_signed_bin& (a As mp_int, Byval b As _Offset) ' as mp_err
162.     Function mp_to_signed_bin_n& (a As mp_int, Byval b As _Offset, outlen As _Unsigned Long) ' as mp_err
163.     Function mp_ubin_size~&& (a As mp_int) ' as uinteger
164.     Function mp_from_ubin& (a As mp_int, Byval buf As _Offset, Byval size As _Unsigned _Integer64) ' as mp_err
165.     Function mp_to_ubin& (a As mp_int, Byval buf As _Offset, Byval maxlen As _Unsigned _Integer64, written As _Unsigned _Integer64) ' as mp_err
166.     Function mp_sbin_size~&& (a As mp_int) ' as uinteger
167.     Function mp_from_sbin& (a As mp_int, Byval buf As _Offset, Byval size As _Unsigned _Integer64) ' as mp_err
168.     Function mp_to_sbin& (a As mp_int, Byval buf As _Offset, Byval maxlen As _Unsigned _Integer64, written As _Unsigned _Integer64) ' as mp_err
169.     Function mp_to_radix& (a As mp_int, str As String, Byval maxlen As _Unsigned _Integer64, written As _Unsigned _Integer64, Byval radix As Long) ' as mp_err
170.     Function mp_radix_size& (a As mp_int, Byval radix As Long, size As Long) ' as mp_err
171.     Function mp_read_radix& (a As mp_int, str As String, Byval radix As Long) ' as mp_err
172.     Function mp_toradix (a As mp_int, str As String, Byval radix As Long) ' as mp_err
173.     Function mp_toradix_n& (a As mp_int, str As String, Byval radix As Long, Byval maxlen As Long) ' as mp_err
174. End Declare
175.  
176. Dim As mp_int n, m, q, r
177. Dim As Long ok, k
178. Dim As Double t
179. t = Timer(.0001)
180. If mp_init(n) <> 0 Then Print "failed to initialize"
181. If mp_init(m) <> 0 Then Print "failed to initiali"
182. If mp_init(q) <> 0 Then Print "failed to initiali"
183. If mp_init(r) <> 0 Then Print "failed to initiali"
184. 'mp_val "2" + String$(100, "0"), n, 10
185. 'ok = mp_n_root&(n, 2, r)
186. 'Print mp_str(r, 10)
187. 'ok = mp_sqrt&(n, m)
188. 'Print mp_str(m, 10)
189. mp_set_i32 n, 2 '4784969
190. ok = mp_expt_u32(n, 1024, n)
191. 'mp_set_i32 m, 10
192. 'ok = mp_expt_u32(m, 1000000, m)
193. 'fib_mod n, m, r
194. 'mp_set_i32 m, 10
195. 'ok = mp_expt_u32(m, 999960, m)
196. 'ok = mp_div(r, m, q, r)
197. 'Print mp_str(q, 10)
198. mp_set_i32 m, 10
199. ok = mp_expt_u32(m, 40, m)
200. fib_mod n, m, r
201. Print mp_str(r, 10)
202. t = Timer(.0001) - t
203. Print t
204. mp_clear r
205. mp_clear q
206. mp_clear m
207. mp_clear n
208.  
209. Sub fib_mod (k As mp_int, m As mp_int, result As mp_int)
210.     Dim As mp_int x, y, xx, temp, tmp
211.     Dim As Long ok, cmp, i, bit_length
212.  
213.     If mp_init(x) <> 0 Then
214.         Print "failed to initialize"
215.         Exit Sub
216.     End If
217.     If mp_init(y) <> 0 Then
218.         Print "failed to initialize"
219.         Exit Sub
220.     End If
221.     If mp_init(xx) <> 0 Then
222.         Print "failed to initialize"
223.         Exit Sub
224.     End If
225.     If mp_init(temp) <> 0 Then
226.         Print "failed to initialize"
227.         Exit Sub
228.     End If
229.     If mp_init(tmp) <> 0 Then
230.         Print "failed to initialize"
231.         Exit Sub
232.     End If
233.  
234.     mp_set_i32 tmp, 1
235.     If mp_cmp(k, tmp) <= 0 Then
236.         Print "mp_cmp(k, tmp) <= 0"
237.         ok = mp_copy(k, result)
238.         Exit Sub
239.     End If
240.  
241.     mp_set_i32 x, 1
242.     mp_set_i32 y, 0
243.     ok = mp_log_u32(k, 2, bit_length)
244.     Print "bit_length="; bit_length
245.     For i = bit_length - 1 To 0 Step -1
246.         ok = mp_sqr(x, tmp)
247.         ok = mp_mod(tmp, m, xx)
248.         ok = mp_mul(x, y, tmp)
249.         ok = mp_add(tmp, tmp, tmp)
250.         ok = mp_add(xx, tmp, tmp)
251.         ok = mp_mod(tmp, m, x)
252.         ok = mp_sqr(y, tmp)
253.         ok = mp_add(xx, tmp, tmp)
254.         ok = mp_mod(tmp, m, y)
255.         If mp_get_bit(k, i) Then
256.             ok = mp_copy(x, temp)
257.             ok = mp_add(x, y, tmp)
258.             ok = mp_mod(tmp, m, x)
259.             ok = mp_copy(temp, y)
260.         End If
261.     Next
262.     ok = mp_copy(x, result)
263.     mp_clear tmp
264.     mp_clear temp
265.     mp_clear xx
266.     mp_clear y
267.     mp_clear x
268. End Sub
269.  
270. Function mp_str$ (n As mp_int, radix As Long)
271.     Dim sresult As String
272.     Dim As Long status, size
273.     status = mp_radix_size&(n, radix, size)
274.     sresult = Space$(size) + Chr$(0)
275.     status = mp_toradix_n(n, sresult, radix, size)
276.     If status = 0 Then
277.         mp_str$ = _Trim$(sresult)
278.     Else
279.         mp_str$ = "error in mp_toradix"
280.     End If
281. End Function
282.  
283. Sub mp_val (s As String, n As mp_int, radix As Long)
284.     Dim value As String
285.     Dim status As Long
286.     Dim As Long ok
287.     value = s + Chr$(0)
288.     status = mp_read_radix(n, value, radix)
289.     If status <> 0 Then Print "could not read number"
290. End Sub
291.