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
|
#include <sollya.h>
sollya_obj_t stupid_wrapper(sollya_obj_t arg1, sollya_obj_t arg2, sollya_obj_t arg3, sollya_obj_t arg4, ...) {
va_list va;
sollya_obj_t a;
va_start(va, arg4);
a = sollya_lib_v_annotatefunction(arg1, arg2, arg3, arg4, va);
va_end(va);
return a;
}
int main(void) {
sollya_obj_t f, g, h, dom, delta;
mpfr_t x, y;
sollya_lib_init();
f = sollya_lib_parse_string("exp(_x_)");
g = sollya_lib_parse_string("function(proc(X,n,p) { var res, oldPrec;\n"
" oldPrec = prec;\n"
" prec = p!;\n"
" \"Using procedure function exponential\";\n"
" res = exp(X);\n"
" prec = oldPrec!;\n"
" return res;\n"
" })");
dom = sollya_lib_parse_string("[-5;5]");
delta = sollya_lib_parse_string("[0]");
h = stupid_wrapper(f, g, dom, delta, NULL);
mpfr_init2(x, 24);
mpfr_init2(y, 24);
mpfr_set_si(x, 2, GMP_RNDN);
sollya_lib_evaluate_function_at_point(y, h, x, NULL);
mpfr_set_si(x, 6, GMP_RNDN);
sollya_lib_evaluate_function_at_point(y, h, x, NULL);
mpfr_clear(y);
mpfr_clear(x);
sollya_lib_clear_obj(f);
sollya_lib_clear_obj(g);
sollya_lib_clear_obj(dom);
sollya_lib_clear_obj(delta);
sollya_lib_clear_obj(h);
sollya_lib_close();
return 0;
}
|
