mpfr

[jack]

here's an almost complete include for the mpfr library https://www.mpfr.org, this is for Windows, datatypes are likely to be different on linux.
[linux version below near end of thread]
warning: not exhaustively tested
mpfr-include.bi

Code: QB64: [Select]

1. Const MPFR_RNDN = 0
2. Const MPFR_RNDZ = 1
3. Const MPFR_RNDU = 2
4. Const MPFR_RNDD = 3
5. Const MPFR_RNDA = 4
6. Const MPFR_RNDF = 5
7. Const MPFR_RNDNA = -1
8.  
9. Type mpz_struct
10.     mp_alloc As Long
11.     mp_size As Long
12.     mp_d As _Offset
13. End Type
14.  
15. Type mpq_struct
16.     mp_num As mpz_struct
17.     mp_den As mpz_struct
18. End Type
19.  
20. Type mpf_struct
21.     mp_prec As Long
22.     mp_size As Long
23.     mp_exp As Long
24.     mp_d As _Offset
25. End Type
26.  
27. Type mpfr_struct
28.     mpfr_prec As _Offset
29.     mpfr_sign As Long
30.     mpfr_exp As _Offset
31.     mpfr_d As _Offset
32. End Type
33.  
34. '    DECLARE DYNAMIC LIBRARY "libgmp-10"
35. '    END DECLARE
36.  
37. Declare Dynamic Library "libmpfr-6"
38.     Function mpfr_get_version$
39.     Function mpfr_get_patches$
40.     Function mpfr_buildopt_tune_case$
41.     Function mpfr_buildopt_tls_p&
42.     Function mpfr_buildopt_decimal_p&
43.     Function mpfr_buildopt_gmpinternals_p&
44.     Function mpfr_get_emin&
45.     Function mpfr_set_emin& (ByVal emin As Long)
46.     Function mpfr_get_emin_min&
47.     Function mpfr_get_emin_max&
48.     Function mpfr_get_emax&
49.     Function mpfr_set_emax& (ByVal emax As Long)
50.     Function mpfr_get_emax_min&
51.     Function mpfr_get_emax_max&
52.     Sub mpfr_set_default_rounding_mode (ByVal rm As Long)
53.     Function mpfr_get_default_rounding_mode&
54.     Function mpfr_print_rnd_mode$ (ByVal rndm As Long)
55.     Sub mpfr_clear_flags
56.     Sub mpfr_clear_underflow
57.     Sub mpfr_clear_overflow
58.     Sub mpfr_clear_divby0
59.     Sub mpfr_clear_nanflag
60.     Sub mpfr_clear_inexflag
61.     Sub mpfr_clear_erangeflag
62.     Sub mpfr_set_underflow
63.     Sub mpfr_set_overflow
64.     Sub mpfr_set_divby0
65.     Sub mpfr_set_nanflag
66.     Sub mpfr_set_inexflag
67.     Sub mpfr_set_erangeflag
68.     Function mpfr_underflow_p&
69.     Function mpfr_overflow_p&
70.     Function mpfr_divby0_p&
71.     Function mpfr_nanflag_p&
72.     Function mpfr_inexflag_p&
73.     Function mpfr_erangeflag_p&
74.     Function mpfr_check_range& (rop As mpfr_struct, Byval op2 As Long, Byval op3 As Long)
75.     Sub mpfr_init2 (rop As mpfr_struct, Byval op1 As Long)
76.     Sub mpfr_init (rop As mpfr_struct)
77.     Sub mpfr_clear (rop As mpfr_struct)
78.     Function mpfr_prec_round& (rop As mpfr_struct, Byval op1 As Long, Byval op2 As Long)
79.     Function mpfr_can_round& (rop As mpfr_struct, Byval op1 As Long, Byval op2 As Long, Byval op3 As Long, Byval op4 As Long)
80.     Function mpfr_min_prec& (rop As mpfr_struct)
81.     Function mpfr_get_exp& (rop As mpfr_struct)
82.     Function mpfr_set_exp& (rop As mpfr_struct, Byval op1 As Long)
83.     Function mpfr_get_prec& (rop As mpfr_struct)
84.     Sub mpfr_set_prec (rop As mpfr_struct, Byval mpfr_prec As Long)
85.     Sub mpfr_set_prec_raw (rop As mpfr_struct, Byval mpfr_prec As Long)
86.     Sub mpfr_set_default_prec (ByVal mpfr_prec As Long)
87.     Function mpfr_get_default_prec& ()
88.     Function mpfr_set_d& (rop As mpfr_struct, Byval x As Double, Byval mpfr_rnd As Long)
89.     Function mpfr_set_flt& (rop As mpfr_struct, Byval x As Single, Byval mpfr_rnd As Long)
90.     Function mpfr_set_ld& (rop As mpfr_struct, Byval x As _Float, Byval mpfr_rnd As Long)
91.     Function mpfr_set_z& (rop As mpfr_struct, opz As mpz_struct, Byval mpfr_rnd As Long)
92.     Function mpfr_set_z_2exp& (rop As mpfr_struct, opz As mpz_struct, Byval mpfr_exp_t As Long, Byval mpfr_rnd As Long)
93.     Sub mpfr_set_nan (rop As mpfr_struct)
94.     Sub mpfr_set_inf (rop As mpfr_struct, Byval n As Long)
95.     Sub mpfr_set_zero (rop As mpfr_struct, Byval n As Long)
96.     Function mpfr_set_f& (rop As mpfr_struct, mpf As mpf_struct, Byval mpfr_rnd As Long)
97.     Function mpfr_get_f& (rop As mpf_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
98.     Function mpfr_set_si& (rop As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
99.     Function mpfr_set_ui& (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
100.     Function mpfr_set_si_2exp& (rop As mpfr_struct, Byval n As Long, Byval mpfr_exp_t As Long, Byval mpfr_rnd As Long)
101.     Function mpfr_set_ui_2exp& (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_exp_t As Long, Byval mpfr_rnd As Long)
102.     Function mpfr_set_q& (rop As mpfr_struct, opq As mpq_struct, Byval mpfr_rnd As Long)
103.     Function mpfr_set_str& (rop As mpfr_struct, char As String, Byval nbase As Long, Byval rd As Long)
104.     Function mpfr_init_set_str& (rop As mpfr_struct, s As String, Byval n As Long, Byval mpfr_rnd As Long)
105.     Function mpfr_set4& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long, Byval n As Long)
106.     Function mpfr_abs& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
107.     Function mpfr_set& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
108.     Function mpfr_neg& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
109.     Function mpfr_signbit& (op As mpfr_struct)
110.     Function mpfr_setsign& (rop As mpfr_struct, op As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
111.     Function mpfr_copysign& (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
112.     Function mpfr_get_z_2exp& (z As mpz_struct, op As mpfr_struct)
113.     Function mpfr_get_flt! (op As mpfr_struct, Byval mpfr_rnd As Long)
114.     Function mpfr_get_d# (op As mpfr_struct, Byval mpfr_rnd As Long)
115.     Function mpfr_get_ld## (op As mpfr_struct, Byval mpfr_rnd As Long)
116.     Function mpfr_get_d1# (op As mpfr_struct)
117.     Function mpfr_get_d_2exp# (n As Long, op As mpfr_struct, Byval mpfr_rnd As Long)
118.     Function mpfr_get_ld_2exp## (n As Long, op As mpfr_struct, Byval mpfr_rnd As Long)
119.     Function mpfr_frexp& (n As Long, rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
120.     Function mpfr_get_si& (op As mpfr_struct, Byval mpfr_rnd As Long)
121.     Function mpfr_get_ui~& (op As mpfr_struct, Byval mpfr_rnd As Long)
122.     Function mpfr_get_str$ (s As String, n As Long, Byval m As Long, Byval k As _Unsigned Long, op As mpfr_struct, Byval mpfr_rnd As Long)
123.     Function mpfr_get_z& (z As mpz_struct, f As mpfr_struct, Byval mpfr_rnd As Long)
124.     Sub mpfr_nextabove (rop As mpfr_struct)
125.     Sub mpfr_nextbelow (rop As mpfr_struct)
126.     Sub mpfr_nexttoward (rop As mpfr_struct, op As mpfr_struct)
127.     Function mpfr_pow& (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
128.     Function mpfr_pow_si& (rop As mpfr_struct, op As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
129.     Function mpfr_pow_ui& (rop As mpfr_struct, op As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
130.     Function mpfr_ui_pow_ui& (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
131.     Function mpfr_ui_pow& (rop As mpfr_struct, Byval culong As _Unsigned Long, op As mpfr_struct, Byval mpfr_rnd As Long)
132.     Function mpfr_pow_z& (rop As mpfr_struct, op1 As mpfr_struct, z As mpz_struct, Byval mpfr_rnd As Long)
133.     Function mpfr_sqrt& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
134.     Function mpfr_sqrt_ui& (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
135.     Function mpfr_rec_sqrt& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
136.     Function mpfr_add& (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
137.     Function mpfr_sub& (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
138.     Function mpfr_mul& (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
139.     Function mpfr_div& (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
140.     Function mpfr_add_ui& (rop As mpfr_struct, op1 As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
141.     Function mpfr_sub_ui& (rop As mpfr_struct, op1 As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
142.     Function mpfr_ui_sub& (rop As mpfr_struct, Byval culong As _Unsigned Long, op As mpfr_struct, Byval mpfr_rnd As Long)
143.     Function mpfr_mul_ui& (rop As mpfr_struct, op1 As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
144.     Function mpfr_div_ui& (rop As mpfr_struct, op1 As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
145.     Function mpfr_ui_div& (rop As mpfr_struct, Byval culong As _Unsigned Long, op As mpfr_struct, Byval mpfr_rnd As Long)
146.     Function mpfr_add_si& (rop As mpfr_struct, op1 As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
147.     Function mpfr_sub_si& (rop As mpfr_struct, op1 As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
148.     Function mpfr_si_sub& (rop As mpfr_struct, Byval n As Long, op As mpfr_struct, Byval mpfr_rnd As Long)
149.     Function mpfr_mul_si& (rop As mpfr_struct, op1 As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
150.     Function mpfr_div_si& (rop As mpfr_struct, op1 As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
151.     Function mpfr_si_div& (rop As mpfr_struct, Byval n As Long, op As mpfr_struct, Byval mpfr_rnd As Long)
152.     Function mpfr_add_d& (rop As mpfr_struct, op1 As mpfr_struct, Byval x As Double, Byval mpfr_rnd As Long)
153.     Function mpfr_sub_d& (rop As mpfr_struct, op1 As mpfr_struct, Byval x As Double, Byval mpfr_rnd As Long)
154.     Function mpfr_d_sub& (rop As mpfr_struct, Byval x As Double, op As mpfr_struct, Byval mpfr_rnd As Long)
155.     Function mpfr_mul_d& (rop As mpfr_struct, op1 As mpfr_struct, Byval x As Double, Byval mpfr_rnd As Long)
156.     Function mpfr_div_d& (rop As mpfr_struct, op1 As mpfr_struct, Byval x As Double, Byval mpfr_rnd As Long)
157.     Function mpfr_d_div& (rop As mpfr_struct, Byval x As Double, op As mpfr_struct, Byval mpfr_rnd As Long)
158.     Function mpfr_sqr& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
159.     Function mpfr_const_pi& (rop As mpfr_struct, Byval op As Long)
160.     Function mpfr_const_log2& (rop As mpfr_struct, Byval mpfr_rnd As Long)
161.     Function mpfr_const_euler& (rop As mpfr_struct, Byval mpfr_rnd As Long)
162.     Function mpfr_const_catalan& (rop As mpfr_struct, Byval mpfr_rnd As Long)
163.     Function mpfr_agm& (rop As mpfr_struct, op1 As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
164.     Function mpfr_log& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
165.     Function mpfr_log2& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
166.     Function mpfr_log10& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
167.     Function mpfr_log1p& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
168.     Function mpfr_exp& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
169.     Function mpfr_exp2& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
170.     Function mpfr_exp10& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
171.     Function mpfr_expm1& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
172.     Function mpfr_eint& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
173.     Function mpfr_li2& (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
174.     Function mpfr_cmp& (op As mpfr_struct, op1 As mpfr_struct)
175.     Function mpfr_cmp3& (op As mpfr_struct, op1 As mpfr_struct, Byval n As Long)
176.     Function mpfr_cmp_d& (rop As mpfr_struct, Byval x As Double)
177.     Function mpfr_cmp_ld& (rop As mpfr_struct, Byval x As _Float)
178.     Function mpfr_cmpabs& (rop As mpfr_struct, op1 As mpfr_struct)
179.     Function mpfr_cmp_ui& (rop As mpfr_struct, Byval culong As _Unsigned Long)
180.     Function mpfr_cmp_si& (rop As mpfr_struct, Byval n As Long)
181.     Function mpfr_cmp_ui_2exp& (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval n As Long)
182.     Function mpfr_cmp_si_2exp& (rop As mpfr_struct, Byval cl As Long, Byval n As Long)
183.     Sub mpfr_reldiff (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
184.     Function mpfr_eq& (rop As mpfr_struct, op1 As mpfr_struct, Byval culong As _Unsigned Long)
185.     Function mpfr_sgn& (rop As mpfr_struct)
186.     Function mpfr_mul_2exp& (rop As mpfr_struct, op As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
187.     Function mpfr_div_2exp& (rop As mpfr_struct, op As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
188.     Function mpfr_mul_2ui& (rop As mpfr_struct, op As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
189.     Function mpfr_div_2ui& (rop As mpfr_struct, op As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
190.     Function mpfr_mul_2si& (rop As mpfr_struct, op As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
191.     Function mpfr_div_2si& (rop As mpfr_struct, op As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
192.     Function mpfr_rint& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
193.     Function mpfr_round& (rop As mpfr_struct, op As mpfr_struct)
194.     Function mpfr_trunc& (rop As mpfr_struct, op As mpfr_struct)
195.     Function mpfr_ceil& (rop As mpfr_struct, op As mpfr_struct)
196.     Function mpfr_floor& (rop As mpfr_struct, op As mpfr_struct)
197.     Function mpfr_rint_round& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
198.     Function mpfr_rint_trunc& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
199.     Function mpfr_rint_ceil& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
200.     Function mpfr_rint_floor& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
201.     Function mpfr_frac& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
202.     Function mpfr_modf& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
203.     Function mpfr_remquo& (rop As mpfr_struct, n As Long, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
204.     Function mpfr_remainder& (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
205.     Function mpfr_fmod& (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
206.     Function mpfr_fits_ulong_p& (op As mpfr_struct, Byval mpfr_rnd As Long)
207.     Function mpfr_fits_slong_p& (op As mpfr_struct, Byval mpfr_rnd As Long)
208.     Function mpfr_fits_uint_p& (op As mpfr_struct, Byval mpfr_rnd As Long)
209.     Function mpfr_fits_sint_p& (op As mpfr_struct, Byval mpfr_rnd As Long)
210.     Function mpfr_fits_ushort_p& (op As mpfr_struct, Byval mpfr_rnd As Long)
211.     Function mpfr_fits_sshort_p& (op As mpfr_struct, Byval mpfr_rnd As Long)
212.     Function mpfr_fits_uintmax_p& (op As mpfr_struct, Byval mpfr_rnd As Long)
213.     Function mpfr_fits_intmax_p& (op As mpfr_struct, Byval mpfr_rnd As Long)
214.     Sub mpfr_extract (z As mpz_struct, op As mpfr_struct, Byval n As _Unsigned Long)
215.     Sub mpfr_swap (rop As mpfr_struct, op As mpfr_struct)
216.     Sub mpfr_dump (op As mpfr_struct)
217.     Function mpfr_nan_p& (op As mpfr_struct)
218.     Function mpfr_inf_p& (op As mpfr_struct)
219.     Function mpfr_number_p& (op As mpfr_struct)
220.     Function mpfr_integer_p& (op As mpfr_struct)
221.     Function mpfr_zero_p& (op As mpfr_struct)
222.     Function mpfr_regular_p& (op As mpfr_struct)
223.     Function mpfr_greater_p& (op As mpfr_struct, op1 As mpfr_struct)
224.     Function mpfr_greaterequal_p& (op As mpfr_struct, op1 As mpfr_struct)
225.     Function mpfr_less_p& (op As mpfr_struct, op1 As mpfr_struct)
226.     Function mpfr_lessequal_p& (op As mpfr_struct, op1 As mpfr_struct)
227.     Function mpfr_lessgreater_p& (op As mpfr_struct, op1 As mpfr_struct)
228.     Function mpfr_equal_p& (op As mpfr_struct, op1 As mpfr_struct)
229.     Function mpfr_unordered_p& (op As mpfr_struct, op1 As mpfr_struct)
230.     Function mpfr_atanh& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
231.     Function mpfr_acosh& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
232.     Function mpfr_asinh& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
233.     Function mpfr_cosh& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
234.     Function mpfr_sinh& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
235.     Function mpfr_tanh& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
236.     Function mpfr_sinh_cosh& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
237.     Function mpfr_sech& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
238.     Function mpfr_csch& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
239.     Function mpfr_coth& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
240.     Function mpfr_acos& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
241.     Function mpfr_asin& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
242.     Function mpfr_atan& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
243.     Function mpfr_sin& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
244.     Function mpfr_sin_cos& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
245.     Function mpfr_cos& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
246.     Function mpfr_tan& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
247.     Function mpfr_atan2& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
248.     Function mpfr_sec& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
249.     Function mpfr_csc& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
250.     Function mpfr_cot& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
251.     Function mpfr_hypot& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
252.     Function mpfr_erf& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
253.     Function mpfr_erfc& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
254.     Function mpfr_cbrt& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
255.     Function mpfr_root& (rop As mpfr_struct, op As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
256.     Function mpfr_gamma& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
257.     Function mpfr_lngamma& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
258.     Function mpfr_lgamma& (rop As mpfr_struct, n As Long, op As mpfr_struct, Byval mpfr_rnd As Long)
259.     Function mpfr_digamma& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
260.     Function mpfr_zeta& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
261.     Function mpfr_zeta_ui& (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
262.     Function mpfr_fac_ui& (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
263.     Function mpfr_j0& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
264.     Function mpfr_j1& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
265.     Function mpfr_jn& (rop As mpfr_struct, Byval n As Long, op As mpfr_struct, Byval mpfr_rnd As Long)
266.     Function mpfr_y0& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
267.     Function mpfr_y1& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
268.     Function mpfr_yn& (rop As mpfr_struct, Byval n As Long, op As mpfr_struct, Byval mpfr_rnd As Long)
269.     Function mpfr_ai& (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
270.     Function mpfr_min& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
271.     Function mpfr_max& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
272.     Function mpfr_dim& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
273.     Function mpfr_mul_z& (rop As mpfr_struct, op As mpfr_struct, z As mpz_struct, Byval mpfr_rnd As Long)
274.     Function mpfr_div_z& (rop As mpfr_struct, op As mpfr_struct, z As mpz_struct, Byval mpfr_rnd As Long)
275.     Function mpfr_add_z& (rop As mpfr_struct, op As mpfr_struct, z As mpz_struct, Byval mpfr_rnd As Long)
276.     Function mpfr_sub_z& (rop As mpfr_struct, op As mpfr_struct, z As mpz_struct, Byval mpfr_rnd As Long)
277.     Function mpfr_z_sub& (rop As mpfr_struct, z As mpz_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
278.     Function mpfr_cmp_z& (op As mpfr_struct, z As mpz_struct)
279.     Function mpfr_mul_q& (rop As mpfr_struct, op As mpfr_struct, q As mpq_struct, Byval mpfr_rnd As Long)
280.     Function mpfr_div_q& (rop As mpfr_struct, op As mpfr_struct, q As mpq_struct, Byval mpfr_rnd As Long)
281.     Function mpfr_add_q& (rop As mpfr_struct, op As mpfr_struct, q As mpq_struct, Byval mpfr_rnd As Long)
282.     Function mpfr_sub_q& (rop As mpfr_struct, op As mpfr_struct, q As mpq_struct, Byval mpfr_rnd As Long)
283.     Function mpfr_cmp_q& (op As mpfr_struct, q As mpq_struct)
284.     Function mpfr_cmp_f& (op As mpfr_struct, f As mpf_struct)
285.     Function mpfr_fma& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
286.     Function mpfr_fms& (rop As mpfr_struct, op As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
287.     Function mpfr_sum& (rop As mpfr_struct, op1 As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
288.     Sub mpfr_free_cache ()
289.     Function mpfr_subnormalize& (rop As mpfr_struct, Byval n As Long, Byval mpfr_rnd As Long)
290.     Function mpfr_strtofr& (rop As mpfr_struct, s As String, s2 As String, Byval n As Long, Byval mpfr_rnd As Long)
291.     Function mpfr_custom_get_size~& (ByVal mpfr_prec As Long)
292.     Sub mpfr_custom_init (any As _Offset, Byval mpfr_prec As Long)
293.     Function mpfr_custom_get_significand%& (op As mpfr_struct)
294.     Function mpfr_custom_get_exp& (op As mpfr_struct)
295.     Sub mpfr_custom_move (rop As mpfr_struct, any As _Offset)
296.     Sub mpfr_custom_init_set (rop As mpfr_struct, Byval n As Long, Byval n2 As Long, Byval mpfr_prec As Long, any As _Offset)
297.     Function mpfr_custom_get_kind& (op As mpfr_struct)
298.     Function mpfr_snprintf& (rop As String, Byval n As _Unsigned _Offset, frmt As String, src2 As mpfr_struct)
299. End Declare
300.  
301. Sub print_mpfr_g (f As mpfr_struct, n As Long, m As Long)
302.     Dim sresult(0) As String, frmt As String
303.     Erase sresult
304.     Dim sresult(0) As String
305.     Dim ln As Long
306.     frmt = "%" + LTrim$(Str$(n)) + "." + LTrim$(Str$(m)) + "RNg" + Chr$(0)
307.     ln = mpfr_snprintf&("", 0, frmt, f)
308.     If ln > 0 Then
309.         sresult(0) = Space$(ln) + Chr$(0)
310.         ln = mpfr_snprintf&(sresult(0), ln + 4, frmt, f)
311.     End If
312.     Print Left$(sresult(0), ln)
313. End Sub
314.  
315. Sub print_mpfr_f (f As mpfr_struct, n As Long, m As Long)
316.     Dim sresult(0) As String, frmt As String
317.     Erase sresult
318.     Dim sresult(0) As String
319.     Dim ln As Long
320.     frmt = "%" + LTrim$(Str$(n)) + "." + LTrim$(Str$(m)) + "RNf" + Chr$(0)
321.     ln = mpfr_snprintf&("", 0, frmt, f)
322.     If ln > 0 Then
323.         sresult(0) = Space$(ln) + Chr$(0)
324.         ln = mpfr_snprintf&(sresult(0), ln + 4, frmt, f)
325.     End If
326.     Print Left$(sresult(0), ln)
327. End Sub
328.  
329. Sub print_mpfr_e (f As mpfr_struct, n As Long, m As Long)
330.     Dim sresult(0) As String, frmt As String
331.     Erase sresult
332.     Dim sresult(0) As String
333.     Dim ln As Long
334.     frmt = "%" + LTrim$(Str$(n)) + "." + LTrim$(Str$(m)) + "RNe" + Chr$(0)
335.     ln = mpfr_snprintf&("", 0, frmt, f)
336.     If ln > 0 Then
337.         sresult(0) = Space$(ln) + Chr$(0)
338.         ln = mpfr_snprintf&(sresult(0), ln + 4, frmt, f)
339.     End If
340.     Print Left$(sresult(0), ln)
341. End Sub
342.  
343. Sub valmpfr (result As mpfr_struct, value As String, nbase As Long)
344.     Dim c As Long
345.     c = mpfr_set_str(result, value, nbase, 0)
346. End Sub
347.  

simple-test.bas

Code: QB64: [Select]

1. '$CONSOLE:ONLY
2. '_DEST _CONSOLE
3.  
4. WIDTH 90, 43
5. _SCREENMOVE 0, 0
6.  
7. main
8.  
9. REM $include: 'mpfr-include.bi'
10.  
11. SUB main
12.     DIM x AS mpfr_struct, y AS mpfr_struct
13.     DIM c AS LONG
14.  
15.     CALL mpfr_init2(x, 1024)
16.     CALL mpfr_init2(y, 1024)
17.  
18.     FOR i = 1 TO 20
19.         c = mpfr_set_si(x, i, MPFR_RNDN)
20.         c = mpfr_log10(y, x, MPFR_RNDN)
21.         PRINT i;
22.         CALL print_mpfr_f(y, 77, 70)
23.     NEXT
24.  
25.     CALL mpfr_clear(y)
26.     CALL mpfr_clear(x)
27. END SUB
28.  

[jack]

a more complex example
lanczos gamma approximation

Code: QB64: [Select]

1.  
2. $CONSOLE:ONLY
3. _DEST _CONSOLE
4.  
5. DIM SHARED prec AS LONG
6. DIM SHARED error_code AS LONG
7. main
8.  
9. REM $include: 'mpfr-include.bi'
10. SUB main
11.     DIM c AS LONG, n AS LONG, i AS LONG
12.  
13.     prec = 1024 ' bits precision, about 300 decimal digits
14.     n = 6
15.     DIM x(n) AS mpfr_struct, g AS mpfr_struct
16.     FOR i = 0 TO n
17.         CALL mpfr_init2(x(i), prec)
18.     NEXT
19.     CALL mpfr_init2(g, prec)
20.     c = mpfr_set_si(g, 5, MPFR_RNDN)
21.     CALL make_lanczos(x(), g, error_code)
22.     PRINT "lanczos approximation for the gamma function"
23.     PRINT "          n = 6, g = 5"
24.     FOR i = 0 TO n
25.         PRINT USING "C(##) = "; i;
26.         CALL print_mpfr_f(x(i), 27, 20)
27.     NEXT
28.     FOR i = 0 TO n
29.         CALL mpfr_clear(x(i))
30.     NEXT
31.  
32.     PRINT "==============================================="
33.     PRINT "lanczos approximation for the gamma function"
34.     PRINT "          n = 14, g = 4.7421875"
35.     ERASE x
36.     n = 14
37.     DIM x(n) AS mpfr_struct
38.     FOR i = 0 TO n
39.         CALL mpfr_init2(x(i), prec)
40.     NEXT
41.     CALL valmpfr(g, "4.7421875", 10)
42.     CALL make_lanczos(x(), g, error_code)
43.     FOR i = 0 TO n
44.         PRINT USING "C(##) = "; i;
45.         CALL print_mpfr_f(x(i), 27, 20)
46.     NEXT
47.     FOR i = 0 TO n
48.         CALL mpfr_clear(x(i))
49.     NEXT
50.     CALL mpfr_clear(g)
51. END SUB
52.  
53. SUB make_lanczos (x() AS mpfr_struct, g AS mpfr_struct, error_code AS LONG)
54.  
55.     DIM i AS LONG, j AS LONG, n AS LONG, c AS LONG
56.     n = UBOUND(x)
57.     DIM a(n, n) AS mpfr_struct
58.     DIM b(n) AS mpfr_struct
59.     DIM sqr2pi AS mpfr_struct, pi2 AS mpfr_struct, fac AS mpfr_struct
60.     DIM one AS mpfr_struct, tmp AS mpfr_struct, tmp2 AS mpfr_struct
61.     FOR i = 0 TO n
62.         FOR j = 0 TO n
63.             CALL mpfr_init2(a(i, j), prec)
64.         NEXT
65.         CALL mpfr_init2(b(i), prec)
66.     NEXT
67.     CALL mpfr_init2(sqr2pi, prec)
68.     CALL mpfr_init2(pi2, prec)
69.     CALL mpfr_init2(fac, prec)
70.     CALL mpfr_init2(one, prec)
71.     CALL mpfr_init2(tmp, prec)
72.     CALL mpfr_init2(tmp2, prec)
73.     c = mpfr_set_si(one, 1, MPFR_RNDN)
74.     c = mpfr_const_pi(pi2, MPFR_RNDN)
75.     c = mpfr_mul_si(pi2, pi2, 2, MPFR_RNDN)
76.     c = mpfr_sqrt(sqr2pi, pi2, MPFR_RNDN)
77.  
78.     'setup the matrix, for n=6 the matrix looks like this
79.  
80.     ' [[ 1, 1,   1/2, 1/3, 1/4,  1/5,  1/6 ]
81.     ' [  1, 1/2, 1/3, 1/4, 1/5,  1/6,  1/7 ]
82.     ' [  1, 1/3, 1/4, 1/5, 1/6,  1/7,  1/8 ]
83.     ' [  1, 1/4, 1/5, 1/6, 1/7,  1/8,  1/9 ]
84.     ' [  1. 1/5, 1/6, 1/7, 1/8,  1/9,  1/10 ]
85.     ' [  1, 1/6, 1/7, 1/8, 1/9,  1/10, 1/11 ]
86.     ' [  1, 1/7, 1/8, 1/9, 1/10, 1/11, 1/12 ]]
87.  
88.     FOR i = 0 TO n
89.         c = mpfr_set_si(a(i, 0), 1, MPFR_RNDN)
90.         FOR j = 1 TO n
91.             c = mpfr_div_si(a(i, j), one, i + j, MPFR_RNDN)
92.         NEXT
93.     NEXT
94.  
95.     'setup the constant terms, for n=6 and g=5 they are
96.     'b( 0) = 41.62443691643906820752
97.     'b( 1) = 16.01231640525168135809
98.     'b( 2) =  9.36473553710404851170
99.     'b( 3) =  6.57049625931606149987
100.     'b( 4) =  5.09521672929646781312
101.     'b( 5) =  4.20380175313001102196
102.     'b( 6) =  3.61487381446333490266
103.  
104.     c = mpfr_set_si(fac, 1, MPFR_RNDN)
105.  
106.     FOR i = 0 TO n
107.         IF i > 1 THEN
108.             c = mpfr_mul_si(fac, fac, i, MPFR_RNDN)
109.         END IF
110.         'b(i)=exp(i+half+g)/(i+half+g)^(i+half)*fac/sqr2pi
111.         c = mpfr_add_d(tmp, g, 0.5, MPFR_RNDN)
112.         c = mpfr_add_si(tmp, tmp, i, MPFR_RNDN)
113.         c = mpfr_exp(b(i), tmp, MPFR_RNDN)
114.         c = mpfr_set_d(tmp2, i + 0.5, MPFR_RNDN)
115.         c = mpfr_pow(tmp, tmp, tmp2, MPFR_RNDN)
116.         c = mpfr_div(b(i), b(i), tmp, MPFR_RNDN)
117.         c = mpfr_mul(b(i), b(i), fac, MPFR_RNDN)
118.         c = mpfr_div(b(i), b(i), sqr2pi, MPFR_RNDN)
119.         'CALL print_mpfr_f(b(i), 27, 20)
120.     NEXT
121.  
122.     'now solve the system of linear equations, for n=6 and g=5 they are
123.     ' [[ 1, 1,   1/2, 1/3, 1/4,  1/5,  1/6 ]    = 41.62443691643906820752
124.     ' [  1, 1/2, 1/3, 1/4, 1/5,  1/6,  1/7 ]    = 16.01231640525168135809
125.     ' [  1, 1/3, 1/4, 1/5, 1/6,  1/7,  1/8 ]    =  9.36473553710404851170
126.     ' [  1, 1/4, 1/5, 1/6, 1/7,  1/8,  1/9 ]    =  6.57049625931606149987
127.     ' [  1. 1/5, 1/6, 1/7, 1/8,  1/9,  1/10 ]   =  5.09521672929646781312
128.     ' [  1, 1/6, 1/7, 1/8, 1/9,  1/10, 1/11 ]   =  4.20380175313001102196
129.     ' [  1, 1/7, 1/8, 1/9, 1/10, 1/11, 1/12 ]]  =  3.61487381446333490266
130.  
131.     'and the solution for the above example is
132.     'C(0) =   1.00000000019001482399
133.     'C(1) =  76.18009172947146348309
134.     'C(2) = -86.50532032941676765250
135.     'C(3) =  24.01409824083091049018
136.     'C(4) =  -1.23173957245015538752
137.     'C(5) =   0.00120865097386617851
138.     'C(6) =  -0.00000539523938495313
139.  
140.     'which is the Lanczos gamma function approximation
141.     'as published in Numerical recipes
142.  
143.     CALL gauss_jordan(a(), b(), x(), n, error_code)
144.  
145.     FOR i = 0 TO n
146.         FOR j = 0 TO n
147.             CALL mpfr_clear(a(i, j))
148.         NEXT
149.         CALL mpfr_clear(b(i))
150.     NEXT
151.     CALL mpfr_clear(sqr2pi)
152.     CALL mpfr_clear(pi2)
153.     CALL mpfr_clear(fac)
154.     CALL mpfr_clear(one)
155.     CALL mpfr_clear(tmp)
156.     CALL mpfr_clear(tmp2)
157. END SUB
158.  
159.  
160.  
161. SUB gauss_jordan (a() AS mpfr_struct, y() AS mpfr_struct, coef() AS mpfr_struct, ncol AS LONG, error_code AS LONG)
162.  
163.     ' matrix solution by Gaussian Elimination
164.  
165.     DIM b(ncol, ncol) AS mpfr_struct, w(ncol) AS mpfr_struct ' work array, ncol long
166.  
167.     DIM i AS LONG, j AS LONG, i1 AS LONG, k AS LONG, l AS LONG, n AS LONG, c AS LONG
168.     DIM hold AS mpfr_struct, sm AS mpfr_struct, tmp AS mpfr_struct
169.     DIM t AS mpfr_struct, ab AS mpfr_struct, big AS mpfr_struct
170.     CONST TRUE = -1
171.     CONST FALSE = NOT TRUE
172.  
173.     error_code = FALSE
174.     n = ncol
175.  
176.     FOR i = 0 TO ncol
177.         FOR j = 0 TO ncol
178.             CALL mpfr_init2(b(i, j), prec)
179.         NEXT
180.         CALL mpfr_init2(w(i), prec)
181.     NEXT
182.     CALL mpfr_init2(hold, prec)
183.     CALL mpfr_init2(sm, prec)
184.     CALL mpfr_init2(t, prec)
185.     CALL mpfr_init2(ab, prec)
186.     CALL mpfr_init2(big, prec)
187.     CALL mpfr_init2(tmp, prec)
188.  
189.     FOR i = 0 TO n
190.         ' copy to work arrays
191.         FOR j = 0 TO n
192.             'b(i,j)=a(i,j)
193.             c = mpfr_set(b(i, j), a(i, j), MPFR_RNDN)
194.         NEXT j
195.         'w(i)=y(i)
196.         c = mpfr_set(w(i), y(i), MPFR_RNDN)
197.         'coef(i)=zero
198.         c = mpfr_set_si(coef(i), 0, MPFR_RNDN)
199.     NEXT
200.     FOR i = 0 TO n - 1
201.         c = mpfr_abs(big, b(i, i), MPFR_RNDN)
202.  
203.         l = i
204.         i1 = i + 1
205.         FOR j = i1 TO n
206.             ' search for largest element
207.             c = mpfr_abs(ab, b(j, i), MPFR_RNDN)
208.             IF mpfr_cmp(ab, big) > 0 THEN
209.                 c = mpfr_set(big, ab, MPFR_RNDN)
210.                 l = j
211.             END IF
212.         NEXT
213.         IF mpfr_cmp_si(big, 0) = 0 THEN
214.             error_code = TRUE
215.         ELSE
216.             IF l <> i THEN
217.                 ' interchange rows to put largest element on diagonal
218.                 FOR j = 0 TO n
219.                     c = mpfr_set(hold, b(l, j), MPFR_RNDN)
220.                     c = mpfr_set(b(l, j), b(i, j), MPFR_RNDN)
221.                     c = mpfr_set(b(i, j), hold, MPFR_RNDN)
222.                 NEXT
223.                 c = mpfr_set(hold, w(l), MPFR_RNDN)
224.                 c = mpfr_set(w(l), w(i), MPFR_RNDN)
225.                 c = mpfr_set(w(i), hold, MPFR_RNDN)
226.             END IF
227.             FOR j = i1 TO n
228.                 c = mpfr_div(t, b(j, i), b(i, i), MPFR_RNDN)
229.                 FOR k = i1 TO n
230.                     c = mpfr_mul(tmp, t, b(i, k), MPFR_RNDN)
231.                     c = mpfr_sub(b(j, k), b(j, k), tmp, MPFR_RNDN)
232.                 NEXT
233.                 c = mpfr_mul(tmp, t, w(i), MPFR_RNDN)
234.                 c = mpfr_sub(w(j), w(j), tmp, MPFR_RNDN)
235.             NEXT ' j-loop
236.         END IF ' if big - else
237.     NEXT ' i-loop
238.     IF mpfr_cmp_si(b(n, n), 0) = 0 THEN
239.         error_code = TRUE
240.     ELSE
241.         c = mpfr_div(coef(n), w(n), b(n, n), MPFR_RNDN)
242.         i = n - 1
243.         ' back substitution
244.         DO
245.             c = mpfr_set_si(sm, 0, MPFR_RNDN)
246.             FOR j = i TO n
247.                 c = mpfr_mul(tmp, b(i, j), coef(j), MPFR_RNDN)
248.                 c = mpfr_add(sm, sm, tmp, MPFR_RNDN)
249.             NEXT
250.             c = mpfr_sub(tmp, w(i), sm, MPFR_RNDN)
251.             c = mpfr_div(coef(i), tmp, b(i, i), MPFR_RNDN)
252.             i = i - 1
253.         LOOP UNTIL i < 0
254.     END IF
255.     IF error_code THEN PRINT "ERROR: Matrix is singular"
256.    
257.     FOR i = 0 TO ncol
258.         FOR j = 0 TO ncol
259.             CALL mpfr_clear(b(i, j))
260.         NEXT
261.         CALL mpfr_clear(w(i))
262.     NEXT
263.     CALL mpfr_clear(hold)
264.     CALL mpfr_clear(sm)
265.     CALL mpfr_clear(t)
266.     CALL mpfr_clear(ab)
267.     CALL mpfr_clear(big)
268.     CALL mpfr_clear(tmp)
269. END SUB
270.  

the output

Code: QB64: [Select]

1. C( 0 ) =      1.00000000019001482399
2. C( 1 ) =     76.18009172947146348309
3. C( 2 ) =    -86.50532032941676765250
4. C( 3 ) =     24.01409824083091049018
5. C( 4 ) =     -1.23173957245015538752
6. C( 5 ) =      0.00120865097386617851
7. C( 6 ) =     -0.00000539523938495313
8.  

example of using the approximation, see https://en.wikipedia.org/wiki/Lanczos_approximation

Code: QB64: [Select]

1. DIM c(6) AS DOUBLE
2. DIM y AS DOUBLE, z AS DOUBLE, sqrt2pi AS DOUBLE, g AS DOUBLE
3. DIM i AS LONG
4.  
5. sqrt2pi = 2.506628274631000502
6. g = 5
7. c(0) = 1.000000000190014824
8. c(1) = 76.180091729471463483
9. c(2) = -86.505320329416767652
10. c(3) = 24.01409824083091049
11. c(4) = -1.2317395724501553875
12. c(5) = 0.0012086509738661785061
13. c(6) = -5.3952393849531283785E-06
14.  
15. INPUT "enter z "; z
16. y = c(0)
17. FOR i = 1 TO 6
18.     y = y + c(i) / (z + i)
19. NEXT
20. y = y * sqrt2pi * (z + g + 0.5) ^ (z + 0.5) * EXP(-(z + g + 0.5))
21. PRINT "gamma( 1 + "; z; ") = "; y
22.  

[jack]

I found the Spouge algorithm for the Gamma function interesting https://en.wikipedia.org/wiki/Spouge%27s_approximation
the coefficients are easy to calculate but don't give very accurate results

let's take a = 7, we have the coefficients ( after dividing the coefficients by Sqr(2*Pi) )

Code: QB64: [Select]

1. C( 0) =  1
2. C( 1) =  394.2326437490294715
3. C( 2) = -661.96874077331269907
4. C( 3) =  348.50416767875391816
5. C( 4) = -62.454919463566450249
6. C( 5) =  2.7792189888444614869
7. C( 6) = -0.0090369795951602295551
8.  

some tests

Code: QB64: [Select]

1. DIM c(6) AS double
2. DIM y AS double, z AS double, a AS double
3. a = 7
4. DIM i AS LONG
5. dim as double fac
6. fac =1
7. dim as double sqrt2pi
8. sqrt2pi = 2.5066282746310005024157652848110452530
9.  
10. C( 0) =  1
11. C( 1) =  394.2326437490294715
12. C( 2) = -661.96874077331269907
13. C( 3) =  348.50416767875391816
14. C( 4) = -62.454919463566450249
15. C( 5) =  2.7792189888444614869
16. C( 6) = -0.0090369795951602295551
17.  
18. for z=1 to 10
19.         fac=fac*z
20.         y = c(0)
21.         FOR i = 1 TO a-1
22.                 y = y + c(i) / (z + i)
23.         NEXT
24.         y =  y * (z+a)^(z+0.5)*exp(-(z+a))*sqrt2pi
25.         PRINT "spouge-gamma("; z; ") = "; y, (y-fac)/fac
26. next
27.  
28. z=0.5
29. fac=0.88622692545275801364908374167057259140
30.  
31. y = c(0)
32. FOR i = 1 TO a-1
33.         y = y + c(i) / (z + i)
34. NEXT
35. y = y * (z + a)^(z+0.5)*exp(-(z+a))*sqrt2pi
36. PRINT "spouge-gamma("; z; ") = "; y, (y-fac)/fac
37.  

the results

spouge-gamma( 1) =  0.9999999996137302    -3.862697939283066e-10
spouge-gamma( 2) =  1.999999998666582     -6.667089103018498e-10
spouge-gamma( 3) =  5.999999994337436     -9.437606252049591e-10
spouge-gamma( 4) =  23.99999997127829     -1.19673811601236e-09
spouge-gamma( 5) =  119.999999829899      -1.417508125693227e-09
spouge-gamma( 6) =  719.9999988446054     -1.604714662385555e-09
spouge-gamma( 7) =  5039.999991128157     -1.760286392839459e-09
spouge-gamma( 8) =  40319.99992389513     -1.887521453422128e-09
spouge-gamma( 9) =  362879.9992778363     -1.99008957806045e-09
spouge-gamma( 10) =  3628799.992482762    -2.071549184825699e-09
spouge-gamma( 0.5) =  0.8862269252253806  -2.565679778746878e-10

now, let's use the lanczos scheme to calculate the coefficients

Code: QB64: [Select]

1. 1, 1/1,1/2,1/3,1/4, 1/5, 1/6    =       Factorial(0)/(7^(1/2)*Exp(-7)*Sqr(2*Pi))
2. 1, 1/2,1/3,1/4,1/5, 1/6, 1/7    =       Factorial(1)/(( 1+7)^( 1+(1/2))*Exp(-( 1+7))*Sqr(2*Pi))
3. 1, 1/3,1/4,1/5,1/6, 1/7, 1/8    =       Factorial(2)/(( 2+7)^( 2+(1/2))*Exp(-( 2+7))*Sqr(2*Pi))
4. 1, 1/4,1/5,1/6,1/7, 1/8, 1/9    =       Factorial(3)/(( 3+7)^( 3+(1/2))*Exp(-( 3+7))*Sqr(2*Pi))
5. 1, 1/5,1/6,1/7,1/8, 1/9, 1/10   =       Factorial(4)/(( 4+7)^( 4+(1/2))*Exp(-( 4+7))*Sqr(2*Pi))
6. 1, 1/6,1/7,1/8,1/9, 1/10,1/11   =       Factorial(5)/(( 5+7)^( 5+(1/2))*Exp(-( 5+7))*Sqr(2*Pi))
7. 1, 1/7,1/8,1/9,1/10,1/11,1/12   =       Factorial(6)/(( 6+7)^( 6+(1/2))*Exp(-( 6+7))*Sqr(2*Pi))
8.  

the solution

Code: QB64: [Select]

1. C( 0) =  0.99999999995273358926
2. C( 1) =  394.23264374929666793
3. C( 2) = -661.96874078550137586
4. C( 3) =  348.50416780528271853
5. C( 4) = -62.454919950924562754
6. C( 5) =  2.7792197239908234111
7. C( 6) = -0.0090372064455749738392
8.  

the results after replacing the coeffients

spouge-lanczos-gamma( 1) =  0.9999999999999999    -1.110223024625157e-16
spouge-lanczos-gamma( 2) =  1.999999999999999     -5.551115123125783e-16
spouge-lanczos-gamma( 3) =  5.999999999999997     -4.440892098500626e-16
spouge-lanczos-gamma( 4) =  24                             -1.480297366166875e-16
spouge-lanczos-gamma( 5) =  119.9999999999999     -4.736951571734001e-16
spouge-lanczos-gamma( 6) =  719.9999999999998     -3.157967714489334e-16
spouge-lanczos-gamma( 7) =  5039.999999999991     -1.804552979708191e-15
spouge-lanczos-gamma( 8) =  40319.99999999953     -1.154913907013242e-14
spouge-lanczos-gamma( 9) =  362879.9999999893     -2.951446651256063e-14
spouge-lanczos-gamma( 10) =  3628799.999999764    -6.506015009507931e-14
spouge-lanczos-gamma( 0.5) =  0.8862269254527551  -3.382431835905346e-15

[_vince]

Need help getting this to work on linux

I'm getting

Code: QB64: [Select]

1. DYNAMIC LIBRARY not found
2. Caused by (or after):DECLARE DYNAMIC LIBRARY "libmpfr",7
3. LINE 37:DECLARE DYNAMIC LIBRARY "libmpfr"

I have libmpfr installed here:

/usr/lib/libmpfr.so.6
/usr/lib/libmpfr.so
/usr/lib/libmpfr.so.6.0.1

none of those work when put in DECLARE LIBRARY

I was hoping to demonstrate an interesting example with fractals since I've used GMP library in C for rendering fractals

[SMcNeill]

Need help getting this to work on linux

I'm getting

Code: QB64: [Select]

1. DYNAMIC LIBRARY not found
2. Caused by (or after):DECLARE DYNAMIC LIBRARY "libmpfr",7
3. LINE 37:DECLARE DYNAMIC LIBRARY "libmpfr"

I have libmpfr installed here:

/usr/lib/libmpfr.so.6
/usr/lib/libmpfr.so
/usr/lib/libmpfr.so.6.0.1

none of those work when put in DECLARE LIBRARY

I was hoping to demonstrate an interesting example with fractals since I've used GMP library in C for rendering fractals

I wish you the best of luck in sorting DECLARE LIBRARY out with Linux.  Try all I might, I never could get it to work with libpng.so, no matter where I put it, or how I declared it.  If somebody can help you get yours working, maybe their advice can help me get mine working as well, so we can do PNG exports on Linux systems.  ;)

[jack]

@_vince
what linux distribution do you use?
32 or 64-bit ?
the problem in my windows translation is that I replaced types like mpfr_rnd_t with long (4 bytes), some of the type's may have different size in linux, will check
too bad that QB64 doesn't have a typedef statement and I am not sure if it could be accomplished with meta commands

[jack]

as I suspected, there are different type siszes in Win64 and Linux64
Win64

SizeOf int 4
SizeOf unsigned int 4
SizeOf long int 4
SizeOf unsigned long int 4
SizeOf long 4
SizeOf unsigned long 4
SizeOf long long int 8
SizeOf unsigned long long int 8
=====================
=========== mpz_struct ==========
SizeOf __mpz_struct 16
SizeOf int 4
SizeOf int 4
SizeOf mp_limb_t 8
=====================
=========== mpq_struct ==========
SizeOf __mpq_struct 32
SizeOf __mpz_struct 16
SizeOf __mpz_struct 16
=====================
=========== mpf_struct ==========
SizeOf __mpf_struct 24
SizeOf int 4
SizeOf int 4
SizeOf mp_exp_t 4
SizeOf mp_limb_t 8
=====================
=========== gmp_randstate_struc ==========
SizeOf __gmp_randstate_struct 32
SizeOf mpz_t 16
SizeOf gmp_randalg_t 4
=====================
=========== mpfr_struct ==========
SizeOf __mpfr_struct 24
SizeOf mpfr_prec_t 4
SizeOf mpfr_sign_t 4
SizeOf mpfr_exp_t 4
SizeOf mp_limb_t 8
=====================
SizeOf mp_bitcnt_t 4
SizeOf size_t 8
SizeOf mp_limb_t 8
SizeOf mp_exp_t 4
SizeOf mp_size_t 4
SizeOf mpfr_t 24
SizeOf mpfr_rnd_t 4
SizeOf mpfr_void 1
SizeOf mpfr_int 4
SizeOf mpfr_long 4
SizeOf mpfr_sign_t 4
SizeOf long double 16

Ubuntu64

SizeOf int 4
SizeOf unsigned int 4
SizeOf long int 8
SizeOf unsigned long int 8
SizeOf long 8
SizeOf unsigned long 8
SizeOf long long int 8
SizeOf unsigned long long int 8
=====================
=========== mpz_struct ==========
SizeOf __mpz_struct 16
SizeOf int 4
SizeOf int 4
SizeOf mp_limb_t 8
=====================
=========== mpq_struct ==========
SizeOf __mpq_struct 32
SizeOf __mpz_struct 16
SizeOf __mpz_struct 16
=====================
=========== mpf_struct ==========
SizeOf __mpf_struct 24
SizeOf int 4
SizeOf int 4
SizeOf mp_exp_t 8
SizeOf mp_limb_t 8
=====================
=========== gmp_randstate_struc ==========
SizeOf __gmp_randstate_struct 32
SizeOf mpz_t 16
SizeOf gmp_randalg_t 4
=====================
=========== mpfr_struct ==========
SizeOf __mpfr_struct 32
SizeOf mpfr_prec_t 8
SizeOf mpfr_sign_t 4
SizeOf mpfr_exp_t 8
SizeOf mp_limb_t 8
=====================
SizeOf mp_bitcnt_t 8
SizeOf size_t 8
SizeOf mp_limb_t 8
SizeOf mp_exp_t 8
SizeOf mp_size_t 8
SizeOf mpfr_t 32
SizeOf mpfr_rnd_t 4
SizeOf mpfr_void 1
SizeOf mpfr_int 4
SizeOf mpfr_long 8
SizeOf mpfr_sign_t 4
SizeOf long double 16

the size of mpfr_struc is 24-bytes in Win64 but 32-bytes in Linux64
the SizeOf long int is 4 in Win64 but 8 in Linux64
the SizeOf mpfr_prec_t is 4 in Win64 but 8 in Linux64
the SizeOf mpfr_exp_t is 4 in Win64 but 8 in Linux64

so I will have to check each declare to make sure the proper size integer is used, that will take time

[jack]

Need help getting this to work on linux

I'm getting

Code: QB64: [Select]

1. DYNAMIC LIBRARY not found
2. Caused by (or after):DECLARE DYNAMIC LIBRARY "libmpfr",7
3. LINE 37:DECLARE DYNAMIC LIBRARY "libmpfr"

I have libmpfr installed here:

/usr/lib/libmpfr.so.6
/usr/lib/libmpfr.so
/usr/lib/libmpfr.so.6.0.1

none of those work when put in DECLARE LIBRARY

I was hoping to demonstrate an interesting example with fractals since I've used GMP library in C for rendering fractals

in your case, since the libraries are in /usr/lib, all you need to do is remove "lib" from the lib name, that is DECLARE DYNAMIC LIBRARY "mpfr" but my include needs fixing before you can use it in linux

@SMcNeill
as you may have guessed, try using png instead of libpng, if libpng is not in /usr/lib then I suggest that you make a symlink, for example sudo ln -s /usr/lib/x86_64-linux-gnu/libpng.so /usr/lib/libpng.so

[jack]

tested on Ubuntu both 32-bit and 64 bit
note: you can either hard-code the library path in the include or make symlinks of libgmp.so and libmpfr.so to /usr/lib
in my case, the libraries in question reside in /usr/lib/i386-linux-gnu/gmp and /usr/lib/i386-linux-gnu/mpfr in Ubuntu-i386
and in /usr/lib/x86_64-linux-gnu/gmp and /usr/lib/x86_64-linux-gnu/mpfr
so I would declare the library as DECLARE DYNAMIC LIBRARY "i386-linux-gnu/gmp" and DECLARE DYNAMIC LIBRARY "i386-linux-gnu/mpfr"
or DECLARE DYNAMIC LIBRARY "x86_64-linux-gnu/gmp" and DECLARE DYNAMIC LIBRARY "x86_64-linux-gnu/mpfr"

mpfr-include.bi

Code: QB64: [Select]

1. CONST MPFR_RNDN = 0
2. CONST MPFR_RNDZ = 1
3. CONST MPFR_RNDU = 2
4. CONST MPFR_RNDD = 3
5. CONST MPFR_RNDA = 4
6. CONST MPFR_RNDF = 5
7. CONST MPFR_RNDNA = -1
8.  
9. TYPE mpz_struct
10.     mp_alloc AS LONG
11.     mp_size AS LONG
12.     mp_d AS _OFFSET
13. END TYPE
14.  
15. TYPE mpq_struct
16.     mp_num AS mpz_struct
17.     mp_den AS mpz_struct
18. END TYPE
19.  
20. TYPE mpf_struct
21.     mp_prec AS LONG
22.     mp_size AS LONG
23.     mp_exp AS _OFFSET
24.     mp_d AS _OFFSET
25. END TYPE
26.  
27. TYPE mpfr_struct
28.     mpfr_prec AS _OFFSET
29.     mpfr_sign AS _OFFSET 'LONG
30.     mpfr_exp AS _OFFSET
31.     mpfr_d AS _OFFSET
32. END TYPE
33.  
34. DECLARE DYNAMIC LIBRARY "gmp"
35. END DECLARE
36.  
37. DECLARE DYNAMIC LIBRARY "mpfr"
38.     SUB mpfr_init2 (x AS mpfr_struct, BYVAL prec AS _OFFSET)
39.     'Sub mpfr_inits2 (byval prec As _offset, x AS mpfr_struct, ...)
40.     SUB mpfr_clear (x AS mpfr_struct)
41.     'Sub mpfr_clears (x AS mpfr_struct, ...)
42.     SUB mpfr_init (x AS mpfr_struct)
43.     'Sub mpfr_inits (x AS mpfr_struct, ...)
44.     SUB mpfr_set_default_prec (BYVAL prec AS _OFFSET)
45.     FUNCTION mpfr_get_default_prec%&
46.     SUB mpfr_set_prec (x AS mpfr_struct, BYVAL prec AS _OFFSET)
47.     FUNCTION mpfr_get_prec%& (x AS mpfr_struct)
48.     FUNCTION mpfr_set& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
49.     FUNCTION mpfr_set_ui& (rop AS mpfr_struct, BYVAL op AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
50.     FUNCTION mpfr_set_si& (rop AS mpfr_struct, BYVAL op AS _OFFSET, BYVAL rndm AS LONG)
51.     ''FUNCTION mpfr_set_uj& (rop AS mpfr_struct, BYVAL op AS _UNSIGNED _INTEGER64, BYVAL rndm AS LONG)
52.     ''FUNCTION mpfr_set_sj& (rop AS mpfr_struct, BYVAL op AS _INTEGER64, BYVAL rndm AS LONG)
53.     FUNCTION mpfr_set_flt& (rop AS mpfr_struct, BYVAL op AS SINGLE, BYVAL rndm AS LONG)
54.     FUNCTION mpfr_set_d& (rop AS mpfr_struct, BYVAL op AS DOUBLE, BYVAL rndm AS LONG)
55.     FUNCTION mpfr_set_ld& (rop AS mpfr_struct, BYVAL op AS _FLOAT, BYVAL rndm AS LONG)
56.     'Function mpfr_set_float128& (rop AS mpfr_struct, __float128 op, Byval rndm as long )
57.     'Function mpfr_set_decimal64& (rop AS mpfr_struct, _Decimal64 op, Byval rndm as long )
58.     FUNCTION mpfr_set_z& (rop AS mpfr_struct, op AS mpz_struct, BYVAL rndm AS LONG)
59.     FUNCTION mpfr_set_q& (rop AS mpfr_struct, op AS mpq_struct, BYVAL rndm AS LONG)
60.     FUNCTION mpfr_set_f& (rop AS mpfr_struct, op AS mpf_struct, BYVAL rndm AS LONG)
61.     FUNCTION mpfr_set_ui_2exp& (rop AS mpfr_struct, BYVAL op AS _UNSIGNED _OFFSET, BYVAL e AS _OFFSET, BYVAL rndm AS LONG)
62.     FUNCTION mpfr_set_si_2exp& (rop AS mpfr_struct, BYVAL op AS _OFFSET, BYVAL e AS _OFFSET, BYVAL rndm AS LONG)
63.     ''FUNCTION mpfr_set_uj_2exp& (rop AS mpfr_struct, BYVAL op AS _UNSIGNED _INTEGER64, BYVAL e AS _INTEGER64, BYVAL rndm AS LONG)
64.     ''FUNCTION mpfr_set_sj_2exp& (rop AS mpfr_struct, BYVAL op AS _INTEGER64, BYVAL e AS _INTEGER64, BYVAL rndm AS LONG)
65.     FUNCTION mpfr_set_z_2exp& (rop AS mpfr_struct, op AS mpz_struct, BYVAL e AS _OFFSET, BYVAL rndm AS LONG)
66.     FUNCTION mpfr_set_str& (rop AS mpfr_struct, s AS STRING, BYVAL base AS LONG, BYVAL rndm AS LONG)
67.     FUNCTION mpfr_strtofr& (rop AS mpfr_struct, nptr AS STRING, endptr AS _OFFSET, BYVAL base AS LONG, BYVAL rndm AS LONG)
68.     SUB mpfr_set_nan (x AS mpfr_struct)
69.     SUB mpfr_set_inf (x AS mpfr_struct, BYVAL sign AS LONG)
70.     SUB mpfr_set_zero (x AS mpfr_struct, BYVAL sign AS LONG)
71.     SUB mpfr_swap (x AS mpfr_struct, y AS mpfr_struct)
72.     FUNCTION mpfr_init_set_str& (x AS mpfr_struct, s AS STRING, BYVAL base AS LONG, BYVAL rndm AS LONG)
73.     FUNCTION mpfr_get_flt! (op AS mpfr_struct, BYVAL rndm AS LONG)
74.     FUNCTION mpfr_get_d# (op AS mpfr_struct, BYVAL rndm AS LONG)
75.     FUNCTION mpfr_get_ld## (op AS mpfr_struct, BYVAL rndm AS LONG)
76.     'Function mpfr_get_float128@__float128 (op AS mpfr_struct, Byval rndm as long )
77.     'Function mpfr_get_decimal64@_Decimal64 (op AS mpfr_struct, Byval rndm as long )
78.     FUNCTION mpfr_get_si%& (op AS mpfr_struct, BYVAL rndm AS LONG)
79.     FUNCTION mpfr_get_ui~%& (op AS mpfr_struct, BYVAL rndm AS LONG)
80.     ''FUNCTION mpfr_get_sj&& (op AS mpfr_struct, BYVAL rndm AS LONG)
81.     ''FUNCTION mpfr_get_uj~&& (op AS mpfr_struct, BYVAL rndm AS LONG)
82.     FUNCTION mpfr_get_d_2exp# (exp1 AS _OFFSET, op AS mpfr_struct, BYVAL rndm AS LONG)
83.     FUNCTION mpfr_get_ld_2exp## (exp1 AS _OFFSET, op AS mpfr_struct, BYVAL rndm AS LONG)
84.     FUNCTION mpfr_frexp& (exp1 AS _OFFSET, y AS mpfr_struct, x AS mpfr_struct, BYVAL rndm AS LONG)
85.     FUNCTION mpfr_get_z_2exp%& (rop AS mpz_struct, op AS mpfr_struct)
86.     FUNCTION mpfr_get_z& (rop AS mpz_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
87.     SUB mpfr_get_q (rop AS mpq_struct, op AS mpfr_struct)
88.     FUNCTION mpfr_get_f& (rop AS mpf_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
89.     FUNCTION mpfr_get_str$ (str1 AS STRING, expptr AS _OFFSET, BYVAL b AS LONG, BYVAL n AS _OFFSET, op AS mpfr_struct, BYVAL rndm AS LONG)
90.     '' SUB mpfr_free_str (str1 AS _OFFSET,) '' not sure about these, mpfr_get_str and mpfr_free_str
91.     FUNCTION mpfr_fits_ulong_p& (op AS mpfr_struct, BYVAL rndm AS LONG)
92.     FUNCTION mpfr_fits_slong_p& (op AS mpfr_struct, BYVAL rndm AS LONG)
93.     FUNCTION mpfr_fits_uint_p& (op AS mpfr_struct, BYVAL rndm AS LONG)
94.     FUNCTION mpfr_fits_sint_p& (op AS mpfr_struct, BYVAL rndm AS LONG)
95.     FUNCTION mpfr_fits_ushort_p& (op AS mpfr_struct, BYVAL rndm AS LONG)
96.     FUNCTION mpfr_fits_sshort_p& (op AS mpfr_struct, BYVAL rndm AS LONG)
97.     FUNCTION mpfr_fits_uintmax_p& (op AS mpfr_struct, BYVAL rndm AS LONG)
98.     FUNCTION mpfr_fits_intmax_p& (op AS mpfr_struct, BYVAL rndm AS LONG)
99.     FUNCTION mpfr_add& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
100.     FUNCTION mpfr_add_ui& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
101.     FUNCTION mpfr_add_si& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _OFFSET, BYVAL rndm AS LONG)
102.     FUNCTION mpfr_add_d& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS DOUBLE, BYVAL rndm AS LONG)
103.     FUNCTION mpfr_add_z& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpz_struct, BYVAL rndm AS LONG)
104.     FUNCTION mpfr_add_q& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpq_struct, BYVAL rndm AS LONG)
105.     FUNCTION mpfr_sub& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
106.     FUNCTION mpfr_ui_sub& (rop AS mpfr_struct, BYVAL op1 AS _UNSIGNED _OFFSET, op2 AS mpfr_struct, BYVAL rndm AS LONG)
107.     FUNCTION mpfr_sub_ui& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
108.     FUNCTION mpfr_si_sub& (rop AS mpfr_struct, BYVAL op1 AS _OFFSET, op2 AS mpfr_struct, BYVAL rndm AS LONG)
109.     FUNCTION mpfr_sub_si& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _OFFSET, BYVAL rndm AS LONG)
110.     FUNCTION mpfr_d_sub& (rop AS mpfr_struct, BYVAL op1 AS DOUBLE, op2 AS mpfr_struct, BYVAL rndm AS LONG)
111.     FUNCTION mpfr_sub_d& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS DOUBLE, BYVAL rndm AS LONG)
112.     FUNCTION mpfr_z_sub& (rop AS mpfr_struct, op1 AS mpz_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
113.     FUNCTION mpfr_sub_z& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpz_struct, BYVAL rndm AS LONG)
114.     FUNCTION mpfr_sub_q& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpq_struct, BYVAL rndm AS LONG)
115.     FUNCTION mpfr_mul& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
116.     FUNCTION mpfr_mul_ui& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
117.     FUNCTION mpfr_mul_si& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _OFFSET, BYVAL rndm AS LONG)
118.     FUNCTION mpfr_mul_d& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS DOUBLE, BYVAL rndm AS LONG)
119.     FUNCTION mpfr_mul_z& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpz_struct, BYVAL rndm AS LONG)
120.     FUNCTION mpfr_mul_q& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpq_struct, BYVAL rndm AS LONG)
121.     FUNCTION mpfr_sqr& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
122.     FUNCTION mpfr_div& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
123.     FUNCTION mpfr_ui_div& (rop AS mpfr_struct, BYVAL op1 AS _UNSIGNED _OFFSET, op2 AS mpfr_struct, BYVAL rndm AS LONG)
124.     FUNCTION mpfr_div_ui& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
125.     FUNCTION mpfr_si_div& (rop AS mpfr_struct, BYVAL op1 AS _OFFSET, op2 AS mpfr_struct, BYVAL rndm AS LONG)
126.     FUNCTION mpfr_div_si& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _OFFSET, BYVAL rndm AS LONG)
127.     FUNCTION mpfr_d_div& (rop AS mpfr_struct, BYVAL op1 AS DOUBLE, op2 AS mpfr_struct, BYVAL rndm AS LONG)
128.     FUNCTION mpfr_div_d& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS DOUBLE, BYVAL rndm AS LONG)
129.     FUNCTION mpfr_div_z& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpz_struct, BYVAL rndm AS LONG)
130.     FUNCTION mpfr_div_q& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpq_struct, BYVAL rndm AS LONG)
131.     FUNCTION mpfr_sqrt& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
132.     FUNCTION mpfr_sqrt_ui& (rop AS mpfr_struct, BYVAL op AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
133.     FUNCTION mpfr_rec_sqrt& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
134.     FUNCTION mpfr_cbrt& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
135.     FUNCTION mpfr_rootn_ui& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL k AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
136.     FUNCTION mpfr_root& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL k AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
137.     FUNCTION mpfr_pow& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
138.     FUNCTION mpfr_pow_ui& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
139.     FUNCTION mpfr_pow_si& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _OFFSET, BYVAL rndm AS LONG)
140.     FUNCTION mpfr_pow_z& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpz_struct, BYVAL rndm AS LONG)
141.     FUNCTION mpfr_ui_pow_ui& (rop AS mpfr_struct, BYVAL op1 AS _UNSIGNED _OFFSET, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
142.     FUNCTION mpfr_ui_pow& (rop AS mpfr_struct, BYVAL op1 AS _UNSIGNED _OFFSET, op2 AS mpfr_struct, BYVAL rndm AS LONG)
143.     FUNCTION mpfr_neg& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
144.     FUNCTION mpfr_abs& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
145.     FUNCTION mpfr_dim& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
146.     FUNCTION mpfr_mul_2ui& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
147.     FUNCTION mpfr_mul_2si& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _OFFSET, BYVAL rndm AS LONG)
148.     FUNCTION mpfr_div_2ui& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
149.     FUNCTION mpfr_div_2si& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _OFFSET, BYVAL rndm AS LONG)
150.     FUNCTION mpfr_cmp& (op1 AS mpfr_struct, op2 AS mpfr_struct)
151.     FUNCTION mpfr_cmp_ui& (op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET)
152.     FUNCTION mpfr_cmp_si& (op1 AS mpfr_struct, BYVAL op2 AS _OFFSET)
153.     FUNCTION mpfr_cmp_d& (op1 AS mpfr_struct, BYVAL op2 AS DOUBLE)
154.     FUNCTION mpfr_cmp_ld& (op1 AS mpfr_struct, BYVAL p2 AS _FLOAT)
155.     FUNCTION mpfr_cmp_z& (op1 AS mpfr_struct, op2 AS mpz_struct)
156.     FUNCTION mpfr_cmp_q& (op1 AS mpfr_struct, op2 AS mpq_struct)
157.     FUNCTION mpfr_cmp_f& (op1 AS mpfr_struct, op2 AS mpf_struct)
158.     FUNCTION mpfr_cmp_ui_2exp& (op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL e AS _OFFSET)
159.     FUNCTION mpfr_cmp_si_2exp& (op1 AS mpfr_struct, BYVAL op2 AS _OFFSET, BYVAL e AS _OFFSET)
160.     FUNCTION mpfr_cmpabs& (op1 AS mpfr_struct, op2 AS mpfr_struct)
161.     FUNCTION mpfr_nan_p& (op AS mpfr_struct)
162.     FUNCTION mpfr_inf_p& (op AS mpfr_struct)
163.     FUNCTION mpfr_number_p& (op AS mpfr_struct)
164.     FUNCTION mpfr_zero_p& (op AS mpfr_struct)
165.     FUNCTION mpfr_regular_p& (op AS mpfr_struct)
166.     FUNCTION mpfr_greater_p& (op1 AS mpfr_struct, op2 AS mpfr_struct)
167.     FUNCTION mpfr_greaterequal_p& (op1 AS mpfr_struct, op2 AS mpfr_struct)
168.     FUNCTION mpfr_less_p& (op1 AS mpfr_struct, op2 AS mpfr_struct)
169.     FUNCTION mpfr_lessequal_p& (op1 AS mpfr_struct, op2 AS mpfr_struct)
170.     FUNCTION mpfr_equal_p& (op1 AS mpfr_struct, op2 AS mpfr_struct)
171.     FUNCTION mpfr_lessgreater_p& (op1 AS mpfr_struct, op2 AS mpfr_struct)
172.     FUNCTION mpfr_unordered_p& (op1 AS mpfr_struct, op2 AS mpfr_struct)
173.     FUNCTION mpfr_log& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
174.     FUNCTION mpfr_log_ui& (rop AS mpfr_struct, BYVAL op AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
175.     FUNCTION mpfr_log2& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
176.     FUNCTION mpfr_log10& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
177.     FUNCTION mpfr_log1p& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
178.     FUNCTION mpfr_exp& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
179.     FUNCTION mpfr_exp2& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
180.     FUNCTION mpfr_exp10& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
181.     FUNCTION mpfr_expm1& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
182.     FUNCTION mpfr_cos& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
183.     FUNCTION mpfr_sin& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
184.     FUNCTION mpfr_tan& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
185.     FUNCTION mpfr_sin_cos& (sop AS mpfr_struct, cop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
186.     FUNCTION mpfr_sec& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
187.     FUNCTION mpfr_csc& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
188.     FUNCTION mpfr_cot& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
189.     FUNCTION mpfr_acos& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
190.     FUNCTION mpfr_asin& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
191.     FUNCTION mpfr_atan& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
192.     FUNCTION mpfr_atan2& (rop AS mpfr_struct, y AS mpfr_struct, x AS mpfr_struct, BYVAL rndm AS LONG)
193.     FUNCTION mpfr_cosh& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
194.     FUNCTION mpfr_sinh& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
195.     FUNCTION mpfr_tanh& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
196.     FUNCTION mpfr_sinh_cosh& (sop AS mpfr_struct, cop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
197.     FUNCTION mpfr_sech& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
198.     FUNCTION mpfr_csch& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
199.     FUNCTION mpfr_coth& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
200.     FUNCTION mpfr_acosh& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
201.     FUNCTION mpfr_asinh& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
202.     FUNCTION mpfr_atanh& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
203.     FUNCTION mpfr_fac_ui& (rop AS mpfr_struct, BYVAL op AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
204.     FUNCTION mpfr_eint& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
205.     FUNCTION mpfr_li2& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
206.     FUNCTION mpfr_gamma& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
207.     FUNCTION mpfr_gamma_inc& (rop AS mpfr_struct, op AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
208.     FUNCTION mpfr_lngamma& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
209.     FUNCTION mpfr_lgamma& (rop AS mpfr_struct, signp AS LONG, op AS mpfr_struct, BYVAL rndm AS LONG)
210.     FUNCTION mpfr_digamma& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
211.     FUNCTION mpfr_beta& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
212.     FUNCTION mpfr_zeta& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
213.     FUNCTION mpfr_zeta_ui& (rop AS mpfr_struct, BYVAL op AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
214.     FUNCTION mpfr_erf& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
215.     FUNCTION mpfr_erfc& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
216.     FUNCTION mpfr_j0& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
217.     FUNCTION mpfr_j1& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
218.     FUNCTION mpfr_jn& (rop AS mpfr_struct, BYVAL n AS _OFFSET, op AS mpfr_struct, BYVAL rndm AS LONG)
219.     FUNCTION mpfr_y0& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
220.     FUNCTION mpfr_y1& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
221.     FUNCTION mpfr_yn& (rop AS mpfr_struct, BYVAL n AS _OFFSET, op AS mpfr_struct, BYVAL rndm AS LONG)
222.     FUNCTION mpfr_fma& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, op3 AS mpfr_struct, BYVAL rndm AS LONG)
223.     FUNCTION mpfr_fms& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, op3 AS mpfr_struct, BYVAL rndm AS LONG)
224.     FUNCTION mpfr_fmma& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, op3 AS mpfr_struct, op4 AS mpfr_struct, BYVAL rndm AS LONG)
225.     FUNCTION mpfr_fmms& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, op3 AS mpfr_struct, op4 AS mpfr_struct, BYVAL rndm AS LONG)
226.     FUNCTION mpfr_agm& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
227.     FUNCTION mpfr_hypot& (rop AS mpfr_struct, x AS mpfr_struct, y AS mpfr_struct, BYVAL rndm AS LONG)
228.     FUNCTION mpfr_ai& (rop AS mpfr_struct, x AS mpfr_struct, BYVAL rndm AS LONG)
229.     FUNCTION mpfr_const_log2& (rop AS mpfr_struct, BYVAL rndm AS LONG)
230.     FUNCTION mpfr_const_pi& (rop AS mpfr_struct, BYVAL rndm AS LONG)
231.     FUNCTION mpfr_const_euler& (rop AS mpfr_struct, BYVAL rndm AS LONG)
232.     FUNCTION mpfr_const_catalan& (rop AS mpfr_struct, BYVAL rndm AS LONG)
233.     SUB mpfr_free_cache
234.     'Sub mpfr_free_cache2 (mpfr_free_cache_t way)
235.     SUB mpfr_free_pool
236.     FUNCTION mpfr_mp_memory_cleanup&
237.     ''Function mpfr_sum& (rop AS mpfr_struct, const mpfr_ptr tab[], byval n as _unsigned _offset, Byval rndm as long )
238.     ''Function mpfr_out_str%&(stream as _offset, byval base as long, byval n as _offset, op AS mpfr_struct, Byval rndm as long )
239.     ''Function mpfr_inp_str%&(rop AS mpfr_struct, stream as _offset, byval base as long, Byval rndm as long )
240.     ''Function mpfr_fpif_export& (stream as _offset, op AS mpfr_struct )
241.     ''Function mpfr_fpif_import& (op AS mpfr_struct, stream as _offset )
242.     SUB mpfr_dump (op AS mpfr_struct)
243.     ''Function mpfr_fprintf& (FILE *stream, const char *template, … )
244.     ''Function mpfr_vfprintf& (FILE *stream, const char *template, va_list ap )
245.     ''Function mpfr_printf& (const char *template, … )
246.     ''Function mpfr_vprintf& (const char *template, va_list ap )
247.     ''Function mpfr_sprintf& (char *buf, const char *template, … )
248.     ''Function mpfr_vsprintf& (char *buf, const char *template, va_list ap )
249.     ''Function mpfr_snprintf& (char *buf, size_t n, const char *template, … )
250.     ''Function mpfr_vsnprintf& (char *buf, size_t n, const char *template, va_list ap )
251.     ''Function mpfr_asprintf& (char **str, const char *template, … )
252.     ''Function mpfr_vasprintf& (char **str, const char *template, va_list ap )
253.     FUNCTION mpfr_rint& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
254.     FUNCTION mpfr_ceil& (rop AS mpfr_struct, op AS mpfr_struct)
255.     FUNCTION mpfr_floor& (rop AS mpfr_struct, op AS mpfr_struct)
256.     FUNCTION mpfr_round& (rop AS mpfr_struct, op AS mpfr_struct)
257.     FUNCTION mpfr_roundeven& (rop AS mpfr_struct, op AS mpfr_struct)
258.     FUNCTION mpfr_trunc& (rop AS mpfr_struct, op AS mpfr_struct)
259.     FUNCTION mpfr_rint_ceil& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
260.     FUNCTION mpfr_rint_floor& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
261.     FUNCTION mpfr_rint_round& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
262.     FUNCTION mpfr_rint_roundeven& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
263.     FUNCTION mpfr_rint_trunc& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
264.     FUNCTION mpfr_frac& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
265.     FUNCTION mpfr_modf& (iop AS mpfr_struct, fop AS mpfr_struct, op AS mpfr_struct, BYVAL rndm AS LONG)
266.     FUNCTION mpfr_fmod& (r AS mpfr_struct, x AS mpfr_struct, y AS mpfr_struct, BYVAL rndm AS LONG)
267.     FUNCTION mpfr_fmodquo& (r AS mpfr_struct, q AS _OFFSET, x AS mpfr_struct, y AS mpfr_struct, BYVAL rndm AS LONG)
268.     FUNCTION mpfr_remainder& (r AS mpfr_struct, x AS mpfr_struct, y AS mpfr_struct, BYVAL rndm AS LONG)
269.     FUNCTION mpfr_remquo& (r AS mpfr_struct, q AS _OFFSET, x AS mpfr_struct, y AS mpfr_struct, BYVAL rndm AS LONG)
270.     FUNCTION mpfr_integer_p& (op AS mpfr_struct)
271.     SUB mpfr_set_default_rounding_mode (BYVAL rndm AS LONG)
272.     FUNCTION mpfr_get_default_rounding_mode&
273.     FUNCTION mpfr_prec_round& (x AS mpfr_struct, BYVAL prec AS _OFFSET, BYVAL rndm AS LONG)
274.     FUNCTION mpfr_can_round& (b AS mpfr_struct, BYVAL err1 AS _OFFSET, BYVAL rndm1 AS LONG, BYVAL rndm2 AS LONG, BYVAL prec AS _OFFSET)
275.     FUNCTION mpfr_min_prec%& (x AS mpfr_struct)
276.     FUNCTION mpfr_print_rnd_mode$ (BYVAL rndm AS LONG)
277.     SUB mpfr_nexttoward (x AS mpfr_struct, y AS mpfr_struct)
278.     SUB mpfr_nextabove (x AS mpfr_struct)
279.     SUB mpfr_nextbelow (x AS mpfr_struct)
280.     FUNCTION mpfr_min& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
281.     FUNCTION mpfr_max& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
282.     ''Function mpfr_urandomb& (rop AS mpfr_struct, state as gmp_randstate_struct )
283.     ''Function mpfr_urandom& (rop AS mpfr_struct, state as gmp_randstate_struct, Byval rndm as long )
284.     ''Function mpfr_nrandom& (rop AS mpfr_struct1, state as gmp_randstate_struct, Byval rndm as long )
285.     ''Function mpfr_grandom& (rop AS mpfr_struct1, rop2 AS mpfr_struct, state as gmp_randstate_struct, Byval rndm as long )
286.     ''Function mpfr_erandom& (rop AS mpfr_struct1, state as gmp_randstate_struct, Byval rndm as long )
287.     FUNCTION mpfr_get_exp%& (x AS mpfr_struct)
288.     FUNCTION mpfr_set_exp& (x AS mpfr_struct, BYVAL e AS _OFFSET)
289.     FUNCTION mpfr_signbit& (op AS mpfr_struct)
290.     FUNCTION mpfr_setsign& (rop AS mpfr_struct, op AS mpfr_struct, BYVAL s AS LONG, BYVAL rndm AS LONG)
291.     FUNCTION mpfr_copysign& (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
292.     FUNCTION mpfr_get_version$
293.     FUNCTION mpfr_get_patches$
294.     FUNCTION mpfr_buildopt_tls_p&
295.     FUNCTION mpfr_buildopt_float128_p&
296.     FUNCTION mpfr_buildopt_decimal_p&
297.     FUNCTION mpfr_buildopt_gmpinternals_p&
298.     FUNCTION mpfr_buildopt_sharedcache_p&
299.     FUNCTION mpfr_buildopt_tune_case$
300.     FUNCTION mpfr_get_emin%&
301.     FUNCTION mpfr_get_emax%&
302.     FUNCTION mpfr_set_emin& (BYVAL exp1 AS _OFFSET)
303.     FUNCTION mpfr_set_emax& (BYVAL exp1 AS _OFFSET)
304.     FUNCTION mpfr_get_emin_min%&
305.     FUNCTION mpfr_get_emin_max%&
306.     FUNCTION mpfr_get_emax_min%&
307.     FUNCTION mpfr_get_emax_max%&
308.     FUNCTION mpfr_check_range& (x AS mpfr_struct, BYVAL t AS LONG, BYVAL rndm AS LONG)
309.     FUNCTION mpfr_subnormalize& (x AS mpfr_struct, BYVAL t AS LONG, BYVAL rndm AS LONG)
310.     SUB mpfr_clear_underflow
311.     SUB mpfr_clear_overflow
312.     SUB mpfr_clear_divby0
313.     SUB mpfr_clear_nanflag
314.     SUB mpfr_clear_inexflag
315.     SUB mpfr_clear_erangeflag
316.     SUB mpfr_clear_flags
317.     SUB mpfr_set_underflow
318.     SUB mpfr_set_overflow
319.     SUB mpfr_set_divby0
320.     SUB mpfr_set_nanflag
321.     SUB mpfr_set_inexflag
322.     SUB mpfr_set_erangeflag
323.     FUNCTION mpfr_underflow_p&
324.     FUNCTION mpfr_overflow_p&
325.     FUNCTION mpfr_divby0_p&
326.     FUNCTION mpfr_nanflag_p&
327.     FUNCTION mpfr_inexflag_p&
328.     FUNCTION mpfr_erangeflag_p&
329.     SUB mpfr_flags_clear (BYVAL mask AS LONG)
330.     SUB mpfr_flags_set (BYVAL mask AS LONG)
331.     FUNCTION mpfr_flags_test& (BYVAL mask AS LONG)
332.     FUNCTION mpfr_flags_save&
333.     SUB mpfr_flags_restore (BYVAL flags AS LONG, BYVAL mask AS LONG)
334.     SUB mpfr_set_prec_raw (x AS mpfr_struct, BYVAL prec AS _OFFSET)
335.     FUNCTION mpfr_eq& (op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL op3 AS _UNSIGNED _OFFSET)
336.     SUB mpfr_reldiff (rop AS mpfr_struct, op1 AS mpfr_struct, op2 AS mpfr_struct, BYVAL rndm AS LONG)
337.     FUNCTION mpfr_mul_2exp& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
338.     FUNCTION mpfr_div_2exp& (rop AS mpfr_struct, op1 AS mpfr_struct, BYVAL op2 AS _UNSIGNED _OFFSET, BYVAL rndm AS LONG)
339.     FUNCTION mpfr_custom_get_size%& (BYVAL prec AS _OFFSET)
340.     SUB mpfr_custom_init (significand AS _OFFSET, BYVAL prec AS _OFFSET)
341.     SUB mpfr_custom_init_set (x AS mpfr_struct, BYVAL kind AS LONG, BYVAL exp1 AS _OFFSET, BYVAL prec AS _OFFSET, significand AS _OFFSET)
342.     FUNCTION mpfr_custom_get_kind& (x AS mpfr_struct)
343.     FUNCTION mpfr_custom_get_significand%& (x AS mpfr_struct)
344.     FUNCTION mpfr_custom_get_exp%& (x AS mpfr_struct)
345.     SUB mpfr_custom_move (x AS mpfr_struct, new_position AS _OFFSET)
346.     FUNCTION mpfr_snprintf& (rop AS STRING, BYVAL n AS _UNSIGNED _OFFSET, frmt AS STRING, src2 AS mpfr_struct)
347. END DECLARE
348.  
349. SUB print_mpfr_g (f AS mpfr_struct, n AS LONG, m AS LONG)
350.     DIM sresult(0) AS STRING, frmt AS STRING
351.     DIM ln AS LONG
352.     frmt = "%" + LTRIM$(STR$(n)) + "." + LTRIM$(STR$(m)) + "RNg"
353.     ln = mpfr_snprintf&("", 0, frmt, f)
354.     IF ln > 0 THEN
355.         sresult(0) = SPACE$(ln)
356.         ln = mpfr_snprintf&(sresult(0), ln + 4, frmt, f)
357.     END IF
358.     PRINT LEFT$(sresult(0), ln)
359. END SUB
360.  
361. SUB print_mpfr_f (f AS mpfr_struct, n AS LONG, m AS LONG)
362.     DIM sresult(0) AS STRING, frmt AS STRING
363.     DIM ln AS LONG
364.     frmt = "%" + LTRIM$(STR$(n)) + "." + LTRIM$(STR$(m)) + "RNf"
365.     ln = mpfr_snprintf&("", 0, frmt, f)
366.     IF ln > 0 THEN
367.         sresult(0) = SPACE$(ln)
368.         ln = mpfr_snprintf&(sresult(0), ln + 4, frmt, f)
369.     END IF
370.     PRINT LEFT$(sresult(0), ln)
371. END SUB
372.  
373. SUB print_mpfr_e (f AS mpfr_struct, n AS LONG, m AS LONG)
374.     DIM sresult(0) AS STRING, frmt AS STRING
375.     DIM ln AS LONG
376.     frmt = "%" + LTRIM$(STR$(n)) + "." + LTRIM$(STR$(m)) + "RNe"
377.     ln = mpfr_snprintf&("", 0, frmt, f)
378.     IF ln > 0 THEN
379.         sresult(0) = SPACE$(ln)
380.         ln = mpfr_snprintf&(sresult(0), ln + 4, frmt, f)
381.     END IF
382.     PRINT LEFT$(sresult(0), ln)
383. END SUB
384.  
385. SUB valmpfr (result AS mpfr_struct, value AS STRING, nbase AS LONG)
386.     DIM c AS LONG
387.     c = mpfr_set_str(result, value, nbase, MPFR_RNDN)
388. END SUB
389.  

[_vince]

I can confirm that removing "lib" makes DECLARE LIBRARY work on linux.  Thanks, Jack

[jack]

@_vince
I am looking forward to your fractal program :)

[jack]

updated first post and also attached the 64-bit mpfr dll
the dll can be in the same place where your program resides or you could put the dll in system32 although some people cringe at the thought of putting third party dll's in system32 so it's up to you where you put it
just in case you didn't know, in Windows x64, system32 is where 64-bit applications/dll's reside and  SysWOW64 is where 32-bit applications/dll's reside

[jack]

the first and last 20 digits of Fibonacci ( 18446744073709551616 ) https://rosettacode.org/wiki/Fibonacci_matrix-exponentiation#Fibmod_method in 0 seconds
get the mpfr dll from opening post

Code: QB64: [Select]

1. 'translated the fib code from https://www.uni-math.gwdg.de/tschinkel/gauss/Fibon.pdf page 14
2. '------------------------------------------------------------------'/
3. '
4. '** Schnelle Berechnung der Fibonacci-Zahlen mittels der Formeln
5. '** fib(2*k-1) = fib(k)**2 + fib(k-1)**2
6. '** fib(2*k) = fib(k)**2 + 2*fib(k)*fib(k-1)
7. '**
8. '** Dabei werden alle Berechnungen mod m durchgef¨uhrt
9. '
10. $Console:Only
11. _Dest _Console
12.  
13. Const MPFR_RNDN = 0
14. Const MPFR_RNDZ = 1
15.  
16. Type mpz_struct
17.     mp_alloc As Long
18.     mp_size As Long
19.     mp_d As _Offset
20. End Type
21.  
22. Type mpfr_struct
23.     mpfr_prec As _Offset
24.     mpfr_sign As Long
25.     mpfr_exp As _Offset
26.     mpfr_d As _Offset
27. End Type
28.  
29. Declare Dynamic Library "libmpfr-6"
30.     Sub mpfr_init2 (rop As mpfr_struct, Byval op1 As Long)
31.     Sub mpfr_clear (rop As mpfr_struct)
32.     Sub mpfr_set_z (rop As mpfr_struct, opz As mpz_struct, Byval mpfr_rnd As Long)
33.     Sub mpfr_set_ui (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
34.     Sub mpfr_ui_pow_ui (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
35.     Sub mpfr_sqrt_ui (rop As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
36.     Sub mpfr_add_ui (rop As mpfr_struct, op1 As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
37.     Sub mpfr_div_ui (rop As mpfr_struct, op1 As mpfr_struct, Byval culong As _Unsigned Long, Byval mpfr_rnd As Long)
38.     Sub mpfr_log10 (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
39.     Sub mpfr_log2 (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
40.     Sub mpfr_exp10 (rop As mpfr_struct, op1 As mpfr_struct, Byval mpfr_rnd As Long)
41.     Sub mpfr_trunc (rop As mpfr_struct, op As mpfr_struct)
42.     Sub mpfr_sub (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
43.     Sub mpfr_mul (rop As mpfr_struct, op1 As mpfr_struct, op2 As mpfr_struct, Byval mpfr_rnd As Long)
44.     Sub mpfr_frac (rop As mpfr_struct, op As mpfr_struct, Byval mpfr_rnd As Long)
45.     Function mpfr_snprintf& (rop As String, Byval n As _Unsigned _Offset, frmt As String, src2 As mpfr_struct)
46.     Function mpfr_get_ui& (mpfr_srcptr As mpfr_struct, Byval mpfr_rnd As Long)
47.     Sub mpfr_get_z (z As mpz_struct, f As mpfr_struct, Byval mpfr_rnd As Long)
48. End Declare
49.  
50. Declare Dynamic Library "libgmp-10"
51.     Sub __gmpz_init (result As mpz_struct)
52.     Sub __gmpz_init_set_ui (rop As mpz_struct, Byval op As _Unsigned Long)
53.     Sub __gmpz_set (rop As mpz_struct, op As mpz_struct)
54.     Function __gmpz_cmp_ui& (mpz_srcptr As mpz_struct, Byval ui As _Unsigned Long)
55.     Sub __gmpz_mul (result As mpz_struct, src1 As mpz_struct, src2 As mpz_struct)
56.     Sub __gmpz_mod (mpz_ptr As mpz_struct, mpz_srcptr As mpz_struct, mpz_srcptr2 As mpz_struct)
57.     Sub __gmpz_mul_ui (result As mpz_struct, src As mpz_struct, Byval ui As _Unsigned Long)
58.     Sub __gmpz_add (result As mpz_struct, src1 As mpz_struct, src2 As mpz_struct)
59.     Function __gmpz_tstbit& (mpz_srcptr As mpz_struct, Byval ui As _Unsigned Long)
60.     Sub __gmpz_clear (result As mpz_struct)
61.     Sub __gmpz_set_ui (mpz_ptr As mpz_struct, Byval ui As _Unsigned Long)
62.     Sub __gmpz_ui_pow_ui (mpz_ptr As mpz_struct, Byval ui1 As _Unsigned Long, Byval ui2 As _Unsigned Long)
63.     Function __gmp_snprintf& (rop As String, Byval n As _Unsigned _Integer64, frmt As String, src2 As mpz_struct)
64. End Declare
65.  
66. 'the first digits are approximated by
67. '10^Frac((log10((1+sqrt(5))/2)*n)-log10(sqrt(5)))
68. 'the number of digits in the Fibonacci number
69. 'is Trunc((log10((1+sqrt(5))/2)*n)-log10(sqrt(5)))+1
70. Dim As Double t
71. Dim As Long c, length, fnp, number
72. length = 21
73. fnp = 64
74. number = 10000000 '4784969
75. t = Timer
76. Dim As String s
77. Dim As mpfr_struct first, five, tmp, nf
78. Dim As mpz_struct flast, n, m
79.  
80. __gmpz_init flast
81. __gmpz_init n
82. __gmpz_init m
83. __gmpz_ui_pow_ui n, 2, fnp
84.  
85. mpfr_init2 first, 256 + fnp
86. mpfr_init2 five, 256 + fnp
87. mpfr_init2 tmp, 256 + fnp
88. mpfr_init2 nf, 256 + fnp
89.  
90. mpfr_ui_pow_ui nf, 2, fnp, MPFR_RNDN
91. 'mpfr_set_si(nf, number, MPFR_RNDN)
92. mpfr_sqrt_ui five, 5, MPFR_RNDN
93. mpfr_add_ui tmp, five, 1, MPFR_RNDN
94. mpfr_div_ui first, tmp, 2, MPFR_RNDN
95. mpfr_log10 tmp, first, MPFR_RNDN
96. mpfr_mul first, tmp, nf, MPFR_RNDN
97. mpfr_log10 tmp, five, MPFR_RNDN
98. mpfr_sub nf, first, tmp, MPFR_RNDN
99. mpfr_frac tmp, nf, MPFR_RNDN
100. mpfr_exp10 first, tmp, MPFR_RNDN
101. mpfr_trunc tmp, nf
102. mpfr_add_ui tmp, tmp, 1, MPFR_RNDN
103. Print "fibonacci "; gmpz$(n)
104. s = print_mpfr_f$(first, length + 4, length - 2)
105.  
106. c = InStr(s, ".")
107. s = Left$(s, c - 1) + Mid$(s, c + 1)
108. c = mpfr_get_ui(tmp, MPFR_RNDN)
109. If c > 0 And c < length Then
110.     Print "the first"; c; " digits are";
111.     Print Left$(s, c)
112. Else
113.     Print "the first"; length - 1; " digits are ";
114.     Print s
115. End If
116.  
117. mpfr_clear nf
118. mpfr_clear five
119. mpfr_clear first
120.  
121. '__gmpz_set_si(n, number)
122. __gmpz_ui_pow_ui m, 10, length
123. fib flast, n, m
124. t = Timer - t
125. If c > 0 And c < length Then
126.     Print "the last"; c; " digits are ";
127. Else
128.     Print "the last"; length - 1; " digits are ";
129. End If
130. Print gmpz$(flast)
131. mpfr_get_z n, tmp, MPFR_RNDN
132. Print "the Fibonacci number has ";
133. Print gmpz$(n);
134. Print " digits"
135. __gmpz_clear m
136. __gmpz_clear n
137. __gmpz_clear flast
138. mpfr_clear tmp
139. Print "elapsed time "; t; " seconds"
140.  
141. Function print_mpfr_f$ (f As mpfr_struct, n As Long, m As Long)
142.     Dim sresult(0) As String, frmt As String
143.     Erase sresult
144.     Dim sresult(0) As String
145.     Dim ln As Long
146.     frmt = "%" + LTrim$(Str$(n)) + "." + LTrim$(Str$(m)) + "RZf" + Chr$(0)
147.     ln = mpfr_snprintf&("" + Chr$(0), 0, frmt, f)
148.     If ln > 0 Then
149.         sresult(0) = Space$(ln) + Chr$(0)
150.         ln = mpfr_snprintf&(sresult(0), ln + 4, frmt, f)
151.     End If
152.     print_mpfr_f$ = Left$(sresult(0), ln)
153. End Function
154.  
155. Function gmpz$ (n As mpz_struct)
156.     Dim sresult As String
157.     Dim ln As _Unsigned _Integer64
158.     sresult = ""
159.     ln = __gmp_snprintf&("" + Chr$(0), 0, "%Zd" + Chr$(0), n)
160.     If ln > 0 Then
161.         sresult = Space$(ln + 5)
162.         ln = __gmp_snprintf&(sresult, ln + 4, "%Zd" + Chr$(0), n)
163.     End If
164.     gmpz$ = sresult
165. End Function
166.  
167. Sub fib (x As mpz_struct, k As mpz_struct, m As mpz_struct)
168.     Dim As mpz_struct b, y, xx, tmp, temp
169.     Dim As mpfr_struct lngth
170.     Dim As Long i, bit_length
171.     mpfr_init2 lngth, 1024
172.     __gmpz_init b
173.     __gmpz_init y
174.     __gmpz_init xx
175.     __gmpz_init tmp
176.     __gmpz_init temp
177.  
178.     If __gmpz_cmp_ui(k, c) <= 0 Then
179.         __gmpz_set x, k
180.         mpfr_clear lngth
181.         __gmpz_clear temp
182.         __gmpz_clear tmp
183.         __gmpz_clear xx
184.         __gmpz_clear y
185.         __gmpz_clear b
186.         Exit Sub
187.     End If
188.     mpfr_set_z lngth, k, MPFR_RNDN
189.     mpfr_log2 lngth, lngth, MPFR_RNDN
190.     mpfr_trunc lngth, lngth
191.     bit_length = mpfr_get_ui(lngth, MPFR_RNDN) 'Fix(Log(__gmpz_get_d(k))/Log(2))
192.     __gmpz_set_ui x, 1
193.     __gmpz_set_ui y, 0
194.     For i = bit_length - 1 To 0 Step -1
195.         __gmpz_mul temp, x, x
196.         __gmpz_mod xx, temp, m
197.         __gmpz_mul temp, x, y
198.         __gmpz_mul_ui tmp, temp, 2
199.         __gmpz_add temp, xx, tmp
200.         __gmpz_mod x, temp, m
201.         'x := (xx + 2*x*y) mod m;
202.         __gmpz_mul tmp, y, y
203.         __gmpz_add temp, xx, tmp
204.         __gmpz_mod y, temp, m
205.         'y := (xx + y*y) mod m;
206.         If __gmpz_tstbit(k, i) Then
207.             __gmpz_set temp, x
208.             __gmpz_add tmp, x, y
209.             __gmpz_mod x, tmp, m
210.             'x := (x + y) mod m;
211.             __gmpz_set y, temp
212.         End If
213.     Next
214.     mpfr_clear lngth
215.     __gmpz_clear temp
216.     __gmpz_clear tmp
217.     __gmpz_clear xx
218.     __gmpz_clear y
219.     __gmpz_clear b
220. End Sub
221.