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
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
|
#include <sollya.h>
int callback(sollya_msg_t msg, void *data) {
(void)data; /* Avoiding "unused parameter" warning */
int message = sollya_lib_get_msg_id(msg);
if (message==SOLLYA_MSG_ROUNDING_OCCURRED_WHILE_READING_A_CONSTANT)
sollya_lib_printf("Caught the message about the fact that faithful rounding has occurred.\n");
else
sollya_lib_printf("Unexpected message caught.\n");
return 0;
}
int main(void) {
sollya_obj_t a,b;
sollya_lib_init();
sollya_lib_install_msg_callback(callback, NULL);
a = SOLLYA_SQRT(SOLLYA_X_); b = sollya_lib_build_function_sqrt(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_EXP(SOLLYA_X_); b = sollya_lib_build_function_exp(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_LOG(SOLLYA_X_); b = sollya_lib_build_function_log(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_LOG2(SOLLYA_X_); b = sollya_lib_build_function_log2(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_LOG10(SOLLYA_X_); b = sollya_lib_build_function_log10(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_SIN(SOLLYA_X_); b = sollya_lib_build_function_sin(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_COS(SOLLYA_X_); b = sollya_lib_build_function_cos(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_TAN(SOLLYA_X_); b = sollya_lib_build_function_tan(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ASIN(SOLLYA_X_); b = sollya_lib_build_function_asin(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ACOS(SOLLYA_X_); b = sollya_lib_build_function_acos(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ATAN(SOLLYA_X_); b = sollya_lib_build_function_atan(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_NEG(SOLLYA_X_); b = sollya_lib_build_function_neg(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ABS(SOLLYA_X_); b = sollya_lib_build_function_abs(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_D(SOLLYA_X_); b = sollya_lib_build_function_double(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_SG(SOLLYA_X_); b = sollya_lib_build_function_single(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_QD(SOLLYA_X_); b = sollya_lib_build_function_quad(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_HP(SOLLYA_X_); b = sollya_lib_build_function_halfprecision(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_DD(SOLLYA_X_); b = sollya_lib_build_function_double_double(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_TD(SOLLYA_X_); b = sollya_lib_build_function_triple_double(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ERF(SOLLYA_X_); b = sollya_lib_build_function_erf(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ERFC(SOLLYA_X_); b = sollya_lib_build_function_erfc(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_LOG1P(SOLLYA_X_); b = sollya_lib_build_function_log1p(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_EXPM1(SOLLYA_X_); b = sollya_lib_build_function_expm1(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_DE(SOLLYA_X_); b = sollya_lib_build_function_doubleextended(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_CEIL(SOLLYA_X_); b = sollya_lib_build_function_ceil(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_FLOOR(SOLLYA_X_); b = sollya_lib_build_function_floor(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_NEARESTINT(SOLLYA_X_); b = sollya_lib_build_function_nearestint(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_SINH(SOLLYA_X_); b = sollya_lib_build_function_sinh(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_COSH(SOLLYA_X_); b = sollya_lib_build_function_cosh(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_TANH(SOLLYA_X_); b = sollya_lib_build_function_tanh(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ASINH(SOLLYA_X_); b = sollya_lib_build_function_asinh(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ACOSH(SOLLYA_X_); b = sollya_lib_build_function_acosh(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ATANH(SOLLYA_X_); b = sollya_lib_build_function_atanh(SOLLYA_X_);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_ADD(SOLLYA_X_, sollya_lib_constant_from_int(1));
b = sollya_lib_build_function_add((SOLLYA_X_), sollya_lib_constant_from_int(1));
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_SUB(SOLLYA_X_, sollya_lib_constant_from_int(1)); b = sollya_lib_build_function_sub((SOLLYA_X_), ( sollya_lib_constant_from_int(1)));
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_MUL(SOLLYA_X_, sollya_lib_constant_from_int(2)); b = sollya_lib_build_function_mul((SOLLYA_X_), ( sollya_lib_constant_from_int(2)));
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_DIV(SOLLYA_X_, sollya_lib_constant_from_int(2)); b = sollya_lib_build_function_div((SOLLYA_X_), ( sollya_lib_constant_from_int(2)));
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_POW(SOLLYA_X_, sollya_lib_constant_from_int(2)); b = sollya_lib_build_function_pow((SOLLYA_X_), ( sollya_lib_constant_from_int(2)));
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_PI; b = sollya_lib_build_function_pi();
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_CONST(43); b = sollya_lib_constant_from_int(43);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_CONST(43.); b = sollya_lib_constant_from_int(43);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_CONST(1180591620717411303424.); b = sollya_lib_constant_from_double(1180591620717411303424.);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_CONST(0.1); b = sollya_lib_parse_string("0.1");
sollya_lib_printf("%b is not %b (and it is expected)\n",a,b);
sollya_lib_clear_obj(b);
b = sollya_lib_parse_string("D(1/10)");
sollya_lib_printf("... but %b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_CONST_UI64(UINT64_C(281474976710656)); b = sollya_lib_constant_from_double(281474976710656.);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_CONST_UI64(UINT64_C(-281474976710656)); b = sollya_lib_constant_from_uint64(UINT64_C(-281474976710656));
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
a = SOLLYA_CONST_SI64(INT64_C(-281474976710656)); b = sollya_lib_constant_from_double(-281474976710656.);
sollya_lib_printf("%b is %b\n",a,b);
sollya_lib_clear_obj(a);
sollya_lib_clear_obj(b);
sollya_lib_close();
return 0;
}
|
