GCC Code Coverage Report

Directory: ./
File: lib/geogram/numerics/predicates/side2.h
Date: 2026-09-07 02:36:43
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 5090168 inline int side2_3d_filter( const double* p0, const double* p1, const double* p2, const double* q0, const double* q1) {
5 5090168 double p1_0_p0_0 = (p1[0] - p0[0]);
6 5090168 double p1_1_p0_1 = (p1[1] - p0[1]);
7 5090168 double p1_2_p0_2 = (p1[2] - p0[2]);
8 double l1;
9 5090168 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 5090168 double p2_0_p0_0 = (p2[0] - p0[0]);
11 5090168 double p2_1_p0_1 = (p2[1] - p0[1]);
12 5090168 double p2_2_p0_2 = (p2[2] - p0[2]);
13 double l2;
14 5090168 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 5090168 double q0_0_p0_0 = (q0[0] - p0[0]);
16 5090168 double q0_1_p0_1 = (q0[1] - p0[1]);
17 5090168 double q0_2_p0_2 = (q0[2] - p0[2]);
18 double a10;
19 5090168 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 5090168 double q1_0_p0_0 = (q1[0] - p0[0]);
21 5090168 double q1_1_p0_1 = (q1[1] - p0[1]);
22 5090168 double q1_2_p0_2 = (q1[2] - p0[2]);
23 double a11;
24 5090168 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 5090168 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 5090168 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 5090168 Delta = (a11 - a10);
31 double DeltaLambda0;
32 5090168 DeltaLambda0 = (a11 - l1);
33 double DeltaLambda1;
34 5090168 DeltaLambda1 = (l1 - a10);
35 double r;
36 5090168 r = (((Delta * l2) - (a20 * DeltaLambda0)) - (a21 * DeltaLambda1));
37 double eps;
38 5090168 double max1 = fabs(p1_0_p0_0);
39
2/2
✓ Branch 0 taken 2518484 times.
✓ Branch 1 taken 2571684 times.
 5090168 if( (max1 < fabs(p1_1_p0_1)) )
40 {
41 2518484 max1 = fabs(p1_1_p0_1);
42 }
43
2/2
✓ Branch 0 taken 1238918 times.
✓ Branch 1 taken 3851250 times.
 5090168 if( (max1 < fabs(p1_2_p0_2)) )
44 {
45 1238918 max1 = fabs(p1_2_p0_2);
46 }
47 5090168 double max2 = fabs(q0_0_p0_0);
48
2/2
✓ Branch 0 taken 1683561 times.
✓ Branch 1 taken 3406607 times.
 5090168 if( (max2 < fabs(q0_1_p0_1)) )
49 {
50 1683561 max2 = fabs(q0_1_p0_1);
51 }
52
2/2
✓ Branch 0 taken 787550 times.
✓ Branch 1 taken 4302618 times.
 5090168 if( (max2 < fabs(q0_2_p0_2)) )
53 {
54 787550 max2 = fabs(q0_2_p0_2);
55 }
56
2/2
✓ Branch 0 taken 1615483 times.
✓ Branch 1 taken 3474685 times.
 5090168 if( (max2 < fabs(q1_0_p0_0)) )
57 {
58 1615483 max2 = fabs(q1_0_p0_0);
59 }
60
2/2
✓ Branch 0 taken 1181794 times.
✓ Branch 1 taken 3908374 times.
 5090168 if( (max2 < fabs(q1_1_p0_1)) )
61 {
62 1181794 max2 = fabs(q1_1_p0_1);
63 }
64
2/2
✓ Branch 0 taken 583735 times.
✓ Branch 1 taken 4506433 times.
 5090168 if( (max2 < fabs(q1_2_p0_2)) )
65 {
66 583735 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 5090168 lower_bound_1 = max1;
73 5090168 upper_bound_1 = max1;
74
2/2
✓ Branch 0 taken 342137 times.
✓ Branch 1 taken 4748031 times.
 5090168 if( (max2 < lower_bound_1) )
75 {
76 342137 lower_bound_1 = max2;
77 }
78 else
79 {
80
2/2
✓ Branch 0 taken 2861698 times.
✓ Branch 1 taken 1886333 times.
 4748031 if( (max2 > upper_bound_1) )
81 {
82 2861698 upper_bound_1 = max2;
83 }
84 }
85
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5090168 times.
 5090168 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 5090168 times.
 5090168 if( (upper_bound_1 > 3.74144419156711063983e+50) )
92 {
93 ✗ return FPG_UNCERTAIN_VALUE;
94 }
95 5090168 eps = (4.44425370757048798480e-15 * (max1 * max2));
96
2/2
✓ Branch 0 taken 2995970 times.
✓ Branch 1 taken 2094198 times.
 5090168 if( (Delta > eps) )
97 {
98 2995970 int_tmp_result = 1;
99 }
100 else
101 {
102
2/2
✓ Branch 0 taken 2093756 times.
✓ Branch 1 taken 442 times.
 2094198 if( (Delta < -eps) )
103 {
104 2093756 int_tmp_result = -1;
105 }
106 else
107 {
108 442 return FPG_UNCERTAIN_VALUE;
109 }
110 }
111 }
112 5089726 Delta_sign = int_tmp_result;
113 5089726 double max3 = max1;
114
2/2
✓ Branch 0 taken 2861256 times.
✓ Branch 1 taken 2228470 times.
 5089726 if( (max3 < max2) )
115 {
116 2861256 max3 = max2;
117 }
118 5089726 double max4 = max2;
119
2/2
✓ Branch 0 taken 772985 times.
✓ Branch 1 taken 4316741 times.
 5089726 if( (max4 < fabs(p2_0_p0_0)) )
120 {
121 772985 max4 = fabs(p2_0_p0_0);
122 }
123
2/2
✓ Branch 0 taken 459873 times.
✓ Branch 1 taken 4629853 times.
 5089726 if( (max4 < fabs(p2_1_p0_1)) )
124 {
125 459873 max4 = fabs(p2_1_p0_1);
126 }
127
2/2
✓ Branch 0 taken 403550 times.
✓ Branch 1 taken 4686176 times.
 5089726 if( (max4 < fabs(p2_2_p0_2)) )
128 {
129 403550 max4 = fabs(p2_2_p0_2);
130 }
131
2/2
✓ Branch 0 taken 1149262 times.
✓ Branch 1 taken 3940464 times.
 5089726 if( (max3 < max4) )
132 {
133 1149262 max3 = max4;
134 }
135 int r_sign;
136 int int_tmp_result_FFWKCAA;
137 5089726 lower_bound_1 = max1;
138 5089726 upper_bound_1 = max1;
139
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5089726 times.
 5089726 if( (max3 < lower_bound_1) )
140 {
141 ✗ lower_bound_1 = max3;
142 }
143 else
144 {
145
2/2
✓ Branch 0 taken 3740763 times.
✓ Branch 1 taken 1348963 times.
 5089726 if( (max3 > upper_bound_1) )
146 {
147 3740763 upper_bound_1 = max3;
148 }
149 }
150
2/2
✓ Branch 0 taken 49723 times.
✓ Branch 1 taken 5040003 times.
 5089726 if( (max4 < lower_bound_1) )
151 {
152 49723 lower_bound_1 = max4;
153 }
154
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 5089726 times.
 5089726 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 5089726 times.
 5089726 if( (upper_bound_1 > 3.74144419156711063983e+50) )
161 {
162 ✗ return FPG_UNCERTAIN_VALUE;
163 }
164 5089726 eps = (8.99983341597279045654e-14 * (((max1 * max4) * max4) * max3));
165
2/2
✓ Branch 0 taken 2783207 times.
✓ Branch 1 taken 2306519 times.
 5089726 if( (r > eps) )
166 {
167 2783207 int_tmp_result_FFWKCAA = 1;
168 }
169 else
170 {
171
2/2
✓ Branch 0 taken 2027441 times.
✓ Branch 1 taken 279078 times.
 2306519 if( (r < -eps) )
172 {
173 2027441 int_tmp_result_FFWKCAA = -1;
174 }
175 else
176 {
177 279078 return FPG_UNCERTAIN_VALUE;
178 }
179 }
180 }
181 4810648 r_sign = int_tmp_result_FFWKCAA;
182 4810648 return (Delta_sign * r_sign);
183 }
184
185
186 189360 inline int side2_4d_filter( const double* p0, const double* p1, const double* p2, const double* q0, const double* q1) {
187 189360 double p1_0_p0_0 = (p1[0] - p0[0]);
188 189360 double p1_1_p0_1 = (p1[1] - p0[1]);
189 189360 double p1_2_p0_2 = (p1[2] - p0[2]);
190 189360 double p1_3_p0_3 = (p1[3] - p0[3]);
191 double l1;
192 189360 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 189360 double p2_0_p0_0 = (p2[0] - p0[0]);
194 189360 double p2_1_p0_1 = (p2[1] - p0[1]);
195 189360 double p2_2_p0_2 = (p2[2] - p0[2]);
196 189360 double p2_3_p0_3 = (p2[3] - p0[3]);
197 double l2;
198 189360 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 189360 double q0_0_p0_0 = (q0[0] - p0[0]);
200 189360 double q0_1_p0_1 = (q0[1] - p0[1]);
201 189360 double q0_2_p0_2 = (q0[2] - p0[2]);
202 189360 double q0_3_p0_3 = (q0[3] - p0[3]);
203 double a10;
204 189360 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 189360 double q1_0_p0_0 = (q1[0] - p0[0]);
206 189360 double q1_1_p0_1 = (q1[1] - p0[1]);
207 189360 double q1_2_p0_2 = (q1[2] - p0[2]);
208 189360 double q1_3_p0_3 = (q1[3] - p0[3]);
209 double a11;
210 189360 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 189360 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 189360 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 189360 Delta = (a11 - a10);
217 double DeltaLambda0;
218 189360 DeltaLambda0 = (a11 - l1);
219 double DeltaLambda1;
220 189360 DeltaLambda1 = (l1 - a10);
221 double r;
222 189360 r = (((Delta * l2) - (a20 * DeltaLambda0)) - (a21 * DeltaLambda1));
223 double eps;
224 189360 double max1 = fabs(p1_2_p0_2);
225
2/2
✓ Branch 0 taken 79455 times.
✓ Branch 1 taken 109905 times.
 189360 if( (max1 < fabs(p1_0_p0_0)) )
226 {
227 79455 max1 = fabs(p1_0_p0_0);
228 }
229
2/2
✓ Branch 0 taken 74903 times.
✓ Branch 1 taken 114457 times.
 189360 if( (max1 < fabs(p1_1_p0_1)) )
230 {
231 74903 max1 = fabs(p1_1_p0_1);
232 }
233
2/2
✓ Branch 0 taken 85809 times.
✓ Branch 1 taken 103551 times.
 189360 if( (max1 < fabs(p1_3_p0_3)) )
234 {
235 85809 max1 = fabs(p1_3_p0_3);
236 }
237 189360 double max2 = fabs(q0_0_p0_0);
238
2/2
✓ Branch 0 taken 94795 times.
✓ Branch 1 taken 94565 times.
 189360 if( (max2 < fabs(q0_1_p0_1)) )
239 {
240 94795 max2 = fabs(q0_1_p0_1);
241 }
242
2/2
✓ Branch 0 taken 66115 times.
✓ Branch 1 taken 123245 times.
 189360 if( (max2 < fabs(q0_2_p0_2)) )
243 {
244 66115 max2 = fabs(q0_2_p0_2);
245 }
246
2/2
✓ Branch 0 taken 150401 times.
✓ Branch 1 taken 38959 times.
 189360 if( (max2 < fabs(q0_3_p0_3)) )
247 {
248 150401 max2 = fabs(q0_3_p0_3);
249 }
250
2/2
✓ Branch 0 taken 15590 times.
✓ Branch 1 taken 173770 times.
 189360 if( (max2 < fabs(q1_0_p0_0)) )
251 {
252 15590 max2 = fabs(q1_0_p0_0);
253 }
254
2/2
✓ Branch 0 taken 7411 times.
✓ Branch 1 taken 181949 times.
 189360 if( (max2 < fabs(q1_1_p0_1)) )
255 {
256 7411 max2 = fabs(q1_1_p0_1);
257 }
258
2/2
✓ Branch 0 taken 5456 times.
✓ Branch 1 taken 183904 times.
 189360 if( (max2 < fabs(q1_2_p0_2)) )
259 {
260 5456 max2 = fabs(q1_2_p0_2);
261 }
262
2/2
✓ Branch 0 taken 73663 times.
✓ Branch 1 taken 115697 times.
 189360 if( (max2 < fabs(q1_3_p0_3)) )
263 {
264 73663 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 189360 lower_bound_1 = max1;
271 189360 upper_bound_1 = max1;
272
2/2
✓ Branch 0 taken 8528 times.
✓ Branch 1 taken 180832 times.
 189360 if( (max2 < lower_bound_1) )
273 {
274 8528 lower_bound_1 = max2;
275 }
276 else
277 {
278
1/2
✓ Branch 0 taken 180832 times.
✗ Branch 1 not taken.
 180832 if( (max2 > upper_bound_1) )
279 {
280 180832 upper_bound_1 = max2;
281 }
282 }
283
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 189360 times.
 189360 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 189360 times.
 189360 if( (upper_bound_1 > 1.87072209578355531992e+50) )
290 {
291 ✗ return FPG_UNCERTAIN_VALUE;
292 }
293 189360 eps = (6.44428177279185717888e-15 * (max1 * max2));
294
2/2
✓ Branch 0 taken 61511 times.
✓ Branch 1 taken 127849 times.
 189360 if( (Delta > eps) )
295 {
296 61511 int_tmp_result = 1;
297 }
298 else
299 {
300
1/2
✓ Branch 0 taken 127849 times.
✗ Branch 1 not taken.
 127849 if( (Delta < -eps) )
301 {
302 127849 int_tmp_result = -1;
303 }
304 else
305 {
306 ✗ return FPG_UNCERTAIN_VALUE;
307 }
308 }
309 }
310 189360 Delta_sign = int_tmp_result;
311 189360 double max3 = max1;
312
2/2
✓ Branch 0 taken 180832 times.
✓ Branch 1 taken 8528 times.
 189360 if( (max3 < max2) )
313 {
314 180832 max3 = max2;
315 }
316 189360 double max4 = max2;
317
2/2
✓ Branch 0 taken 167 times.
✓ Branch 1 taken 189193 times.
 189360 if( (max4 < fabs(p2_0_p0_0)) )
318 {
319 167 max4 = fabs(p2_0_p0_0);
320 }
321
2/2
✓ Branch 0 taken 5962 times.
✓ Branch 1 taken 183398 times.
 189360 if( (max4 < fabs(p2_2_p0_2)) )
322 {
323 5962 max4 = fabs(p2_2_p0_2);
324 }
325
2/2
✓ Branch 0 taken 5440 times.
✓ Branch 1 taken 183920 times.
 189360 if( (max4 < fabs(p2_1_p0_1)) )
326 {
327 5440 max4 = fabs(p2_1_p0_1);
328 }
329
2/2
✓ Branch 0 taken 3983 times.
✓ Branch 1 taken 185377 times.
 189360 if( (max4 < fabs(p2_3_p0_3)) )
330 {
331 3983 max4 = fabs(p2_3_p0_3);
332 }
333
2/2
✓ Branch 0 taken 12188 times.
✓ Branch 1 taken 177172 times.
 189360 if( (max3 < max4) )
334 {
335 12188 max3 = max4;
336 }
337 int r_sign;
338 int int_tmp_result_FFWKCAA;
339 189360 lower_bound_1 = max3;
340 189360 upper_bound_1 = max3;
341
2/2
✓ Branch 0 taken 187009 times.
✓ Branch 1 taken 2351 times.
 189360 if( (max1 < lower_bound_1) )
342 {
343 187009 lower_bound_1 = max1;
344 }
345
2/2
✓ Branch 0 taken 2351 times.
✓ Branch 1 taken 187009 times.
 189360 if( (max4 < lower_bound_1) )
346 {
347 2351 lower_bound_1 = max4;
348 }
349
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 189360 times.
 189360 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 189360 times.
 189360 if( (upper_bound_1 > 1.87072209578355531992e+50) )
356 {
357 ✗ return FPG_UNCERTAIN_VALUE;
358 }
359 189360 eps = (1.72443682410932010423e-13 * (((max1 * max4) * max4) * max3));
360
2/2
✓ Branch 0 taken 72432 times.
✓ Branch 1 taken 116928 times.
 189360 if( (r > eps) )
361 {
362 72432 int_tmp_result_FFWKCAA = 1;
363 }
364 else
365 {
366
1/2
✓ Branch 0 taken 116928 times.
✗ Branch 1 not taken.
 116928 if( (r < -eps) )
367 {
368 116928 int_tmp_result_FFWKCAA = -1;
369 }
370 else
371 {
372 ✗ return FPG_UNCERTAIN_VALUE;
373 }
374 }
375 }
376 189360 r_sign = int_tmp_result_FFWKCAA;
377 189360 return (Delta_sign * r_sign);
378 }
379
380
381 260704 inline int side2_6d_filter( const double* p0, const double* p1, const double* p2, const double* q0, const double* q1) {
382 260704 double p1_0_p0_0 = (p1[0] - p0[0]);
383 260704 double p1_1_p0_1 = (p1[1] - p0[1]);
384 260704 double p1_2_p0_2 = (p1[2] - p0[2]);
385 260704 double p1_3_p0_3 = (p1[3] - p0[3]);
386 260704 double p1_4_p0_4 = (p1[4] - p0[4]);
387 260704 double p1_5_p0_5 = (p1[5] - p0[5]);
388 double l1;
389 260704 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 260704 double p2_0_p0_0 = (p2[0] - p0[0]);
391 260704 double p2_1_p0_1 = (p2[1] - p0[1]);
392 260704 double p2_2_p0_2 = (p2[2] - p0[2]);
393 260704 double p2_3_p0_3 = (p2[3] - p0[3]);
394 260704 double p2_4_p0_4 = (p2[4] - p0[4]);
395 260704 double p2_5_p0_5 = (p2[5] - p0[5]);
396 double l2;
397 260704 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 260704 double q0_0_p0_0 = (q0[0] - p0[0]);
399 260704 double q0_1_p0_1 = (q0[1] - p0[1]);
400 260704 double q0_2_p0_2 = (q0[2] - p0[2]);
401 260704 double q0_3_p0_3 = (q0[3] - p0[3]);
402 260704 double q0_4_p0_4 = (q0[4] - p0[4]);
403 260704 double q0_5_p0_5 = (q0[5] - p0[5]);
404 double a10;
405 260704 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 260704 double q1_0_p0_0 = (q1[0] - p0[0]);
407 260704 double q1_1_p0_1 = (q1[1] - p0[1]);
408 260704 double q1_2_p0_2 = (q1[2] - p0[2]);
409 260704 double q1_3_p0_3 = (q1[3] - p0[3]);
410 260704 double q1_4_p0_4 = (q1[4] - p0[4]);
411 260704 double q1_5_p0_5 = (q1[5] - p0[5]);
412 double a11;
413 260704 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 260704 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 260704 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 260704 Delta = (a11 - a10);
420 double DeltaLambda0;
421 260704 DeltaLambda0 = (a11 - l1);
422 double DeltaLambda1;
423 260704 DeltaLambda1 = (l1 - a10);
424 double r;
425 260704 r = (((Delta * l2) - (a20 * DeltaLambda0)) - (a21 * DeltaLambda1));
426 double eps;
427 260704 double max1 = fabs(p1_1_p0_1);
428
2/2
✓ Branch 0 taken 133392 times.
✓ Branch 1 taken 127312 times.
 260704 if( (max1 < fabs(p1_2_p0_2)) )
429 {
430 133392 max1 = fabs(p1_2_p0_2);
431 }
432
2/2
✓ Branch 0 taken 74682 times.
✓ Branch 1 taken 186022 times.
 260704 if( (max1 < fabs(p1_3_p0_3)) )
433 {
434 74682 max1 = fabs(p1_3_p0_3);
435 }
436
2/2
✓ Branch 0 taken 1913 times.
✓ Branch 1 taken 258791 times.
 260704 if( (max1 < fabs(p1_0_p0_0)) )
437 {
438 1913 max1 = fabs(p1_0_p0_0);
439 }
440
2/2
✓ Branch 0 taken 71674 times.
✓ Branch 1 taken 189030 times.
 260704 if( (max1 < fabs(p1_4_p0_4)) )
441 {
442 71674 max1 = fabs(p1_4_p0_4);
443 }
444
2/2
✓ Branch 0 taken 16895 times.
✓ Branch 1 taken 243809 times.
 260704 if( (max1 < fabs(p1_5_p0_5)) )
445 {
446 16895 max1 = fabs(p1_5_p0_5);
447 }
448 260704 double max2 = fabs(q0_0_p0_0);
449
2/2
✓ Branch 0 taken 144383 times.
✓ Branch 1 taken 116321 times.
 260704 if( (max2 < fabs(q0_1_p0_1)) )
450 {
451 144383 max2 = fabs(q0_1_p0_1);
452 }
453
2/2
✓ Branch 0 taken 95125 times.
✓ Branch 1 taken 165579 times.
 260704 if( (max2 < fabs(q0_2_p0_2)) )
454 {
455 95125 max2 = fabs(q0_2_p0_2);
456 }
457
2/2
✓ Branch 0 taken 207927 times.
✓ Branch 1 taken 52777 times.
 260704 if( (max2 < fabs(q0_3_p0_3)) )
458 {
459 207927 max2 = fabs(q0_3_p0_3);
460 }
461
2/2
✓ Branch 0 taken 199175 times.
✓ Branch 1 taken 61529 times.
 260704 if( (max2 < fabs(q0_4_p0_4)) )
462 {
463 199175 max2 = fabs(q0_4_p0_4);
464 }
465
2/2
✓ Branch 0 taken 42533 times.
✓ Branch 1 taken 218171 times.
 260704 if( (max2 < fabs(q0_5_p0_5)) )
466 {
467 42533 max2 = fabs(q0_5_p0_5);
468 }
469
2/2
✓ Branch 0 taken 25200 times.
✓ Branch 1 taken 235504 times.
 260704 if( (max2 < fabs(q1_0_p0_0)) )
470 {
471 25200 max2 = fabs(q1_0_p0_0);
472 }
473
2/2
✓ Branch 0 taken 6437 times.
✓ Branch 1 taken 254267 times.
 260704 if( (max2 < fabs(q1_1_p0_1)) )
474 {
475 6437 max2 = fabs(q1_1_p0_1);
476 }
477
2/2
✓ Branch 0 taken 5445 times.
✓ Branch 1 taken 255259 times.
 260704 if( (max2 < fabs(q1_2_p0_2)) )
478 {
479 5445 max2 = fabs(q1_2_p0_2);
480 }
481
2/2
✓ Branch 0 taken 87353 times.
✓ Branch 1 taken 173351 times.
 260704 if( (max2 < fabs(q1_3_p0_3)) )
482 {
483 87353 max2 = fabs(q1_3_p0_3);
484 }
485
2/2
✓ Branch 0 taken 103690 times.
✓ Branch 1 taken 157014 times.
 260704 if( (max2 < fabs(q1_4_p0_4)) )
486 {
487 103690 max2 = fabs(q1_4_p0_4);
488 }
489
2/2
✓ Branch 0 taken 28959 times.
✓ Branch 1 taken 231745 times.
 260704 if( (max2 < fabs(q1_5_p0_5)) )
490 {
491 28959 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 260704 lower_bound_1 = max1;
498 260704 upper_bound_1 = max1;
499
2/2
✓ Branch 0 taken 8933 times.
✓ Branch 1 taken 251771 times.
 260704 if( (max2 < lower_bound_1) )
500 {
501 8933 lower_bound_1 = max2;
502 }
503 else
504 {
505
1/2
✓ Branch 0 taken 251771 times.
✗ Branch 1 not taken.
 251771 if( (max2 > upper_bound_1) )
506 {
507 251771 upper_bound_1 = max2;
508 }
509 }
510
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 260704 times.
 260704 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 260704 times.
 260704 if( (upper_bound_1 > 1.87072209578355531992e+50) )
517 {
518 ✗ return FPG_UNCERTAIN_VALUE;
519 }
520 260704 eps = (1.11111223981318615596e-14 * (max1 * max2));
521
2/2
✓ Branch 0 taken 96236 times.
✓ Branch 1 taken 164468 times.
 260704 if( (Delta > eps) )
522 {
523 96236 int_tmp_result = 1;
524 }
525 else
526 {
527
1/2
✓ Branch 0 taken 164468 times.
✗ Branch 1 not taken.
 164468 if( (Delta < -eps) )
528 {
529 164468 int_tmp_result = -1;
530 }
531 else
532 {
533 ✗ return FPG_UNCERTAIN_VALUE;
534 }
535 }
536 }
537 260704 Delta_sign = int_tmp_result;
538 260704 double max3 = max1;
539
2/2
✓ Branch 0 taken 251771 times.
✓ Branch 1 taken 8933 times.
 260704 if( (max3 < max2) )
540 {
541 251771 max3 = max2;
542 }
543 260704 double max4 = max2;
544
2/2
✓ Branch 0 taken 5356 times.
✓ Branch 1 taken 255348 times.
 260704 if( (max4 < fabs(p2_2_p0_2)) )
545 {
546 5356 max4 = fabs(p2_2_p0_2);
547 }
548
2/2
✓ Branch 0 taken 2788 times.
✓ Branch 1 taken 257916 times.
 260704 if( (max4 < fabs(p2_4_p0_4)) )
549 {
550 2788 max4 = fabs(p2_4_p0_4);
551 }
552
2/2
✓ Branch 0 taken 953 times.
✓ Branch 1 taken 259751 times.
 260704 if( (max4 < fabs(p2_3_p0_3)) )
553 {
554 953 max4 = fabs(p2_3_p0_3);
555 }
556
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 260704 times.
 260704 if( (max4 < fabs(p2_0_p0_0)) )
557 {
558 ✗ max4 = fabs(p2_0_p0_0);
559 }
560
2/2
✓ Branch 0 taken 3929 times.
✓ Branch 1 taken 256775 times.
 260704 if( (max4 < fabs(p2_1_p0_1)) )
561 {
562 3929 max4 = fabs(p2_1_p0_1);
563 }
564
2/2
✓ Branch 0 taken 197 times.
✓ Branch 1 taken 260507 times.
 260704 if( (max4 < fabs(p2_5_p0_5)) )
565 {
566 197 max4 = fabs(p2_5_p0_5);
567 }
568
2/2
✓ Branch 0 taken 7639 times.
✓ Branch 1 taken 253065 times.
 260704 if( (max3 < max4) )
569 {
570 7639 max3 = max4;
571 }
572 int r_sign;
573 int int_tmp_result_FFWKCAA;
574 260704 lower_bound_1 = max1;
575 260704 upper_bound_1 = max1;
576
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 260704 times.
 260704 if( (max3 < lower_bound_1) )
577 {
578 ✗ lower_bound_1 = max3;
579 }
580 else
581 {
582
2/2
✓ Branch 0 taken 257530 times.
✓ Branch 1 taken 3174 times.
 260704 if( (max3 > upper_bound_1) )
583 {
584 257530 upper_bound_1 = max3;
585 }
586 }
587
2/2
✓ Branch 0 taken 3174 times.
✓ Branch 1 taken 257530 times.
 260704 if( (max4 < lower_bound_1) )
588 {
589 3174 lower_bound_1 = max4;
590 }
591
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 260704 times.
 260704 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 260704 times.
 260704 if( (upper_bound_1 > 1.87072209578355531992e+50) )
598 {
599 ✗ return FPG_UNCERTAIN_VALUE;
600 }
601 260704 eps = (4.40007476026584016994e-13 * (((max1 * max4) * max4) * max3));
602
2/2
✓ Branch 0 taken 106023 times.
✓ Branch 1 taken 154681 times.
 260704 if( (r > eps) )
603 {
604 106023 int_tmp_result_FFWKCAA = 1;
605 }
606 else
607 {
608
1/2
✓ Branch 0 taken 154681 times.
✗ Branch 1 not taken.
 154681 if( (r < -eps) )
609 {
610 154681 int_tmp_result_FFWKCAA = -1;
611 }
612 else
613 {
614 ✗ return FPG_UNCERTAIN_VALUE;
615 }
616 }
617 }
618 260704 r_sign = int_tmp_result_FFWKCAA;
619 260704 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 258047 inline int side2_8d_filter( const double* p0, const double* p1, const double* p2, const double* q0, const double* q1) {
886 258047 double p1_0_p0_0 = (p1[0] - p0[0]);
887 258047 double p1_1_p0_1 = (p1[1] - p0[1]);
888 258047 double p1_2_p0_2 = (p1[2] - p0[2]);
889 258047 double p1_3_p0_3 = (p1[3] - p0[3]);
890 258047 double p1_4_p0_4 = (p1[4] - p0[4]);
891 258047 double p1_5_p0_5 = (p1[5] - p0[5]);
892 258047 double p1_6_p0_6 = (p1[6] - p0[6]);
893 258047 double p1_7_p0_7 = (p1[7] - p0[7]);
894 double l1;
895 258047 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 258047 double p2_0_p0_0 = (p2[0] - p0[0]);
897 258047 double p2_1_p0_1 = (p2[1] - p0[1]);
898 258047 double p2_2_p0_2 = (p2[2] - p0[2]);
899 258047 double p2_3_p0_3 = (p2[3] - p0[3]);
900 258047 double p2_4_p0_4 = (p2[4] - p0[4]);
901 258047 double p2_5_p0_5 = (p2[5] - p0[5]);
902 258047 double p2_6_p0_6 = (p2[6] - p0[6]);
903 258047 double p2_7_p0_7 = (p2[7] - p0[7]);
904 double l2;
905 258047 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 258047 double q0_0_p0_0 = (q0[0] - p0[0]);
907 258047 double q0_1_p0_1 = (q0[1] - p0[1]);
908 258047 double q0_2_p0_2 = (q0[2] - p0[2]);
909 258047 double q0_3_p0_3 = (q0[3] - p0[3]);
910 258047 double q0_4_p0_4 = (q0[4] - p0[4]);
911 258047 double q0_5_p0_5 = (q0[5] - p0[5]);
912 258047 double q0_6_p0_6 = (q0[6] - p0[6]);
913 258047 double q0_7_p0_7 = (q0[7] - p0[7]);
914 double a10;
915 258047 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 258047 double q1_0_p0_0 = (q1[0] - p0[0]);
917 258047 double q1_1_p0_1 = (q1[1] - p0[1]);
918 258047 double q1_2_p0_2 = (q1[2] - p0[2]);
919 258047 double q1_3_p0_3 = (q1[3] - p0[3]);
920 258047 double q1_4_p0_4 = (q1[4] - p0[4]);
921 258047 double q1_5_p0_5 = (q1[5] - p0[5]);
922 258047 double q1_6_p0_6 = (q1[6] - p0[6]);
923 258047 double q1_7_p0_7 = (q1[7] - p0[7]);
924 double a11;
925 258047 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 258047 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 258047 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 258047 Delta = (a11 - a10);
932 double DeltaLambda0;
933 258047 DeltaLambda0 = (a11 - l1);
934 double DeltaLambda1;
935 258047 DeltaLambda1 = (l1 - a10);
936 double r;
937 258047 r = (((Delta * l2) - (a20 * DeltaLambda0)) - (a21 * DeltaLambda1));
938 double eps;
939 258047 double max1 = fabs(p1_4_p0_4);
940
2/2
✓ Branch 0 taken 79299 times.
✓ Branch 1 taken 178748 times.
 258047 if( (max1 < fabs(p1_3_p0_3)) )
941 {
942 79299 max1 = fabs(p1_3_p0_3);
943 }
944
2/2
✓ Branch 0 taken 87483 times.
✓ Branch 1 taken 170564 times.
 258047 if( (max1 < fabs(p1_7_p0_7)) )
945 {
946 87483 max1 = fabs(p1_7_p0_7);
947 }
948
2/2
✓ Branch 0 taken 2520 times.
✓ Branch 1 taken 255527 times.
 258047 if( (max1 < fabs(p1_0_p0_0)) )
949 {
950 2520 max1 = fabs(p1_0_p0_0);
951 }
952
2/2
✓ Branch 0 taken 73045 times.
✓ Branch 1 taken 185002 times.
 258047 if( (max1 < fabs(p1_6_p0_6)) )
953 {
954 73045 max1 = fabs(p1_6_p0_6);
955 }
956
2/2
✓ Branch 0 taken 103511 times.
✓ Branch 1 taken 154536 times.
 258047 if( (max1 < fabs(p1_2_p0_2)) )
957 {
958 103511 max1 = fabs(p1_2_p0_2);
959 }
960
2/2
✓ Branch 0 taken 63862 times.
✓ Branch 1 taken 194185 times.
 258047 if( (max1 < fabs(p1_1_p0_1)) )
961 {
962 63862 max1 = fabs(p1_1_p0_1);
963 }
964
2/2
✓ Branch 0 taken 16631 times.
✓ Branch 1 taken 241416 times.
 258047 if( (max1 < fabs(p1_5_p0_5)) )
965 {
966 16631 max1 = fabs(p1_5_p0_5);
967 }
968 258047 double max2 = fabs(q0_0_p0_0);
969
2/2
✓ Branch 0 taken 136251 times.
✓ Branch 1 taken 121796 times.
 258047 if( (max2 < fabs(q0_1_p0_1)) )
970 {
971 136251 max2 = fabs(q0_1_p0_1);
972 }
973
2/2
✓ Branch 0 taken 92756 times.
✓ Branch 1 taken 165291 times.
 258047 if( (max2 < fabs(q0_2_p0_2)) )
974 {
975 92756 max2 = fabs(q0_2_p0_2);
976 }
977
2/2
✓ Branch 0 taken 202366 times.
✓ Branch 1 taken 55681 times.
 258047 if( (max2 < fabs(q0_3_p0_3)) )
978 {
979 202366 max2 = fabs(q0_3_p0_3);
980 }
981
2/2
✓ Branch 0 taken 200969 times.
✓ Branch 1 taken 57078 times.
 258047 if( (max2 < fabs(q0_4_p0_4)) )
982 {
983 200969 max2 = fabs(q0_4_p0_4);
984 }
985
2/2
✓ Branch 0 taken 37496 times.
✓ Branch 1 taken 220551 times.
 258047 if( (max2 < fabs(q0_5_p0_5)) )
986 {
987 37496 max2 = fabs(q0_5_p0_5);
988 }
989
2/2
✓ Branch 0 taken 49869 times.
✓ Branch 1 taken 208178 times.
 258047 if( (max2 < fabs(q0_6_p0_6)) )
990 {
991 49869 max2 = fabs(q0_6_p0_6);
992 }
993
2/2
✓ Branch 0 taken 61768 times.
✓ Branch 1 taken 196279 times.
 258047 if( (max2 < fabs(q0_7_p0_7)) )
994 {
995 61768 max2 = fabs(q0_7_p0_7);
996 }
997
2/2
✓ Branch 0 taken 19514 times.
✓ Branch 1 taken 238533 times.
 258047 if( (max2 < fabs(q1_0_p0_0)) )
998 {
999 19514 max2 = fabs(q1_0_p0_0);
1000 }
1001
2/2
✓ Branch 0 taken 4195 times.
✓ Branch 1 taken 253852 times.
 258047 if( (max2 < fabs(q1_1_p0_1)) )
1002 {
1003 4195 max2 = fabs(q1_1_p0_1);
1004 }
1005
2/2
✓ Branch 0 taken 5609 times.
✓ Branch 1 taken 252438 times.
 258047 if( (max2 < fabs(q1_2_p0_2)) )
1006 {
1007 5609 max2 = fabs(q1_2_p0_2);
1008 }
1009
2/2
✓ Branch 0 taken 78928 times.
✓ Branch 1 taken 179119 times.
 258047 if( (max2 < fabs(q1_3_p0_3)) )
1010 {
1011 78928 max2 = fabs(q1_3_p0_3);
1012 }
1013
2/2
✓ Branch 0 taken 92985 times.
✓ Branch 1 taken 165062 times.
 258047 if( (max2 < fabs(q1_4_p0_4)) )
1014 {
1015 92985 max2 = fabs(q1_4_p0_4);
1016 }
1017
2/2
✓ Branch 0 taken 28037 times.
✓ Branch 1 taken 230010 times.
 258047 if( (max2 < fabs(q1_5_p0_5)) )
1018 {
1019 28037 max2 = fabs(q1_5_p0_5);
1020 }
1021
2/2
✓ Branch 0 taken 26844 times.
✓ Branch 1 taken 231203 times.
 258047 if( (max2 < fabs(q1_6_p0_6)) )
1022 {
1023 26844 max2 = fabs(q1_6_p0_6);
1024 }
1025
2/2
✓ Branch 0 taken 30698 times.
✓ Branch 1 taken 227349 times.
 258047 if( (max2 < fabs(q1_7_p0_7)) )
1026 {
1027 30698 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 258047 lower_bound_1 = max2;
1034 258047 upper_bound_1 = max2;
1035
2/2
✓ Branch 0 taken 247995 times.
✓ Branch 1 taken 10052 times.
 258047 if( (max1 < lower_bound_1) )
1036 {
1037 247995 lower_bound_1 = max1;
1038 }
1039 else
1040 {
1041
1/2
✓ Branch 0 taken 10052 times.
✗ Branch 1 not taken.
 10052 if( (max1 > upper_bound_1) )
1042 {
1043 10052 upper_bound_1 = max1;
1044 }
1045 }
1046
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 258047 times.
 258047 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 258047 times.
 258047 if( (upper_bound_1 > 1.87072209578355531992e+50) )
1053 {
1054 ✗ return FPG_UNCERTAIN_VALUE;
1055 }
1056 258047 eps = (1.66670090166682227006e-14 * (max1 * max2));
1057
2/2
✓ Branch 0 taken 98816 times.
✓ Branch 1 taken 159231 times.
 258047 if( (Delta > eps) )
1058 {
1059 98816 int_tmp_result = 1;
1060 }
1061 else
1062 {
1063
1/2
✓ Branch 0 taken 159231 times.
✗ Branch 1 not taken.
 159231 if( (Delta < -eps) )
1064 {
1065 159231 int_tmp_result = -1;
1066 }
1067 else
1068 {
1069 ✗ return FPG_UNCERTAIN_VALUE;
1070 }
1071 }
1072 }
1073 258047 Delta_sign = int_tmp_result;
1074 258047 double max3 = max2;
1075
2/2
✓ Branch 0 taken 10052 times.
✓ Branch 1 taken 247995 times.
 258047 if( (max3 < max1) )
1076 {
1077 10052 max3 = max1;
1078 }
1079 258047 double max4 = max2;
1080
2/2
✓ Branch 0 taken 2125 times.
✓ Branch 1 taken 255922 times.
 258047 if( (max4 < fabs(p2_4_p0_4)) )
1081 {
1082 2125 max4 = fabs(p2_4_p0_4);
1083 }
1084
2/2
✓ Branch 0 taken 5444 times.
✓ Branch 1 taken 252603 times.
 258047 if( (max4 < fabs(p2_2_p0_2)) )
1085 {
1086 5444 max4 = fabs(p2_2_p0_2);
1087 }
1088
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 258047 times.
 258047 if( (max4 < fabs(p2_0_p0_0)) )
1089 {
1090 ✗ max4 = fabs(p2_0_p0_0);
1091 }
1092
2/2
✓ Branch 0 taken 4071 times.
✓ Branch 1 taken 253976 times.
 258047 if( (max4 < fabs(p2_1_p0_1)) )
1093 {
1094 4071 max4 = fabs(p2_1_p0_1);
1095 }
1096
2/2
✓ Branch 0 taken 304 times.
✓ Branch 1 taken 257743 times.
 258047 if( (max4 < fabs(p2_3_p0_3)) )
1097 {
1098 304 max4 = fabs(p2_3_p0_3);
1099 }
1100
2/2
✓ Branch 0 taken 133 times.
✓ Branch 1 taken 257914 times.
 258047 if( (max4 < fabs(p2_5_p0_5)) )
1101 {
1102 133 max4 = fabs(p2_5_p0_5);
1103 }
1104
2/2
✓ Branch 0 taken 884 times.
✓ Branch 1 taken 257163 times.
 258047 if( (max4 < fabs(p2_6_p0_6)) )
1105 {
1106 884 max4 = fabs(p2_6_p0_6);
1107 }
1108
2/2
✓ Branch 0 taken 380 times.
✓ Branch 1 taken 257667 times.
 258047 if( (max4 < fabs(p2_7_p0_7)) )
1109 {
1110 380 max4 = fabs(p2_7_p0_7);
1111 }
1112
2/2
✓ Branch 0 taken 6245 times.
✓ Branch 1 taken 251802 times.
 258047 if( (max3 < max4) )
1113 {
1114 6245 max3 = max4;
1115 }
1116 int r_sign;
1117 int int_tmp_result_FFWKCAA;
1118 258047 lower_bound_1 = max1;
1119 258047 upper_bound_1 = max1;
1120
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 258047 times.
 258047 if( (max3 < lower_bound_1) )
1121 {
1122 ✗ lower_bound_1 = max3;
1123 }
1124 else
1125 {
1126
2/2
✓ Branch 0 taken 253478 times.
✓ Branch 1 taken 4569 times.
 258047 if( (max3 > upper_bound_1) )
1127 {
1128 253478 upper_bound_1 = max3;
1129 }
1130 }
1131
2/2
✓ Branch 0 taken 4569 times.
✓ Branch 1 taken 253478 times.
 258047 if( (max4 < lower_bound_1) )
1132 {
1133 4569 lower_bound_1 = max4;
1134 }
1135
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 258047 times.
 258047 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 258047 times.
 258047 if( (upper_bound_1 > 1.87072209578355531992e+50) )
1142 {
1143 ✗ return FPG_UNCERTAIN_VALUE;
1144 }
1145 258047 eps = (8.71140112255785451890e-13 * (((max1 * max4) * max4) * max3));
1146
2/2
✓ Branch 0 taken 107760 times.
✓ Branch 1 taken 150287 times.
 258047 if( (r > eps) )
1147 {
1148 107760 int_tmp_result_FFWKCAA = 1;
1149 }
1150 else
1151 {
1152
1/2
✓ Branch 0 taken 150287 times.
✗ Branch 1 not taken.
 150287 if( (r < -eps) )
1153 {
1154 150287 int_tmp_result_FFWKCAA = -1;
1155 }
1156 else
1157 {
1158 ✗ return FPG_UNCERTAIN_VALUE;
1159 }
1160 }
1161 }
1162 258047 r_sign = int_tmp_result_FFWKCAA;
1163 258047 return (Delta_sign * r_sign);
1164 }
1165