GCC Code Coverage Report

Directory: ./
File: lib/geogram/numerics/predicates/side2.h
Date: 2026-09-07 02:37:58
Exec Total Coverage
Lines: 431 594 72.6%
Functions: 4 5 80.0%
Branches: 240 352 68.2%
Line Branch Exec Source
1 /* Automatically generated code, do not edit */
2 /* Generated from source file: side2.pck */
3
4 5087500 inline int side2_3d_filter( const double* p0, const double* p1, const double* p2, const double* q0, const double* q1) {
5 5087500 double p1_0_p0_0 = (p1[0] - p0[0]);
6 5087500 double p1_1_p0_1 = (p1[1] - p0[1]);
7 5087500 double p1_2_p0_2 = (p1[2] - p0[2]);
8 double l1;
9 5087500 l1 = (1 * (((p1_0_p0_0 * p1_0_p0_0) + (p1_1_p0_1 * p1_1_p0_1)) + (p1_2_p0_2 * p1_2_p0_2)));
10 5087500 double p2_0_p0_0 = (p2[0] - p0[0]);
11 5087500 double p2_1_p0_1 = (p2[1] - p0[1]);
12 5087500 double p2_2_p0_2 = (p2[2] - p0[2]);
13 double l2;
14 5087500 l2 = (1 * (((p2_0_p0_0 * p2_0_p0_0) + (p2_1_p0_1 * p2_1_p0_1)) + (p2_2_p0_2 * p2_2_p0_2)));
15 5087500 double q0_0_p0_0 = (q0[0] - p0[0]);
16 5087500 double q0_1_p0_1 = (q0[1] - p0[1]);
17 5087500 double q0_2_p0_2 = (q0[2] - p0[2]);
18 double a10;
19 5087500 a10 = (2 * (((p1_0_p0_0 * q0_0_p0_0) + (p1_1_p0_1 * q0_1_p0_1)) + (p1_2_p0_2 * q0_2_p0_2)));
20 5087500 double q1_0_p0_0 = (q1[0] - p0[0]);
21 5087500 double q1_1_p0_1 = (q1[1] - p0[1]);
22 5087500 double q1_2_p0_2 = (q1[2] - p0[2]);
23 double a11;
24 5087500 a11 = (2 * (((p1_0_p0_0 * q1_0_p0_0) + (p1_1_p0_1 * q1_1_p0_1)) + (p1_2_p0_2 * q1_2_p0_2)));
25 double a20;
26 5087500 a20 = (2 * (((p2_0_p0_0 * q0_0_p0_0) + (p2_1_p0_1 * q0_1_p0_1)) + (p2_2_p0_2 * q0_2_p0_2)));
27 double a21;
28 5087500 a21 = (2 * (((p2_0_p0_0 * q1_0_p0_0) + (p2_1_p0_1 * q1_1_p0_1)) + (p2_2_p0_2 * q1_2_p0_2)));
29 double Delta;
30 5087500 Delta = (a11 - a10);
31 double DeltaLambda0;
32 5087500 DeltaLambda0 = (a11 - l1);
33 double DeltaLambda1;
34 5087500 DeltaLambda1 = (l1 - a10);
35 double r;
36 5087500 r = (((Delta * l2) - (a20 * DeltaLambda0)) - (a21 * DeltaLambda1));
37 double eps;
38 5087500 double max1 = fabs(p1_0_p0_0);
39
2/2
✓ Branch 0 taken 2521999 times.
✓ Branch 1 taken 2565501 times.
 5087500 if( (max1 < fabs(p1_1_p0_1)) )
40 {
41 2521999 max1 = fabs(p1_1_p0_1);
42 }
43
2/2
✓ Branch 0 taken 1240740 times.
✓ Branch 1 taken 3846760 times.
 5087500 if( (max1 < fabs(p1_2_p0_2)) )
44 {
45 1240740 max1 = fabs(p1_2_p0_2);
46 }
47 5087500 double max2 = fabs(q0_0_p0_0);
48
2/2
✓ Branch 0 taken 1684220 times.
✓ Branch 1 taken 3403280 times.
 5087500 if( (max2 < fabs(q0_1_p0_1)) )
49 {
50 1684220 max2 = fabs(q0_1_p0_1);
51 }
52
2/2
✓ Branch 0 taken 787965 times.
✓ Branch 1 taken 4299535 times.
 5087500 if( (max2 < fabs(q0_2_p0_2)) )
53 {
54 787965 max2 = fabs(q0_2_p0_2);
55 }
56
2/2
✓ Branch 0 taken 1616768 times.
✓ Branch 1 taken 3470732 times.
 5087500 if( (max2 < fabs(q1_0_p0_0)) )
57 {
58 1616768 max2 = fabs(q1_0_p0_0);
59 }
60
2/2
✓ Branch 0 taken 1180793 times.
✓ Branch 1 taken 3906707 times.
 5087500 if( (max2 < fabs(q1_1_p0_1)) )
61 {
62 1180793 max2 = fabs(q1_1_p0_1);
63 }
64
2/2
✓ Branch 0 taken 583144 times.
✓ Branch 1 taken 4504356 times.
 5087500 if( (max2 < fabs(q1_2_p0_2)) )
65 {
66 583144 max2 = fabs(q1_2_p0_2);
67 }
68 double lower_bound_1;
69 double upper_bound_1;
70 int Delta_sign;
71 int int_tmp_result;
72 5087500 lower_bound_1 = max1;
73 5087500 upper_bound_1 = max1;
74
2/2
✓ Branch 0 taken 341591 times.
✓ Branch 1 taken 4745909 times.
 5087500 if( (max2 < lower_bound_1) )
75 {
76 341591 lower_bound_1 = max2;
77 }
78 else
79 {
80
2/2
✓ Branch 0 taken 2859576 times.
✓ Branch 1 taken 1886333 times.
 4745909 if( (max2 > upper_bound_1) )
81 {
82 2859576 upper_bound_1 = max2;
83 }
84 }
85
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5087500 times.
 5087500 if( (lower_bound_1 < 2.23755023300058943229e-147) )
86 {
87 ✗ return FPG_UNCERTAIN_VALUE;
88 }
89 else
90 {
91
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5087500 times.
 5087500 if( (upper_bound_1 > 3.74144419156711063983e+50) )
92 {
93 ✗ return FPG_UNCERTAIN_VALUE;
94 }
95 5087500 eps = (4.44425370757048798480e-15 * (max1 * max2));
96
2/2
✓ Branch 0 taken 2994940 times.
✓ Branch 1 taken 2092560 times.
 5087500 if( (Delta > eps) )
97 {
98 2994940 int_tmp_result = 1;
99 }
100 else
101 {
102
2/2
✓ Branch 0 taken 2092118 times.
✓ Branch 1 taken 442 times.
 2092560 if( (Delta < -eps) )
103 {
104 2092118 int_tmp_result = -1;
105 }
106 else
107 {
108 442 return FPG_UNCERTAIN_VALUE;
109 }
110 }
111 }
112 5087058 Delta_sign = int_tmp_result;
113 5087058 double max3 = max1;
114
2/2
✓ Branch 0 taken 2859134 times.
✓ Branch 1 taken 2227924 times.
 5087058 if( (max3 < max2) )
115 {
116 2859134 max3 = max2;
117 }
118 5087058 double max4 = max2;
119
2/2
✓ Branch 0 taken 772031 times.
✓ Branch 1 taken 4315027 times.
 5087058 if( (max4 < fabs(p2_0_p0_0)) )
120 {
121 772031 max4 = fabs(p2_0_p0_0);
122 }
123
2/2
✓ Branch 0 taken 458891 times.
✓ Branch 1 taken 4628167 times.
 5087058 if( (max4 < fabs(p2_1_p0_1)) )
124 {
125 458891 max4 = fabs(p2_1_p0_1);
126 }
127
2/2
✓ Branch 0 taken 402928 times.
✓ Branch 1 taken 4684130 times.
 5087058 if( (max4 < fabs(p2_2_p0_2)) )
128 {
129 402928 max4 = fabs(p2_2_p0_2);
130 }
131
2/2
✓ Branch 0 taken 1146682 times.
✓ Branch 1 taken 3940376 times.
 5087058 if( (max3 < max4) )
132 {
133 1146682 max3 = max4;
134 }
135 int r_sign;
136 int int_tmp_result_FFWKCAA;
137 5087058 lower_bound_1 = max1;
138 5087058 upper_bound_1 = max1;
139
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5087058 times.
 5087058 if( (max3 < lower_bound_1) )
140 {
141 ✗ lower_bound_1 = max3;
142 }
143 else
144 {
145
2/2
✓ Branch 0 taken 3738078 times.
✓ Branch 1 taken 1348980 times.
 5087058 if( (max3 > upper_bound_1) )
146 {
147 3738078 upper_bound_1 = max3;
148 }
149 }
150
2/2
✓ Branch 0 taken 49740 times.
✓ Branch 1 taken 5037318 times.
 5087058 if( (max4 < lower_bound_1) )
151 {
152 49740 lower_bound_1 = max4;
153 }
154
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5087058 times.
 5087058 if( (lower_bound_1 < 2.22985945097100191780e-74) )
155 {
156 ✗ return FPG_UNCERTAIN_VALUE;
157 }
158 else
159 {
160
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5087058 times.
 5087058 if( (upper_bound_1 > 3.74144419156711063983e+50) )
161 {
162 ✗ return FPG_UNCERTAIN_VALUE;
163 }
164 5087058 eps = (8.99983341597279045654e-14 * (((max1 * max4) * max4) * max3));
165
2/2
✓ Branch 0 taken 2781865 times.
✓ Branch 1 taken 2305193 times.
 5087058 if( (r > eps) )
166 {
167 2781865 int_tmp_result_FFWKCAA = 1;
168 }
169 else
170 {
171
2/2
✓ Branch 0 taken 2026115 times.
✓ Branch 1 taken 279078 times.
 2305193 if( (r < -eps) )
172 {
173 2026115 int_tmp_result_FFWKCAA = -1;
174 }
175 else
176 {
177 279078 return FPG_UNCERTAIN_VALUE;
178 }
179 }
180 }
181 4807980 r_sign = int_tmp_result_FFWKCAA;
182 4807980 return (Delta_sign * r_sign);
183 }
184
185
186 187020 inline int side2_4d_filter( const double* p0, const double* p1, const double* p2, const double* q0, const double* q1) {
187 187020 double p1_0_p0_0 = (p1[0] - p0[0]);
188 187020 double p1_1_p0_1 = (p1[1] - p0[1]);
189 187020 double p1_2_p0_2 = (p1[2] - p0[2]);
190 187020 double p1_3_p0_3 = (p1[3] - p0[3]);
191 double l1;
192 187020 l1 = (1 * ((((p1_0_p0_0 * p1_0_p0_0) + (p1_1_p0_1 * p1_1_p0_1)) + (p1_2_p0_2 * p1_2_p0_2)) + (p1_3_p0_3 * p1_3_p0_3)));
193 187020 double p2_0_p0_0 = (p2[0] - p0[0]);
194 187020 double p2_1_p0_1 = (p2[1] - p0[1]);
195 187020 double p2_2_p0_2 = (p2[2] - p0[2]);
196 187020 double p2_3_p0_3 = (p2[3] - p0[3]);
197 double l2;
198 187020 l2 = (1 * ((((p2_0_p0_0 * p2_0_p0_0) + (p2_1_p0_1 * p2_1_p0_1)) + (p2_2_p0_2 * p2_2_p0_2)) + (p2_3_p0_3 * p2_3_p0_3)));
199 187020 double q0_0_p0_0 = (q0[0] - p0[0]);
200 187020 double q0_1_p0_1 = (q0[1] - p0[1]);
201 187020 double q0_2_p0_2 = (q0[2] - p0[2]);
202 187020 double q0_3_p0_3 = (q0[3] - p0[3]);
203 double a10;
204 187020 a10 = (2 * ((((p1_0_p0_0 * q0_0_p0_0) + (p1_1_p0_1 * q0_1_p0_1)) + (p1_2_p0_2 * q0_2_p0_2)) + (p1_3_p0_3 * q0_3_p0_3)));
205 187020 double q1_0_p0_0 = (q1[0] - p0[0]);
206 187020 double q1_1_p0_1 = (q1[1] - p0[1]);
207 187020 double q1_2_p0_2 = (q1[2] - p0[2]);
208 187020 double q1_3_p0_3 = (q1[3] - p0[3]);
209 double a11;
210 187020 a11 = (2 * ((((p1_0_p0_0 * q1_0_p0_0) + (p1_1_p0_1 * q1_1_p0_1)) + (p1_2_p0_2 * q1_2_p0_2)) + (p1_3_p0_3 * q1_3_p0_3)));
211 double a20;
212 187020 a20 = (2 * ((((p2_0_p0_0 * q0_0_p0_0) + (p2_1_p0_1 * q0_1_p0_1)) + (p2_2_p0_2 * q0_2_p0_2)) + (p2_3_p0_3 * q0_3_p0_3)));
213 double a21;
214 187020 a21 = (2 * ((((p2_0_p0_0 * q1_0_p0_0) + (p2_1_p0_1 * q1_1_p0_1)) + (p2_2_p0_2 * q1_2_p0_2)) + (p2_3_p0_3 * q1_3_p0_3)));
215 double Delta;
216 187020 Delta = (a11 - a10);
217 double DeltaLambda0;
218 187020 DeltaLambda0 = (a11 - l1);
219 double DeltaLambda1;
220 187020 DeltaLambda1 = (l1 - a10);
221 double r;
222 187020 r = (((Delta * l2) - (a20 * DeltaLambda0)) - (a21 * DeltaLambda1));
223 double eps;
224 187020 double max1 = fabs(p1_2_p0_2);
225
2/2
✓ Branch 0 taken 77333 times.
✓ Branch 1 taken 109687 times.
 187020 if( (max1 < fabs(p1_0_p0_0)) )
226 {
227 77333 max1 = fabs(p1_0_p0_0);
228 }
229
2/2
✓ Branch 0 taken 75428 times.
✓ Branch 1 taken 111592 times.
 187020 if( (max1 < fabs(p1_1_p0_1)) )
230 {
231 75428 max1 = fabs(p1_1_p0_1);
232 }
233
2/2
✓ Branch 0 taken 62405 times.
✓ Branch 1 taken 124615 times.
 187020 if( (max1 < fabs(p1_3_p0_3)) )
234 {
235 62405 max1 = fabs(p1_3_p0_3);
236 }
237 187020 double max2 = fabs(q0_0_p0_0);
238
2/2
✓ Branch 0 taken 91158 times.
✓ Branch 1 taken 95862 times.
 187020 if( (max2 < fabs(q0_1_p0_1)) )
239 {
240 91158 max2 = fabs(q0_1_p0_1);
241 }
242
2/2
✓ Branch 0 taken 60516 times.
✓ Branch 1 taken 126504 times.
 187020 if( (max2 < fabs(q0_2_p0_2)) )
243 {
244 60516 max2 = fabs(q0_2_p0_2);
245 }
246
2/2
✓ Branch 0 taken 150696 times.
✓ Branch 1 taken 36324 times.
 187020 if( (max2 < fabs(q0_3_p0_3)) )
247 {
248 150696 max2 = fabs(q0_3_p0_3);
249 }
250
2/2
✓ Branch 0 taken 8291 times.
✓ Branch 1 taken 178729 times.
 187020 if( (max2 < fabs(q1_0_p0_0)) )
251 {
252 8291 max2 = fabs(q1_0_p0_0);
253 }
254
2/2
✓ Branch 0 taken 4445 times.
✓ Branch 1 taken 182575 times.
 187020 if( (max2 < fabs(q1_1_p0_1)) )
255 {
256 4445 max2 = fabs(q1_1_p0_1);
257 }
258
2/2
✓ Branch 0 taken 3760 times.
✓ Branch 1 taken 183260 times.
 187020 if( (max2 < fabs(q1_2_p0_2)) )
259 {
260 3760 max2 = fabs(q1_2_p0_2);
261 }
262
2/2
✓ Branch 0 taken 68261 times.
✓ Branch 1 taken 118759 times.
 187020 if( (max2 < fabs(q1_3_p0_3)) )
263 {
264 68261 max2 = fabs(q1_3_p0_3);
265 }
266 double lower_bound_1;
267 double upper_bound_1;
268 int Delta_sign;
269 int int_tmp_result;
270 187020 lower_bound_1 = max1;
271 187020 upper_bound_1 = max1;
272
2/2
✓ Branch 0 taken 7667 times.
✓ Branch 1 taken 179353 times.
 187020 if( (max2 < lower_bound_1) )
273 {
274 7667 lower_bound_1 = max2;
275 }
276 else
277 {
278
1/2
✓ Branch 0 taken 179353 times.
✗ Branch 1 not taken.
 179353 if( (max2 > upper_bound_1) )
279 {
280 179353 upper_bound_1 = max2;
281 }
282 }
283
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 187020 times.
 187020 if( (lower_bound_1 < 1.85816790703293534018e-147) )
284 {
285 ✗ return FPG_UNCERTAIN_VALUE;
286 }
287 else
288 {
289
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 187020 times.
 187020 if( (upper_bound_1 > 1.87072209578355531992e+50) )
290 {
291 ✗ return FPG_UNCERTAIN_VALUE;
292 }
293 187020 eps = (6.44428177279185717888e-15 * (max1 * max2));
294
2/2
✓ Branch 0 taken 56245 times.
✓ Branch 1 taken 130775 times.
 187020 if( (Delta > eps) )
295 {
296 56245 int_tmp_result = 1;
297 }
298 else
299 {
300
1/2
✓ Branch 0 taken 130775 times.
✗ Branch 1 not taken.
 130775 if( (Delta < -eps) )
301 {
302 130775 int_tmp_result = -1;
303 }
304 else
305 {
306 ✗ return FPG_UNCERTAIN_VALUE;
307 }
308 }
309 }
310 187020 Delta_sign = int_tmp_result;
311 187020 double max3 = max1;
312
2/2
✓ Branch 0 taken 179353 times.
✓ Branch 1 taken 7667 times.
 187020 if( (max3 < max2) )
313 {
314 179353 max3 = max2;
315 }
316 187020 double max4 = max2;
317
2/2
✓ Branch 0 taken 194 times.
✓ Branch 1 taken 186826 times.
 187020 if( (max4 < fabs(p2_0_p0_0)) )
318 {
319 194 max4 = fabs(p2_0_p0_0);
320 }
321
2/2
✓ Branch 0 taken 6207 times.
✓ Branch 1 taken 180813 times.
 187020 if( (max4 < fabs(p2_2_p0_2)) )
322 {
323 6207 max4 = fabs(p2_2_p0_2);
324 }
325
2/2
✓ Branch 0 taken 5172 times.
✓ Branch 1 taken 181848 times.
 187020 if( (max4 < fabs(p2_1_p0_1)) )
326 {
327 5172 max4 = fabs(p2_1_p0_1);
328 }
329
2/2
✓ Branch 0 taken 4179 times.
✓ Branch 1 taken 182841 times.
 187020 if( (max4 < fabs(p2_3_p0_3)) )
330 {
331 4179 max4 = fabs(p2_3_p0_3);
332 }
333
2/2
✓ Branch 0 taken 12289 times.
✓ Branch 1 taken 174731 times.
 187020 if( (max3 < max4) )
334 {
335 12289 max3 = max4;
336 }
337 int r_sign;
338 int int_tmp_result_FFWKCAA;
339 187020 lower_bound_1 = max3;
340 187020 upper_bound_1 = max3;
341
2/2
✓ Branch 0 taken 184699 times.
✓ Branch 1 taken 2321 times.
 187020 if( (max1 < lower_bound_1) )
342 {
343 184699 lower_bound_1 = max1;
344 }
345
2/2
✓ Branch 0 taken 2321 times.
✓ Branch 1 taken 184699 times.
 187020 if( (max4 < lower_bound_1) )
346 {
347 2321 lower_bound_1 = max4;
348 }
349
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 187020 times.
 187020 if( (lower_bound_1 < 1.89528395402941802921e-74) )
350 {
351 ✗ return FPG_UNCERTAIN_VALUE;
352 }
353 else
354 {
355
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 187020 times.
 187020 if( (upper_bound_1 > 1.87072209578355531992e+50) )
356 {
357 ✗ return FPG_UNCERTAIN_VALUE;
358 }
359 187020 eps = (1.72443682410932010423e-13 * (((max1 * max4) * max4) * max3));
360
2/2
✓ Branch 0 taken 67387 times.
✓ Branch 1 taken 119633 times.
 187020 if( (r > eps) )
361 {
362 67387 int_tmp_result_FFWKCAA = 1;
363 }
364 else
365 {
366
1/2
✓ Branch 0 taken 119633 times.
✗ Branch 1 not taken.
 119633 if( (r < -eps) )
367 {
368 119633 int_tmp_result_FFWKCAA = -1;
369 }
370 else
371 {
372 ✗ return FPG_UNCERTAIN_VALUE;
373 }
374 }
375 }
376 187020 r_sign = int_tmp_result_FFWKCAA;
377 187020 return (Delta_sign * r_sign);
378 }
379
380
381 269465 inline int side2_6d_filter( const double* p0, const double* p1, const double* p2, const double* q0, const double* q1) {
382 269465 double p1_0_p0_0 = (p1[0] - p0[0]);
383 269465 double p1_1_p0_1 = (p1[1] - p0[1]);
384 269465 double p1_2_p0_2 = (p1[2] - p0[2]);
385 269465 double p1_3_p0_3 = (p1[3] - p0[3]);
386 269465 double p1_4_p0_4 = (p1[4] - p0[4]);
387 269465 double p1_5_p0_5 = (p1[5] - p0[5]);
388 double l1;
389 269465 l1 = (1 * ((((((p1_0_p0_0 * p1_0_p0_0) + (p1_1_p0_1 * p1_1_p0_1)) + (p1_2_p0_2 * p1_2_p0_2)) + (p1_3_p0_3 * p1_3_p0_3)) + (p1_4_p0_4 * p1_4_p0_4)) + (p1_5_p0_5 * p1_5_p0_5)));
390 269465 double p2_0_p0_0 = (p2[0] - p0[0]);
391 269465 double p2_1_p0_1 = (p2[1] - p0[1]);
392 269465 double p2_2_p0_2 = (p2[2] - p0[2]);
393 269465 double p2_3_p0_3 = (p2[3] - p0[3]);
394 269465 double p2_4_p0_4 = (p2[4] - p0[4]);
395 269465 double p2_5_p0_5 = (p2[5] - p0[5]);
396 double l2;
397 269465 l2 = (1 * ((((((p2_0_p0_0 * p2_0_p0_0) + (p2_1_p0_1 * p2_1_p0_1)) + (p2_2_p0_2 * p2_2_p0_2)) + (p2_3_p0_3 * p2_3_p0_3)) + (p2_4_p0_4 * p2_4_p0_4)) + (p2_5_p0_5 * p2_5_p0_5)));
398 269465 double q0_0_p0_0 = (q0[0] - p0[0]);
399 269465 double q0_1_p0_1 = (q0[1] - p0[1]);
400 269465 double q0_2_p0_2 = (q0[2] - p0[2]);
401 269465 double q0_3_p0_3 = (q0[3] - p0[3]);
402 269465 double q0_4_p0_4 = (q0[4] - p0[4]);
403 269465 double q0_5_p0_5 = (q0[5] - p0[5]);
404 double a10;
405 269465 a10 = (2 * ((((((p1_0_p0_0 * q0_0_p0_0) + (p1_1_p0_1 * q0_1_p0_1)) + (p1_2_p0_2 * q0_2_p0_2)) + (p1_3_p0_3 * q0_3_p0_3)) + (p1_4_p0_4 * q0_4_p0_4)) + (p1_5_p0_5 * q0_5_p0_5)));
406 269465 double q1_0_p0_0 = (q1[0] - p0[0]);
407 269465 double q1_1_p0_1 = (q1[1] - p0[1]);
408 269465 double q1_2_p0_2 = (q1[2] - p0[2]);
409 269465 double q1_3_p0_3 = (q1[3] - p0[3]);
410 269465 double q1_4_p0_4 = (q1[4] - p0[4]);
411 269465 double q1_5_p0_5 = (q1[5] - p0[5]);
412 double a11;
413 269465 a11 = (2 * ((((((p1_0_p0_0 * q1_0_p0_0) + (p1_1_p0_1 * q1_1_p0_1)) + (p1_2_p0_2 * q1_2_p0_2)) + (p1_3_p0_3 * q1_3_p0_3)) + (p1_4_p0_4 * q1_4_p0_4)) + (p1_5_p0_5 * q1_5_p0_5)));
414 double a20;
415 269465 a20 = (2 * ((((((p2_0_p0_0 * q0_0_p0_0) + (p2_1_p0_1 * q0_1_p0_1)) + (p2_2_p0_2 * q0_2_p0_2)) + (p2_3_p0_3 * q0_3_p0_3)) + (p2_4_p0_4 * q0_4_p0_4)) + (p2_5_p0_5 * q0_5_p0_5)));
416 double a21;
417 269465 a21 = (2 * ((((((p2_0_p0_0 * q1_0_p0_0) + (p2_1_p0_1 * q1_1_p0_1)) + (p2_2_p0_2 * q1_2_p0_2)) + (p2_3_p0_3 * q1_3_p0_3)) + (p2_4_p0_4 * q1_4_p0_4)) + (p2_5_p0_5 * q1_5_p0_5)));
418 double Delta;
419 269465 Delta = (a11 - a10);
420 double DeltaLambda0;
421 269465 DeltaLambda0 = (a11 - l1);
422 double DeltaLambda1;
423 269465 DeltaLambda1 = (l1 - a10);
424 double r;
425 269465 r = (((Delta * l2) - (a20 * DeltaLambda0)) - (a21 * DeltaLambda1));
426 double eps;
427 269465 double max1 = fabs(p1_1_p0_1);
428
2/2
✓ Branch 0 taken 129567 times.
✓ Branch 1 taken 139898 times.
 269465 if( (max1 < fabs(p1_2_p0_2)) )
429 {
430 129567 max1 = fabs(p1_2_p0_2);
431 }
432
2/2
✓ Branch 0 taken 59504 times.
✓ Branch 1 taken 209961 times.
 269465 if( (max1 < fabs(p1_3_p0_3)) )
433 {
434 59504 max1 = fabs(p1_3_p0_3);
435 }
436
2/2
✓ Branch 0 taken 1543 times.
✓ Branch 1 taken 267922 times.
 269465 if( (max1 < fabs(p1_0_p0_0)) )
437 {
438 1543 max1 = fabs(p1_0_p0_0);
439 }
440
2/2
✓ Branch 0 taken 51282 times.
✓ Branch 1 taken 218183 times.
 269465 if( (max1 < fabs(p1_4_p0_4)) )
441 {
442 51282 max1 = fabs(p1_4_p0_4);
443 }
444
2/2
✓ Branch 0 taken 10269 times.
✓ Branch 1 taken 259196 times.
 269465 if( (max1 < fabs(p1_5_p0_5)) )
445 {
446 10269 max1 = fabs(p1_5_p0_5);
447 }
448 269465 double max2 = fabs(q0_0_p0_0);
449
2/2
✓ Branch 0 taken 173101 times.
✓ Branch 1 taken 96364 times.
 269465 if( (max2 < fabs(q0_1_p0_1)) )
450 {
451 173101 max2 = fabs(q0_1_p0_1);
452 }
453
2/2
✓ Branch 0 taken 98196 times.
✓ Branch 1 taken 171269 times.
 269465 if( (max2 < fabs(q0_2_p0_2)) )
454 {
455 98196 max2 = fabs(q0_2_p0_2);
456 }
457
2/2
✓ Branch 0 taken 211749 times.
✓ Branch 1 taken 57716 times.
 269465 if( (max2 < fabs(q0_3_p0_3)) )
458 {
459 211749 max2 = fabs(q0_3_p0_3);
460 }
461
2/2
✓ Branch 0 taken 204240 times.
✓ Branch 1 taken 65225 times.
 269465 if( (max2 < fabs(q0_4_p0_4)) )
462 {
463 204240 max2 = fabs(q0_4_p0_4);
464 }
465
2/2
✓ Branch 0 taken 43692 times.
✓ Branch 1 taken 225773 times.
 269465 if( (max2 < fabs(q0_5_p0_5)) )
466 {
467 43692 max2 = fabs(q0_5_p0_5);
468 }
469
2/2
✓ Branch 0 taken 23472 times.
✓ Branch 1 taken 245993 times.
 269465 if( (max2 < fabs(q1_0_p0_0)) )
470 {
471 23472 max2 = fabs(q1_0_p0_0);
472 }
473
2/2
✓ Branch 0 taken 7711 times.
✓ Branch 1 taken 261754 times.
 269465 if( (max2 < fabs(q1_1_p0_1)) )
474 {
475 7711 max2 = fabs(q1_1_p0_1);
476 }
477
2/2
✓ Branch 0 taken 5785 times.
✓ Branch 1 taken 263680 times.
 269465 if( (max2 < fabs(q1_2_p0_2)) )
478 {
479 5785 max2 = fabs(q1_2_p0_2);
480 }
481
2/2
✓ Branch 0 taken 96729 times.
✓ Branch 1 taken 172736 times.
 269465 if( (max2 < fabs(q1_3_p0_3)) )
482 {
483 96729 max2 = fabs(q1_3_p0_3);
484 }
485
2/2
✓ Branch 0 taken 118330 times.
✓ Branch 1 taken 151135 times.
 269465 if( (max2 < fabs(q1_4_p0_4)) )
486 {
487 118330 max2 = fabs(q1_4_p0_4);
488 }
489
2/2
✓ Branch 0 taken 31567 times.
✓ Branch 1 taken 237898 times.
 269465 if( (max2 < fabs(q1_5_p0_5)) )
490 {
491 31567 max2 = fabs(q1_5_p0_5);
492 }
493 double lower_bound_1;
494 double upper_bound_1;
495 int Delta_sign;
496 int int_tmp_result;
497 269465 lower_bound_1 = max1;
498 269465 upper_bound_1 = max1;
499
2/2
✓ Branch 0 taken 10899 times.
✓ Branch 1 taken 258566 times.
 269465 if( (max2 < lower_bound_1) )
500 {
501 10899 lower_bound_1 = max2;
502 }
503 else
504 {
505
1/2
✓ Branch 0 taken 258566 times.
✗ Branch 1 not taken.
 258566 if( (max2 > upper_bound_1) )
506 {
507 258566 upper_bound_1 = max2;
508 }
509 }
510
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269465 times.
 269465 if( (lower_bound_1 < 1.41511993781011659868e-147) )
511 {
512 ✗ return FPG_UNCERTAIN_VALUE;
513 }
514 else
515 {
516
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269465 times.
 269465 if( (upper_bound_1 > 1.87072209578355531992e+50) )
517 {
518 ✗ return FPG_UNCERTAIN_VALUE;
519 }
520 269465 eps = (1.11111223981318615596e-14 * (max1 * max2));
521
2/2
✓ Branch 0 taken 98963 times.
✓ Branch 1 taken 170502 times.
 269465 if( (Delta > eps) )
522 {
523 98963 int_tmp_result = 1;
524 }
525 else
526 {
527
1/2
✓ Branch 0 taken 170502 times.
✗ Branch 1 not taken.
 170502 if( (Delta < -eps) )
528 {
529 170502 int_tmp_result = -1;
530 }
531 else
532 {
533 ✗ return FPG_UNCERTAIN_VALUE;
534 }
535 }
536 }
537 269465 Delta_sign = int_tmp_result;
538 269465 double max3 = max1;
539
2/2
✓ Branch 0 taken 258566 times.
✓ Branch 1 taken 10899 times.
 269465 if( (max3 < max2) )
540 {
541 258566 max3 = max2;
542 }
543 269465 double max4 = max2;
544
2/2
✓ Branch 0 taken 6836 times.
✓ Branch 1 taken 262629 times.
 269465 if( (max4 < fabs(p2_2_p0_2)) )
545 {
546 6836 max4 = fabs(p2_2_p0_2);
547 }
548
2/2
✓ Branch 0 taken 2924 times.
✓ Branch 1 taken 266541 times.
 269465 if( (max4 < fabs(p2_4_p0_4)) )
549 {
550 2924 max4 = fabs(p2_4_p0_4);
551 }
552
2/2
✓ Branch 0 taken 1198 times.
✓ Branch 1 taken 268267 times.
 269465 if( (max4 < fabs(p2_3_p0_3)) )
553 {
554 1198 max4 = fabs(p2_3_p0_3);
555 }
556
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269465 times.
 269465 if( (max4 < fabs(p2_0_p0_0)) )
557 {
558 ✗ max4 = fabs(p2_0_p0_0);
559 }
560
2/2
✓ Branch 0 taken 4455 times.
✓ Branch 1 taken 265010 times.
 269465 if( (max4 < fabs(p2_1_p0_1)) )
561 {
562 4455 max4 = fabs(p2_1_p0_1);
563 }
564
2/2
✓ Branch 0 taken 220 times.
✓ Branch 1 taken 269245 times.
 269465 if( (max4 < fabs(p2_5_p0_5)) )
565 {
566 220 max4 = fabs(p2_5_p0_5);
567 }
568
2/2
✓ Branch 0 taken 8656 times.
✓ Branch 1 taken 260809 times.
 269465 if( (max3 < max4) )
569 {
570 8656 max3 = max4;
571 }
572 int r_sign;
573 int int_tmp_result_FFWKCAA;
574 269465 lower_bound_1 = max1;
575 269465 upper_bound_1 = max1;
576
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269465 times.
 269465 if( (max3 < lower_bound_1) )
577 {
578 ✗ lower_bound_1 = max3;
579 }
580 else
581 {
582
2/2
✓ Branch 0 taken 265575 times.
✓ Branch 1 taken 3890 times.
 269465 if( (max3 > upper_bound_1) )
583 {
584 265575 upper_bound_1 = max3;
585 }
586 }
587
2/2
✓ Branch 0 taken 3890 times.
✓ Branch 1 taken 265575 times.
 269465 if( (max4 < lower_bound_1) )
588 {
589 3890 lower_bound_1 = max4;
590 }
591
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269465 times.
 269465 if( (lower_bound_1 < 1.49958502193059513986e-74) )
592 {
593 ✗ return FPG_UNCERTAIN_VALUE;
594 }
595 else
596 {
597
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269465 times.
 269465 if( (upper_bound_1 > 1.87072209578355531992e+50) )
598 {
599 ✗ return FPG_UNCERTAIN_VALUE;
600 }
601 269465 eps = (4.40007476026584016994e-13 * (((max1 * max4) * max4) * max3));
602
2/2
✓ Branch 0 taken 107475 times.
✓ Branch 1 taken 161990 times.
 269465 if( (r > eps) )
603 {
604 107475 int_tmp_result_FFWKCAA = 1;
605 }
606 else
607 {
608
1/2
✓ Branch 0 taken 161990 times.
✗ Branch 1 not taken.
 161990 if( (r < -eps) )
609 {
610 161990 int_tmp_result_FFWKCAA = -1;
611 }
612 else
613 {
614 ✗ return FPG_UNCERTAIN_VALUE;
615 }
616 }
617 }
618 269465 r_sign = int_tmp_result_FFWKCAA;
619 269465 return (Delta_sign * r_sign);
620 }
621
622
623 ✗ inline int side2_7d_filter( const double* p0, const double* p1, const double* p2, const double* q0, const double* q1) {
624 ✗ double p1_0_p0_0 = (p1[0] - p0[0]);
625 ✗ double p1_1_p0_1 = (p1[1] - p0[1]);
626 ✗ double p1_2_p0_2 = (p1[2] - p0[2]);
627 ✗ double p1_3_p0_3 = (p1[3] - p0[3]);
628 ✗ double p1_4_p0_4 = (p1[4] - p0[4]);
629 ✗ double p1_5_p0_5 = (p1[5] - p0[5]);
630 ✗ double p1_6_p0_6 = (p1[6] - p0[6]);
631 double l1;
632 ✗ l1 = (1 * (((((((p1_0_p0_0 * p1_0_p0_0) + (p1_1_p0_1 * p1_1_p0_1)) + (p1_2_p0_2 * p1_2_p0_2)) + (p1_3_p0_3 * p1_3_p0_3)) + (p1_4_p0_4 * p1_4_p0_4)) + (p1_5_p0_5 * p1_5_p0_5)) + (p1_6_p0_6 * p1_6_p0_6)));
633 ✗ double p2_0_p0_0 = (p2[0] - p0[0]);
634 ✗ double p2_1_p0_1 = (p2[1] - p0[1]);
635 ✗ double p2_2_p0_2 = (p2[2] - p0[2]);
636 ✗ double p2_3_p0_3 = (p2[3] - p0[3]);
637 ✗ double p2_4_p0_4 = (p2[4] - p0[4]);
638 ✗ double p2_5_p0_5 = (p2[5] - p0[5]);
639 ✗ double p2_6_p0_6 = (p2[6] - p0[6]);
640 double l2;
641 ✗ l2 = (1 * (((((((p2_0_p0_0 * p2_0_p0_0) + (p2_1_p0_1 * p2_1_p0_1)) + (p2_2_p0_2 * p2_2_p0_2)) + (p2_3_p0_3 * p2_3_p0_3)) + (p2_4_p0_4 * p2_4_p0_4)) + (p2_5_p0_5 * p2_5_p0_5)) + (p2_6_p0_6 * p2_6_p0_6)));
642 ✗ double q0_0_p0_0 = (q0[0] - p0[0]);
643 ✗ double q0_1_p0_1 = (q0[1] - p0[1]);
644 ✗ double q0_2_p0_2 = (q0[2] - p0[2]);
645 ✗ double q0_3_p0_3 = (q0[3] - p0[3]);
646 ✗ double q0_4_p0_4 = (q0[4] - p0[4]);
647 ✗ double q0_5_p0_5 = (q0[5] - p0[5]);
648 ✗ double q0_6_p0_6 = (q0[6] - p0[6]);
649 double a10;
650 ✗ a10 = (2 * (((((((p1_0_p0_0 * q0_0_p0_0) + (p1_1_p0_1 * q0_1_p0_1)) + (p1_2_p0_2 * q0_2_p0_2)) + (p1_3_p0_3 * q0_3_p0_3)) + (p1_4_p0_4 * q0_4_p0_4)) + (p1_5_p0_5 * q0_5_p0_5)) + (p1_6_p0_6 * q0_6_p0_6)));
651 ✗ double q1_0_p0_0 = (q1[0] - p0[0]);
652 ✗ double q1_1_p0_1 = (q1[1] - p0[1]);
653 ✗ double q1_2_p0_2 = (q1[2] - p0[2]);
654 ✗ double q1_3_p0_3 = (q1[3] - p0[3]);
655 ✗ double q1_4_p0_4 = (q1[4] - p0[4]);
656 ✗ double q1_5_p0_5 = (q1[5] - p0[5]);
657 ✗ double q1_6_p0_6 = (q1[6] - p0[6]);
658 double a11;
659 ✗ a11 = (2 * (((((((p1_0_p0_0 * q1_0_p0_0) + (p1_1_p0_1 * q1_1_p0_1)) + (p1_2_p0_2 * q1_2_p0_2)) + (p1_3_p0_3 * q1_3_p0_3)) + (p1_4_p0_4 * q1_4_p0_4)) + (p1_5_p0_5 * q1_5_p0_5)) + (p1_6_p0_6 * q1_6_p0_6)));
660 double a20;
661 ✗ a20 = (2 * (((((((p2_0_p0_0 * q0_0_p0_0) + (p2_1_p0_1 * q0_1_p0_1)) + (p2_2_p0_2 * q0_2_p0_2)) + (p2_3_p0_3 * q0_3_p0_3)) + (p2_4_p0_4 * q0_4_p0_4)) + (p2_5_p0_5 * q0_5_p0_5)) + (p2_6_p0_6 * q0_6_p0_6)));
662 double a21;
663 ✗ a21 = (2 * (((((((p2_0_p0_0 * q1_0_p0_0) + (p2_1_p0_1 * q1_1_p0_1)) + (p2_2_p0_2 * q1_2_p0_2)) + (p2_3_p0_3 * q1_3_p0_3)) + (p2_4_p0_4 * q1_4_p0_4)) + (p2_5_p0_5 * q1_5_p0_5)) + (p2_6_p0_6 * q1_6_p0_6)));
664 double Delta;
665 ✗ Delta = (a11 - a10);
666 double DeltaLambda0;
667 ✗ DeltaLambda0 = (a11 - l1);
668 double DeltaLambda1;
669 ✗ DeltaLambda1 = (l1 - a10);
670 double r;
671 ✗ r = (((Delta * l2) - (a20 * DeltaLambda0)) - (a21 * DeltaLambda1));
672 double eps;
673 ✗ double max1 = fabs(p1_2_p0_2);
674 ✗ if( (max1 < fabs(p1_1_p0_1)) )
675 {
676 ✗ max1 = fabs(p1_1_p0_1);
677 }
678 ✗ if( (max1 < fabs(p1_0_p0_0)) )
679 {
680 ✗ max1 = fabs(p1_0_p0_0);
681 }
682 ✗ if( (max1 < fabs(p1_4_p0_4)) )
683 {
684 ✗ max1 = fabs(p1_4_p0_4);
685 }
686 ✗ if( (max1 < fabs(p1_3_p0_3)) )
687 {
688 ✗ max1 = fabs(p1_3_p0_3);
689 }
690 ✗ if( (max1 < fabs(p1_5_p0_5)) )
691 {
692 ✗ max1 = fabs(p1_5_p0_5);
693 }
694 ✗ if( (max1 < fabs(p1_6_p0_6)) )
695 {
696 ✗ max1 = fabs(p1_6_p0_6);
697 }
698 ✗ double max2 = fabs(q0_0_p0_0);
699 ✗ if( (max2 < fabs(q0_1_p0_1)) )
700 {
701 ✗ max2 = fabs(q0_1_p0_1);
702 }
703 ✗ if( (max2 < fabs(q0_2_p0_2)) )
704 {
705 ✗ max2 = fabs(q0_2_p0_2);
706 }
707 ✗ if( (max2 < fabs(q0_3_p0_3)) )
708 {
709 ✗ max2 = fabs(q0_3_p0_3);
710 }
711 ✗ if( (max2 < fabs(q0_4_p0_4)) )
712 {
713 ✗ max2 = fabs(q0_4_p0_4);
714 }
715 ✗ if( (max2 < fabs(q0_5_p0_5)) )
716 {
717 ✗ max2 = fabs(q0_5_p0_5);
718 }
719 ✗ if( (max2 < fabs(q0_6_p0_6)) )
720 {
721 ✗ max2 = fabs(q0_6_p0_6);
722 }
723 ✗ if( (max2 < fabs(q1_0_p0_0)) )
724 {
725 ✗ max2 = fabs(q1_0_p0_0);
726 }
727 ✗ if( (max2 < fabs(q1_1_p0_1)) )
728 {
729 ✗ max2 = fabs(q1_1_p0_1);
730 }
731 ✗ if( (max2 < fabs(q1_2_p0_2)) )
732 {
733 ✗ max2 = fabs(q1_2_p0_2);
734 }
735 ✗ if( (max2 < fabs(q1_3_p0_3)) )
736 {
737 ✗ max2 = fabs(q1_3_p0_3);
738 }
739 ✗ if( (max2 < fabs(q1_4_p0_4)) )
740 {
741 ✗ max2 = fabs(q1_4_p0_4);
742 }
743 ✗ if( (max2 < fabs(q1_5_p0_5)) )
744 {
745 ✗ max2 = fabs(q1_5_p0_5);
746 }
747 ✗ if( (max2 < fabs(q1_6_p0_6)) )
748 {
749 ✗ max2 = fabs(q1_6_p0_6);
750 }
751 double lower_bound_1;
752 double upper_bound_1;
753 int Delta_sign;
754 int int_tmp_result;
755 ✗ lower_bound_1 = max1;
756 ✗ upper_bound_1 = max1;
757 ✗ if( (max2 < lower_bound_1) )
758 {
759 ✗ lower_bound_1 = max2;
760 }
761 else
762 {
763 ✗ if( (max2 > upper_bound_1) )
764 {
765 ✗ upper_bound_1 = max2;
766 }
767 }
768 ✗ if( (lower_bound_1 < 1.27080861580266953580e-147) )
769 {
770 ✗ return FPG_UNCERTAIN_VALUE;
771 }
772 else
773 {
774 ✗ if( (upper_bound_1 > 1.87072209578355531992e+50) )
775 {
776 ✗ return FPG_UNCERTAIN_VALUE;
777 }
778 ✗ eps = (1.37779349582504943796e-14 * (max1 * max2));
779 ✗ if( (Delta > eps) )
780 {
781 ✗ int_tmp_result = 1;
782 }
783 else
784 {
785 ✗ if( (Delta < -eps) )
786 {
787 ✗ int_tmp_result = -1;
788 }
789 else
790 {
791 ✗ return FPG_UNCERTAIN_VALUE;
792 }
793 }
794 }
795 ✗ Delta_sign = int_tmp_result;
796 ✗ double max3 = max1;
797 ✗ if( (max3 < max2) )
798 {
799 ✗ max3 = max2;
800 }
801 ✗ double max4 = max2;
802 ✗ if( (max4 < fabs(p2_0_p0_0)) )
803 {
804 ✗ max4 = fabs(p2_0_p0_0);
805 }
806 ✗ if( (max4 < fabs(p2_1_p0_1)) )
807 {
808 ✗ max4 = fabs(p2_1_p0_1);
809 }
810 ✗ if( (max4 < fabs(p2_2_p0_2)) )
811 {
812 ✗ max4 = fabs(p2_2_p0_2);
813 }
814 ✗ if( (max4 < fabs(p2_3_p0_3)) )
815 {
816 ✗ max4 = fabs(p2_3_p0_3);
817 }
818 ✗ if( (max4 < fabs(p2_4_p0_4)) )
819 {
820 ✗ max4 = fabs(p2_4_p0_4);
821 }
822 ✗ if( (max4 < fabs(p2_5_p0_5)) )
823 {
824 ✗ max4 = fabs(p2_5_p0_5);
825 }
826 ✗ if( (max4 < fabs(p2_6_p0_6)) )
827 {
828 ✗ max4 = fabs(p2_6_p0_6);
829 }
830 ✗ if( (max3 < max4) )
831 {
832 ✗ max3 = max4;
833 }
834 int r_sign;
835 int int_tmp_result_FFWKCAA;
836 ✗ lower_bound_1 = max1;
837 ✗ upper_bound_1 = max1;
838 ✗ if( (max3 < lower_bound_1) )
839 {
840 ✗ lower_bound_1 = max3;
841 }
842 else
843 {
844 ✗ if( (max3 > upper_bound_1) )
845 {
846 ✗ upper_bound_1 = max3;
847 }
848 }
849 ✗ if( (max4 < lower_bound_1) )
850 {
851 ✗ lower_bound_1 = max4;
852 }
853 ✗ if( (lower_bound_1 < 1.36918881183883509035e-74) )
854 {
855 ✗ return FPG_UNCERTAIN_VALUE;
856 }
857 else
858 {
859 ✗ if( (upper_bound_1 > 1.87072209578355531992e+50) )
860 {
861 ✗ return FPG_UNCERTAIN_VALUE;
862 }
863 ✗ eps = (6.33127335329798996022e-13 * (((max1 * max4) * max4) * max3));
864 ✗ if( (r > eps) )
865 {
866 ✗ int_tmp_result_FFWKCAA = 1;
867 }
868 else
869 {
870 ✗ if( (r < -eps) )
871 {
872 ✗ int_tmp_result_FFWKCAA = -1;
873 }
874 else
875 {
876 ✗ return FPG_UNCERTAIN_VALUE;
877 }
878 }
879 }
880 ✗ r_sign = int_tmp_result_FFWKCAA;
881 ✗ return (Delta_sign * r_sign);
882 }
883
884
885 269311 inline int side2_8d_filter( const double* p0, const double* p1, const double* p2, const double* q0, const double* q1) {
886 269311 double p1_0_p0_0 = (p1[0] - p0[0]);
887 269311 double p1_1_p0_1 = (p1[1] - p0[1]);
888 269311 double p1_2_p0_2 = (p1[2] - p0[2]);
889 269311 double p1_3_p0_3 = (p1[3] - p0[3]);
890 269311 double p1_4_p0_4 = (p1[4] - p0[4]);
891 269311 double p1_5_p0_5 = (p1[5] - p0[5]);
892 269311 double p1_6_p0_6 = (p1[6] - p0[6]);
893 269311 double p1_7_p0_7 = (p1[7] - p0[7]);
894 double l1;
895 269311 l1 = (1 * ((((((((p1_0_p0_0 * p1_0_p0_0) + (p1_1_p0_1 * p1_1_p0_1)) + (p1_2_p0_2 * p1_2_p0_2)) + (p1_3_p0_3 * p1_3_p0_3)) + (p1_4_p0_4 * p1_4_p0_4)) + (p1_5_p0_5 * p1_5_p0_5)) + (p1_6_p0_6 * p1_6_p0_6)) + (p1_7_p0_7 * p1_7_p0_7)));
896 269311 double p2_0_p0_0 = (p2[0] - p0[0]);
897 269311 double p2_1_p0_1 = (p2[1] - p0[1]);
898 269311 double p2_2_p0_2 = (p2[2] - p0[2]);
899 269311 double p2_3_p0_3 = (p2[3] - p0[3]);
900 269311 double p2_4_p0_4 = (p2[4] - p0[4]);
901 269311 double p2_5_p0_5 = (p2[5] - p0[5]);
902 269311 double p2_6_p0_6 = (p2[6] - p0[6]);
903 269311 double p2_7_p0_7 = (p2[7] - p0[7]);
904 double l2;
905 269311 l2 = (1 * ((((((((p2_0_p0_0 * p2_0_p0_0) + (p2_1_p0_1 * p2_1_p0_1)) + (p2_2_p0_2 * p2_2_p0_2)) + (p2_3_p0_3 * p2_3_p0_3)) + (p2_4_p0_4 * p2_4_p0_4)) + (p2_5_p0_5 * p2_5_p0_5)) + (p2_6_p0_6 * p2_6_p0_6)) + (p2_7_p0_7 * p2_7_p0_7)));
906 269311 double q0_0_p0_0 = (q0[0] - p0[0]);
907 269311 double q0_1_p0_1 = (q0[1] - p0[1]);
908 269311 double q0_2_p0_2 = (q0[2] - p0[2]);
909 269311 double q0_3_p0_3 = (q0[3] - p0[3]);
910 269311 double q0_4_p0_4 = (q0[4] - p0[4]);
911 269311 double q0_5_p0_5 = (q0[5] - p0[5]);
912 269311 double q0_6_p0_6 = (q0[6] - p0[6]);
913 269311 double q0_7_p0_7 = (q0[7] - p0[7]);
914 double a10;
915 269311 a10 = (2 * ((((((((p1_0_p0_0 * q0_0_p0_0) + (p1_1_p0_1 * q0_1_p0_1)) + (p1_2_p0_2 * q0_2_p0_2)) + (p1_3_p0_3 * q0_3_p0_3)) + (p1_4_p0_4 * q0_4_p0_4)) + (p1_5_p0_5 * q0_5_p0_5)) + (p1_6_p0_6 * q0_6_p0_6)) + (p1_7_p0_7 * q0_7_p0_7)));
916 269311 double q1_0_p0_0 = (q1[0] - p0[0]);
917 269311 double q1_1_p0_1 = (q1[1] - p0[1]);
918 269311 double q1_2_p0_2 = (q1[2] - p0[2]);
919 269311 double q1_3_p0_3 = (q1[3] - p0[3]);
920 269311 double q1_4_p0_4 = (q1[4] - p0[4]);
921 269311 double q1_5_p0_5 = (q1[5] - p0[5]);
922 269311 double q1_6_p0_6 = (q1[6] - p0[6]);
923 269311 double q1_7_p0_7 = (q1[7] - p0[7]);
924 double a11;
925 269311 a11 = (2 * ((((((((p1_0_p0_0 * q1_0_p0_0) + (p1_1_p0_1 * q1_1_p0_1)) + (p1_2_p0_2 * q1_2_p0_2)) + (p1_3_p0_3 * q1_3_p0_3)) + (p1_4_p0_4 * q1_4_p0_4)) + (p1_5_p0_5 * q1_5_p0_5)) + (p1_6_p0_6 * q1_6_p0_6)) + (p1_7_p0_7 * q1_7_p0_7)));
926 double a20;
927 269311 a20 = (2 * ((((((((p2_0_p0_0 * q0_0_p0_0) + (p2_1_p0_1 * q0_1_p0_1)) + (p2_2_p0_2 * q0_2_p0_2)) + (p2_3_p0_3 * q0_3_p0_3)) + (p2_4_p0_4 * q0_4_p0_4)) + (p2_5_p0_5 * q0_5_p0_5)) + (p2_6_p0_6 * q0_6_p0_6)) + (p2_7_p0_7 * q0_7_p0_7)));
928 double a21;
929 269311 a21 = (2 * ((((((((p2_0_p0_0 * q1_0_p0_0) + (p2_1_p0_1 * q1_1_p0_1)) + (p2_2_p0_2 * q1_2_p0_2)) + (p2_3_p0_3 * q1_3_p0_3)) + (p2_4_p0_4 * q1_4_p0_4)) + (p2_5_p0_5 * q1_5_p0_5)) + (p2_6_p0_6 * q1_6_p0_6)) + (p2_7_p0_7 * q1_7_p0_7)));
930 double Delta;
931 269311 Delta = (a11 - a10);
932 double DeltaLambda0;
933 269311 DeltaLambda0 = (a11 - l1);
934 double DeltaLambda1;
935 269311 DeltaLambda1 = (l1 - a10);
936 double r;
937 269311 r = (((Delta * l2) - (a20 * DeltaLambda0)) - (a21 * DeltaLambda1));
938 double eps;
939 269311 double max1 = fabs(p1_4_p0_4);
940
2/2
✓ Branch 0 taken 86236 times.
✓ Branch 1 taken 183075 times.
 269311 if( (max1 < fabs(p1_3_p0_3)) )
941 {
942 86236 max1 = fabs(p1_3_p0_3);
943 }
944
2/2
✓ Branch 0 taken 91497 times.
✓ Branch 1 taken 177814 times.
 269311 if( (max1 < fabs(p1_7_p0_7)) )
945 {
946 91497 max1 = fabs(p1_7_p0_7);
947 }
948
2/2
✓ Branch 0 taken 879 times.
✓ Branch 1 taken 268432 times.
 269311 if( (max1 < fabs(p1_0_p0_0)) )
949 {
950 879 max1 = fabs(p1_0_p0_0);
951 }
952
2/2
✓ Branch 0 taken 72871 times.
✓ Branch 1 taken 196440 times.
 269311 if( (max1 < fabs(p1_6_p0_6)) )
953 {
954 72871 max1 = fabs(p1_6_p0_6);
955 }
956
2/2
✓ Branch 0 taken 112495 times.
✓ Branch 1 taken 156816 times.
 269311 if( (max1 < fabs(p1_2_p0_2)) )
957 {
958 112495 max1 = fabs(p1_2_p0_2);
959 }
960
2/2
✓ Branch 0 taken 77946 times.
✓ Branch 1 taken 191365 times.
 269311 if( (max1 < fabs(p1_1_p0_1)) )
961 {
962 77946 max1 = fabs(p1_1_p0_1);
963 }
964
2/2
✓ Branch 0 taken 12745 times.
✓ Branch 1 taken 256566 times.
 269311 if( (max1 < fabs(p1_5_p0_5)) )
965 {
966 12745 max1 = fabs(p1_5_p0_5);
967 }
968 269311 double max2 = fabs(q0_0_p0_0);
969
2/2
✓ Branch 0 taken 144757 times.
✓ Branch 1 taken 124554 times.
 269311 if( (max2 < fabs(q0_1_p0_1)) )
970 {
971 144757 max2 = fabs(q0_1_p0_1);
972 }
973
2/2
✓ Branch 0 taken 91576 times.
✓ Branch 1 taken 177735 times.
 269311 if( (max2 < fabs(q0_2_p0_2)) )
974 {
975 91576 max2 = fabs(q0_2_p0_2);
976 }
977
2/2
✓ Branch 0 taken 208860 times.
✓ Branch 1 taken 60451 times.
 269311 if( (max2 < fabs(q0_3_p0_3)) )
978 {
979 208860 max2 = fabs(q0_3_p0_3);
980 }
981
2/2
✓ Branch 0 taken 207621 times.
✓ Branch 1 taken 61690 times.
 269311 if( (max2 < fabs(q0_4_p0_4)) )
982 {
983 207621 max2 = fabs(q0_4_p0_4);
984 }
985
2/2
✓ Branch 0 taken 43713 times.
✓ Branch 1 taken 225598 times.
 269311 if( (max2 < fabs(q0_5_p0_5)) )
986 {
987 43713 max2 = fabs(q0_5_p0_5);
988 }
989
2/2
✓ Branch 0 taken 51270 times.
✓ Branch 1 taken 218041 times.
 269311 if( (max2 < fabs(q0_6_p0_6)) )
990 {
991 51270 max2 = fabs(q0_6_p0_6);
992 }
993
2/2
✓ Branch 0 taken 61842 times.
✓ Branch 1 taken 207469 times.
 269311 if( (max2 < fabs(q0_7_p0_7)) )
994 {
995 61842 max2 = fabs(q0_7_p0_7);
996 }
997
2/2
✓ Branch 0 taken 21798 times.
✓ Branch 1 taken 247513 times.
 269311 if( (max2 < fabs(q1_0_p0_0)) )
998 {
999 21798 max2 = fabs(q1_0_p0_0);
1000 }
1001
2/2
✓ Branch 0 taken 5223 times.
✓ Branch 1 taken 264088 times.
 269311 if( (max2 < fabs(q1_1_p0_1)) )
1002 {
1003 5223 max2 = fabs(q1_1_p0_1);
1004 }
1005
2/2
✓ Branch 0 taken 6636 times.
✓ Branch 1 taken 262675 times.
 269311 if( (max2 < fabs(q1_2_p0_2)) )
1006 {
1007 6636 max2 = fabs(q1_2_p0_2);
1008 }
1009
2/2
✓ Branch 0 taken 79623 times.
✓ Branch 1 taken 189688 times.
 269311 if( (max2 < fabs(q1_3_p0_3)) )
1010 {
1011 79623 max2 = fabs(q1_3_p0_3);
1012 }
1013
2/2
✓ Branch 0 taken 89454 times.
✓ Branch 1 taken 179857 times.
 269311 if( (max2 < fabs(q1_4_p0_4)) )
1014 {
1015 89454 max2 = fabs(q1_4_p0_4);
1016 }
1017
2/2
✓ Branch 0 taken 25193 times.
✓ Branch 1 taken 244118 times.
 269311 if( (max2 < fabs(q1_5_p0_5)) )
1018 {
1019 25193 max2 = fabs(q1_5_p0_5);
1020 }
1021
2/2
✓ Branch 0 taken 17399 times.
✓ Branch 1 taken 251912 times.
 269311 if( (max2 < fabs(q1_6_p0_6)) )
1022 {
1023 17399 max2 = fabs(q1_6_p0_6);
1024 }
1025
2/2
✓ Branch 0 taken 34441 times.
✓ Branch 1 taken 234870 times.
 269311 if( (max2 < fabs(q1_7_p0_7)) )
1026 {
1027 34441 max2 = fabs(q1_7_p0_7);
1028 }
1029 double lower_bound_1;
1030 double upper_bound_1;
1031 int Delta_sign;
1032 int int_tmp_result;
1033 269311 lower_bound_1 = max2;
1034 269311 upper_bound_1 = max2;
1035
2/2
✓ Branch 0 taken 258554 times.
✓ Branch 1 taken 10757 times.
 269311 if( (max1 < lower_bound_1) )
1036 {
1037 258554 lower_bound_1 = max1;
1038 }
1039 else
1040 {
1041
1/2
✓ Branch 0 taken 10757 times.
✗ Branch 1 not taken.
 10757 if( (max1 > upper_bound_1) )
1042 {
1043 10757 upper_bound_1 = max1;
1044 }
1045 }
1046
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269311 times.
 269311 if( (lower_bound_1 < 1.15542931091530087067e-147) )
1047 {
1048 ✗ return FPG_UNCERTAIN_VALUE;
1049 }
1050 else
1051 {
1052
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269311 times.
 269311 if( (upper_bound_1 > 1.87072209578355531992e+50) )
1053 {
1054 ✗ return FPG_UNCERTAIN_VALUE;
1055 }
1056 269311 eps = (1.66670090166682227006e-14 * (max1 * max2));
1057
2/2
✓ Branch 0 taken 101649 times.
✓ Branch 1 taken 167662 times.
 269311 if( (Delta > eps) )
1058 {
1059 101649 int_tmp_result = 1;
1060 }
1061 else
1062 {
1063
1/2
✓ Branch 0 taken 167662 times.
✗ Branch 1 not taken.
 167662 if( (Delta < -eps) )
1064 {
1065 167662 int_tmp_result = -1;
1066 }
1067 else
1068 {
1069 ✗ return FPG_UNCERTAIN_VALUE;
1070 }
1071 }
1072 }
1073 269311 Delta_sign = int_tmp_result;
1074 269311 double max3 = max2;
1075
2/2
✓ Branch 0 taken 10757 times.
✓ Branch 1 taken 258554 times.
 269311 if( (max3 < max1) )
1076 {
1077 10757 max3 = max1;
1078 }
1079 269311 double max4 = max2;
1080
2/2
✓ Branch 0 taken 2364 times.
✓ Branch 1 taken 266947 times.
 269311 if( (max4 < fabs(p2_4_p0_4)) )
1081 {
1082 2364 max4 = fabs(p2_4_p0_4);
1083 }
1084
2/2
✓ Branch 0 taken 6060 times.
✓ Branch 1 taken 263251 times.
 269311 if( (max4 < fabs(p2_2_p0_2)) )
1085 {
1086 6060 max4 = fabs(p2_2_p0_2);
1087 }
1088
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269311 times.
 269311 if( (max4 < fabs(p2_0_p0_0)) )
1089 {
1090 ✗ max4 = fabs(p2_0_p0_0);
1091 }
1092
2/2
✓ Branch 0 taken 4919 times.
✓ Branch 1 taken 264392 times.
 269311 if( (max4 < fabs(p2_1_p0_1)) )
1093 {
1094 4919 max4 = fabs(p2_1_p0_1);
1095 }
1096
2/2
✓ Branch 0 taken 381 times.
✓ Branch 1 taken 268930 times.
 269311 if( (max4 < fabs(p2_3_p0_3)) )
1097 {
1098 381 max4 = fabs(p2_3_p0_3);
1099 }
1100
2/2
✓ Branch 0 taken 162 times.
✓ Branch 1 taken 269149 times.
 269311 if( (max4 < fabs(p2_5_p0_5)) )
1101 {
1102 162 max4 = fabs(p2_5_p0_5);
1103 }
1104
2/2
✓ Branch 0 taken 771 times.
✓ Branch 1 taken 268540 times.
 269311 if( (max4 < fabs(p2_6_p0_6)) )
1105 {
1106 771 max4 = fabs(p2_6_p0_6);
1107 }
1108
2/2
✓ Branch 0 taken 397 times.
✓ Branch 1 taken 268914 times.
 269311 if( (max4 < fabs(p2_7_p0_7)) )
1109 {
1110 397 max4 = fabs(p2_7_p0_7);
1111 }
1112
2/2
✓ Branch 0 taken 7571 times.
✓ Branch 1 taken 261740 times.
 269311 if( (max3 < max4) )
1113 {
1114 7571 max3 = max4;
1115 }
1116 int r_sign;
1117 int int_tmp_result_FFWKCAA;
1118 269311 lower_bound_1 = max1;
1119 269311 upper_bound_1 = max1;
1120
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269311 times.
 269311 if( (max3 < lower_bound_1) )
1121 {
1122 ✗ lower_bound_1 = max3;
1123 }
1124 else
1125 {
1126
2/2
✓ Branch 0 taken 264808 times.
✓ Branch 1 taken 4503 times.
 269311 if( (max3 > upper_bound_1) )
1127 {
1128 264808 upper_bound_1 = max3;
1129 }
1130 }
1131
2/2
✓ Branch 0 taken 4503 times.
✓ Branch 1 taken 264808 times.
 269311 if( (max4 < lower_bound_1) )
1132 {
1133 4503 lower_bound_1 = max4;
1134 }
1135
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269311 times.
 269311 if( (lower_bound_1 < 1.26419510663115923609e-74) )
1136 {
1137 ✗ return FPG_UNCERTAIN_VALUE;
1138 }
1139 else
1140 {
1141
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 269311 times.
 269311 if( (upper_bound_1 > 1.87072209578355531992e+50) )
1142 {
1143 ✗ return FPG_UNCERTAIN_VALUE;
1144 }
1145 269311 eps = (8.71140112255785451890e-13 * (((max1 * max4) * max4) * max3));
1146
2/2
✓ Branch 0 taken 110310 times.
✓ Branch 1 taken 159001 times.
 269311 if( (r > eps) )
1147 {
1148 110310 int_tmp_result_FFWKCAA = 1;
1149 }
1150 else
1151 {
1152
1/2
✓ Branch 0 taken 159001 times.
✗ Branch 1 not taken.
 159001 if( (r < -eps) )
1153 {
1154 159001 int_tmp_result_FFWKCAA = -1;
1155 }
1156 else
1157 {
1158 ✗ return FPG_UNCERTAIN_VALUE;
1159 }
1160 }
1161 }
1162 269311 r_sign = int_tmp_result_FFWKCAA;
1163 269311 return (Delta_sign * r_sign);
1164 }
1165