GCC Code Coverage Report


Directory: ./
File: lib/geogram/numerics/predicates/side4.h
Date: 2026-09-07 02:37:58
Exec Total Coverage
Lines: 980 1349 72.6%
Functions: 4 5 80.0%
Branches: 589 842 70.0%

Line Branch Exec Source
1 /* Automatically generated code, do not edit */
2 /* Generated from source file: side4.pck */
3
4 30447796 inline int side4_3d_filter( const double* p0, const double* p1, const double* p2, const double* p3, const double* p4) {
5 double a11;
6 30447796 a11 = (p1[0] - p0[0]);
7 double a12;
8 30447796 a12 = (p1[1] - p0[1]);
9 double a13;
10 30447796 a13 = (p1[2] - p0[2]);
11 30447796 double p1_0_p0_0 = (p1[0] - p0[0]);
12 30447796 double p1_1_p0_1 = (p1[1] - p0[1]);
13 30447796 double p1_2_p0_2 = (p1[2] - p0[2]);
14 double a14;
15 30447796 a14 = -(((p1_0_p0_0 * p1_0_p0_0) + (p1_1_p0_1 * p1_1_p0_1)) + (p1_2_p0_2 * p1_2_p0_2));
16 double a21;
17 30447796 a21 = (p2[0] - p0[0]);
18 double a22;
19 30447796 a22 = (p2[1] - p0[1]);
20 double a23;
21 30447796 a23 = (p2[2] - p0[2]);
22 30447796 double p2_0_p0_0 = (p2[0] - p0[0]);
23 30447796 double p2_1_p0_1 = (p2[1] - p0[1]);
24 30447796 double p2_2_p0_2 = (p2[2] - p0[2]);
25 double a24;
26 30447796 a24 = -(((p2_0_p0_0 * p2_0_p0_0) + (p2_1_p0_1 * p2_1_p0_1)) + (p2_2_p0_2 * p2_2_p0_2));
27 double a31;
28 30447796 a31 = (p3[0] - p0[0]);
29 double a32;
30 30447796 a32 = (p3[1] - p0[1]);
31 double a33;
32 30447796 a33 = (p3[2] - p0[2]);
33 30447796 double p3_0_p0_0 = (p3[0] - p0[0]);
34 30447796 double p3_1_p0_1 = (p3[1] - p0[1]);
35 30447796 double p3_2_p0_2 = (p3[2] - p0[2]);
36 double a34;
37 30447796 a34 = -(((p3_0_p0_0 * p3_0_p0_0) + (p3_1_p0_1 * p3_1_p0_1)) + (p3_2_p0_2 * p3_2_p0_2));
38 double a41;
39 30447796 a41 = (p4[0] - p0[0]);
40 double a42;
41 30447796 a42 = (p4[1] - p0[1]);
42 double a43;
43 30447796 a43 = (p4[2] - p0[2]);
44 30447796 double p4_0_p0_0 = (p4[0] - p0[0]);
45 30447796 double p4_1_p0_1 = (p4[1] - p0[1]);
46 30447796 double p4_2_p0_2 = (p4[2] - p0[2]);
47 double a44;
48 30447796 a44 = -(((p4_0_p0_0 * p4_0_p0_0) + (p4_1_p0_1 * p4_1_p0_1)) + (p4_2_p0_2 * p4_2_p0_2));
49 double Delta1;
50 30447796 Delta1 = (((a21 * ((a32 * a43) - (a33 * a42))) - (a31 * ((a22 * a43) - (a23 * a42)))) + (a41 * ((a22 * a33) - (a23 * a32))));
51 double Delta2;
52 30447796 Delta2 = (((a11 * ((a32 * a43) - (a33 * a42))) - (a31 * ((a12 * a43) - (a13 * a42)))) + (a41 * ((a12 * a33) - (a13 * a32))));
53 double Delta3;
54 30447796 Delta3 = (((a11 * ((a22 * a43) - (a23 * a42))) - (a21 * ((a12 * a43) - (a13 * a42)))) + (a41 * ((a12 * a23) - (a13 * a22))));
55 double Delta4;
56 30447796 Delta4 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
57 double r;
58 30447796 r = ((((Delta1 * a14) - (Delta2 * a24)) + (Delta3 * a34)) - (Delta4 * a44));
59 double eps;
60 30447796 double max1 = fabs(a11);
61
2/2
✓ Branch 0 taken 15182604 times.
✓ Branch 1 taken 15265192 times.
30447796 if( (max1 < fabs(a21)) )
62 {
63 15182604 max1 = fabs(a21);
64 }
65
2/2
✓ Branch 0 taken 9877754 times.
✓ Branch 1 taken 20570042 times.
30447796 if( (max1 < fabs(a31)) )
66 {
67 9877754 max1 = fabs(a31);
68 }
69 30447796 double max2 = fabs(a12);
70
2/2
✓ Branch 0 taken 19874634 times.
✓ Branch 1 taken 10573162 times.
30447796 if( (max2 < fabs(a13)) )
71 {
72 19874634 max2 = fabs(a13);
73 }
74
2/2
✓ Branch 0 taken 5270329 times.
✓ Branch 1 taken 25177467 times.
30447796 if( (max2 < fabs(a22)) )
75 {
76 5270329 max2 = fabs(a22);
77 }
78
2/2
✓ Branch 0 taken 4892717 times.
✓ Branch 1 taken 25555079 times.
30447796 if( (max2 < fabs(a23)) )
79 {
80 4892717 max2 = fabs(a23);
81 }
82 30447796 double max3 = fabs(a22);
83
2/2
✓ Branch 0 taken 19957568 times.
✓ Branch 1 taken 10490228 times.
30447796 if( (max3 < fabs(a23)) )
84 {
85 19957568 max3 = fabs(a23);
86 }
87
2/2
✓ Branch 0 taken 5223506 times.
✓ Branch 1 taken 25224290 times.
30447796 if( (max3 < fabs(a32)) )
88 {
89 5223506 max3 = fabs(a32);
90 }
91
2/2
✓ Branch 0 taken 5277367 times.
✓ Branch 1 taken 25170429 times.
30447796 if( (max3 < fabs(a33)) )
92 {
93 5277367 max3 = fabs(a33);
94 }
95 double lower_bound_1;
96 double upper_bound_1;
97 int Delta4_sign;
98 int int_tmp_result;
99 30447796 lower_bound_1 = max1;
100 30447796 upper_bound_1 = max1;
101
2/2
✓ Branch 0 taken 7615186 times.
✓ Branch 1 taken 22832610 times.
30447796 if( (max3 < lower_bound_1) )
102 {
103 7615186 lower_bound_1 = max3;
104 }
105 else
106 {
107
1/2
✓ Branch 0 taken 22832610 times.
✗ Branch 1 not taken.
22832610 if( (max3 > upper_bound_1) )
108 {
109 22832610 upper_bound_1 = max3;
110 }
111 }
112
2/2
✓ Branch 0 taken 2679394 times.
✓ Branch 1 taken 27768402 times.
30447796 if( (max2 < lower_bound_1) )
113 {
114 2679394 lower_bound_1 = max2;
115 }
116 else
117 {
118
2/2
✓ Branch 0 taken 4409375 times.
✓ Branch 1 taken 23359027 times.
27768402 if( (max2 > upper_bound_1) )
119 {
120 4409375 upper_bound_1 = max2;
121 }
122 }
123
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 30447796 times.
30447796 if( (lower_bound_1 < 1.63288018496748314939e-98) )
124 {
125 return FPG_UNCERTAIN_VALUE;
126 }
127 else
128 {
129
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 30447796 times.
30447796 if( (upper_bound_1 > 3.21387608851797948065e+60) )
130 {
131 return FPG_UNCERTAIN_VALUE;
132 }
133 30447796 eps = (5.11071278299732992696e-15 * ((max2 * max3) * max1));
134
2/2
✓ Branch 0 taken 15158455 times.
✓ Branch 1 taken 15289341 times.
30447796 if( (Delta4 > eps) )
135 {
136 15158455 int_tmp_result = 1;
137 }
138 else
139 {
140
1/2
✓ Branch 0 taken 15289341 times.
✗ Branch 1 not taken.
15289341 if( (Delta4 < -eps) )
141 {
142 15289341 int_tmp_result = -1;
143 }
144 else
145 {
146 return FPG_UNCERTAIN_VALUE;
147 }
148 }
149 }
150 30447796 Delta4_sign = int_tmp_result;
151 int int_tmp_result_FFWKCAA;
152 30447796 double max4 = max1;
153
2/2
✓ Branch 0 taken 14034577 times.
✓ Branch 1 taken 16413219 times.
30447796 if( (max4 < fabs(a41)) )
154 {
155 14034577 max4 = fabs(a41);
156 }
157 30447796 double max5 = max3;
158
2/2
✓ Branch 0 taken 7009523 times.
✓ Branch 1 taken 23438273 times.
30447796 if( (max5 < max2) )
159 {
160 7009523 max5 = max2;
161 }
162 30447796 double max6 = max3;
163
2/2
✓ Branch 0 taken 5547173 times.
✓ Branch 1 taken 24900623 times.
30447796 if( (max6 < fabs(a42)) )
164 {
165 5547173 max6 = fabs(a42);
166 }
167
2/2
✓ Branch 0 taken 4519564 times.
✓ Branch 1 taken 25928232 times.
30447796 if( (max6 < fabs(a43)) )
168 {
169 4519564 max6 = fabs(a43);
170 }
171 30447796 double max7 = fabs(p1_0_p0_0);
172
2/2
✓ Branch 0 taken 14820328 times.
✓ Branch 1 taken 15627468 times.
30447796 if( (max7 < fabs(p1_1_p0_1)) )
173 {
174 14820328 max7 = fabs(p1_1_p0_1);
175 }
176
2/2
✓ Branch 0 taken 16329464 times.
✓ Branch 1 taken 14118332 times.
30447796 if( (max7 < fabs(p1_2_p0_2)) )
177 {
178 16329464 max7 = fabs(p1_2_p0_2);
179 }
180
2/2
✓ Branch 0 taken 3715133 times.
✓ Branch 1 taken 26732663 times.
30447796 if( (max7 < fabs(p2_0_p0_0)) )
181 {
182 3715133 max7 = fabs(p2_0_p0_0);
183 }
184
2/2
✓ Branch 0 taken 3079639 times.
✓ Branch 1 taken 27368157 times.
30447796 if( (max7 < fabs(p2_2_p0_2)) )
185 {
186 3079639 max7 = fabs(p2_2_p0_2);
187 }
188
2/2
✓ Branch 0 taken 3328780 times.
✓ Branch 1 taken 27119016 times.
30447796 if( (max7 < fabs(p2_1_p0_1)) )
189 {
190 3328780 max7 = fabs(p2_1_p0_1);
191 }
192
2/2
✓ Branch 0 taken 2551130 times.
✓ Branch 1 taken 27896666 times.
30447796 if( (max7 < fabs(p3_0_p0_0)) )
193 {
194 2551130 max7 = fabs(p3_0_p0_0);
195 }
196
2/2
✓ Branch 0 taken 2023689 times.
✓ Branch 1 taken 28424107 times.
30447796 if( (max7 < fabs(p3_1_p0_1)) )
197 {
198 2023689 max7 = fabs(p3_1_p0_1);
199 }
200
2/2
✓ Branch 0 taken 2186860 times.
✓ Branch 1 taken 28260936 times.
30447796 if( (max7 < fabs(p3_2_p0_2)) )
201 {
202 2186860 max7 = fabs(p3_2_p0_2);
203 }
204
2/2
✓ Branch 0 taken 3088310 times.
✓ Branch 1 taken 27359486 times.
30447796 if( (max7 < fabs(p4_0_p0_0)) )
205 {
206 3088310 max7 = fabs(p4_0_p0_0);
207 }
208
2/2
✓ Branch 0 taken 2781729 times.
✓ Branch 1 taken 27666067 times.
30447796 if( (max7 < fabs(p4_1_p0_1)) )
209 {
210 2781729 max7 = fabs(p4_1_p0_1);
211 }
212
2/2
✓ Branch 0 taken 2450583 times.
✓ Branch 1 taken 27997213 times.
30447796 if( (max7 < fabs(p4_2_p0_2)) )
213 {
214 2450583 max7 = fabs(p4_2_p0_2);
215 }
216 30447796 lower_bound_1 = max7;
217 30447796 upper_bound_1 = max7;
218
2/2
✓ Branch 0 taken 23013229 times.
✓ Branch 1 taken 7434567 times.
30447796 if( (max4 < lower_bound_1) )
219 {
220 23013229 lower_bound_1 = max4;
221 }
222 else
223 {
224
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 7434567 times.
7434567 if( (max4 > upper_bound_1) )
225 {
226 upper_bound_1 = max4;
227 }
228 }
229
2/2
✓ Branch 0 taken 9400746 times.
✓ Branch 1 taken 21047050 times.
30447796 if( (max5 < lower_bound_1) )
230 {
231 9400746 lower_bound_1 = max5;
232 }
233 else
234 {
235
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 21047050 times.
21047050 if( (max5 > upper_bound_1) )
236 {
237 upper_bound_1 = max5;
238 }
239 }
240
2/2
✓ Branch 0 taken 1497418 times.
✓ Branch 1 taken 28950378 times.
30447796 if( (max6 < lower_bound_1) )
241 {
242 1497418 lower_bound_1 = max6;
243 }
244 else
245 {
246
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28950378 times.
28950378 if( (max6 > upper_bound_1) )
247 {
248 upper_bound_1 = max6;
249 }
250 }
251
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 30447796 times.
30447796 if( (lower_bound_1 < 1.12285198342304832993e-59) )
252 {
253 return FPG_UNCERTAIN_VALUE;
254 }
255 else
256 {
257
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 30447796 times.
30447796 if( (upper_bound_1 > 3.21387608851797948065e+60) )
258 {
259 return FPG_UNCERTAIN_VALUE;
260 }
261 30447796 eps = (1.24661365310273025710e-13 * ((((max5 * max6) * max4) * max7) * max7));
262
2/2
✓ Branch 0 taken 15036910 times.
✓ Branch 1 taken 15410886 times.
30447796 if( (r > eps) )
263 {
264 15036910 int_tmp_result_FFWKCAA = 1;
265 }
266 else
267 {
268
2/2
✓ Branch 0 taken 15312501 times.
✓ Branch 1 taken 98385 times.
15410886 if( (r < -eps) )
269 {
270 15312501 int_tmp_result_FFWKCAA = -1;
271 }
272 else
273 {
274 98385 return FPG_UNCERTAIN_VALUE;
275 }
276 }
277 }
278 30349411 return (Delta4_sign * int_tmp_result_FFWKCAA);
279 }
280
281
282 28747 inline int side4_4d_filter( const double* p0, const double* p1, const double* p2, const double* p3, const double* p4, const double* q0, const double* q1, const double* q2, const double* q3) {
283 28747 double p1_0_p0_0 = (p1[0] - p0[0]);
284 28747 double p1_1_p0_1 = (p1[1] - p0[1]);
285 28747 double p1_2_p0_2 = (p1[2] - p0[2]);
286 28747 double p1_3_p0_3 = (p1[3] - p0[3]);
287 double l1;
288 28747 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)));
289 28747 double p2_0_p0_0 = (p2[0] - p0[0]);
290 28747 double p2_1_p0_1 = (p2[1] - p0[1]);
291 28747 double p2_2_p0_2 = (p2[2] - p0[2]);
292 28747 double p2_3_p0_3 = (p2[3] - p0[3]);
293 double l2;
294 28747 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)));
295 28747 double p3_0_p0_0 = (p3[0] - p0[0]);
296 28747 double p3_1_p0_1 = (p3[1] - p0[1]);
297 28747 double p3_2_p0_2 = (p3[2] - p0[2]);
298 28747 double p3_3_p0_3 = (p3[3] - p0[3]);
299 double l3;
300 28747 l3 = (1 * ((((p3_0_p0_0 * p3_0_p0_0) + (p3_1_p0_1 * p3_1_p0_1)) + (p3_2_p0_2 * p3_2_p0_2)) + (p3_3_p0_3 * p3_3_p0_3)));
301 28747 double p4_0_p0_0 = (p4[0] - p0[0]);
302 28747 double p4_1_p0_1 = (p4[1] - p0[1]);
303 28747 double p4_2_p0_2 = (p4[2] - p0[2]);
304 28747 double p4_3_p0_3 = (p4[3] - p0[3]);
305 double l4;
306 28747 l4 = (1 * ((((p4_0_p0_0 * p4_0_p0_0) + (p4_1_p0_1 * p4_1_p0_1)) + (p4_2_p0_2 * p4_2_p0_2)) + (p4_3_p0_3 * p4_3_p0_3)));
307 28747 double q0_0_p0_0 = (q0[0] - p0[0]);
308 28747 double q0_1_p0_1 = (q0[1] - p0[1]);
309 28747 double q0_2_p0_2 = (q0[2] - p0[2]);
310 28747 double q0_3_p0_3 = (q0[3] - p0[3]);
311 double a10;
312 28747 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)));
313 28747 double q1_0_p0_0 = (q1[0] - p0[0]);
314 28747 double q1_1_p0_1 = (q1[1] - p0[1]);
315 28747 double q1_2_p0_2 = (q1[2] - p0[2]);
316 28747 double q1_3_p0_3 = (q1[3] - p0[3]);
317 double a11;
318 28747 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)));
319 28747 double q2_0_p0_0 = (q2[0] - p0[0]);
320 28747 double q2_1_p0_1 = (q2[1] - p0[1]);
321 28747 double q2_2_p0_2 = (q2[2] - p0[2]);
322 28747 double q2_3_p0_3 = (q2[3] - p0[3]);
323 double a12;
324 28747 a12 = (2 * ((((p1_0_p0_0 * q2_0_p0_0) + (p1_1_p0_1 * q2_1_p0_1)) + (p1_2_p0_2 * q2_2_p0_2)) + (p1_3_p0_3 * q2_3_p0_3)));
325 28747 double q3_0_p0_0 = (q3[0] - p0[0]);
326 28747 double q3_1_p0_1 = (q3[1] - p0[1]);
327 28747 double q3_2_p0_2 = (q3[2] - p0[2]);
328 28747 double q3_3_p0_3 = (q3[3] - p0[3]);
329 double a13;
330 28747 a13 = (2 * ((((p1_0_p0_0 * q3_0_p0_0) + (p1_1_p0_1 * q3_1_p0_1)) + (p1_2_p0_2 * q3_2_p0_2)) + (p1_3_p0_3 * q3_3_p0_3)));
331 double a20;
332 28747 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)));
333 double a21;
334 28747 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)));
335 double a22;
336 28747 a22 = (2 * ((((p2_0_p0_0 * q2_0_p0_0) + (p2_1_p0_1 * q2_1_p0_1)) + (p2_2_p0_2 * q2_2_p0_2)) + (p2_3_p0_3 * q2_3_p0_3)));
337 double a23;
338 28747 a23 = (2 * ((((p2_0_p0_0 * q3_0_p0_0) + (p2_1_p0_1 * q3_1_p0_1)) + (p2_2_p0_2 * q3_2_p0_2)) + (p2_3_p0_3 * q3_3_p0_3)));
339 double a30;
340 28747 a30 = (2 * ((((p3_0_p0_0 * q0_0_p0_0) + (p3_1_p0_1 * q0_1_p0_1)) + (p3_2_p0_2 * q0_2_p0_2)) + (p3_3_p0_3 * q0_3_p0_3)));
341 double a31;
342 28747 a31 = (2 * ((((p3_0_p0_0 * q1_0_p0_0) + (p3_1_p0_1 * q1_1_p0_1)) + (p3_2_p0_2 * q1_2_p0_2)) + (p3_3_p0_3 * q1_3_p0_3)));
343 double a32;
344 28747 a32 = (2 * ((((p3_0_p0_0 * q2_0_p0_0) + (p3_1_p0_1 * q2_1_p0_1)) + (p3_2_p0_2 * q2_2_p0_2)) + (p3_3_p0_3 * q2_3_p0_3)));
345 double a33;
346 28747 a33 = (2 * ((((p3_0_p0_0 * q3_0_p0_0) + (p3_1_p0_1 * q3_1_p0_1)) + (p3_2_p0_2 * q3_2_p0_2)) + (p3_3_p0_3 * q3_3_p0_3)));
347 double a40;
348 28747 a40 = (2 * ((((p4_0_p0_0 * q0_0_p0_0) + (p4_1_p0_1 * q0_1_p0_1)) + (p4_2_p0_2 * q0_2_p0_2)) + (p4_3_p0_3 * q0_3_p0_3)));
349 double a41;
350 28747 a41 = (2 * ((((p4_0_p0_0 * q1_0_p0_0) + (p4_1_p0_1 * q1_1_p0_1)) + (p4_2_p0_2 * q1_2_p0_2)) + (p4_3_p0_3 * q1_3_p0_3)));
351 double a42;
352 28747 a42 = (2 * ((((p4_0_p0_0 * q2_0_p0_0) + (p4_1_p0_1 * q2_1_p0_1)) + (p4_2_p0_2 * q2_2_p0_2)) + (p4_3_p0_3 * q2_3_p0_3)));
353 double a43;
354 28747 a43 = (2 * ((((p4_0_p0_0 * q3_0_p0_0) + (p4_1_p0_1 * q3_1_p0_1)) + (p4_2_p0_2 * q3_2_p0_2)) + (p4_3_p0_3 * q3_3_p0_3)));
355 double b00;
356 28747 b00 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
357 double b01;
358 28747 b01 = -((((a22 * a33) - (a23 * a32)) + ((a23 * a31) - (a21 * a33))) + ((a21 * a32) - (a22 * a31)));
359 double b02;
360 28747 b02 = ((((a12 * a33) - (a13 * a32)) + ((a13 * a31) - (a11 * a33))) + ((a11 * a32) - (a12 * a31)));
361 double b03;
362 28747 b03 = -((((a12 * a23) - (a13 * a22)) + ((a13 * a21) - (a11 * a23))) + ((a11 * a22) - (a12 * a21)));
363 double b10;
364 28747 b10 = -(((a10 * ((a22 * a33) - (a23 * a32))) - (a20 * ((a12 * a33) - (a13 * a32)))) + (a30 * ((a12 * a23) - (a13 * a22))));
365 double b11;
366 28747 b11 = ((((a22 * a33) - (a23 * a32)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a32) - (a22 * a30)));
367 double b12;
368 28747 b12 = -((((a12 * a33) - (a13 * a32)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a32) - (a12 * a30)));
369 double b13;
370 28747 b13 = ((((a12 * a23) - (a13 * a22)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a22) - (a12 * a20)));
371 double b20;
372 28747 b20 = (((a10 * ((a21 * a33) - (a23 * a31))) - (a20 * ((a11 * a33) - (a13 * a31)))) + (a30 * ((a11 * a23) - (a13 * a21))));
373 double b21;
374 28747 b21 = -((((a21 * a33) - (a23 * a31)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a31) - (a21 * a30)));
375 double b22;
376 28747 b22 = ((((a11 * a33) - (a13 * a31)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a31) - (a11 * a30)));
377 double b23;
378 28747 b23 = -((((a11 * a23) - (a13 * a21)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a21) - (a11 * a20)));
379 double b30;
380 28747 b30 = -(((a10 * ((a21 * a32) - (a22 * a31))) - (a20 * ((a11 * a32) - (a12 * a31)))) + (a30 * ((a11 * a22) - (a12 * a21))));
381 double b31;
382 28747 b31 = ((((a21 * a32) - (a22 * a31)) + ((a22 * a30) - (a20 * a32))) + ((a20 * a31) - (a21 * a30)));
383 double b32;
384 28747 b32 = -((((a11 * a32) - (a12 * a31)) + ((a12 * a30) - (a10 * a32))) + ((a10 * a31) - (a11 * a30)));
385 double b33;
386 28747 b33 = ((((a11 * a22) - (a12 * a21)) + ((a12 * a20) - (a10 * a22))) + ((a10 * a21) - (a11 * a20)));
387 double Delta;
388 28747 Delta = (((b00 + b10) + b20) + b30);
389 double DeltaLambda0;
390 28747 DeltaLambda0 = ((((b01 * l1) + (b02 * l2)) + (b03 * l3)) + b00);
391 double DeltaLambda1;
392 28747 DeltaLambda1 = ((((b11 * l1) + (b12 * l2)) + (b13 * l3)) + b10);
393 double DeltaLambda2;
394 28747 DeltaLambda2 = ((((b21 * l1) + (b22 * l2)) + (b23 * l3)) + b20);
395 double DeltaLambda3;
396 28747 DeltaLambda3 = ((((b31 * l1) + (b32 * l2)) + (b33 * l3)) + b30);
397 double r;
398 28747 r = ((Delta * l4) - ((((a40 * DeltaLambda0) + (a41 * DeltaLambda1)) + (a42 * DeltaLambda2)) + (a43 * DeltaLambda3)));
399 double eps;
400 28747 double max1 = fabs(p1_3_p0_3);
401
2/2
✓ Branch 0 taken 2929 times.
✓ Branch 1 taken 25818 times.
28747 if( (max1 < fabs(p1_0_p0_0)) )
402 {
403 2929 max1 = fabs(p1_0_p0_0);
404 }
405
2/2
✓ Branch 0 taken 14323 times.
✓ Branch 1 taken 14424 times.
28747 if( (max1 < fabs(p1_1_p0_1)) )
406 {
407 14323 max1 = fabs(p1_1_p0_1);
408 }
409
2/2
✓ Branch 0 taken 9950 times.
✓ Branch 1 taken 18797 times.
28747 if( (max1 < fabs(p1_2_p0_2)) )
410 {
411 9950 max1 = fabs(p1_2_p0_2);
412 }
413 28747 double max2 = fabs(p2_2_p0_2);
414
2/2
✓ Branch 0 taken 14352 times.
✓ Branch 1 taken 14395 times.
28747 if( (max2 < fabs(p2_1_p0_1)) )
415 {
416 14352 max2 = fabs(p2_1_p0_1);
417 }
418
2/2
✓ Branch 0 taken 7536 times.
✓ Branch 1 taken 21211 times.
28747 if( (max2 < fabs(p2_3_p0_3)) )
419 {
420 7536 max2 = fabs(p2_3_p0_3);
421 }
422
2/2
✓ Branch 0 taken 148 times.
✓ Branch 1 taken 28599 times.
28747 if( (max2 < fabs(p2_0_p0_0)) )
423 {
424 148 max2 = fabs(p2_0_p0_0);
425 }
426 28747 double max3 = fabs(p3_0_p0_0);
427
2/2
✓ Branch 0 taken 21100 times.
✓ Branch 1 taken 7647 times.
28747 if( (max3 < fabs(p3_1_p0_1)) )
428 {
429 21100 max3 = fabs(p3_1_p0_1);
430 }
431
2/2
✓ Branch 0 taken 13937 times.
✓ Branch 1 taken 14810 times.
28747 if( (max3 < fabs(p3_3_p0_3)) )
432 {
433 13937 max3 = fabs(p3_3_p0_3);
434 }
435
2/2
✓ Branch 0 taken 10987 times.
✓ Branch 1 taken 17760 times.
28747 if( (max3 < fabs(p3_2_p0_2)) )
436 {
437 10987 max3 = fabs(p3_2_p0_2);
438 }
439 28747 double max4 = fabs(q0_3_p0_3);
440
2/2
✓ Branch 0 taken 181 times.
✓ Branch 1 taken 28566 times.
28747 if( (max4 < fabs(q0_0_p0_0)) )
441 {
442 181 max4 = fabs(q0_0_p0_0);
443 }
444
2/2
✓ Branch 0 taken 3380 times.
✓ Branch 1 taken 25367 times.
28747 if( (max4 < fabs(q0_1_p0_1)) )
445 {
446 3380 max4 = fabs(q0_1_p0_1);
447 }
448
2/2
✓ Branch 0 taken 3246 times.
✓ Branch 1 taken 25501 times.
28747 if( (max4 < fabs(q0_2_p0_2)) )
449 {
450 3246 max4 = fabs(q0_2_p0_2);
451 }
452
2/2
✓ Branch 0 taken 2195 times.
✓ Branch 1 taken 26552 times.
28747 if( (max4 < fabs(q1_0_p0_0)) )
453 {
454 2195 max4 = fabs(q1_0_p0_0);
455 }
456
2/2
✓ Branch 0 taken 1781 times.
✓ Branch 1 taken 26966 times.
28747 if( (max4 < fabs(q1_1_p0_1)) )
457 {
458 1781 max4 = fabs(q1_1_p0_1);
459 }
460
2/2
✓ Branch 0 taken 1648 times.
✓ Branch 1 taken 27099 times.
28747 if( (max4 < fabs(q1_2_p0_2)) )
461 {
462 1648 max4 = fabs(q1_2_p0_2);
463 }
464
2/2
✓ Branch 0 taken 11103 times.
✓ Branch 1 taken 17644 times.
28747 if( (max4 < fabs(q1_3_p0_3)) )
465 {
466 11103 max4 = fabs(q1_3_p0_3);
467 }
468 28747 double max5 = fabs(q1_0_p0_0);
469
2/2
✓ Branch 0 taken 16782 times.
✓ Branch 1 taken 11965 times.
28747 if( (max5 < fabs(q1_1_p0_1)) )
470 {
471 16782 max5 = fabs(q1_1_p0_1);
472 }
473
2/2
✓ Branch 0 taken 10026 times.
✓ Branch 1 taken 18721 times.
28747 if( (max5 < fabs(q1_2_p0_2)) )
474 {
475 10026 max5 = fabs(q1_2_p0_2);
476 }
477
2/2
✓ Branch 0 taken 21815 times.
✓ Branch 1 taken 6932 times.
28747 if( (max5 < fabs(q1_3_p0_3)) )
478 {
479 21815 max5 = fabs(q1_3_p0_3);
480 }
481
2/2
✓ Branch 0 taken 2572 times.
✓ Branch 1 taken 26175 times.
28747 if( (max5 < fabs(q2_0_p0_0)) )
482 {
483 2572 max5 = fabs(q2_0_p0_0);
484 }
485
2/2
✓ Branch 0 taken 1877 times.
✓ Branch 1 taken 26870 times.
28747 if( (max5 < fabs(q2_1_p0_1)) )
486 {
487 1877 max5 = fabs(q2_1_p0_1);
488 }
489
2/2
✓ Branch 0 taken 1726 times.
✓ Branch 1 taken 27021 times.
28747 if( (max5 < fabs(q2_2_p0_2)) )
490 {
491 1726 max5 = fabs(q2_2_p0_2);
492 }
493
2/2
✓ Branch 0 taken 12293 times.
✓ Branch 1 taken 16454 times.
28747 if( (max5 < fabs(q2_3_p0_3)) )
494 {
495 12293 max5 = fabs(q2_3_p0_3);
496 }
497 28747 double max6 = fabs(q2_0_p0_0);
498
2/2
✓ Branch 0 taken 16608 times.
✓ Branch 1 taken 12139 times.
28747 if( (max6 < fabs(q2_1_p0_1)) )
499 {
500 16608 max6 = fabs(q2_1_p0_1);
501 }
502
2/2
✓ Branch 0 taken 10748 times.
✓ Branch 1 taken 17999 times.
28747 if( (max6 < fabs(q2_2_p0_2)) )
503 {
504 10748 max6 = fabs(q2_2_p0_2);
505 }
506
2/2
✓ Branch 0 taken 22939 times.
✓ Branch 1 taken 5808 times.
28747 if( (max6 < fabs(q2_3_p0_3)) )
507 {
508 22939 max6 = fabs(q2_3_p0_3);
509 }
510
2/2
✓ Branch 0 taken 3131 times.
✓ Branch 1 taken 25616 times.
28747 if( (max6 < fabs(q3_0_p0_0)) )
511 {
512 3131 max6 = fabs(q3_0_p0_0);
513 }
514
2/2
✓ Branch 0 taken 2393 times.
✓ Branch 1 taken 26354 times.
28747 if( (max6 < fabs(q3_1_p0_1)) )
515 {
516 2393 max6 = fabs(q3_1_p0_1);
517 }
518
2/2
✓ Branch 0 taken 1483 times.
✓ Branch 1 taken 27264 times.
28747 if( (max6 < fabs(q3_2_p0_2)) )
519 {
520 1483 max6 = fabs(q3_2_p0_2);
521 }
522
2/2
✓ Branch 0 taken 13269 times.
✓ Branch 1 taken 15478 times.
28747 if( (max6 < fabs(q3_3_p0_3)) )
523 {
524 13269 max6 = fabs(q3_3_p0_3);
525 }
526 double lower_bound_1;
527 double upper_bound_1;
528 int Delta_sign;
529 int int_tmp_result;
530 28747 lower_bound_1 = max1;
531 28747 upper_bound_1 = max1;
532
2/2
✓ Branch 0 taken 14322 times.
✓ Branch 1 taken 14425 times.
28747 if( (max2 < lower_bound_1) )
533 {
534 14322 lower_bound_1 = max2;
535 }
536 else
537 {
538
1/2
✓ Branch 0 taken 14425 times.
✗ Branch 1 not taken.
14425 if( (max2 > upper_bound_1) )
539 {
540 14425 upper_bound_1 = max2;
541 }
542 }
543
2/2
✓ Branch 0 taken 9932 times.
✓ Branch 1 taken 18815 times.
28747 if( (max3 < lower_bound_1) )
544 {
545 9932 lower_bound_1 = max3;
546 }
547 else
548 {
549
2/2
✓ Branch 0 taken 9080 times.
✓ Branch 1 taken 9735 times.
18815 if( (max3 > upper_bound_1) )
550 {
551 9080 upper_bound_1 = max3;
552 }
553 }
554
2/2
✓ Branch 0 taken 672 times.
✓ Branch 1 taken 28075 times.
28747 if( (max4 < lower_bound_1) )
555 {
556 672 lower_bound_1 = max4;
557 }
558 else
559 {
560
2/2
✓ Branch 0 taken 26884 times.
✓ Branch 1 taken 1191 times.
28075 if( (max4 > upper_bound_1) )
561 {
562 26884 upper_bound_1 = max4;
563 }
564 }
565
2/2
✓ Branch 0 taken 538 times.
✓ Branch 1 taken 28209 times.
28747 if( (max5 < lower_bound_1) )
566 {
567 538 lower_bound_1 = max5;
568 }
569 else
570 {
571
2/2
✓ Branch 0 taken 9047 times.
✓ Branch 1 taken 19162 times.
28209 if( (max5 > upper_bound_1) )
572 {
573 9047 upper_bound_1 = max5;
574 }
575 }
576
2/2
✓ Branch 0 taken 241 times.
✓ Branch 1 taken 28506 times.
28747 if( (max6 < lower_bound_1) )
577 {
578 241 lower_bound_1 = max6;
579 }
580 else
581 {
582
2/2
✓ Branch 0 taken 7780 times.
✓ Branch 1 taken 20726 times.
28506 if( (max6 > upper_bound_1) )
583 {
584 7780 upper_bound_1 = max6;
585 }
586 }
587
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28747 times.
28747 if( (lower_bound_1 < 4.14607644401726239868e-50) )
588 {
589 return FPG_UNCERTAIN_VALUE;
590 }
591 else
592 {
593
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28747 times.
28747 if( (upper_bound_1 > 4.83570327845851562508e+24) )
594 {
595 return FPG_UNCERTAIN_VALUE;
596 }
597 28747 eps = (4.38046888801178809320e-12 * (((((max1 * max4) * max2) * max5) * max3) * max6));
598
2/2
✓ Branch 0 taken 14316 times.
✓ Branch 1 taken 14431 times.
28747 if( (Delta > eps) )
599 {
600 14316 int_tmp_result = 1;
601 }
602 else
603 {
604
1/2
✓ Branch 0 taken 14431 times.
✗ Branch 1 not taken.
14431 if( (Delta < -eps) )
605 {
606 14431 int_tmp_result = -1;
607 }
608 else
609 {
610 return FPG_UNCERTAIN_VALUE;
611 }
612 }
613 }
614 28747 Delta_sign = int_tmp_result;
615 int int_tmp_result_FFWKCAA;
616 28747 double max7 = max1;
617
2/2
✓ Branch 0 taken 14425 times.
✓ Branch 1 taken 14322 times.
28747 if( (max7 < max2) )
618 {
619 14425 max7 = max2;
620 }
621
2/2
✓ Branch 0 taken 9080 times.
✓ Branch 1 taken 19667 times.
28747 if( (max7 < max3) )
622 {
623 9080 max7 = max3;
624 }
625
2/2
✓ Branch 0 taken 27362 times.
✓ Branch 1 taken 1385 times.
28747 if( (max7 < max6) )
626 {
627 27362 max7 = max6;
628 }
629 28747 double max8 = max1;
630
2/2
✓ Branch 0 taken 9560 times.
✓ Branch 1 taken 19187 times.
28747 if( (max8 < fabs(p4_1_p0_1)) )
631 {
632 9560 max8 = fabs(p4_1_p0_1);
633 }
634
2/2
✓ Branch 0 taken 7713 times.
✓ Branch 1 taken 21034 times.
28747 if( (max8 < fabs(p4_2_p0_2)) )
635 {
636 7713 max8 = fabs(p4_2_p0_2);
637 }
638
2/2
✓ Branch 0 taken 1719 times.
✓ Branch 1 taken 27028 times.
28747 if( (max8 < fabs(p4_0_p0_0)) )
639 {
640 1719 max8 = fabs(p4_0_p0_0);
641 }
642
2/2
✓ Branch 0 taken 7059 times.
✓ Branch 1 taken 21688 times.
28747 if( (max8 < fabs(p4_3_p0_3)) )
643 {
644 7059 max8 = fabs(p4_3_p0_3);
645 }
646
2/2
✓ Branch 0 taken 1545 times.
✓ Branch 1 taken 27202 times.
28747 if( (max7 < max8) )
647 {
648 1545 max7 = max8;
649 }
650 28747 double max9 = max1;
651
2/2
✓ Branch 0 taken 27366 times.
✓ Branch 1 taken 1381 times.
28747 if( (max9 < max5) )
652 {
653 27366 max9 = max5;
654 }
655
2/2
✓ Branch 0 taken 2620 times.
✓ Branch 1 taken 26127 times.
28747 if( (max9 < max8) )
656 {
657 2620 max9 = max8;
658 }
659 double max10;
660 28747 double max11 = max4;
661
2/2
✓ Branch 0 taken 9276 times.
✓ Branch 1 taken 19471 times.
28747 if( (max11 < max5) )
662 {
663 9276 max11 = max5;
664 }
665 28747 max10 = max11;
666
2/2
✓ Branch 0 taken 337 times.
✓ Branch 1 taken 28410 times.
28747 if( (max10 < max1) )
667 {
668 337 max10 = max1;
669 }
670
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28747 times.
28747 if( (max10 < max4) )
671 {
672 max10 = max4;
673 }
674
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28747 times.
28747 if( (max10 < max5) )
675 {
676 max10 = max5;
677 }
678
2/2
✓ Branch 0 taken 7781 times.
✓ Branch 1 taken 20966 times.
28747 if( (max10 < max6) )
679 {
680 7781 max10 = max6;
681 }
682 28747 lower_bound_1 = max10;
683 28747 upper_bound_1 = max10;
684
2/2
✓ Branch 0 taken 7792 times.
✓ Branch 1 taken 20955 times.
28747 if( (max11 < lower_bound_1) )
685 {
686 7792 lower_bound_1 = max11;
687 }
688 else
689 {
690
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20955 times.
20955 if( (max11 > upper_bound_1) )
691 {
692 upper_bound_1 = max11;
693 }
694 }
695
2/2
✓ Branch 0 taken 28476 times.
✓ Branch 1 taken 271 times.
28747 if( (max2 < lower_bound_1) )
696 {
697 28476 lower_bound_1 = max2;
698 }
699
2/2
✓ Branch 0 taken 14671 times.
✓ Branch 1 taken 14076 times.
28747 if( (max3 < lower_bound_1) )
700 {
701 14671 lower_bound_1 = max3;
702 }
703
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28747 times.
28747 if( (max7 < lower_bound_1) )
704 {
705 lower_bound_1 = max7;
706 }
707 else
708 {
709
2/2
✓ Branch 0 taken 56 times.
✓ Branch 1 taken 28691 times.
28747 if( (max7 > upper_bound_1) )
710 {
711 56 upper_bound_1 = max7;
712 }
713 }
714
2/2
✓ Branch 0 taken 1493 times.
✓ Branch 1 taken 27254 times.
28747 if( (max8 < lower_bound_1) )
715 {
716 1493 lower_bound_1 = max8;
717 }
718
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28747 times.
28747 if( (max9 < lower_bound_1) )
719 {
720 lower_bound_1 = max9;
721 }
722 else
723 {
724
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28747 times.
28747 if( (max9 > upper_bound_1) )
725 {
726 upper_bound_1 = max9;
727 }
728 }
729
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28747 times.
28747 if( (lower_bound_1 < 6.06263132863556750071e-38) )
730 {
731 return FPG_UNCERTAIN_VALUE;
732 }
733 else
734 {
735
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 28747 times.
28747 if( (upper_bound_1 > 4.83570327845851562508e+24) )
736 {
737 return FPG_UNCERTAIN_VALUE;
738 }
739 28747 eps = (1.21914442286604163181e-10 * (((((((max8 * max11) * max2) * max10) * max3) * max10) * max9) * max7));
740
2/2
✓ Branch 0 taken 14336 times.
✓ Branch 1 taken 14411 times.
28747 if( (r > eps) )
741 {
742 14336 int_tmp_result_FFWKCAA = 1;
743 }
744 else
745 {
746
1/2
✓ Branch 0 taken 14411 times.
✗ Branch 1 not taken.
14411 if( (r < -eps) )
747 {
748 14411 int_tmp_result_FFWKCAA = -1;
749 }
750 else
751 {
752 return FPG_UNCERTAIN_VALUE;
753 }
754 }
755 }
756 28747 return (Delta_sign * int_tmp_result_FFWKCAA);
757 }
758
759
760 20219 inline int side4_6d_filter( const double* p0, const double* p1, const double* p2, const double* p3, const double* p4, const double* q0, const double* q1, const double* q2, const double* q3) {
761 20219 double p1_0_p0_0 = (p1[0] - p0[0]);
762 20219 double p1_1_p0_1 = (p1[1] - p0[1]);
763 20219 double p1_2_p0_2 = (p1[2] - p0[2]);
764 20219 double p1_3_p0_3 = (p1[3] - p0[3]);
765 20219 double p1_4_p0_4 = (p1[4] - p0[4]);
766 20219 double p1_5_p0_5 = (p1[5] - p0[5]);
767 double l1;
768 20219 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)));
769 20219 double p2_0_p0_0 = (p2[0] - p0[0]);
770 20219 double p2_1_p0_1 = (p2[1] - p0[1]);
771 20219 double p2_2_p0_2 = (p2[2] - p0[2]);
772 20219 double p2_3_p0_3 = (p2[3] - p0[3]);
773 20219 double p2_4_p0_4 = (p2[4] - p0[4]);
774 20219 double p2_5_p0_5 = (p2[5] - p0[5]);
775 double l2;
776 20219 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)));
777 20219 double p3_0_p0_0 = (p3[0] - p0[0]);
778 20219 double p3_1_p0_1 = (p3[1] - p0[1]);
779 20219 double p3_2_p0_2 = (p3[2] - p0[2]);
780 20219 double p3_3_p0_3 = (p3[3] - p0[3]);
781 20219 double p3_4_p0_4 = (p3[4] - p0[4]);
782 20219 double p3_5_p0_5 = (p3[5] - p0[5]);
783 double l3;
784 20219 l3 = (1 * ((((((p3_0_p0_0 * p3_0_p0_0) + (p3_1_p0_1 * p3_1_p0_1)) + (p3_2_p0_2 * p3_2_p0_2)) + (p3_3_p0_3 * p3_3_p0_3)) + (p3_4_p0_4 * p3_4_p0_4)) + (p3_5_p0_5 * p3_5_p0_5)));
785 20219 double p4_0_p0_0 = (p4[0] - p0[0]);
786 20219 double p4_1_p0_1 = (p4[1] - p0[1]);
787 20219 double p4_2_p0_2 = (p4[2] - p0[2]);
788 20219 double p4_3_p0_3 = (p4[3] - p0[3]);
789 20219 double p4_4_p0_4 = (p4[4] - p0[4]);
790 20219 double p4_5_p0_5 = (p4[5] - p0[5]);
791 double l4;
792 20219 l4 = (1 * ((((((p4_0_p0_0 * p4_0_p0_0) + (p4_1_p0_1 * p4_1_p0_1)) + (p4_2_p0_2 * p4_2_p0_2)) + (p4_3_p0_3 * p4_3_p0_3)) + (p4_4_p0_4 * p4_4_p0_4)) + (p4_5_p0_5 * p4_5_p0_5)));
793 20219 double q0_0_p0_0 = (q0[0] - p0[0]);
794 20219 double q0_1_p0_1 = (q0[1] - p0[1]);
795 20219 double q0_2_p0_2 = (q0[2] - p0[2]);
796 20219 double q0_3_p0_3 = (q0[3] - p0[3]);
797 20219 double q0_4_p0_4 = (q0[4] - p0[4]);
798 20219 double q0_5_p0_5 = (q0[5] - p0[5]);
799 double a10;
800 20219 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)));
801 20219 double q1_0_p0_0 = (q1[0] - p0[0]);
802 20219 double q1_1_p0_1 = (q1[1] - p0[1]);
803 20219 double q1_2_p0_2 = (q1[2] - p0[2]);
804 20219 double q1_3_p0_3 = (q1[3] - p0[3]);
805 20219 double q1_4_p0_4 = (q1[4] - p0[4]);
806 20219 double q1_5_p0_5 = (q1[5] - p0[5]);
807 double a11;
808 20219 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)));
809 20219 double q2_0_p0_0 = (q2[0] - p0[0]);
810 20219 double q2_1_p0_1 = (q2[1] - p0[1]);
811 20219 double q2_2_p0_2 = (q2[2] - p0[2]);
812 20219 double q2_3_p0_3 = (q2[3] - p0[3]);
813 20219 double q2_4_p0_4 = (q2[4] - p0[4]);
814 20219 double q2_5_p0_5 = (q2[5] - p0[5]);
815 double a12;
816 20219 a12 = (2 * ((((((p1_0_p0_0 * q2_0_p0_0) + (p1_1_p0_1 * q2_1_p0_1)) + (p1_2_p0_2 * q2_2_p0_2)) + (p1_3_p0_3 * q2_3_p0_3)) + (p1_4_p0_4 * q2_4_p0_4)) + (p1_5_p0_5 * q2_5_p0_5)));
817 20219 double q3_0_p0_0 = (q3[0] - p0[0]);
818 20219 double q3_1_p0_1 = (q3[1] - p0[1]);
819 20219 double q3_2_p0_2 = (q3[2] - p0[2]);
820 20219 double q3_3_p0_3 = (q3[3] - p0[3]);
821 20219 double q3_4_p0_4 = (q3[4] - p0[4]);
822 20219 double q3_5_p0_5 = (q3[5] - p0[5]);
823 double a13;
824 20219 a13 = (2 * ((((((p1_0_p0_0 * q3_0_p0_0) + (p1_1_p0_1 * q3_1_p0_1)) + (p1_2_p0_2 * q3_2_p0_2)) + (p1_3_p0_3 * q3_3_p0_3)) + (p1_4_p0_4 * q3_4_p0_4)) + (p1_5_p0_5 * q3_5_p0_5)));
825 double a20;
826 20219 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)));
827 double a21;
828 20219 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)));
829 double a22;
830 20219 a22 = (2 * ((((((p2_0_p0_0 * q2_0_p0_0) + (p2_1_p0_1 * q2_1_p0_1)) + (p2_2_p0_2 * q2_2_p0_2)) + (p2_3_p0_3 * q2_3_p0_3)) + (p2_4_p0_4 * q2_4_p0_4)) + (p2_5_p0_5 * q2_5_p0_5)));
831 double a23;
832 20219 a23 = (2 * ((((((p2_0_p0_0 * q3_0_p0_0) + (p2_1_p0_1 * q3_1_p0_1)) + (p2_2_p0_2 * q3_2_p0_2)) + (p2_3_p0_3 * q3_3_p0_3)) + (p2_4_p0_4 * q3_4_p0_4)) + (p2_5_p0_5 * q3_5_p0_5)));
833 double a30;
834 20219 a30 = (2 * ((((((p3_0_p0_0 * q0_0_p0_0) + (p3_1_p0_1 * q0_1_p0_1)) + (p3_2_p0_2 * q0_2_p0_2)) + (p3_3_p0_3 * q0_3_p0_3)) + (p3_4_p0_4 * q0_4_p0_4)) + (p3_5_p0_5 * q0_5_p0_5)));
835 double a31;
836 20219 a31 = (2 * ((((((p3_0_p0_0 * q1_0_p0_0) + (p3_1_p0_1 * q1_1_p0_1)) + (p3_2_p0_2 * q1_2_p0_2)) + (p3_3_p0_3 * q1_3_p0_3)) + (p3_4_p0_4 * q1_4_p0_4)) + (p3_5_p0_5 * q1_5_p0_5)));
837 double a32;
838 20219 a32 = (2 * ((((((p3_0_p0_0 * q2_0_p0_0) + (p3_1_p0_1 * q2_1_p0_1)) + (p3_2_p0_2 * q2_2_p0_2)) + (p3_3_p0_3 * q2_3_p0_3)) + (p3_4_p0_4 * q2_4_p0_4)) + (p3_5_p0_5 * q2_5_p0_5)));
839 double a33;
840 20219 a33 = (2 * ((((((p3_0_p0_0 * q3_0_p0_0) + (p3_1_p0_1 * q3_1_p0_1)) + (p3_2_p0_2 * q3_2_p0_2)) + (p3_3_p0_3 * q3_3_p0_3)) + (p3_4_p0_4 * q3_4_p0_4)) + (p3_5_p0_5 * q3_5_p0_5)));
841 double a40;
842 20219 a40 = (2 * ((((((p4_0_p0_0 * q0_0_p0_0) + (p4_1_p0_1 * q0_1_p0_1)) + (p4_2_p0_2 * q0_2_p0_2)) + (p4_3_p0_3 * q0_3_p0_3)) + (p4_4_p0_4 * q0_4_p0_4)) + (p4_5_p0_5 * q0_5_p0_5)));
843 double a41;
844 20219 a41 = (2 * ((((((p4_0_p0_0 * q1_0_p0_0) + (p4_1_p0_1 * q1_1_p0_1)) + (p4_2_p0_2 * q1_2_p0_2)) + (p4_3_p0_3 * q1_3_p0_3)) + (p4_4_p0_4 * q1_4_p0_4)) + (p4_5_p0_5 * q1_5_p0_5)));
845 double a42;
846 20219 a42 = (2 * ((((((p4_0_p0_0 * q2_0_p0_0) + (p4_1_p0_1 * q2_1_p0_1)) + (p4_2_p0_2 * q2_2_p0_2)) + (p4_3_p0_3 * q2_3_p0_3)) + (p4_4_p0_4 * q2_4_p0_4)) + (p4_5_p0_5 * q2_5_p0_5)));
847 double a43;
848 20219 a43 = (2 * ((((((p4_0_p0_0 * q3_0_p0_0) + (p4_1_p0_1 * q3_1_p0_1)) + (p4_2_p0_2 * q3_2_p0_2)) + (p4_3_p0_3 * q3_3_p0_3)) + (p4_4_p0_4 * q3_4_p0_4)) + (p4_5_p0_5 * q3_5_p0_5)));
849 double b00;
850 20219 b00 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
851 double b01;
852 20219 b01 = -((((a22 * a33) - (a23 * a32)) + ((a23 * a31) - (a21 * a33))) + ((a21 * a32) - (a22 * a31)));
853 double b02;
854 20219 b02 = ((((a12 * a33) - (a13 * a32)) + ((a13 * a31) - (a11 * a33))) + ((a11 * a32) - (a12 * a31)));
855 double b03;
856 20219 b03 = -((((a12 * a23) - (a13 * a22)) + ((a13 * a21) - (a11 * a23))) + ((a11 * a22) - (a12 * a21)));
857 double b10;
858 20219 b10 = -(((a10 * ((a22 * a33) - (a23 * a32))) - (a20 * ((a12 * a33) - (a13 * a32)))) + (a30 * ((a12 * a23) - (a13 * a22))));
859 double b11;
860 20219 b11 = ((((a22 * a33) - (a23 * a32)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a32) - (a22 * a30)));
861 double b12;
862 20219 b12 = -((((a12 * a33) - (a13 * a32)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a32) - (a12 * a30)));
863 double b13;
864 20219 b13 = ((((a12 * a23) - (a13 * a22)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a22) - (a12 * a20)));
865 double b20;
866 20219 b20 = (((a10 * ((a21 * a33) - (a23 * a31))) - (a20 * ((a11 * a33) - (a13 * a31)))) + (a30 * ((a11 * a23) - (a13 * a21))));
867 double b21;
868 20219 b21 = -((((a21 * a33) - (a23 * a31)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a31) - (a21 * a30)));
869 double b22;
870 20219 b22 = ((((a11 * a33) - (a13 * a31)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a31) - (a11 * a30)));
871 double b23;
872 20219 b23 = -((((a11 * a23) - (a13 * a21)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a21) - (a11 * a20)));
873 double b30;
874 20219 b30 = -(((a10 * ((a21 * a32) - (a22 * a31))) - (a20 * ((a11 * a32) - (a12 * a31)))) + (a30 * ((a11 * a22) - (a12 * a21))));
875 double b31;
876 20219 b31 = ((((a21 * a32) - (a22 * a31)) + ((a22 * a30) - (a20 * a32))) + ((a20 * a31) - (a21 * a30)));
877 double b32;
878 20219 b32 = -((((a11 * a32) - (a12 * a31)) + ((a12 * a30) - (a10 * a32))) + ((a10 * a31) - (a11 * a30)));
879 double b33;
880 20219 b33 = ((((a11 * a22) - (a12 * a21)) + ((a12 * a20) - (a10 * a22))) + ((a10 * a21) - (a11 * a20)));
881 double Delta;
882 20219 Delta = (((b00 + b10) + b20) + b30);
883 double DeltaLambda0;
884 20219 DeltaLambda0 = ((((b01 * l1) + (b02 * l2)) + (b03 * l3)) + b00);
885 double DeltaLambda1;
886 20219 DeltaLambda1 = ((((b11 * l1) + (b12 * l2)) + (b13 * l3)) + b10);
887 double DeltaLambda2;
888 20219 DeltaLambda2 = ((((b21 * l1) + (b22 * l2)) + (b23 * l3)) + b20);
889 double DeltaLambda3;
890 20219 DeltaLambda3 = ((((b31 * l1) + (b32 * l2)) + (b33 * l3)) + b30);
891 double r;
892 20219 r = ((Delta * l4) - ((((a40 * DeltaLambda0) + (a41 * DeltaLambda1)) + (a42 * DeltaLambda2)) + (a43 * DeltaLambda3)));
893 double eps;
894 20219 double max1 = fabs(p3_2_p0_2);
895
2/2
✓ Branch 0 taken 4879 times.
✓ Branch 1 taken 15340 times.
20219 if( (max1 < fabs(p3_0_p0_0)) )
896 {
897 4879 max1 = fabs(p3_0_p0_0);
898 }
899
2/2
✓ Branch 0 taken 7767 times.
✓ Branch 1 taken 12452 times.
20219 if( (max1 < fabs(p3_3_p0_3)) )
900 {
901 7767 max1 = fabs(p3_3_p0_3);
902 }
903
2/2
✓ Branch 0 taken 6442 times.
✓ Branch 1 taken 13777 times.
20219 if( (max1 < fabs(p3_4_p0_4)) )
904 {
905 6442 max1 = fabs(p3_4_p0_4);
906 }
907
2/2
✓ Branch 0 taken 8343 times.
✓ Branch 1 taken 11876 times.
20219 if( (max1 < fabs(p3_1_p0_1)) )
908 {
909 8343 max1 = fabs(p3_1_p0_1);
910 }
911
2/2
✓ Branch 0 taken 699 times.
✓ Branch 1 taken 19520 times.
20219 if( (max1 < fabs(p3_5_p0_5)) )
912 {
913 699 max1 = fabs(p3_5_p0_5);
914 }
915 20219 double max2 = fabs(p2_1_p0_1);
916
2/2
✓ Branch 0 taken 9114 times.
✓ Branch 1 taken 11105 times.
20219 if( (max2 < fabs(p2_4_p0_4)) )
917 {
918 9114 max2 = fabs(p2_4_p0_4);
919 }
920
2/2
✓ Branch 0 taken 8302 times.
✓ Branch 1 taken 11917 times.
20219 if( (max2 < fabs(p2_2_p0_2)) )
921 {
922 8302 max2 = fabs(p2_2_p0_2);
923 }
924
2/2
✓ Branch 0 taken 79 times.
✓ Branch 1 taken 20140 times.
20219 if( (max2 < fabs(p2_0_p0_0)) )
925 {
926 79 max2 = fabs(p2_0_p0_0);
927 }
928
2/2
✓ Branch 0 taken 1022 times.
✓ Branch 1 taken 19197 times.
20219 if( (max2 < fabs(p2_3_p0_3)) )
929 {
930 1022 max2 = fabs(p2_3_p0_3);
931 }
932
2/2
✓ Branch 0 taken 738 times.
✓ Branch 1 taken 19481 times.
20219 if( (max2 < fabs(p2_5_p0_5)) )
933 {
934 738 max2 = fabs(p2_5_p0_5);
935 }
936 20219 double max3 = fabs(p1_0_p0_0);
937
2/2
✓ Branch 0 taken 15006 times.
✓ Branch 1 taken 5213 times.
20219 if( (max3 < fabs(p1_1_p0_1)) )
938 {
939 15006 max3 = fabs(p1_1_p0_1);
940 }
941
2/2
✓ Branch 0 taken 9005 times.
✓ Branch 1 taken 11214 times.
20219 if( (max3 < fabs(p1_2_p0_2)) )
942 {
943 9005 max3 = fabs(p1_2_p0_2);
944 }
945
2/2
✓ Branch 0 taken 2670 times.
✓ Branch 1 taken 17549 times.
20219 if( (max3 < fabs(p1_3_p0_3)) )
946 {
947 2670 max3 = fabs(p1_3_p0_3);
948 }
949
2/2
✓ Branch 0 taken 2701 times.
✓ Branch 1 taken 17518 times.
20219 if( (max3 < fabs(p1_4_p0_4)) )
950 {
951 2701 max3 = fabs(p1_4_p0_4);
952 }
953
2/2
✓ Branch 0 taken 700 times.
✓ Branch 1 taken 19519 times.
20219 if( (max3 < fabs(p1_5_p0_5)) )
954 {
955 700 max3 = fabs(p1_5_p0_5);
956 }
957 20219 double max4 = fabs(q0_0_p0_0);
958
2/2
✓ Branch 0 taken 13324 times.
✓ Branch 1 taken 6895 times.
20219 if( (max4 < fabs(q0_1_p0_1)) )
959 {
960 13324 max4 = fabs(q0_1_p0_1);
961 }
962
2/2
✓ Branch 0 taken 7658 times.
✓ Branch 1 taken 12561 times.
20219 if( (max4 < fabs(q0_2_p0_2)) )
963 {
964 7658 max4 = fabs(q0_2_p0_2);
965 }
966
2/2
✓ Branch 0 taken 14465 times.
✓ Branch 1 taken 5754 times.
20219 if( (max4 < fabs(q0_3_p0_3)) )
967 {
968 14465 max4 = fabs(q0_3_p0_3);
969 }
970
2/2
✓ Branch 0 taken 14703 times.
✓ Branch 1 taken 5516 times.
20219 if( (max4 < fabs(q0_4_p0_4)) )
971 {
972 14703 max4 = fabs(q0_4_p0_4);
973 }
974
2/2
✓ Branch 0 taken 3123 times.
✓ Branch 1 taken 17096 times.
20219 if( (max4 < fabs(q0_5_p0_5)) )
975 {
976 3123 max4 = fabs(q0_5_p0_5);
977 }
978
2/2
✓ Branch 0 taken 1704 times.
✓ Branch 1 taken 18515 times.
20219 if( (max4 < fabs(q1_0_p0_0)) )
979 {
980 1704 max4 = fabs(q1_0_p0_0);
981 }
982
2/2
✓ Branch 0 taken 1419 times.
✓ Branch 1 taken 18800 times.
20219 if( (max4 < fabs(q1_1_p0_1)) )
983 {
984 1419 max4 = fabs(q1_1_p0_1);
985 }
986
2/2
✓ Branch 0 taken 1044 times.
✓ Branch 1 taken 19175 times.
20219 if( (max4 < fabs(q1_2_p0_2)) )
987 {
988 1044 max4 = fabs(q1_2_p0_2);
989 }
990
2/2
✓ Branch 0 taken 5046 times.
✓ Branch 1 taken 15173 times.
20219 if( (max4 < fabs(q1_3_p0_3)) )
991 {
992 5046 max4 = fabs(q1_3_p0_3);
993 }
994
2/2
✓ Branch 0 taken 5930 times.
✓ Branch 1 taken 14289 times.
20219 if( (max4 < fabs(q1_4_p0_4)) )
995 {
996 5930 max4 = fabs(q1_4_p0_4);
997 }
998
2/2
✓ Branch 0 taken 2358 times.
✓ Branch 1 taken 17861 times.
20219 if( (max4 < fabs(q1_5_p0_5)) )
999 {
1000 2358 max4 = fabs(q1_5_p0_5);
1001 }
1002 20219 double max5 = fabs(q1_0_p0_0);
1003
2/2
✓ Branch 0 taken 13116 times.
✓ Branch 1 taken 7103 times.
20219 if( (max5 < fabs(q1_1_p0_1)) )
1004 {
1005 13116 max5 = fabs(q1_1_p0_1);
1006 }
1007
2/2
✓ Branch 0 taken 7415 times.
✓ Branch 1 taken 12804 times.
20219 if( (max5 < fabs(q1_2_p0_2)) )
1008 {
1009 7415 max5 = fabs(q1_2_p0_2);
1010 }
1011
2/2
✓ Branch 0 taken 14176 times.
✓ Branch 1 taken 6043 times.
20219 if( (max5 < fabs(q1_3_p0_3)) )
1012 {
1013 14176 max5 = fabs(q1_3_p0_3);
1014 }
1015
2/2
✓ Branch 0 taken 13754 times.
✓ Branch 1 taken 6465 times.
20219 if( (max5 < fabs(q1_4_p0_4)) )
1016 {
1017 13754 max5 = fabs(q1_4_p0_4);
1018 }
1019
2/2
✓ Branch 0 taken 3231 times.
✓ Branch 1 taken 16988 times.
20219 if( (max5 < fabs(q1_5_p0_5)) )
1020 {
1021 3231 max5 = fabs(q1_5_p0_5);
1022 }
1023
2/2
✓ Branch 0 taken 1659 times.
✓ Branch 1 taken 18560 times.
20219 if( (max5 < fabs(q2_0_p0_0)) )
1024 {
1025 1659 max5 = fabs(q2_0_p0_0);
1026 }
1027
2/2
✓ Branch 0 taken 1087 times.
✓ Branch 1 taken 19132 times.
20219 if( (max5 < fabs(q2_1_p0_1)) )
1028 {
1029 1087 max5 = fabs(q2_1_p0_1);
1030 }
1031
2/2
✓ Branch 0 taken 728 times.
✓ Branch 1 taken 19491 times.
20219 if( (max5 < fabs(q2_2_p0_2)) )
1032 {
1033 728 max5 = fabs(q2_2_p0_2);
1034 }
1035
2/2
✓ Branch 0 taken 5384 times.
✓ Branch 1 taken 14835 times.
20219 if( (max5 < fabs(q2_3_p0_3)) )
1036 {
1037 5384 max5 = fabs(q2_3_p0_3);
1038 }
1039
2/2
✓ Branch 0 taken 7591 times.
✓ Branch 1 taken 12628 times.
20219 if( (max5 < fabs(q2_4_p0_4)) )
1040 {
1041 7591 max5 = fabs(q2_4_p0_4);
1042 }
1043
2/2
✓ Branch 0 taken 1423 times.
✓ Branch 1 taken 18796 times.
20219 if( (max5 < fabs(q2_5_p0_5)) )
1044 {
1045 1423 max5 = fabs(q2_5_p0_5);
1046 }
1047 20219 double max6 = fabs(q2_0_p0_0);
1048
2/2
✓ Branch 0 taken 13254 times.
✓ Branch 1 taken 6965 times.
20219 if( (max6 < fabs(q2_1_p0_1)) )
1049 {
1050 13254 max6 = fabs(q2_1_p0_1);
1051 }
1052
2/2
✓ Branch 0 taken 7621 times.
✓ Branch 1 taken 12598 times.
20219 if( (max6 < fabs(q2_2_p0_2)) )
1053 {
1054 7621 max6 = fabs(q2_2_p0_2);
1055 }
1056
2/2
✓ Branch 0 taken 14378 times.
✓ Branch 1 taken 5841 times.
20219 if( (max6 < fabs(q2_3_p0_3)) )
1057 {
1058 14378 max6 = fabs(q2_3_p0_3);
1059 }
1060
2/2
✓ Branch 0 taken 14733 times.
✓ Branch 1 taken 5486 times.
20219 if( (max6 < fabs(q2_4_p0_4)) )
1061 {
1062 14733 max6 = fabs(q2_4_p0_4);
1063 }
1064
2/2
✓ Branch 0 taken 2995 times.
✓ Branch 1 taken 17224 times.
20219 if( (max6 < fabs(q2_5_p0_5)) )
1065 {
1066 2995 max6 = fabs(q2_5_p0_5);
1067 }
1068
2/2
✓ Branch 0 taken 2537 times.
✓ Branch 1 taken 17682 times.
20219 if( (max6 < fabs(q3_0_p0_0)) )
1069 {
1070 2537 max6 = fabs(q3_0_p0_0);
1071 }
1072
2/2
✓ Branch 0 taken 1502 times.
✓ Branch 1 taken 18717 times.
20219 if( (max6 < fabs(q3_1_p0_1)) )
1073 {
1074 1502 max6 = fabs(q3_1_p0_1);
1075 }
1076
2/2
✓ Branch 0 taken 801 times.
✓ Branch 1 taken 19418 times.
20219 if( (max6 < fabs(q3_2_p0_2)) )
1077 {
1078 801 max6 = fabs(q3_2_p0_2);
1079 }
1080
2/2
✓ Branch 0 taken 6940 times.
✓ Branch 1 taken 13279 times.
20219 if( (max6 < fabs(q3_3_p0_3)) )
1081 {
1082 6940 max6 = fabs(q3_3_p0_3);
1083 }
1084
2/2
✓ Branch 0 taken 7100 times.
✓ Branch 1 taken 13119 times.
20219 if( (max6 < fabs(q3_4_p0_4)) )
1085 {
1086 7100 max6 = fabs(q3_4_p0_4);
1087 }
1088
2/2
✓ Branch 0 taken 2514 times.
✓ Branch 1 taken 17705 times.
20219 if( (max6 < fabs(q3_5_p0_5)) )
1089 {
1090 2514 max6 = fabs(q3_5_p0_5);
1091 }
1092 double lower_bound_1;
1093 double upper_bound_1;
1094 int Delta_sign;
1095 int int_tmp_result;
1096 20219 lower_bound_1 = max4;
1097 20219 upper_bound_1 = max4;
1098
2/2
✓ Branch 0 taken 19103 times.
✓ Branch 1 taken 1116 times.
20219 if( (max1 < lower_bound_1) )
1099 {
1100 19103 lower_bound_1 = max1;
1101 }
1102 else
1103 {
1104
1/2
✓ Branch 0 taken 1116 times.
✗ Branch 1 not taken.
1116 if( (max1 > upper_bound_1) )
1105 {
1106 1116 upper_bound_1 = max1;
1107 }
1108 }
1109
2/2
✓ Branch 0 taken 9606 times.
✓ Branch 1 taken 10613 times.
20219 if( (max2 < lower_bound_1) )
1110 {
1111 9606 lower_bound_1 = max2;
1112 }
1113 else
1114 {
1115
2/2
✓ Branch 0 taken 714 times.
✓ Branch 1 taken 9899 times.
10613 if( (max2 > upper_bound_1) )
1116 {
1117 714 upper_bound_1 = max2;
1118 }
1119 }
1120
2/2
✓ Branch 0 taken 6408 times.
✓ Branch 1 taken 13811 times.
20219 if( (max3 < lower_bound_1) )
1121 {
1122 6408 lower_bound_1 = max3;
1123 }
1124 else
1125 {
1126
2/2
✓ Branch 0 taken 536 times.
✓ Branch 1 taken 13275 times.
13811 if( (max3 > upper_bound_1) )
1127 {
1128 536 upper_bound_1 = max3;
1129 }
1130 }
1131
2/2
✓ Branch 0 taken 485 times.
✓ Branch 1 taken 19734 times.
20219 if( (max5 < lower_bound_1) )
1132 {
1133 485 lower_bound_1 = max5;
1134 }
1135 else
1136 {
1137
2/2
✓ Branch 0 taken 6541 times.
✓ Branch 1 taken 13193 times.
19734 if( (max5 > upper_bound_1) )
1138 {
1139 6541 upper_bound_1 = max5;
1140 }
1141 }
1142
2/2
✓ Branch 0 taken 390 times.
✓ Branch 1 taken 19829 times.
20219 if( (max6 < lower_bound_1) )
1143 {
1144 390 lower_bound_1 = max6;
1145 }
1146 else
1147 {
1148
2/2
✓ Branch 0 taken 5277 times.
✓ Branch 1 taken 14552 times.
19829 if( (max6 > upper_bound_1) )
1149 {
1150 5277 upper_bound_1 = max6;
1151 }
1152 }
1153
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20219 times.
20219 if( (lower_bound_1 < 3.31864264949884013629e-50) )
1154 {
1155 return FPG_UNCERTAIN_VALUE;
1156 }
1157 else
1158 {
1159
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20219 times.
20219 if( (upper_bound_1 > 4.83570327845851562508e+24) )
1160 {
1161 return FPG_UNCERTAIN_VALUE;
1162 }
1163 20219 eps = (1.66564133587113165316e-11 * (((((max3 * max4) * max2) * max5) * max1) * max6));
1164
2/2
✓ Branch 0 taken 10105 times.
✓ Branch 1 taken 10114 times.
20219 if( (Delta > eps) )
1165 {
1166 10105 int_tmp_result = 1;
1167 }
1168 else
1169 {
1170
1/2
✓ Branch 0 taken 10114 times.
✗ Branch 1 not taken.
10114 if( (Delta < -eps) )
1171 {
1172 10114 int_tmp_result = -1;
1173 }
1174 else
1175 {
1176 return FPG_UNCERTAIN_VALUE;
1177 }
1178 }
1179 }
1180 20219 Delta_sign = int_tmp_result;
1181 int int_tmp_result_FFWKCAA;
1182 20219 double max7 = max1;
1183
2/2
✓ Branch 0 taken 10198 times.
✓ Branch 1 taken 10021 times.
20219 if( (max7 < max2) )
1184 {
1185 10198 max7 = max2;
1186 }
1187
2/2
✓ Branch 0 taken 6898 times.
✓ Branch 1 taken 13321 times.
20219 if( (max7 < max3) )
1188 {
1189 6898 max7 = max3;
1190 }
1191 20219 double max8 = max3;
1192
2/2
✓ Branch 0 taken 1013 times.
✓ Branch 1 taken 19206 times.
20219 if( (max8 < fabs(p4_5_p0_5)) )
1193 {
1194 1013 max8 = fabs(p4_5_p0_5);
1195 }
1196
2/2
✓ Branch 0 taken 43 times.
✓ Branch 1 taken 20176 times.
20219 if( (max8 < fabs(p4_0_p0_0)) )
1197 {
1198 43 max8 = fabs(p4_0_p0_0);
1199 }
1200
2/2
✓ Branch 0 taken 6897 times.
✓ Branch 1 taken 13322 times.
20219 if( (max8 < fabs(p4_1_p0_1)) )
1201 {
1202 6897 max8 = fabs(p4_1_p0_1);
1203 }
1204
2/2
✓ Branch 0 taken 6221 times.
✓ Branch 1 taken 13998 times.
20219 if( (max8 < fabs(p4_2_p0_2)) )
1205 {
1206 6221 max8 = fabs(p4_2_p0_2);
1207 }
1208
2/2
✓ Branch 0 taken 759 times.
✓ Branch 1 taken 19460 times.
20219 if( (max8 < fabs(p4_3_p0_3)) )
1209 {
1210 759 max8 = fabs(p4_3_p0_3);
1211 }
1212
2/2
✓ Branch 0 taken 2382 times.
✓ Branch 1 taken 17837 times.
20219 if( (max8 < fabs(p4_4_p0_4)) )
1213 {
1214 2382 max8 = fabs(p4_4_p0_4);
1215 }
1216
2/2
✓ Branch 0 taken 10652 times.
✓ Branch 1 taken 9567 times.
20219 if( (max7 < max8) )
1217 {
1218 10652 max7 = max8;
1219 }
1220
2/2
✓ Branch 0 taken 18775 times.
✓ Branch 1 taken 1444 times.
20219 if( (max7 < max6) )
1221 {
1222 18775 max7 = max6;
1223 }
1224 20219 double max9 = max3;
1225
2/2
✓ Branch 0 taken 14831 times.
✓ Branch 1 taken 5388 times.
20219 if( (max9 < max8) )
1226 {
1227 14831 max9 = max8;
1228 }
1229
2/2
✓ Branch 0 taken 18081 times.
✓ Branch 1 taken 2138 times.
20219 if( (max9 < max5) )
1230 {
1231 18081 max9 = max5;
1232 }
1233 20219 double max10 = max4;
1234
2/2
✓ Branch 0 taken 1117 times.
✓ Branch 1 taken 19102 times.
20219 if( (max10 < max3) )
1235 {
1236 1117 max10 = max3;
1237 }
1238
2/2
✓ Branch 0 taken 6609 times.
✓ Branch 1 taken 13610 times.
20219 if( (max10 < max5) )
1239 {
1240 6609 max10 = max5;
1241 }
1242 20219 double max11 = max4;
1243
2/2
✓ Branch 0 taken 6690 times.
✓ Branch 1 taken 13529 times.
20219 if( (max11 < max5) )
1244 {
1245 6690 max11 = max5;
1246 }
1247
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20219 times.
20219 if( (max10 < max11) )
1248 {
1249 max10 = max11;
1250 }
1251
2/2
✓ Branch 0 taken 5277 times.
✓ Branch 1 taken 14942 times.
20219 if( (max10 < max6) )
1252 {
1253 5277 max10 = max6;
1254 }
1255 20219 lower_bound_1 = max1;
1256 20219 upper_bound_1 = max1;
1257
2/2
✓ Branch 0 taken 10021 times.
✓ Branch 1 taken 10198 times.
20219 if( (max2 < lower_bound_1) )
1258 {
1259 10021 lower_bound_1 = max2;
1260 }
1261
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20219 times.
20219 if( (max7 < lower_bound_1) )
1262 {
1263 lower_bound_1 = max7;
1264 }
1265 else
1266 {
1267
2/2
✓ Branch 0 taken 20037 times.
✓ Branch 1 taken 182 times.
20219 if( (max7 > upper_bound_1) )
1268 {
1269 20037 upper_bound_1 = max7;
1270 }
1271 }
1272
2/2
✓ Branch 0 taken 1674 times.
✓ Branch 1 taken 18545 times.
20219 if( (max8 < lower_bound_1) )
1273 {
1274 1674 lower_bound_1 = max8;
1275 }
1276
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20219 times.
20219 if( (max9 < lower_bound_1) )
1277 {
1278 lower_bound_1 = max9;
1279 }
1280 else
1281 {
1282
2/2
✓ Branch 0 taken 5975 times.
✓ Branch 1 taken 14244 times.
20219 if( (max9 > upper_bound_1) )
1283 {
1284 5975 upper_bound_1 = max9;
1285 }
1286 }
1287
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20219 times.
20219 if( (max10 < lower_bound_1) )
1288 {
1289 lower_bound_1 = max10;
1290 }
1291 else
1292 {
1293
2/2
✓ Branch 0 taken 4864 times.
✓ Branch 1 taken 15355 times.
20219 if( (max10 > upper_bound_1) )
1294 {
1295 4864 upper_bound_1 = max10;
1296 }
1297 }
1298
2/2
✓ Branch 0 taken 229 times.
✓ Branch 1 taken 19990 times.
20219 if( (max11 < lower_bound_1) )
1299 {
1300 229 lower_bound_1 = max11;
1301 }
1302 else
1303 {
1304
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19990 times.
19990 if( (max11 > upper_bound_1) )
1305 {
1306 upper_bound_1 = max11;
1307 }
1308 }
1309
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20219 times.
20219 if( (lower_bound_1 < 4.87975611107819181771e-38) )
1310 {
1311 return FPG_UNCERTAIN_VALUE;
1312 }
1313 else
1314 {
1315
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 20219 times.
20219 if( (upper_bound_1 > 4.83570327845851562508e+24) )
1316 {
1317 return FPG_UNCERTAIN_VALUE;
1318 }
1319 20219 eps = (6.92085098542795335117e-10 * (((((((max8 * max11) * max2) * max10) * max1) * max10) * max9) * max7));
1320
2/2
✓ Branch 0 taken 10105 times.
✓ Branch 1 taken 10114 times.
20219 if( (r > eps) )
1321 {
1322 10105 int_tmp_result_FFWKCAA = 1;
1323 }
1324 else
1325 {
1326
1/2
✓ Branch 0 taken 10114 times.
✗ Branch 1 not taken.
10114 if( (r < -eps) )
1327 {
1328 10114 int_tmp_result_FFWKCAA = -1;
1329 }
1330 else
1331 {
1332 return FPG_UNCERTAIN_VALUE;
1333 }
1334 }
1335 }
1336 20219 return (Delta_sign * int_tmp_result_FFWKCAA);
1337 }
1338
1339
1340 inline int side4_7d_filter( const double* p0, const double* p1, const double* p2, const double* p3, const double* p4, const double* q0, const double* q1, const double* q2, const double* q3) {
1341 double p1_0_p0_0 = (p1[0] - p0[0]);
1342 double p1_1_p0_1 = (p1[1] - p0[1]);
1343 double p1_2_p0_2 = (p1[2] - p0[2]);
1344 double p1_3_p0_3 = (p1[3] - p0[3]);
1345 double p1_4_p0_4 = (p1[4] - p0[4]);
1346 double p1_5_p0_5 = (p1[5] - p0[5]);
1347 double p1_6_p0_6 = (p1[6] - p0[6]);
1348 double l1;
1349 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)));
1350 double p2_0_p0_0 = (p2[0] - p0[0]);
1351 double p2_1_p0_1 = (p2[1] - p0[1]);
1352 double p2_2_p0_2 = (p2[2] - p0[2]);
1353 double p2_3_p0_3 = (p2[3] - p0[3]);
1354 double p2_4_p0_4 = (p2[4] - p0[4]);
1355 double p2_5_p0_5 = (p2[5] - p0[5]);
1356 double p2_6_p0_6 = (p2[6] - p0[6]);
1357 double l2;
1358 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)));
1359 double p3_0_p0_0 = (p3[0] - p0[0]);
1360 double p3_1_p0_1 = (p3[1] - p0[1]);
1361 double p3_2_p0_2 = (p3[2] - p0[2]);
1362 double p3_3_p0_3 = (p3[3] - p0[3]);
1363 double p3_4_p0_4 = (p3[4] - p0[4]);
1364 double p3_5_p0_5 = (p3[5] - p0[5]);
1365 double p3_6_p0_6 = (p3[6] - p0[6]);
1366 double l3;
1367 l3 = (1 * (((((((p3_0_p0_0 * p3_0_p0_0) + (p3_1_p0_1 * p3_1_p0_1)) + (p3_2_p0_2 * p3_2_p0_2)) + (p3_3_p0_3 * p3_3_p0_3)) + (p3_4_p0_4 * p3_4_p0_4)) + (p3_5_p0_5 * p3_5_p0_5)) + (p3_6_p0_6 * p3_6_p0_6)));
1368 double p4_0_p0_0 = (p4[0] - p0[0]);
1369 double p4_1_p0_1 = (p4[1] - p0[1]);
1370 double p4_2_p0_2 = (p4[2] - p0[2]);
1371 double p4_3_p0_3 = (p4[3] - p0[3]);
1372 double p4_4_p0_4 = (p4[4] - p0[4]);
1373 double p4_5_p0_5 = (p4[5] - p0[5]);
1374 double p4_6_p0_6 = (p4[6] - p0[6]);
1375 double l4;
1376 l4 = (1 * (((((((p4_0_p0_0 * p4_0_p0_0) + (p4_1_p0_1 * p4_1_p0_1)) + (p4_2_p0_2 * p4_2_p0_2)) + (p4_3_p0_3 * p4_3_p0_3)) + (p4_4_p0_4 * p4_4_p0_4)) + (p4_5_p0_5 * p4_5_p0_5)) + (p4_6_p0_6 * p4_6_p0_6)));
1377 double q0_0_p0_0 = (q0[0] - p0[0]);
1378 double q0_1_p0_1 = (q0[1] - p0[1]);
1379 double q0_2_p0_2 = (q0[2] - p0[2]);
1380 double q0_3_p0_3 = (q0[3] - p0[3]);
1381 double q0_4_p0_4 = (q0[4] - p0[4]);
1382 double q0_5_p0_5 = (q0[5] - p0[5]);
1383 double q0_6_p0_6 = (q0[6] - p0[6]);
1384 double a10;
1385 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)));
1386 double q1_0_p0_0 = (q1[0] - p0[0]);
1387 double q1_1_p0_1 = (q1[1] - p0[1]);
1388 double q1_2_p0_2 = (q1[2] - p0[2]);
1389 double q1_3_p0_3 = (q1[3] - p0[3]);
1390 double q1_4_p0_4 = (q1[4] - p0[4]);
1391 double q1_5_p0_5 = (q1[5] - p0[5]);
1392 double q1_6_p0_6 = (q1[6] - p0[6]);
1393 double a11;
1394 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)));
1395 double q2_0_p0_0 = (q2[0] - p0[0]);
1396 double q2_1_p0_1 = (q2[1] - p0[1]);
1397 double q2_2_p0_2 = (q2[2] - p0[2]);
1398 double q2_3_p0_3 = (q2[3] - p0[3]);
1399 double q2_4_p0_4 = (q2[4] - p0[4]);
1400 double q2_5_p0_5 = (q2[5] - p0[5]);
1401 double q2_6_p0_6 = (q2[6] - p0[6]);
1402 double a12;
1403 a12 = (2 * (((((((p1_0_p0_0 * q2_0_p0_0) + (p1_1_p0_1 * q2_1_p0_1)) + (p1_2_p0_2 * q2_2_p0_2)) + (p1_3_p0_3 * q2_3_p0_3)) + (p1_4_p0_4 * q2_4_p0_4)) + (p1_5_p0_5 * q2_5_p0_5)) + (p1_6_p0_6 * q2_6_p0_6)));
1404 double q3_0_p0_0 = (q3[0] - p0[0]);
1405 double q3_1_p0_1 = (q3[1] - p0[1]);
1406 double q3_2_p0_2 = (q3[2] - p0[2]);
1407 double q3_3_p0_3 = (q3[3] - p0[3]);
1408 double q3_4_p0_4 = (q3[4] - p0[4]);
1409 double q3_5_p0_5 = (q3[5] - p0[5]);
1410 double q3_6_p0_6 = (q3[6] - p0[6]);
1411 double a13;
1412 a13 = (2 * (((((((p1_0_p0_0 * q3_0_p0_0) + (p1_1_p0_1 * q3_1_p0_1)) + (p1_2_p0_2 * q3_2_p0_2)) + (p1_3_p0_3 * q3_3_p0_3)) + (p1_4_p0_4 * q3_4_p0_4)) + (p1_5_p0_5 * q3_5_p0_5)) + (p1_6_p0_6 * q3_6_p0_6)));
1413 double a20;
1414 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)));
1415 double a21;
1416 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)));
1417 double a22;
1418 a22 = (2 * (((((((p2_0_p0_0 * q2_0_p0_0) + (p2_1_p0_1 * q2_1_p0_1)) + (p2_2_p0_2 * q2_2_p0_2)) + (p2_3_p0_3 * q2_3_p0_3)) + (p2_4_p0_4 * q2_4_p0_4)) + (p2_5_p0_5 * q2_5_p0_5)) + (p2_6_p0_6 * q2_6_p0_6)));
1419 double a23;
1420 a23 = (2 * (((((((p2_0_p0_0 * q3_0_p0_0) + (p2_1_p0_1 * q3_1_p0_1)) + (p2_2_p0_2 * q3_2_p0_2)) + (p2_3_p0_3 * q3_3_p0_3)) + (p2_4_p0_4 * q3_4_p0_4)) + (p2_5_p0_5 * q3_5_p0_5)) + (p2_6_p0_6 * q3_6_p0_6)));
1421 double a30;
1422 a30 = (2 * (((((((p3_0_p0_0 * q0_0_p0_0) + (p3_1_p0_1 * q0_1_p0_1)) + (p3_2_p0_2 * q0_2_p0_2)) + (p3_3_p0_3 * q0_3_p0_3)) + (p3_4_p0_4 * q0_4_p0_4)) + (p3_5_p0_5 * q0_5_p0_5)) + (p3_6_p0_6 * q0_6_p0_6)));
1423 double a31;
1424 a31 = (2 * (((((((p3_0_p0_0 * q1_0_p0_0) + (p3_1_p0_1 * q1_1_p0_1)) + (p3_2_p0_2 * q1_2_p0_2)) + (p3_3_p0_3 * q1_3_p0_3)) + (p3_4_p0_4 * q1_4_p0_4)) + (p3_5_p0_5 * q1_5_p0_5)) + (p3_6_p0_6 * q1_6_p0_6)));
1425 double a32;
1426 a32 = (2 * (((((((p3_0_p0_0 * q2_0_p0_0) + (p3_1_p0_1 * q2_1_p0_1)) + (p3_2_p0_2 * q2_2_p0_2)) + (p3_3_p0_3 * q2_3_p0_3)) + (p3_4_p0_4 * q2_4_p0_4)) + (p3_5_p0_5 * q2_5_p0_5)) + (p3_6_p0_6 * q2_6_p0_6)));
1427 double a33;
1428 a33 = (2 * (((((((p3_0_p0_0 * q3_0_p0_0) + (p3_1_p0_1 * q3_1_p0_1)) + (p3_2_p0_2 * q3_2_p0_2)) + (p3_3_p0_3 * q3_3_p0_3)) + (p3_4_p0_4 * q3_4_p0_4)) + (p3_5_p0_5 * q3_5_p0_5)) + (p3_6_p0_6 * q3_6_p0_6)));
1429 double a40;
1430 a40 = (2 * (((((((p4_0_p0_0 * q0_0_p0_0) + (p4_1_p0_1 * q0_1_p0_1)) + (p4_2_p0_2 * q0_2_p0_2)) + (p4_3_p0_3 * q0_3_p0_3)) + (p4_4_p0_4 * q0_4_p0_4)) + (p4_5_p0_5 * q0_5_p0_5)) + (p4_6_p0_6 * q0_6_p0_6)));
1431 double a41;
1432 a41 = (2 * (((((((p4_0_p0_0 * q1_0_p0_0) + (p4_1_p0_1 * q1_1_p0_1)) + (p4_2_p0_2 * q1_2_p0_2)) + (p4_3_p0_3 * q1_3_p0_3)) + (p4_4_p0_4 * q1_4_p0_4)) + (p4_5_p0_5 * q1_5_p0_5)) + (p4_6_p0_6 * q1_6_p0_6)));
1433 double a42;
1434 a42 = (2 * (((((((p4_0_p0_0 * q2_0_p0_0) + (p4_1_p0_1 * q2_1_p0_1)) + (p4_2_p0_2 * q2_2_p0_2)) + (p4_3_p0_3 * q2_3_p0_3)) + (p4_4_p0_4 * q2_4_p0_4)) + (p4_5_p0_5 * q2_5_p0_5)) + (p4_6_p0_6 * q2_6_p0_6)));
1435 double a43;
1436 a43 = (2 * (((((((p4_0_p0_0 * q3_0_p0_0) + (p4_1_p0_1 * q3_1_p0_1)) + (p4_2_p0_2 * q3_2_p0_2)) + (p4_3_p0_3 * q3_3_p0_3)) + (p4_4_p0_4 * q3_4_p0_4)) + (p4_5_p0_5 * q3_5_p0_5)) + (p4_6_p0_6 * q3_6_p0_6)));
1437 double b00;
1438 b00 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
1439 double b01;
1440 b01 = -((((a22 * a33) - (a23 * a32)) + ((a23 * a31) - (a21 * a33))) + ((a21 * a32) - (a22 * a31)));
1441 double b02;
1442 b02 = ((((a12 * a33) - (a13 * a32)) + ((a13 * a31) - (a11 * a33))) + ((a11 * a32) - (a12 * a31)));
1443 double b03;
1444 b03 = -((((a12 * a23) - (a13 * a22)) + ((a13 * a21) - (a11 * a23))) + ((a11 * a22) - (a12 * a21)));
1445 double b10;
1446 b10 = -(((a10 * ((a22 * a33) - (a23 * a32))) - (a20 * ((a12 * a33) - (a13 * a32)))) + (a30 * ((a12 * a23) - (a13 * a22))));
1447 double b11;
1448 b11 = ((((a22 * a33) - (a23 * a32)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a32) - (a22 * a30)));
1449 double b12;
1450 b12 = -((((a12 * a33) - (a13 * a32)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a32) - (a12 * a30)));
1451 double b13;
1452 b13 = ((((a12 * a23) - (a13 * a22)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a22) - (a12 * a20)));
1453 double b20;
1454 b20 = (((a10 * ((a21 * a33) - (a23 * a31))) - (a20 * ((a11 * a33) - (a13 * a31)))) + (a30 * ((a11 * a23) - (a13 * a21))));
1455 double b21;
1456 b21 = -((((a21 * a33) - (a23 * a31)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a31) - (a21 * a30)));
1457 double b22;
1458 b22 = ((((a11 * a33) - (a13 * a31)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a31) - (a11 * a30)));
1459 double b23;
1460 b23 = -((((a11 * a23) - (a13 * a21)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a21) - (a11 * a20)));
1461 double b30;
1462 b30 = -(((a10 * ((a21 * a32) - (a22 * a31))) - (a20 * ((a11 * a32) - (a12 * a31)))) + (a30 * ((a11 * a22) - (a12 * a21))));
1463 double b31;
1464 b31 = ((((a21 * a32) - (a22 * a31)) + ((a22 * a30) - (a20 * a32))) + ((a20 * a31) - (a21 * a30)));
1465 double b32;
1466 b32 = -((((a11 * a32) - (a12 * a31)) + ((a12 * a30) - (a10 * a32))) + ((a10 * a31) - (a11 * a30)));
1467 double b33;
1468 b33 = ((((a11 * a22) - (a12 * a21)) + ((a12 * a20) - (a10 * a22))) + ((a10 * a21) - (a11 * a20)));
1469 double Delta;
1470 Delta = (((b00 + b10) + b20) + b30);
1471 double DeltaLambda0;
1472 DeltaLambda0 = ((((b01 * l1) + (b02 * l2)) + (b03 * l3)) + b00);
1473 double DeltaLambda1;
1474 DeltaLambda1 = ((((b11 * l1) + (b12 * l2)) + (b13 * l3)) + b10);
1475 double DeltaLambda2;
1476 DeltaLambda2 = ((((b21 * l1) + (b22 * l2)) + (b23 * l3)) + b20);
1477 double DeltaLambda3;
1478 DeltaLambda3 = ((((b31 * l1) + (b32 * l2)) + (b33 * l3)) + b30);
1479 double r;
1480 r = ((Delta * l4) - ((((a40 * DeltaLambda0) + (a41 * DeltaLambda1)) + (a42 * DeltaLambda2)) + (a43 * DeltaLambda3)));
1481 double eps;
1482 double max1 = fabs(p1_0_p0_0);
1483 if( (max1 < fabs(p1_1_p0_1)) )
1484 {
1485 max1 = fabs(p1_1_p0_1);
1486 }
1487 if( (max1 < fabs(p1_2_p0_2)) )
1488 {
1489 max1 = fabs(p1_2_p0_2);
1490 }
1491 if( (max1 < fabs(p1_3_p0_3)) )
1492 {
1493 max1 = fabs(p1_3_p0_3);
1494 }
1495 if( (max1 < fabs(p1_4_p0_4)) )
1496 {
1497 max1 = fabs(p1_4_p0_4);
1498 }
1499 if( (max1 < fabs(p1_5_p0_5)) )
1500 {
1501 max1 = fabs(p1_5_p0_5);
1502 }
1503 if( (max1 < fabs(p1_6_p0_6)) )
1504 {
1505 max1 = fabs(p1_6_p0_6);
1506 }
1507 double max2 = fabs(p3_0_p0_0);
1508 if( (max2 < fabs(p3_4_p0_4)) )
1509 {
1510 max2 = fabs(p3_4_p0_4);
1511 }
1512 if( (max2 < fabs(p3_2_p0_2)) )
1513 {
1514 max2 = fabs(p3_2_p0_2);
1515 }
1516 if( (max2 < fabs(p3_1_p0_1)) )
1517 {
1518 max2 = fabs(p3_1_p0_1);
1519 }
1520 if( (max2 < fabs(p3_3_p0_3)) )
1521 {
1522 max2 = fabs(p3_3_p0_3);
1523 }
1524 if( (max2 < fabs(p3_5_p0_5)) )
1525 {
1526 max2 = fabs(p3_5_p0_5);
1527 }
1528 if( (max2 < fabs(p3_6_p0_6)) )
1529 {
1530 max2 = fabs(p3_6_p0_6);
1531 }
1532 double max3 = fabs(p2_5_p0_5);
1533 if( (max3 < fabs(p2_2_p0_2)) )
1534 {
1535 max3 = fabs(p2_2_p0_2);
1536 }
1537 if( (max3 < fabs(p2_3_p0_3)) )
1538 {
1539 max3 = fabs(p2_3_p0_3);
1540 }
1541 if( (max3 < fabs(p2_0_p0_0)) )
1542 {
1543 max3 = fabs(p2_0_p0_0);
1544 }
1545 if( (max3 < fabs(p2_1_p0_1)) )
1546 {
1547 max3 = fabs(p2_1_p0_1);
1548 }
1549 if( (max3 < fabs(p2_6_p0_6)) )
1550 {
1551 max3 = fabs(p2_6_p0_6);
1552 }
1553 if( (max3 < fabs(p2_4_p0_4)) )
1554 {
1555 max3 = fabs(p2_4_p0_4);
1556 }
1557 double max4 = fabs(q0_0_p0_0);
1558 if( (max4 < fabs(q0_1_p0_1)) )
1559 {
1560 max4 = fabs(q0_1_p0_1);
1561 }
1562 if( (max4 < fabs(q0_2_p0_2)) )
1563 {
1564 max4 = fabs(q0_2_p0_2);
1565 }
1566 if( (max4 < fabs(q0_3_p0_3)) )
1567 {
1568 max4 = fabs(q0_3_p0_3);
1569 }
1570 if( (max4 < fabs(q0_4_p0_4)) )
1571 {
1572 max4 = fabs(q0_4_p0_4);
1573 }
1574 if( (max4 < fabs(q0_5_p0_5)) )
1575 {
1576 max4 = fabs(q0_5_p0_5);
1577 }
1578 if( (max4 < fabs(q0_6_p0_6)) )
1579 {
1580 max4 = fabs(q0_6_p0_6);
1581 }
1582 if( (max4 < fabs(q1_0_p0_0)) )
1583 {
1584 max4 = fabs(q1_0_p0_0);
1585 }
1586 if( (max4 < fabs(q1_1_p0_1)) )
1587 {
1588 max4 = fabs(q1_1_p0_1);
1589 }
1590 if( (max4 < fabs(q1_2_p0_2)) )
1591 {
1592 max4 = fabs(q1_2_p0_2);
1593 }
1594 if( (max4 < fabs(q1_3_p0_3)) )
1595 {
1596 max4 = fabs(q1_3_p0_3);
1597 }
1598 if( (max4 < fabs(q1_4_p0_4)) )
1599 {
1600 max4 = fabs(q1_4_p0_4);
1601 }
1602 if( (max4 < fabs(q1_5_p0_5)) )
1603 {
1604 max4 = fabs(q1_5_p0_5);
1605 }
1606 if( (max4 < fabs(q1_6_p0_6)) )
1607 {
1608 max4 = fabs(q1_6_p0_6);
1609 }
1610 double max5 = fabs(q1_0_p0_0);
1611 if( (max5 < fabs(q1_1_p0_1)) )
1612 {
1613 max5 = fabs(q1_1_p0_1);
1614 }
1615 if( (max5 < fabs(q1_2_p0_2)) )
1616 {
1617 max5 = fabs(q1_2_p0_2);
1618 }
1619 if( (max5 < fabs(q1_3_p0_3)) )
1620 {
1621 max5 = fabs(q1_3_p0_3);
1622 }
1623 if( (max5 < fabs(q1_4_p0_4)) )
1624 {
1625 max5 = fabs(q1_4_p0_4);
1626 }
1627 if( (max5 < fabs(q1_5_p0_5)) )
1628 {
1629 max5 = fabs(q1_5_p0_5);
1630 }
1631 if( (max5 < fabs(q1_6_p0_6)) )
1632 {
1633 max5 = fabs(q1_6_p0_6);
1634 }
1635 if( (max5 < fabs(q2_0_p0_0)) )
1636 {
1637 max5 = fabs(q2_0_p0_0);
1638 }
1639 if( (max5 < fabs(q2_1_p0_1)) )
1640 {
1641 max5 = fabs(q2_1_p0_1);
1642 }
1643 if( (max5 < fabs(q2_2_p0_2)) )
1644 {
1645 max5 = fabs(q2_2_p0_2);
1646 }
1647 if( (max5 < fabs(q2_3_p0_3)) )
1648 {
1649 max5 = fabs(q2_3_p0_3);
1650 }
1651 if( (max5 < fabs(q2_4_p0_4)) )
1652 {
1653 max5 = fabs(q2_4_p0_4);
1654 }
1655 if( (max5 < fabs(q2_5_p0_5)) )
1656 {
1657 max5 = fabs(q2_5_p0_5);
1658 }
1659 if( (max5 < fabs(q2_6_p0_6)) )
1660 {
1661 max5 = fabs(q2_6_p0_6);
1662 }
1663 double max6 = fabs(q2_0_p0_0);
1664 if( (max6 < fabs(q2_1_p0_1)) )
1665 {
1666 max6 = fabs(q2_1_p0_1);
1667 }
1668 if( (max6 < fabs(q2_2_p0_2)) )
1669 {
1670 max6 = fabs(q2_2_p0_2);
1671 }
1672 if( (max6 < fabs(q2_3_p0_3)) )
1673 {
1674 max6 = fabs(q2_3_p0_3);
1675 }
1676 if( (max6 < fabs(q2_4_p0_4)) )
1677 {
1678 max6 = fabs(q2_4_p0_4);
1679 }
1680 if( (max6 < fabs(q2_5_p0_5)) )
1681 {
1682 max6 = fabs(q2_5_p0_5);
1683 }
1684 if( (max6 < fabs(q2_6_p0_6)) )
1685 {
1686 max6 = fabs(q2_6_p0_6);
1687 }
1688 if( (max6 < fabs(q3_0_p0_0)) )
1689 {
1690 max6 = fabs(q3_0_p0_0);
1691 }
1692 if( (max6 < fabs(q3_1_p0_1)) )
1693 {
1694 max6 = fabs(q3_1_p0_1);
1695 }
1696 if( (max6 < fabs(q3_2_p0_2)) )
1697 {
1698 max6 = fabs(q3_2_p0_2);
1699 }
1700 if( (max6 < fabs(q3_3_p0_3)) )
1701 {
1702 max6 = fabs(q3_3_p0_3);
1703 }
1704 if( (max6 < fabs(q3_4_p0_4)) )
1705 {
1706 max6 = fabs(q3_4_p0_4);
1707 }
1708 if( (max6 < fabs(q3_5_p0_5)) )
1709 {
1710 max6 = fabs(q3_5_p0_5);
1711 }
1712 if( (max6 < fabs(q3_6_p0_6)) )
1713 {
1714 max6 = fabs(q3_6_p0_6);
1715 }
1716 double lower_bound_1;
1717 double upper_bound_1;
1718 int Delta_sign;
1719 int int_tmp_result;
1720 lower_bound_1 = max3;
1721 upper_bound_1 = max3;
1722 if( (max4 < lower_bound_1) )
1723 {
1724 lower_bound_1 = max4;
1725 }
1726 else
1727 {
1728 if( (max4 > upper_bound_1) )
1729 {
1730 upper_bound_1 = max4;
1731 }
1732 }
1733 if( (max1 < lower_bound_1) )
1734 {
1735 lower_bound_1 = max1;
1736 }
1737 else
1738 {
1739 if( (max1 > upper_bound_1) )
1740 {
1741 upper_bound_1 = max1;
1742 }
1743 }
1744 if( (max2 < lower_bound_1) )
1745 {
1746 lower_bound_1 = max2;
1747 }
1748 else
1749 {
1750 if( (max2 > upper_bound_1) )
1751 {
1752 upper_bound_1 = max2;
1753 }
1754 }
1755 if( (max5 < lower_bound_1) )
1756 {
1757 lower_bound_1 = max5;
1758 }
1759 else
1760 {
1761 if( (max5 > upper_bound_1) )
1762 {
1763 upper_bound_1 = max5;
1764 }
1765 }
1766 if( (max6 < lower_bound_1) )
1767 {
1768 lower_bound_1 = max6;
1769 }
1770 else
1771 {
1772 if( (max6 > upper_bound_1) )
1773 {
1774 upper_bound_1 = max6;
1775 }
1776 }
1777 if( (lower_bound_1 < 3.04548303565602498901e-50) )
1778 {
1779 return FPG_UNCERTAIN_VALUE;
1780 }
1781 else
1782 {
1783 if( (upper_bound_1 > 4.83570327845851562508e+24) )
1784 {
1785 return FPG_UNCERTAIN_VALUE;
1786 }
1787 eps = (2.78873548804336160566e-11 * (((((max1 * max4) * max3) * max5) * max2) * max6));
1788 if( (Delta > eps) )
1789 {
1790 int_tmp_result = 1;
1791 }
1792 else
1793 {
1794 if( (Delta < -eps) )
1795 {
1796 int_tmp_result = -1;
1797 }
1798 else
1799 {
1800 return FPG_UNCERTAIN_VALUE;
1801 }
1802 }
1803 }
1804 Delta_sign = int_tmp_result;
1805 int int_tmp_result_FFWKCAA;
1806 double max7 = max3;
1807 if( (max7 < max1) )
1808 {
1809 max7 = max1;
1810 }
1811 if( (max7 < max2) )
1812 {
1813 max7 = max2;
1814 }
1815 if( (max7 < max6) )
1816 {
1817 max7 = max6;
1818 }
1819 double max8 = max1;
1820 if( (max8 < fabs(p4_4_p0_4)) )
1821 {
1822 max8 = fabs(p4_4_p0_4);
1823 }
1824 if( (max8 < fabs(p4_1_p0_1)) )
1825 {
1826 max8 = fabs(p4_1_p0_1);
1827 }
1828 if( (max8 < fabs(p4_5_p0_5)) )
1829 {
1830 max8 = fabs(p4_5_p0_5);
1831 }
1832 if( (max8 < fabs(p4_0_p0_0)) )
1833 {
1834 max8 = fabs(p4_0_p0_0);
1835 }
1836 if( (max8 < fabs(p4_2_p0_2)) )
1837 {
1838 max8 = fabs(p4_2_p0_2);
1839 }
1840 if( (max8 < fabs(p4_3_p0_3)) )
1841 {
1842 max8 = fabs(p4_3_p0_3);
1843 }
1844 if( (max8 < fabs(p4_6_p0_6)) )
1845 {
1846 max8 = fabs(p4_6_p0_6);
1847 }
1848 if( (max7 < max8) )
1849 {
1850 max7 = max8;
1851 }
1852 double max9 = max1;
1853 if( (max9 < max5) )
1854 {
1855 max9 = max5;
1856 }
1857 if( (max9 < max8) )
1858 {
1859 max9 = max8;
1860 }
1861 double max10 = max4;
1862 if( (max10 < max1) )
1863 {
1864 max10 = max1;
1865 }
1866 if( (max10 < max5) )
1867 {
1868 max10 = max5;
1869 }
1870 if( (max10 < max6) )
1871 {
1872 max10 = max6;
1873 }
1874 double max11 = max4;
1875 if( (max11 < max5) )
1876 {
1877 max11 = max5;
1878 }
1879 if( (max10 < max11) )
1880 {
1881 max10 = max11;
1882 }
1883 lower_bound_1 = max3;
1884 upper_bound_1 = max3;
1885 if( (max2 < lower_bound_1) )
1886 {
1887 lower_bound_1 = max2;
1888 }
1889 if( (max7 < lower_bound_1) )
1890 {
1891 lower_bound_1 = max7;
1892 }
1893 else
1894 {
1895 if( (max7 > upper_bound_1) )
1896 {
1897 upper_bound_1 = max7;
1898 }
1899 }
1900 if( (max8 < lower_bound_1) )
1901 {
1902 lower_bound_1 = max8;
1903 }
1904 if( (max9 < lower_bound_1) )
1905 {
1906 lower_bound_1 = max9;
1907 }
1908 else
1909 {
1910 if( (max9 > upper_bound_1) )
1911 {
1912 upper_bound_1 = max9;
1913 }
1914 }
1915 if( (max10 < lower_bound_1) )
1916 {
1917 lower_bound_1 = max10;
1918 }
1919 else
1920 {
1921 if( (max10 > upper_bound_1) )
1922 {
1923 upper_bound_1 = max10;
1924 }
1925 }
1926 if( (max11 < lower_bound_1) )
1927 {
1928 lower_bound_1 = max11;
1929 }
1930 else
1931 {
1932 if( (max11 > upper_bound_1) )
1933 {
1934 upper_bound_1 = max11;
1935 }
1936 }
1937 if( (lower_bound_1 < 4.48906690519700369396e-38) )
1938 {
1939 return FPG_UNCERTAIN_VALUE;
1940 }
1941 else
1942 {
1943 if( (upper_bound_1 > 4.83570327845851562508e+24) )
1944 {
1945 return FPG_UNCERTAIN_VALUE;
1946 }
1947 eps = (1.34926049830188433875e-09 * (((((((max8 * max11) * max3) * max10) * max2) * max10) * max9) * max7));
1948 if( (r > eps) )
1949 {
1950 int_tmp_result_FFWKCAA = 1;
1951 }
1952 else
1953 {
1954 if( (r < -eps) )
1955 {
1956 int_tmp_result_FFWKCAA = -1;
1957 }
1958 else
1959 {
1960 return FPG_UNCERTAIN_VALUE;
1961 }
1962 }
1963 }
1964 return (Delta_sign * int_tmp_result_FFWKCAA);
1965 }
1966
1967
1968 19207 inline int side4_8d_filter( const double* p0, const double* p1, const double* p2, const double* p3, const double* p4, const double* q0, const double* q1, const double* q2, const double* q3) {
1969 19207 double p1_0_p0_0 = (p1[0] - p0[0]);
1970 19207 double p1_1_p0_1 = (p1[1] - p0[1]);
1971 19207 double p1_2_p0_2 = (p1[2] - p0[2]);
1972 19207 double p1_3_p0_3 = (p1[3] - p0[3]);
1973 19207 double p1_4_p0_4 = (p1[4] - p0[4]);
1974 19207 double p1_5_p0_5 = (p1[5] - p0[5]);
1975 19207 double p1_6_p0_6 = (p1[6] - p0[6]);
1976 19207 double p1_7_p0_7 = (p1[7] - p0[7]);
1977 double l1;
1978 19207 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)));
1979 19207 double p2_0_p0_0 = (p2[0] - p0[0]);
1980 19207 double p2_1_p0_1 = (p2[1] - p0[1]);
1981 19207 double p2_2_p0_2 = (p2[2] - p0[2]);
1982 19207 double p2_3_p0_3 = (p2[3] - p0[3]);
1983 19207 double p2_4_p0_4 = (p2[4] - p0[4]);
1984 19207 double p2_5_p0_5 = (p2[5] - p0[5]);
1985 19207 double p2_6_p0_6 = (p2[6] - p0[6]);
1986 19207 double p2_7_p0_7 = (p2[7] - p0[7]);
1987 double l2;
1988 19207 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)));
1989 19207 double p3_0_p0_0 = (p3[0] - p0[0]);
1990 19207 double p3_1_p0_1 = (p3[1] - p0[1]);
1991 19207 double p3_2_p0_2 = (p3[2] - p0[2]);
1992 19207 double p3_3_p0_3 = (p3[3] - p0[3]);
1993 19207 double p3_4_p0_4 = (p3[4] - p0[4]);
1994 19207 double p3_5_p0_5 = (p3[5] - p0[5]);
1995 19207 double p3_6_p0_6 = (p3[6] - p0[6]);
1996 19207 double p3_7_p0_7 = (p3[7] - p0[7]);
1997 double l3;
1998 19207 l3 = (1 * ((((((((p3_0_p0_0 * p3_0_p0_0) + (p3_1_p0_1 * p3_1_p0_1)) + (p3_2_p0_2 * p3_2_p0_2)) + (p3_3_p0_3 * p3_3_p0_3)) + (p3_4_p0_4 * p3_4_p0_4)) + (p3_5_p0_5 * p3_5_p0_5)) + (p3_6_p0_6 * p3_6_p0_6)) + (p3_7_p0_7 * p3_7_p0_7)));
1999 19207 double p4_0_p0_0 = (p4[0] - p0[0]);
2000 19207 double p4_1_p0_1 = (p4[1] - p0[1]);
2001 19207 double p4_2_p0_2 = (p4[2] - p0[2]);
2002 19207 double p4_3_p0_3 = (p4[3] - p0[3]);
2003 19207 double p4_4_p0_4 = (p4[4] - p0[4]);
2004 19207 double p4_5_p0_5 = (p4[5] - p0[5]);
2005 19207 double p4_6_p0_6 = (p4[6] - p0[6]);
2006 19207 double p4_7_p0_7 = (p4[7] - p0[7]);
2007 double l4;
2008 19207 l4 = (1 * ((((((((p4_0_p0_0 * p4_0_p0_0) + (p4_1_p0_1 * p4_1_p0_1)) + (p4_2_p0_2 * p4_2_p0_2)) + (p4_3_p0_3 * p4_3_p0_3)) + (p4_4_p0_4 * p4_4_p0_4)) + (p4_5_p0_5 * p4_5_p0_5)) + (p4_6_p0_6 * p4_6_p0_6)) + (p4_7_p0_7 * p4_7_p0_7)));
2009 19207 double q0_0_p0_0 = (q0[0] - p0[0]);
2010 19207 double q0_1_p0_1 = (q0[1] - p0[1]);
2011 19207 double q0_2_p0_2 = (q0[2] - p0[2]);
2012 19207 double q0_3_p0_3 = (q0[3] - p0[3]);
2013 19207 double q0_4_p0_4 = (q0[4] - p0[4]);
2014 19207 double q0_5_p0_5 = (q0[5] - p0[5]);
2015 19207 double q0_6_p0_6 = (q0[6] - p0[6]);
2016 19207 double q0_7_p0_7 = (q0[7] - p0[7]);
2017 double a10;
2018 19207 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)));
2019 19207 double q1_0_p0_0 = (q1[0] - p0[0]);
2020 19207 double q1_1_p0_1 = (q1[1] - p0[1]);
2021 19207 double q1_2_p0_2 = (q1[2] - p0[2]);
2022 19207 double q1_3_p0_3 = (q1[3] - p0[3]);
2023 19207 double q1_4_p0_4 = (q1[4] - p0[4]);
2024 19207 double q1_5_p0_5 = (q1[5] - p0[5]);
2025 19207 double q1_6_p0_6 = (q1[6] - p0[6]);
2026 19207 double q1_7_p0_7 = (q1[7] - p0[7]);
2027 double a11;
2028 19207 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)));
2029 19207 double q2_0_p0_0 = (q2[0] - p0[0]);
2030 19207 double q2_1_p0_1 = (q2[1] - p0[1]);
2031 19207 double q2_2_p0_2 = (q2[2] - p0[2]);
2032 19207 double q2_3_p0_3 = (q2[3] - p0[3]);
2033 19207 double q2_4_p0_4 = (q2[4] - p0[4]);
2034 19207 double q2_5_p0_5 = (q2[5] - p0[5]);
2035 19207 double q2_6_p0_6 = (q2[6] - p0[6]);
2036 19207 double q2_7_p0_7 = (q2[7] - p0[7]);
2037 double a12;
2038 19207 a12 = (2 * ((((((((p1_0_p0_0 * q2_0_p0_0) + (p1_1_p0_1 * q2_1_p0_1)) + (p1_2_p0_2 * q2_2_p0_2)) + (p1_3_p0_3 * q2_3_p0_3)) + (p1_4_p0_4 * q2_4_p0_4)) + (p1_5_p0_5 * q2_5_p0_5)) + (p1_6_p0_6 * q2_6_p0_6)) + (p1_7_p0_7 * q2_7_p0_7)));
2039 19207 double q3_0_p0_0 = (q3[0] - p0[0]);
2040 19207 double q3_1_p0_1 = (q3[1] - p0[1]);
2041 19207 double q3_2_p0_2 = (q3[2] - p0[2]);
2042 19207 double q3_3_p0_3 = (q3[3] - p0[3]);
2043 19207 double q3_4_p0_4 = (q3[4] - p0[4]);
2044 19207 double q3_5_p0_5 = (q3[5] - p0[5]);
2045 19207 double q3_6_p0_6 = (q3[6] - p0[6]);
2046 19207 double q3_7_p0_7 = (q3[7] - p0[7]);
2047 double a13;
2048 19207 a13 = (2 * ((((((((p1_0_p0_0 * q3_0_p0_0) + (p1_1_p0_1 * q3_1_p0_1)) + (p1_2_p0_2 * q3_2_p0_2)) + (p1_3_p0_3 * q3_3_p0_3)) + (p1_4_p0_4 * q3_4_p0_4)) + (p1_5_p0_5 * q3_5_p0_5)) + (p1_6_p0_6 * q3_6_p0_6)) + (p1_7_p0_7 * q3_7_p0_7)));
2049 double a20;
2050 19207 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)));
2051 double a21;
2052 19207 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)));
2053 double a22;
2054 19207 a22 = (2 * ((((((((p2_0_p0_0 * q2_0_p0_0) + (p2_1_p0_1 * q2_1_p0_1)) + (p2_2_p0_2 * q2_2_p0_2)) + (p2_3_p0_3 * q2_3_p0_3)) + (p2_4_p0_4 * q2_4_p0_4)) + (p2_5_p0_5 * q2_5_p0_5)) + (p2_6_p0_6 * q2_6_p0_6)) + (p2_7_p0_7 * q2_7_p0_7)));
2055 double a23;
2056 19207 a23 = (2 * ((((((((p2_0_p0_0 * q3_0_p0_0) + (p2_1_p0_1 * q3_1_p0_1)) + (p2_2_p0_2 * q3_2_p0_2)) + (p2_3_p0_3 * q3_3_p0_3)) + (p2_4_p0_4 * q3_4_p0_4)) + (p2_5_p0_5 * q3_5_p0_5)) + (p2_6_p0_6 * q3_6_p0_6)) + (p2_7_p0_7 * q3_7_p0_7)));
2057 double a30;
2058 19207 a30 = (2 * ((((((((p3_0_p0_0 * q0_0_p0_0) + (p3_1_p0_1 * q0_1_p0_1)) + (p3_2_p0_2 * q0_2_p0_2)) + (p3_3_p0_3 * q0_3_p0_3)) + (p3_4_p0_4 * q0_4_p0_4)) + (p3_5_p0_5 * q0_5_p0_5)) + (p3_6_p0_6 * q0_6_p0_6)) + (p3_7_p0_7 * q0_7_p0_7)));
2059 double a31;
2060 19207 a31 = (2 * ((((((((p3_0_p0_0 * q1_0_p0_0) + (p3_1_p0_1 * q1_1_p0_1)) + (p3_2_p0_2 * q1_2_p0_2)) + (p3_3_p0_3 * q1_3_p0_3)) + (p3_4_p0_4 * q1_4_p0_4)) + (p3_5_p0_5 * q1_5_p0_5)) + (p3_6_p0_6 * q1_6_p0_6)) + (p3_7_p0_7 * q1_7_p0_7)));
2061 double a32;
2062 19207 a32 = (2 * ((((((((p3_0_p0_0 * q2_0_p0_0) + (p3_1_p0_1 * q2_1_p0_1)) + (p3_2_p0_2 * q2_2_p0_2)) + (p3_3_p0_3 * q2_3_p0_3)) + (p3_4_p0_4 * q2_4_p0_4)) + (p3_5_p0_5 * q2_5_p0_5)) + (p3_6_p0_6 * q2_6_p0_6)) + (p3_7_p0_7 * q2_7_p0_7)));
2063 double a33;
2064 19207 a33 = (2 * ((((((((p3_0_p0_0 * q3_0_p0_0) + (p3_1_p0_1 * q3_1_p0_1)) + (p3_2_p0_2 * q3_2_p0_2)) + (p3_3_p0_3 * q3_3_p0_3)) + (p3_4_p0_4 * q3_4_p0_4)) + (p3_5_p0_5 * q3_5_p0_5)) + (p3_6_p0_6 * q3_6_p0_6)) + (p3_7_p0_7 * q3_7_p0_7)));
2065 double a40;
2066 19207 a40 = (2 * ((((((((p4_0_p0_0 * q0_0_p0_0) + (p4_1_p0_1 * q0_1_p0_1)) + (p4_2_p0_2 * q0_2_p0_2)) + (p4_3_p0_3 * q0_3_p0_3)) + (p4_4_p0_4 * q0_4_p0_4)) + (p4_5_p0_5 * q0_5_p0_5)) + (p4_6_p0_6 * q0_6_p0_6)) + (p4_7_p0_7 * q0_7_p0_7)));
2067 double a41;
2068 19207 a41 = (2 * ((((((((p4_0_p0_0 * q1_0_p0_0) + (p4_1_p0_1 * q1_1_p0_1)) + (p4_2_p0_2 * q1_2_p0_2)) + (p4_3_p0_3 * q1_3_p0_3)) + (p4_4_p0_4 * q1_4_p0_4)) + (p4_5_p0_5 * q1_5_p0_5)) + (p4_6_p0_6 * q1_6_p0_6)) + (p4_7_p0_7 * q1_7_p0_7)));
2069 double a42;
2070 19207 a42 = (2 * ((((((((p4_0_p0_0 * q2_0_p0_0) + (p4_1_p0_1 * q2_1_p0_1)) + (p4_2_p0_2 * q2_2_p0_2)) + (p4_3_p0_3 * q2_3_p0_3)) + (p4_4_p0_4 * q2_4_p0_4)) + (p4_5_p0_5 * q2_5_p0_5)) + (p4_6_p0_6 * q2_6_p0_6)) + (p4_7_p0_7 * q2_7_p0_7)));
2071 double a43;
2072 19207 a43 = (2 * ((((((((p4_0_p0_0 * q3_0_p0_0) + (p4_1_p0_1 * q3_1_p0_1)) + (p4_2_p0_2 * q3_2_p0_2)) + (p4_3_p0_3 * q3_3_p0_3)) + (p4_4_p0_4 * q3_4_p0_4)) + (p4_5_p0_5 * q3_5_p0_5)) + (p4_6_p0_6 * q3_6_p0_6)) + (p4_7_p0_7 * q3_7_p0_7)));
2073 double b00;
2074 19207 b00 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
2075 double b01;
2076 19207 b01 = -((((a22 * a33) - (a23 * a32)) + ((a23 * a31) - (a21 * a33))) + ((a21 * a32) - (a22 * a31)));
2077 double b02;
2078 19207 b02 = ((((a12 * a33) - (a13 * a32)) + ((a13 * a31) - (a11 * a33))) + ((a11 * a32) - (a12 * a31)));
2079 double b03;
2080 19207 b03 = -((((a12 * a23) - (a13 * a22)) + ((a13 * a21) - (a11 * a23))) + ((a11 * a22) - (a12 * a21)));
2081 double b10;
2082 19207 b10 = -(((a10 * ((a22 * a33) - (a23 * a32))) - (a20 * ((a12 * a33) - (a13 * a32)))) + (a30 * ((a12 * a23) - (a13 * a22))));
2083 double b11;
2084 19207 b11 = ((((a22 * a33) - (a23 * a32)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a32) - (a22 * a30)));
2085 double b12;
2086 19207 b12 = -((((a12 * a33) - (a13 * a32)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a32) - (a12 * a30)));
2087 double b13;
2088 19207 b13 = ((((a12 * a23) - (a13 * a22)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a22) - (a12 * a20)));
2089 double b20;
2090 19207 b20 = (((a10 * ((a21 * a33) - (a23 * a31))) - (a20 * ((a11 * a33) - (a13 * a31)))) + (a30 * ((a11 * a23) - (a13 * a21))));
2091 double b21;
2092 19207 b21 = -((((a21 * a33) - (a23 * a31)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a31) - (a21 * a30)));
2093 double b22;
2094 19207 b22 = ((((a11 * a33) - (a13 * a31)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a31) - (a11 * a30)));
2095 double b23;
2096 19207 b23 = -((((a11 * a23) - (a13 * a21)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a21) - (a11 * a20)));
2097 double b30;
2098 19207 b30 = -(((a10 * ((a21 * a32) - (a22 * a31))) - (a20 * ((a11 * a32) - (a12 * a31)))) + (a30 * ((a11 * a22) - (a12 * a21))));
2099 double b31;
2100 19207 b31 = ((((a21 * a32) - (a22 * a31)) + ((a22 * a30) - (a20 * a32))) + ((a20 * a31) - (a21 * a30)));
2101 double b32;
2102 19207 b32 = -((((a11 * a32) - (a12 * a31)) + ((a12 * a30) - (a10 * a32))) + ((a10 * a31) - (a11 * a30)));
2103 double b33;
2104 19207 b33 = ((((a11 * a22) - (a12 * a21)) + ((a12 * a20) - (a10 * a22))) + ((a10 * a21) - (a11 * a20)));
2105 double Delta;
2106 19207 Delta = (((b00 + b10) + b20) + b30);
2107 double DeltaLambda0;
2108 19207 DeltaLambda0 = ((((b01 * l1) + (b02 * l2)) + (b03 * l3)) + b00);
2109 double DeltaLambda1;
2110 19207 DeltaLambda1 = ((((b11 * l1) + (b12 * l2)) + (b13 * l3)) + b10);
2111 double DeltaLambda2;
2112 19207 DeltaLambda2 = ((((b21 * l1) + (b22 * l2)) + (b23 * l3)) + b20);
2113 double DeltaLambda3;
2114 19207 DeltaLambda3 = ((((b31 * l1) + (b32 * l2)) + (b33 * l3)) + b30);
2115 double r;
2116 19207 r = ((Delta * l4) - ((((a40 * DeltaLambda0) + (a41 * DeltaLambda1)) + (a42 * DeltaLambda2)) + (a43 * DeltaLambda3)));
2117 double eps;
2118 19207 double max1 = fabs(p2_5_p0_5);
2119
2/2
✓ Branch 0 taken 13819 times.
✓ Branch 1 taken 5388 times.
19207 if( (max1 < fabs(p2_3_p0_3)) )
2120 {
2121 13819 max1 = fabs(p2_3_p0_3);
2122 }
2123
2/2
✓ Branch 0 taken 581 times.
✓ Branch 1 taken 18626 times.
19207 if( (max1 < fabs(p2_0_p0_0)) )
2124 {
2125 581 max1 = fabs(p2_0_p0_0);
2126 }
2127
2/2
✓ Branch 0 taken 12636 times.
✓ Branch 1 taken 6571 times.
19207 if( (max1 < fabs(p2_1_p0_1)) )
2128 {
2129 12636 max1 = fabs(p2_1_p0_1);
2130 }
2131
2/2
✓ Branch 0 taken 3790 times.
✓ Branch 1 taken 15417 times.
19207 if( (max1 < fabs(p2_6_p0_6)) )
2132 {
2133 3790 max1 = fabs(p2_6_p0_6);
2134 }
2135
2/2
✓ Branch 0 taken 7948 times.
✓ Branch 1 taken 11259 times.
19207 if( (max1 < fabs(p2_2_p0_2)) )
2136 {
2137 7948 max1 = fabs(p2_2_p0_2);
2138 }
2139
2/2
✓ Branch 0 taken 2116 times.
✓ Branch 1 taken 17091 times.
19207 if( (max1 < fabs(p2_4_p0_4)) )
2140 {
2141 2116 max1 = fabs(p2_4_p0_4);
2142 }
2143
2/2
✓ Branch 0 taken 1651 times.
✓ Branch 1 taken 17556 times.
19207 if( (max1 < fabs(p2_7_p0_7)) )
2144 {
2145 1651 max1 = fabs(p2_7_p0_7);
2146 }
2147 19207 double max2 = fabs(p1_4_p0_4);
2148
2/2
✓ Branch 0 taken 7645 times.
✓ Branch 1 taken 11562 times.
19207 if( (max2 < fabs(p1_3_p0_3)) )
2149 {
2150 7645 max2 = fabs(p1_3_p0_3);
2151 }
2152
2/2
✓ Branch 0 taken 7487 times.
✓ Branch 1 taken 11720 times.
19207 if( (max2 < fabs(p1_7_p0_7)) )
2153 {
2154 7487 max2 = fabs(p1_7_p0_7);
2155 }
2156
2/2
✓ Branch 0 taken 48 times.
✓ Branch 1 taken 19159 times.
19207 if( (max2 < fabs(p1_0_p0_0)) )
2157 {
2158 48 max2 = fabs(p1_0_p0_0);
2159 }
2160
2/2
✓ Branch 0 taken 11223 times.
✓ Branch 1 taken 7984 times.
19207 if( (max2 < fabs(p1_2_p0_2)) )
2161 {
2162 11223 max2 = fabs(p1_2_p0_2);
2163 }
2164
2/2
✓ Branch 0 taken 848 times.
✓ Branch 1 taken 18359 times.
19207 if( (max2 < fabs(p1_5_p0_5)) )
2165 {
2166 848 max2 = fabs(p1_5_p0_5);
2167 }
2168
2/2
✓ Branch 0 taken 7355 times.
✓ Branch 1 taken 11852 times.
19207 if( (max2 < fabs(p1_1_p0_1)) )
2169 {
2170 7355 max2 = fabs(p1_1_p0_1);
2171 }
2172
2/2
✓ Branch 0 taken 941 times.
✓ Branch 1 taken 18266 times.
19207 if( (max2 < fabs(p1_6_p0_6)) )
2173 {
2174 941 max2 = fabs(p1_6_p0_6);
2175 }
2176 19207 double max3 = fabs(p3_3_p0_3);
2177
2/2
✓ Branch 0 taken 1271 times.
✓ Branch 1 taken 17936 times.
19207 if( (max3 < fabs(p3_0_p0_0)) )
2178 {
2179 1271 max3 = fabs(p3_0_p0_0);
2180 }
2181
2/2
✓ Branch 0 taken 13700 times.
✓ Branch 1 taken 5507 times.
19207 if( (max3 < fabs(p3_1_p0_1)) )
2182 {
2183 13700 max3 = fabs(p3_1_p0_1);
2184 }
2185
2/2
✓ Branch 0 taken 8418 times.
✓ Branch 1 taken 10789 times.
19207 if( (max3 < fabs(p3_2_p0_2)) )
2186 {
2187 8418 max3 = fabs(p3_2_p0_2);
2188 }
2189
2/2
✓ Branch 0 taken 2624 times.
✓ Branch 1 taken 16583 times.
19207 if( (max3 < fabs(p3_4_p0_4)) )
2190 {
2191 2624 max3 = fabs(p3_4_p0_4);
2192 }
2193
2/2
✓ Branch 0 taken 636 times.
✓ Branch 1 taken 18571 times.
19207 if( (max3 < fabs(p3_5_p0_5)) )
2194 {
2195 636 max3 = fabs(p3_5_p0_5);
2196 }
2197
2/2
✓ Branch 0 taken 1090 times.
✓ Branch 1 taken 18117 times.
19207 if( (max3 < fabs(p3_6_p0_6)) )
2198 {
2199 1090 max3 = fabs(p3_6_p0_6);
2200 }
2201
2/2
✓ Branch 0 taken 1357 times.
✓ Branch 1 taken 17850 times.
19207 if( (max3 < fabs(p3_7_p0_7)) )
2202 {
2203 1357 max3 = fabs(p3_7_p0_7);
2204 }
2205 19207 double max4 = fabs(q0_0_p0_0);
2206
2/2
✓ Branch 0 taken 11890 times.
✓ Branch 1 taken 7317 times.
19207 if( (max4 < fabs(q0_1_p0_1)) )
2207 {
2208 11890 max4 = fabs(q0_1_p0_1);
2209 }
2210
2/2
✓ Branch 0 taken 8008 times.
✓ Branch 1 taken 11199 times.
19207 if( (max4 < fabs(q0_2_p0_2)) )
2211 {
2212 8008 max4 = fabs(q0_2_p0_2);
2213 }
2214
2/2
✓ Branch 0 taken 12373 times.
✓ Branch 1 taken 6834 times.
19207 if( (max4 < fabs(q0_3_p0_3)) )
2215 {
2216 12373 max4 = fabs(q0_3_p0_3);
2217 }
2218
2/2
✓ Branch 0 taken 13435 times.
✓ Branch 1 taken 5772 times.
19207 if( (max4 < fabs(q0_4_p0_4)) )
2219 {
2220 13435 max4 = fabs(q0_4_p0_4);
2221 }
2222
2/2
✓ Branch 0 taken 2820 times.
✓ Branch 1 taken 16387 times.
19207 if( (max4 < fabs(q0_5_p0_5)) )
2223 {
2224 2820 max4 = fabs(q0_5_p0_5);
2225 }
2226
2/2
✓ Branch 0 taken 2969 times.
✓ Branch 1 taken 16238 times.
19207 if( (max4 < fabs(q0_6_p0_6)) )
2227 {
2228 2969 max4 = fabs(q0_6_p0_6);
2229 }
2230
2/2
✓ Branch 0 taken 5302 times.
✓ Branch 1 taken 13905 times.
19207 if( (max4 < fabs(q0_7_p0_7)) )
2231 {
2232 5302 max4 = fabs(q0_7_p0_7);
2233 }
2234
2/2
✓ Branch 0 taken 1310 times.
✓ Branch 1 taken 17897 times.
19207 if( (max4 < fabs(q1_0_p0_0)) )
2235 {
2236 1310 max4 = fabs(q1_0_p0_0);
2237 }
2238
2/2
✓ Branch 0 taken 1196 times.
✓ Branch 1 taken 18011 times.
19207 if( (max4 < fabs(q1_1_p0_1)) )
2239 {
2240 1196 max4 = fabs(q1_1_p0_1);
2241 }
2242
2/2
✓ Branch 0 taken 745 times.
✓ Branch 1 taken 18462 times.
19207 if( (max4 < fabs(q1_2_p0_2)) )
2243 {
2244 745 max4 = fabs(q1_2_p0_2);
2245 }
2246
2/2
✓ Branch 0 taken 4009 times.
✓ Branch 1 taken 15198 times.
19207 if( (max4 < fabs(q1_3_p0_3)) )
2247 {
2248 4009 max4 = fabs(q1_3_p0_3);
2249 }
2250
2/2
✓ Branch 0 taken 3832 times.
✓ Branch 1 taken 15375 times.
19207 if( (max4 < fabs(q1_4_p0_4)) )
2251 {
2252 3832 max4 = fabs(q1_4_p0_4);
2253 }
2254
2/2
✓ Branch 0 taken 2207 times.
✓ Branch 1 taken 17000 times.
19207 if( (max4 < fabs(q1_5_p0_5)) )
2255 {
2256 2207 max4 = fabs(q1_5_p0_5);
2257 }
2258
2/2
✓ Branch 0 taken 1563 times.
✓ Branch 1 taken 17644 times.
19207 if( (max4 < fabs(q1_6_p0_6)) )
2259 {
2260 1563 max4 = fabs(q1_6_p0_6);
2261 }
2262
2/2
✓ Branch 0 taken 2583 times.
✓ Branch 1 taken 16624 times.
19207 if( (max4 < fabs(q1_7_p0_7)) )
2263 {
2264 2583 max4 = fabs(q1_7_p0_7);
2265 }
2266 19207 double max5 = fabs(q1_0_p0_0);
2267
2/2
✓ Branch 0 taken 13138 times.
✓ Branch 1 taken 6069 times.
19207 if( (max5 < fabs(q1_1_p0_1)) )
2268 {
2269 13138 max5 = fabs(q1_1_p0_1);
2270 }
2271
2/2
✓ Branch 0 taken 7127 times.
✓ Branch 1 taken 12080 times.
19207 if( (max5 < fabs(q1_2_p0_2)) )
2272 {
2273 7127 max5 = fabs(q1_2_p0_2);
2274 }
2275
2/2
✓ Branch 0 taken 11706 times.
✓ Branch 1 taken 7501 times.
19207 if( (max5 < fabs(q1_3_p0_3)) )
2276 {
2277 11706 max5 = fabs(q1_3_p0_3);
2278 }
2279
2/2
✓ Branch 0 taken 11994 times.
✓ Branch 1 taken 7213 times.
19207 if( (max5 < fabs(q1_4_p0_4)) )
2280 {
2281 11994 max5 = fabs(q1_4_p0_4);
2282 }
2283
2/2
✓ Branch 0 taken 3000 times.
✓ Branch 1 taken 16207 times.
19207 if( (max5 < fabs(q1_5_p0_5)) )
2284 {
2285 3000 max5 = fabs(q1_5_p0_5);
2286 }
2287
2/2
✓ Branch 0 taken 3005 times.
✓ Branch 1 taken 16202 times.
19207 if( (max5 < fabs(q1_6_p0_6)) )
2288 {
2289 3005 max5 = fabs(q1_6_p0_6);
2290 }
2291
2/2
✓ Branch 0 taken 4676 times.
✓ Branch 1 taken 14531 times.
19207 if( (max5 < fabs(q1_7_p0_7)) )
2292 {
2293 4676 max5 = fabs(q1_7_p0_7);
2294 }
2295
2/2
✓ Branch 0 taken 1170 times.
✓ Branch 1 taken 18037 times.
19207 if( (max5 < fabs(q2_0_p0_0)) )
2296 {
2297 1170 max5 = fabs(q2_0_p0_0);
2298 }
2299
2/2
✓ Branch 0 taken 971 times.
✓ Branch 1 taken 18236 times.
19207 if( (max5 < fabs(q2_1_p0_1)) )
2300 {
2301 971 max5 = fabs(q2_1_p0_1);
2302 }
2303
2/2
✓ Branch 0 taken 944 times.
✓ Branch 1 taken 18263 times.
19207 if( (max5 < fabs(q2_2_p0_2)) )
2304 {
2305 944 max5 = fabs(q2_2_p0_2);
2306 }
2307
2/2
✓ Branch 0 taken 4439 times.
✓ Branch 1 taken 14768 times.
19207 if( (max5 < fabs(q2_3_p0_3)) )
2308 {
2309 4439 max5 = fabs(q2_3_p0_3);
2310 }
2311
2/2
✓ Branch 0 taken 5281 times.
✓ Branch 1 taken 13926 times.
19207 if( (max5 < fabs(q2_4_p0_4)) )
2312 {
2313 5281 max5 = fabs(q2_4_p0_4);
2314 }
2315
2/2
✓ Branch 0 taken 1362 times.
✓ Branch 1 taken 17845 times.
19207 if( (max5 < fabs(q2_5_p0_5)) )
2316 {
2317 1362 max5 = fabs(q2_5_p0_5);
2318 }
2319
2/2
✓ Branch 0 taken 1388 times.
✓ Branch 1 taken 17819 times.
19207 if( (max5 < fabs(q2_6_p0_6)) )
2320 {
2321 1388 max5 = fabs(q2_6_p0_6);
2322 }
2323
2/2
✓ Branch 0 taken 2538 times.
✓ Branch 1 taken 16669 times.
19207 if( (max5 < fabs(q2_7_p0_7)) )
2324 {
2325 2538 max5 = fabs(q2_7_p0_7);
2326 }
2327 19207 double max6 = fabs(q2_0_p0_0);
2328
2/2
✓ Branch 0 taken 12018 times.
✓ Branch 1 taken 7189 times.
19207 if( (max6 < fabs(q2_1_p0_1)) )
2329 {
2330 12018 max6 = fabs(q2_1_p0_1);
2331 }
2332
2/2
✓ Branch 0 taken 7672 times.
✓ Branch 1 taken 11535 times.
19207 if( (max6 < fabs(q2_2_p0_2)) )
2333 {
2334 7672 max6 = fabs(q2_2_p0_2);
2335 }
2336
2/2
✓ Branch 0 taken 12623 times.
✓ Branch 1 taken 6584 times.
19207 if( (max6 < fabs(q2_3_p0_3)) )
2337 {
2338 12623 max6 = fabs(q2_3_p0_3);
2339 }
2340
2/2
✓ Branch 0 taken 13664 times.
✓ Branch 1 taken 5543 times.
19207 if( (max6 < fabs(q2_4_p0_4)) )
2341 {
2342 13664 max6 = fabs(q2_4_p0_4);
2343 }
2344
2/2
✓ Branch 0 taken 2735 times.
✓ Branch 1 taken 16472 times.
19207 if( (max6 < fabs(q2_5_p0_5)) )
2345 {
2346 2735 max6 = fabs(q2_5_p0_5);
2347 }
2348
2/2
✓ Branch 0 taken 3020 times.
✓ Branch 1 taken 16187 times.
19207 if( (max6 < fabs(q2_6_p0_6)) )
2349 {
2350 3020 max6 = fabs(q2_6_p0_6);
2351 }
2352
2/2
✓ Branch 0 taken 5457 times.
✓ Branch 1 taken 13750 times.
19207 if( (max6 < fabs(q2_7_p0_7)) )
2353 {
2354 5457 max6 = fabs(q2_7_p0_7);
2355 }
2356
2/2
✓ Branch 0 taken 1819 times.
✓ Branch 1 taken 17388 times.
19207 if( (max6 < fabs(q3_0_p0_0)) )
2357 {
2358 1819 max6 = fabs(q3_0_p0_0);
2359 }
2360
2/2
✓ Branch 0 taken 1184 times.
✓ Branch 1 taken 18023 times.
19207 if( (max6 < fabs(q3_1_p0_1)) )
2361 {
2362 1184 max6 = fabs(q3_1_p0_1);
2363 }
2364
2/2
✓ Branch 0 taken 662 times.
✓ Branch 1 taken 18545 times.
19207 if( (max6 < fabs(q3_2_p0_2)) )
2365 {
2366 662 max6 = fabs(q3_2_p0_2);
2367 }
2368
2/2
✓ Branch 0 taken 5426 times.
✓ Branch 1 taken 13781 times.
19207 if( (max6 < fabs(q3_3_p0_3)) )
2369 {
2370 5426 max6 = fabs(q3_3_p0_3);
2371 }
2372
2/2
✓ Branch 0 taken 5370 times.
✓ Branch 1 taken 13837 times.
19207 if( (max6 < fabs(q3_4_p0_4)) )
2373 {
2374 5370 max6 = fabs(q3_4_p0_4);
2375 }
2376
2/2
✓ Branch 0 taken 1964 times.
✓ Branch 1 taken 17243 times.
19207 if( (max6 < fabs(q3_5_p0_5)) )
2377 {
2378 1964 max6 = fabs(q3_5_p0_5);
2379 }
2380
2/2
✓ Branch 0 taken 1650 times.
✓ Branch 1 taken 17557 times.
19207 if( (max6 < fabs(q3_6_p0_6)) )
2381 {
2382 1650 max6 = fabs(q3_6_p0_6);
2383 }
2384
2/2
✓ Branch 0 taken 3044 times.
✓ Branch 1 taken 16163 times.
19207 if( (max6 < fabs(q3_7_p0_7)) )
2385 {
2386 3044 max6 = fabs(q3_7_p0_7);
2387 }
2388 double lower_bound_1;
2389 double upper_bound_1;
2390 int Delta_sign;
2391 int int_tmp_result;
2392 19207 lower_bound_1 = max4;
2393 19207 upper_bound_1 = max4;
2394
2/2
✓ Branch 0 taken 18014 times.
✓ Branch 1 taken 1193 times.
19207 if( (max1 < lower_bound_1) )
2395 {
2396 18014 lower_bound_1 = max1;
2397 }
2398 else
2399 {
2400
1/2
✓ Branch 0 taken 1193 times.
✗ Branch 1 not taken.
1193 if( (max1 > upper_bound_1) )
2401 {
2402 1193 upper_bound_1 = max1;
2403 }
2404 }
2405
2/2
✓ Branch 0 taken 9220 times.
✓ Branch 1 taken 9987 times.
19207 if( (max3 < lower_bound_1) )
2406 {
2407 9220 lower_bound_1 = max3;
2408 }
2409 else
2410 {
2411
2/2
✓ Branch 0 taken 904 times.
✓ Branch 1 taken 9083 times.
9987 if( (max3 > upper_bound_1) )
2412 {
2413 904 upper_bound_1 = max3;
2414 }
2415 }
2416
2/2
✓ Branch 0 taken 5853 times.
✓ Branch 1 taken 13354 times.
19207 if( (max2 < lower_bound_1) )
2417 {
2418 5853 lower_bound_1 = max2;
2419 }
2420 else
2421 {
2422
2/2
✓ Branch 0 taken 601 times.
✓ Branch 1 taken 12753 times.
13354 if( (max2 > upper_bound_1) )
2423 {
2424 601 upper_bound_1 = max2;
2425 }
2426 }
2427
2/2
✓ Branch 0 taken 503 times.
✓ Branch 1 taken 18704 times.
19207 if( (max6 < lower_bound_1) )
2428 {
2429 503 lower_bound_1 = max6;
2430 }
2431 else
2432 {
2433
2/2
✓ Branch 0 taken 9753 times.
✓ Branch 1 taken 8951 times.
18704 if( (max6 > upper_bound_1) )
2434 {
2435 9753 upper_bound_1 = max6;
2436 }
2437 }
2438
2/2
✓ Branch 0 taken 325 times.
✓ Branch 1 taken 18882 times.
19207 if( (max5 < lower_bound_1) )
2439 {
2440 325 lower_bound_1 = max5;
2441 }
2442 else
2443 {
2444
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 18882 times.
18882 if( (max5 > upper_bound_1) )
2445 {
2446 upper_bound_1 = max5;
2447 }
2448 }
2449
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19207 times.
19207 if( (lower_bound_1 < 2.82528483194754087282e-50) )
2450 {
2451 return FPG_UNCERTAIN_VALUE;
2452 }
2453 else
2454 {
2455
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19207 times.
19207 if( (upper_bound_1 > 4.83570327845851562508e+24) )
2456 {
2457 return FPG_UNCERTAIN_VALUE;
2458 }
2459 19207 eps = (4.37492894694731040560e-11 * (((((max2 * max4) * max1) * max5) * max3) * max6));
2460
2/2
✓ Branch 0 taken 9673 times.
✓ Branch 1 taken 9534 times.
19207 if( (Delta > eps) )
2461 {
2462 9673 int_tmp_result = 1;
2463 }
2464 else
2465 {
2466
1/2
✓ Branch 0 taken 9534 times.
✗ Branch 1 not taken.
9534 if( (Delta < -eps) )
2467 {
2468 9534 int_tmp_result = -1;
2469 }
2470 else
2471 {
2472 return FPG_UNCERTAIN_VALUE;
2473 }
2474 }
2475 }
2476 19207 Delta_sign = int_tmp_result;
2477 int int_tmp_result_FFWKCAA;
2478 19207 double max7 = max1;
2479
2/2
✓ Branch 0 taken 9601 times.
✓ Branch 1 taken 9606 times.
19207 if( (max7 < max3) )
2480 {
2481 9601 max7 = max3;
2482 }
2483 19207 double max8 = max2;
2484
2/2
✓ Branch 0 taken 682 times.
✓ Branch 1 taken 18525 times.
19207 if( (max8 < fabs(p4_3_p0_3)) )
2485 {
2486 682 max8 = fabs(p4_3_p0_3);
2487 }
2488
2/2
✓ Branch 0 taken 7000 times.
✓ Branch 1 taken 12207 times.
19207 if( (max8 < fabs(p4_1_p0_1)) )
2489 {
2490 7000 max8 = fabs(p4_1_p0_1);
2491 }
2492
2/2
✓ Branch 0 taken 1065 times.
✓ Branch 1 taken 18142 times.
19207 if( (max8 < fabs(p4_4_p0_4)) )
2493 {
2494 1065 max8 = fabs(p4_4_p0_4);
2495 }
2496
2/2
✓ Branch 0 taken 60 times.
✓ Branch 1 taken 19147 times.
19207 if( (max8 < fabs(p4_0_p0_0)) )
2497 {
2498 60 max8 = fabs(p4_0_p0_0);
2499 }
2500
2/2
✓ Branch 0 taken 5590 times.
✓ Branch 1 taken 13617 times.
19207 if( (max8 < fabs(p4_2_p0_2)) )
2501 {
2502 5590 max8 = fabs(p4_2_p0_2);
2503 }
2504
2/2
✓ Branch 0 taken 692 times.
✓ Branch 1 taken 18515 times.
19207 if( (max8 < fabs(p4_5_p0_5)) )
2505 {
2506 692 max8 = fabs(p4_5_p0_5);
2507 }
2508
2/2
✓ Branch 0 taken 633 times.
✓ Branch 1 taken 18574 times.
19207 if( (max8 < fabs(p4_6_p0_6)) )
2509 {
2510 633 max8 = fabs(p4_6_p0_6);
2511 }
2512
2/2
✓ Branch 0 taken 455 times.
✓ Branch 1 taken 18752 times.
19207 if( (max8 < fabs(p4_7_p0_7)) )
2513 {
2514 455 max8 = fabs(p4_7_p0_7);
2515 }
2516
2/2
✓ Branch 0 taken 13073 times.
✓ Branch 1 taken 6134 times.
19207 if( (max7 < max8) )
2517 {
2518 13073 max7 = max8;
2519 }
2520
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19207 times.
19207 if( (max7 < max2) )
2521 {
2522 max7 = max2;
2523 }
2524
2/2
✓ Branch 0 taken 16992 times.
✓ Branch 1 taken 2215 times.
19207 if( (max7 < max6) )
2525 {
2526 16992 max7 = max6;
2527 }
2528 19207 double max9 = max8;
2529
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19207 times.
19207 if( (max9 < max2) )
2530 {
2531 max9 = max2;
2532 }
2533
2/2
✓ Branch 0 taken 16813 times.
✓ Branch 1 taken 2394 times.
19207 if( (max9 < max5) )
2534 {
2535 16813 max9 = max5;
2536 }
2537 19207 double max10 = max4;
2538 19207 double max11 = max4;
2539
2/2
✓ Branch 0 taken 5568 times.
✓ Branch 1 taken 13639 times.
19207 if( (max11 < max5) )
2540 {
2541 5568 max11 = max5;
2542 }
2543
2/2
✓ Branch 0 taken 5568 times.
✓ Branch 1 taken 13639 times.
19207 if( (max10 < max11) )
2544 {
2545 5568 max10 = max11;
2546 }
2547
2/2
✓ Branch 0 taken 561 times.
✓ Branch 1 taken 18646 times.
19207 if( (max10 < max2) )
2548 {
2549 561 max10 = max2;
2550 }
2551
2/2
✓ Branch 0 taken 5372 times.
✓ Branch 1 taken 13835 times.
19207 if( (max10 < max6) )
2552 {
2553 5372 max10 = max6;
2554 }
2555
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19207 times.
19207 if( (max10 < max5) )
2556 {
2557 max10 = max5;
2558 }
2559 19207 lower_bound_1 = max11;
2560 19207 upper_bound_1 = max11;
2561
2/2
✓ Branch 0 taken 6042 times.
✓ Branch 1 taken 13165 times.
19207 if( (max9 < lower_bound_1) )
2562 {
2563 6042 lower_bound_1 = max9;
2564 }
2565 else
2566 {
2567
2/2
✓ Branch 0 taken 1056 times.
✓ Branch 1 taken 12109 times.
13165 if( (max9 > upper_bound_1) )
2568 {
2569 1056 upper_bound_1 = max9;
2570 }
2571 }
2572
2/2
✓ Branch 0 taken 5714 times.
✓ Branch 1 taken 13493 times.
19207 if( (max7 < lower_bound_1) )
2573 {
2574 5714 lower_bound_1 = max7;
2575 }
2576 else
2577 {
2578
2/2
✓ Branch 0 taken 5372 times.
✓ Branch 1 taken 8121 times.
13493 if( (max7 > upper_bound_1) )
2579 {
2580 5372 upper_bound_1 = max7;
2581 }
2582 }
2583
2/2
✓ Branch 0 taken 18250 times.
✓ Branch 1 taken 957 times.
19207 if( (max1 < lower_bound_1) )
2584 {
2585 18250 lower_bound_1 = max1;
2586 }
2587
2/2
✓ Branch 0 taken 9454 times.
✓ Branch 1 taken 9753 times.
19207 if( (max3 < lower_bound_1) )
2588 {
2589 9454 lower_bound_1 = max3;
2590 }
2591
2/2
✓ Branch 0 taken 1474 times.
✓ Branch 1 taken 17733 times.
19207 if( (max8 < lower_bound_1) )
2592 {
2593 1474 lower_bound_1 = max8;
2594 }
2595
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19207 times.
19207 if( (max10 < lower_bound_1) )
2596 {
2597 lower_bound_1 = max10;
2598 }
2599 else
2600 {
2601
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19207 times.
19207 if( (max10 > upper_bound_1) )
2602 {
2603 upper_bound_1 = max10;
2604 }
2605 }
2606
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19207 times.
19207 if( (lower_bound_1 < 4.17402518597284772324e-38) )
2607 {
2608 return FPG_UNCERTAIN_VALUE;
2609 }
2610 else
2611 {
2612
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 19207 times.
19207 if( (upper_bound_1 > 4.83570327845851562508e+24) )
2613 {
2614 return FPG_UNCERTAIN_VALUE;
2615 }
2616 19207 eps = (2.41492645607254025015e-09 * (((((((max8 * max11) * max1) * max10) * max3) * max10) * max9) * max7));
2617
2/2
✓ Branch 0 taken 9614 times.
✓ Branch 1 taken 9593 times.
19207 if( (r > eps) )
2618 {
2619 9614 int_tmp_result_FFWKCAA = 1;
2620 }
2621 else
2622 {
2623
1/2
✓ Branch 0 taken 9593 times.
✗ Branch 1 not taken.
9593 if( (r < -eps) )
2624 {
2625 9593 int_tmp_result_FFWKCAA = -1;
2626 }
2627 else
2628 {
2629 return FPG_UNCERTAIN_VALUE;
2630 }
2631 }
2632 }
2633 19207 return (Delta_sign * int_tmp_result_FFWKCAA);
2634 }
2635