1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
|
#include <sollya.h>
#include <mpfr.h>
#include <mpfi.h>
#include <stdlib.h>
#include <string.h>
typedef struct data_struct_t {
char text[32];
int counter;
} data_t;
void euler_gamma(mpfr_t res, mp_prec_t prec, void *ptr) {
data_t *data;
data = (data_t *) ptr;
sollya_lib_printf(">>>>%s<<<<>>>>%d<<<<\n", data->text, data->counter);
(data->counter)++;
mpfr_set_prec(res, prec);
mpfr_const_euler(res, GMP_RNDN);
return;
}
void euler_gamma2(mpfr_t res, mp_prec_t prec, void *ptr) {
data_t *data;
data = (data_t *) ptr;
sollya_lib_printf(">>>>%s<<<<>>>>%d<<<<\n", data->text, data->counter);
(data->counter)++;
mpfr_set_prec(res, prec);
mpfr_const_euler(res, GMP_RNDN);
return;
}
void dealloc(void *ptr) {
data_t *data;
data = (data_t *) ptr;
sollya_lib_printf("Deallocation function called for the data pointer (%s<<<>>>%d)\n", data->text, data->counter);
return;
}
int main(void) {
sollya_obj_t f,g;
mpfr_t x,y;
data_t data = { "Hello world", 0 };
sollya_lib_init();
f = sollya_lib_build_function_libraryconstant_with_data("superconst", euler_gamma, &data, NULL);
sollya_lib_printf("%b (expecting: superconst)\n", f);
mpfr_init2(x, 30);
mpfr_init2(y, 50);
mpfr_set_ui(x, 2, GMP_RNDN);
sollya_lib_evaluate_function_at_point(y, f, x, NULL);
sollya_lib_printf("%v (expecting: 0.5772...)\n", y);
sollya_lib_clear_obj(f);
f = sollya_lib_parse_string("superconst");
sollya_lib_printf("%b (expecting: superconst)\n", f);
mpfr_set_ui(x, 2, GMP_RNDN);
sollya_lib_evaluate_function_at_point(y, f, x, NULL);
sollya_lib_printf("%v (expecting: 0.5772...)\n", y);
sollya_lib_clear_obj(f);
f = sollya_lib_libraryconstant_with_data("superconst2", euler_gamma2, &data, dealloc);
g = SOLLYA_ADD(SOLLYA_CONST(2), sollya_lib_copy_obj(f));
sollya_lib_printf("%b (expecting: superconst2)\n", f);
sollya_lib_evaluate_function_at_point(y, f, x, NULL);
sollya_lib_printf("%v (expecting: 0.5772...)\n", y);
sollya_lib_clear_obj(f);
sollya_lib_printf("Deallocation function must not yet have been called\n");
sollya_lib_clear_obj(g);
mpfr_clear(x);
mpfr_clear(y);
sollya_lib_close();
return 0;
}
|
