4inline int side4_3d_filter(
const double* p0,
const double* p1,
const double* p2,
const double* p3,
const double* p4) {
10 a13 = (p1[2] - p0[2]);
11 double p1_0_p0_0 = (p1[0] - p0[0]);
12 double p1_1_p0_1 = (p1[1] - p0[1]);
13 double p1_2_p0_2 = (p1[2] - p0[2]);
15 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));
17 a21 = (p2[0] - p0[0]);
19 a22 = (p2[1] - p0[1]);
21 a23 = (p2[2] - p0[2]);
22 double p2_0_p0_0 = (p2[0] - p0[0]);
23 double p2_1_p0_1 = (p2[1] - p0[1]);
24 double p2_2_p0_2 = (p2[2] - p0[2]);
26 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));
28 a31 = (p3[0] - p0[0]);
30 a32 = (p3[1] - p0[1]);
32 a33 = (p3[2] - p0[2]);
33 double p3_0_p0_0 = (p3[0] - p0[0]);
34 double p3_1_p0_1 = (p3[1] - p0[1]);
35 double p3_2_p0_2 = (p3[2] - p0[2]);
37 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));
39 a41 = (p4[0] - p0[0]);
41 a42 = (p4[1] - p0[1]);
43 a43 = (p4[2] - p0[2]);
44 double p4_0_p0_0 = (p4[0] - p0[0]);
45 double p4_1_p0_1 = (p4[1] - p0[1]);
46 double p4_2_p0_2 = (p4[2] - p0[2]);
48 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));
50 Delta1 = (((a21 * ((a32 * a43) - (a33 * a42))) - (a31 * ((a22 * a43) - (a23 * a42)))) + (a41 * ((a22 * a33) - (a23 * a32))));
52 Delta2 = (((a11 * ((a32 * a43) - (a33 * a42))) - (a31 * ((a12 * a43) - (a13 * a42)))) + (a41 * ((a12 * a33) - (a13 * a32))));
54 Delta3 = (((a11 * ((a22 * a43) - (a23 * a42))) - (a21 * ((a12 * a43) - (a13 * a42)))) + (a41 * ((a12 * a23) - (a13 * a22))));
56 Delta4 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
58 r = ((((Delta1 * a14) - (Delta2 * a24)) + (Delta3 * a34)) - (Delta4 * a44));
60 double max1 = fabs(a11);
61 if( (max1 < fabs(a21)) )
65 if( (max1 < fabs(a31)) )
69 double max2 = fabs(a12);
70 if( (max2 < fabs(a13)) )
74 if( (max2 < fabs(a22)) )
78 if( (max2 < fabs(a23)) )
82 double max3 = fabs(a22);
83 if( (max3 < fabs(a23)) )
87 if( (max3 < fabs(a32)) )
91 if( (max3 < fabs(a33)) )
100 upper_bound_1 = max1;
101 if( (max3 < lower_bound_1) )
103 lower_bound_1 = max3;
107 if( (max3 > upper_bound_1) )
109 upper_bound_1 = max3;
112 if( (max2 < lower_bound_1) )
114 lower_bound_1 = max2;
118 if( (max2 > upper_bound_1) )
120 upper_bound_1 = max2;
123 if( (lower_bound_1 < 1.63288018496748314939e-98) )
125 return FPG_UNCERTAIN_VALUE;
129 if( (upper_bound_1 > 3.21387608851797948065e+60) )
131 return FPG_UNCERTAIN_VALUE;
133 eps = (5.11071278299732992696e-15 * ((max2 * max3) * max1));
140 if( (Delta4 < -eps) )
146 return FPG_UNCERTAIN_VALUE;
150 Delta4_sign = int_tmp_result;
151 int int_tmp_result_FFWKCAA;
153 if( (max4 < fabs(a41)) )
163 if( (max6 < fabs(a42)) )
167 if( (max6 < fabs(a43)) )
171 double max7 = fabs(p1_0_p0_0);
172 if( (max7 < fabs(p1_1_p0_1)) )
174 max7 = fabs(p1_1_p0_1);
176 if( (max7 < fabs(p1_2_p0_2)) )
178 max7 = fabs(p1_2_p0_2);
180 if( (max7 < fabs(p2_0_p0_0)) )
182 max7 = fabs(p2_0_p0_0);
184 if( (max7 < fabs(p2_2_p0_2)) )
186 max7 = fabs(p2_2_p0_2);
188 if( (max7 < fabs(p2_1_p0_1)) )
190 max7 = fabs(p2_1_p0_1);
192 if( (max7 < fabs(p3_0_p0_0)) )
194 max7 = fabs(p3_0_p0_0);
196 if( (max7 < fabs(p3_1_p0_1)) )
198 max7 = fabs(p3_1_p0_1);
200 if( (max7 < fabs(p3_2_p0_2)) )
202 max7 = fabs(p3_2_p0_2);
204 if( (max7 < fabs(p4_0_p0_0)) )
206 max7 = fabs(p4_0_p0_0);
208 if( (max7 < fabs(p4_1_p0_1)) )
210 max7 = fabs(p4_1_p0_1);
212 if( (max7 < fabs(p4_2_p0_2)) )
214 max7 = fabs(p4_2_p0_2);
216 lower_bound_1 = max7;
217 upper_bound_1 = max7;
218 if( (max4 < lower_bound_1) )
220 lower_bound_1 = max4;
224 if( (max4 > upper_bound_1) )
226 upper_bound_1 = max4;
229 if( (max5 < lower_bound_1) )
231 lower_bound_1 = max5;
235 if( (max5 > upper_bound_1) )
237 upper_bound_1 = max5;
240 if( (max6 < lower_bound_1) )
242 lower_bound_1 = max6;
246 if( (max6 > upper_bound_1) )
248 upper_bound_1 = max6;
251 if( (lower_bound_1 < 1.12285198342304832993e-59) )
253 return FPG_UNCERTAIN_VALUE;
257 if( (upper_bound_1 > 3.21387608851797948065e+60) )
259 return FPG_UNCERTAIN_VALUE;
261 eps = (1.24661365310273025710e-13 * ((((max5 * max6) * max4) * max7) * max7));
264 int_tmp_result_FFWKCAA = 1;
270 int_tmp_result_FFWKCAA = -1;
274 return FPG_UNCERTAIN_VALUE;
278 return (Delta4_sign * int_tmp_result_FFWKCAA);
282inline 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 double p1_0_p0_0 = (p1[0] - p0[0]);
284 double p1_1_p0_1 = (p1[1] - p0[1]);
285 double p1_2_p0_2 = (p1[2] - p0[2]);
286 double p1_3_p0_3 = (p1[3] - p0[3]);
288 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 double p2_0_p0_0 = (p2[0] - p0[0]);
290 double p2_1_p0_1 = (p2[1] - p0[1]);
291 double p2_2_p0_2 = (p2[2] - p0[2]);
292 double p2_3_p0_3 = (p2[3] - p0[3]);
294 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 double p3_0_p0_0 = (p3[0] - p0[0]);
296 double p3_1_p0_1 = (p3[1] - p0[1]);
297 double p3_2_p0_2 = (p3[2] - p0[2]);
298 double p3_3_p0_3 = (p3[3] - p0[3]);
300 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 double p4_0_p0_0 = (p4[0] - p0[0]);
302 double p4_1_p0_1 = (p4[1] - p0[1]);
303 double p4_2_p0_2 = (p4[2] - p0[2]);
304 double p4_3_p0_3 = (p4[3] - p0[3]);
306 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 double q0_0_p0_0 = (q0[0] - p0[0]);
308 double q0_1_p0_1 = (q0[1] - p0[1]);
309 double q0_2_p0_2 = (q0[2] - p0[2]);
310 double q0_3_p0_3 = (q0[3] - p0[3]);
312 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 double q1_0_p0_0 = (q1[0] - p0[0]);
314 double q1_1_p0_1 = (q1[1] - p0[1]);
315 double q1_2_p0_2 = (q1[2] - p0[2]);
316 double q1_3_p0_3 = (q1[3] - p0[3]);
318 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 double q2_0_p0_0 = (q2[0] - p0[0]);
320 double q2_1_p0_1 = (q2[1] - p0[1]);
321 double q2_2_p0_2 = (q2[2] - p0[2]);
322 double q2_3_p0_3 = (q2[3] - p0[3]);
324 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 double q3_0_p0_0 = (q3[0] - p0[0]);
326 double q3_1_p0_1 = (q3[1] - p0[1]);
327 double q3_2_p0_2 = (q3[2] - p0[2]);
328 double q3_3_p0_3 = (q3[3] - p0[3]);
330 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)));
332 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)));
334 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)));
336 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)));
338 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)));
340 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)));
342 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)));
344 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)));
346 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)));
348 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)));
350 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)));
352 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)));
354 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)));
356 b00 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
358 b01 = -((((a22 * a33) - (a23 * a32)) + ((a23 * a31) - (a21 * a33))) + ((a21 * a32) - (a22 * a31)));
360 b02 = ((((a12 * a33) - (a13 * a32)) + ((a13 * a31) - (a11 * a33))) + ((a11 * a32) - (a12 * a31)));
362 b03 = -((((a12 * a23) - (a13 * a22)) + ((a13 * a21) - (a11 * a23))) + ((a11 * a22) - (a12 * a21)));
364 b10 = -(((a10 * ((a22 * a33) - (a23 * a32))) - (a20 * ((a12 * a33) - (a13 * a32)))) + (a30 * ((a12 * a23) - (a13 * a22))));
366 b11 = ((((a22 * a33) - (a23 * a32)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a32) - (a22 * a30)));
368 b12 = -((((a12 * a33) - (a13 * a32)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a32) - (a12 * a30)));
370 b13 = ((((a12 * a23) - (a13 * a22)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a22) - (a12 * a20)));
372 b20 = (((a10 * ((a21 * a33) - (a23 * a31))) - (a20 * ((a11 * a33) - (a13 * a31)))) + (a30 * ((a11 * a23) - (a13 * a21))));
374 b21 = -((((a21 * a33) - (a23 * a31)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a31) - (a21 * a30)));
376 b22 = ((((a11 * a33) - (a13 * a31)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a31) - (a11 * a30)));
378 b23 = -((((a11 * a23) - (a13 * a21)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a21) - (a11 * a20)));
380 b30 = -(((a10 * ((a21 * a32) - (a22 * a31))) - (a20 * ((a11 * a32) - (a12 * a31)))) + (a30 * ((a11 * a22) - (a12 * a21))));
382 b31 = ((((a21 * a32) - (a22 * a31)) + ((a22 * a30) - (a20 * a32))) + ((a20 * a31) - (a21 * a30)));
384 b32 = -((((a11 * a32) - (a12 * a31)) + ((a12 * a30) - (a10 * a32))) + ((a10 * a31) - (a11 * a30)));
386 b33 = ((((a11 * a22) - (a12 * a21)) + ((a12 * a20) - (a10 * a22))) + ((a10 * a21) - (a11 * a20)));
388 Delta = (((b00 + b10) + b20) + b30);
390 DeltaLambda0 = ((((b01 * l1) + (b02 * l2)) + (b03 * l3)) + b00);
392 DeltaLambda1 = ((((b11 * l1) + (b12 * l2)) + (b13 * l3)) + b10);
394 DeltaLambda2 = ((((b21 * l1) + (b22 * l2)) + (b23 * l3)) + b20);
396 DeltaLambda3 = ((((b31 * l1) + (b32 * l2)) + (b33 * l3)) + b30);
398 r = ((Delta * l4) - ((((a40 * DeltaLambda0) + (a41 * DeltaLambda1)) + (a42 * DeltaLambda2)) + (a43 * DeltaLambda3)));
400 double max1 = fabs(p1_3_p0_3);
401 if( (max1 < fabs(p1_0_p0_0)) )
403 max1 = fabs(p1_0_p0_0);
405 if( (max1 < fabs(p1_1_p0_1)) )
407 max1 = fabs(p1_1_p0_1);
409 if( (max1 < fabs(p1_2_p0_2)) )
411 max1 = fabs(p1_2_p0_2);
413 double max2 = fabs(p2_2_p0_2);
414 if( (max2 < fabs(p2_1_p0_1)) )
416 max2 = fabs(p2_1_p0_1);
418 if( (max2 < fabs(p2_3_p0_3)) )
420 max2 = fabs(p2_3_p0_3);
422 if( (max2 < fabs(p2_0_p0_0)) )
424 max2 = fabs(p2_0_p0_0);
426 double max3 = fabs(p3_0_p0_0);
427 if( (max3 < fabs(p3_1_p0_1)) )
429 max3 = fabs(p3_1_p0_1);
431 if( (max3 < fabs(p3_3_p0_3)) )
433 max3 = fabs(p3_3_p0_3);
435 if( (max3 < fabs(p3_2_p0_2)) )
437 max3 = fabs(p3_2_p0_2);
439 double max4 = fabs(q0_3_p0_3);
440 if( (max4 < fabs(q0_0_p0_0)) )
442 max4 = fabs(q0_0_p0_0);
444 if( (max4 < fabs(q0_1_p0_1)) )
446 max4 = fabs(q0_1_p0_1);
448 if( (max4 < fabs(q0_2_p0_2)) )
450 max4 = fabs(q0_2_p0_2);
452 if( (max4 < fabs(q1_0_p0_0)) )
454 max4 = fabs(q1_0_p0_0);
456 if( (max4 < fabs(q1_1_p0_1)) )
458 max4 = fabs(q1_1_p0_1);
460 if( (max4 < fabs(q1_2_p0_2)) )
462 max4 = fabs(q1_2_p0_2);
464 if( (max4 < fabs(q1_3_p0_3)) )
466 max4 = fabs(q1_3_p0_3);
468 double max5 = fabs(q1_0_p0_0);
469 if( (max5 < fabs(q1_1_p0_1)) )
471 max5 = fabs(q1_1_p0_1);
473 if( (max5 < fabs(q1_2_p0_2)) )
475 max5 = fabs(q1_2_p0_2);
477 if( (max5 < fabs(q1_3_p0_3)) )
479 max5 = fabs(q1_3_p0_3);
481 if( (max5 < fabs(q2_0_p0_0)) )
483 max5 = fabs(q2_0_p0_0);
485 if( (max5 < fabs(q2_1_p0_1)) )
487 max5 = fabs(q2_1_p0_1);
489 if( (max5 < fabs(q2_2_p0_2)) )
491 max5 = fabs(q2_2_p0_2);
493 if( (max5 < fabs(q2_3_p0_3)) )
495 max5 = fabs(q2_3_p0_3);
497 double max6 = fabs(q2_0_p0_0);
498 if( (max6 < fabs(q2_1_p0_1)) )
500 max6 = fabs(q2_1_p0_1);
502 if( (max6 < fabs(q2_2_p0_2)) )
504 max6 = fabs(q2_2_p0_2);
506 if( (max6 < fabs(q2_3_p0_3)) )
508 max6 = fabs(q2_3_p0_3);
510 if( (max6 < fabs(q3_0_p0_0)) )
512 max6 = fabs(q3_0_p0_0);
514 if( (max6 < fabs(q3_1_p0_1)) )
516 max6 = fabs(q3_1_p0_1);
518 if( (max6 < fabs(q3_2_p0_2)) )
520 max6 = fabs(q3_2_p0_2);
522 if( (max6 < fabs(q3_3_p0_3)) )
524 max6 = fabs(q3_3_p0_3);
526 double lower_bound_1;
527 double upper_bound_1;
530 lower_bound_1 = max1;
531 upper_bound_1 = max1;
532 if( (max2 < lower_bound_1) )
534 lower_bound_1 = max2;
538 if( (max2 > upper_bound_1) )
540 upper_bound_1 = max2;
543 if( (max3 < lower_bound_1) )
545 lower_bound_1 = max3;
549 if( (max3 > upper_bound_1) )
551 upper_bound_1 = max3;
554 if( (max4 < lower_bound_1) )
556 lower_bound_1 = max4;
560 if( (max4 > upper_bound_1) )
562 upper_bound_1 = max4;
565 if( (max5 < lower_bound_1) )
567 lower_bound_1 = max5;
571 if( (max5 > upper_bound_1) )
573 upper_bound_1 = max5;
576 if( (max6 < lower_bound_1) )
578 lower_bound_1 = max6;
582 if( (max6 > upper_bound_1) )
584 upper_bound_1 = max6;
587 if( (lower_bound_1 < 4.14607644401726239868e-50) )
589 return FPG_UNCERTAIN_VALUE;
593 if( (upper_bound_1 > 4.83570327845851562508e+24) )
595 return FPG_UNCERTAIN_VALUE;
597 eps = (4.38046888801178809320e-12 * (((((max1 * max4) * max2) * max5) * max3) * max6));
610 return FPG_UNCERTAIN_VALUE;
614 Delta_sign = int_tmp_result;
615 int int_tmp_result_FFWKCAA;
630 if( (max8 < fabs(p4_1_p0_1)) )
632 max8 = fabs(p4_1_p0_1);
634 if( (max8 < fabs(p4_2_p0_2)) )
636 max8 = fabs(p4_2_p0_2);
638 if( (max8 < fabs(p4_0_p0_0)) )
640 max8 = fabs(p4_0_p0_0);
642 if( (max8 < fabs(p4_3_p0_3)) )
644 max8 = fabs(p4_3_p0_3);
682 lower_bound_1 = max10;
683 upper_bound_1 = max10;
684 if( (max11 < lower_bound_1) )
686 lower_bound_1 = max11;
690 if( (max11 > upper_bound_1) )
692 upper_bound_1 = max11;
695 if( (max2 < lower_bound_1) )
697 lower_bound_1 = max2;
699 if( (max3 < lower_bound_1) )
701 lower_bound_1 = max3;
703 if( (max7 < lower_bound_1) )
705 lower_bound_1 = max7;
709 if( (max7 > upper_bound_1) )
711 upper_bound_1 = max7;
714 if( (max8 < lower_bound_1) )
716 lower_bound_1 = max8;
718 if( (max9 < lower_bound_1) )
720 lower_bound_1 = max9;
724 if( (max9 > upper_bound_1) )
726 upper_bound_1 = max9;
729 if( (lower_bound_1 < 6.06263132863556750071e-38) )
731 return FPG_UNCERTAIN_VALUE;
735 if( (upper_bound_1 > 4.83570327845851562508e+24) )
737 return FPG_UNCERTAIN_VALUE;
739 eps = (1.21914442286604163181e-10 * (((((((max8 * max11) * max2) * max10) * max3) * max10) * max9) * max7));
742 int_tmp_result_FFWKCAA = 1;
748 int_tmp_result_FFWKCAA = -1;
752 return FPG_UNCERTAIN_VALUE;
756 return (Delta_sign * int_tmp_result_FFWKCAA);
760inline 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 double p1_0_p0_0 = (p1[0] - p0[0]);
762 double p1_1_p0_1 = (p1[1] - p0[1]);
763 double p1_2_p0_2 = (p1[2] - p0[2]);
764 double p1_3_p0_3 = (p1[3] - p0[3]);
765 double p1_4_p0_4 = (p1[4] - p0[4]);
766 double p1_5_p0_5 = (p1[5] - p0[5]);
768 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 double p2_0_p0_0 = (p2[0] - p0[0]);
770 double p2_1_p0_1 = (p2[1] - p0[1]);
771 double p2_2_p0_2 = (p2[2] - p0[2]);
772 double p2_3_p0_3 = (p2[3] - p0[3]);
773 double p2_4_p0_4 = (p2[4] - p0[4]);
774 double p2_5_p0_5 = (p2[5] - p0[5]);
776 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 double p3_0_p0_0 = (p3[0] - p0[0]);
778 double p3_1_p0_1 = (p3[1] - p0[1]);
779 double p3_2_p0_2 = (p3[2] - p0[2]);
780 double p3_3_p0_3 = (p3[3] - p0[3]);
781 double p3_4_p0_4 = (p3[4] - p0[4]);
782 double p3_5_p0_5 = (p3[5] - p0[5]);
784 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 double p4_0_p0_0 = (p4[0] - p0[0]);
786 double p4_1_p0_1 = (p4[1] - p0[1]);
787 double p4_2_p0_2 = (p4[2] - p0[2]);
788 double p4_3_p0_3 = (p4[3] - p0[3]);
789 double p4_4_p0_4 = (p4[4] - p0[4]);
790 double p4_5_p0_5 = (p4[5] - p0[5]);
792 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 double q0_0_p0_0 = (q0[0] - p0[0]);
794 double q0_1_p0_1 = (q0[1] - p0[1]);
795 double q0_2_p0_2 = (q0[2] - p0[2]);
796 double q0_3_p0_3 = (q0[3] - p0[3]);
797 double q0_4_p0_4 = (q0[4] - p0[4]);
798 double q0_5_p0_5 = (q0[5] - p0[5]);
800 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 double q1_0_p0_0 = (q1[0] - p0[0]);
802 double q1_1_p0_1 = (q1[1] - p0[1]);
803 double q1_2_p0_2 = (q1[2] - p0[2]);
804 double q1_3_p0_3 = (q1[3] - p0[3]);
805 double q1_4_p0_4 = (q1[4] - p0[4]);
806 double q1_5_p0_5 = (q1[5] - p0[5]);
808 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 double q2_0_p0_0 = (q2[0] - p0[0]);
810 double q2_1_p0_1 = (q2[1] - p0[1]);
811 double q2_2_p0_2 = (q2[2] - p0[2]);
812 double q2_3_p0_3 = (q2[3] - p0[3]);
813 double q2_4_p0_4 = (q2[4] - p0[4]);
814 double q2_5_p0_5 = (q2[5] - p0[5]);
816 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 double q3_0_p0_0 = (q3[0] - p0[0]);
818 double q3_1_p0_1 = (q3[1] - p0[1]);
819 double q3_2_p0_2 = (q3[2] - p0[2]);
820 double q3_3_p0_3 = (q3[3] - p0[3]);
821 double q3_4_p0_4 = (q3[4] - p0[4]);
822 double q3_5_p0_5 = (q3[5] - p0[5]);
824 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)));
826 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)));
828 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)));
830 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)));
832 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)));
834 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)));
836 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)));
838 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)));
840 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)));
842 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)));
844 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)));
846 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)));
848 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)));
850 b00 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
852 b01 = -((((a22 * a33) - (a23 * a32)) + ((a23 * a31) - (a21 * a33))) + ((a21 * a32) - (a22 * a31)));
854 b02 = ((((a12 * a33) - (a13 * a32)) + ((a13 * a31) - (a11 * a33))) + ((a11 * a32) - (a12 * a31)));
856 b03 = -((((a12 * a23) - (a13 * a22)) + ((a13 * a21) - (a11 * a23))) + ((a11 * a22) - (a12 * a21)));
858 b10 = -(((a10 * ((a22 * a33) - (a23 * a32))) - (a20 * ((a12 * a33) - (a13 * a32)))) + (a30 * ((a12 * a23) - (a13 * a22))));
860 b11 = ((((a22 * a33) - (a23 * a32)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a32) - (a22 * a30)));
862 b12 = -((((a12 * a33) - (a13 * a32)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a32) - (a12 * a30)));
864 b13 = ((((a12 * a23) - (a13 * a22)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a22) - (a12 * a20)));
866 b20 = (((a10 * ((a21 * a33) - (a23 * a31))) - (a20 * ((a11 * a33) - (a13 * a31)))) + (a30 * ((a11 * a23) - (a13 * a21))));
868 b21 = -((((a21 * a33) - (a23 * a31)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a31) - (a21 * a30)));
870 b22 = ((((a11 * a33) - (a13 * a31)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a31) - (a11 * a30)));
872 b23 = -((((a11 * a23) - (a13 * a21)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a21) - (a11 * a20)));
874 b30 = -(((a10 * ((a21 * a32) - (a22 * a31))) - (a20 * ((a11 * a32) - (a12 * a31)))) + (a30 * ((a11 * a22) - (a12 * a21))));
876 b31 = ((((a21 * a32) - (a22 * a31)) + ((a22 * a30) - (a20 * a32))) + ((a20 * a31) - (a21 * a30)));
878 b32 = -((((a11 * a32) - (a12 * a31)) + ((a12 * a30) - (a10 * a32))) + ((a10 * a31) - (a11 * a30)));
880 b33 = ((((a11 * a22) - (a12 * a21)) + ((a12 * a20) - (a10 * a22))) + ((a10 * a21) - (a11 * a20)));
882 Delta = (((b00 + b10) + b20) + b30);
884 DeltaLambda0 = ((((b01 * l1) + (b02 * l2)) + (b03 * l3)) + b00);
886 DeltaLambda1 = ((((b11 * l1) + (b12 * l2)) + (b13 * l3)) + b10);
888 DeltaLambda2 = ((((b21 * l1) + (b22 * l2)) + (b23 * l3)) + b20);
890 DeltaLambda3 = ((((b31 * l1) + (b32 * l2)) + (b33 * l3)) + b30);
892 r = ((Delta * l4) - ((((a40 * DeltaLambda0) + (a41 * DeltaLambda1)) + (a42 * DeltaLambda2)) + (a43 * DeltaLambda3)));
894 double max1 = fabs(p3_2_p0_2);
895 if( (max1 < fabs(p3_0_p0_0)) )
897 max1 = fabs(p3_0_p0_0);
899 if( (max1 < fabs(p3_3_p0_3)) )
901 max1 = fabs(p3_3_p0_3);
903 if( (max1 < fabs(p3_4_p0_4)) )
905 max1 = fabs(p3_4_p0_4);
907 if( (max1 < fabs(p3_1_p0_1)) )
909 max1 = fabs(p3_1_p0_1);
911 if( (max1 < fabs(p3_5_p0_5)) )
913 max1 = fabs(p3_5_p0_5);
915 double max2 = fabs(p2_1_p0_1);
916 if( (max2 < fabs(p2_4_p0_4)) )
918 max2 = fabs(p2_4_p0_4);
920 if( (max2 < fabs(p2_2_p0_2)) )
922 max2 = fabs(p2_2_p0_2);
924 if( (max2 < fabs(p2_0_p0_0)) )
926 max2 = fabs(p2_0_p0_0);
928 if( (max2 < fabs(p2_3_p0_3)) )
930 max2 = fabs(p2_3_p0_3);
932 if( (max2 < fabs(p2_5_p0_5)) )
934 max2 = fabs(p2_5_p0_5);
936 double max3 = fabs(p1_0_p0_0);
937 if( (max3 < fabs(p1_1_p0_1)) )
939 max3 = fabs(p1_1_p0_1);
941 if( (max3 < fabs(p1_2_p0_2)) )
943 max3 = fabs(p1_2_p0_2);
945 if( (max3 < fabs(p1_3_p0_3)) )
947 max3 = fabs(p1_3_p0_3);
949 if( (max3 < fabs(p1_4_p0_4)) )
951 max3 = fabs(p1_4_p0_4);
953 if( (max3 < fabs(p1_5_p0_5)) )
955 max3 = fabs(p1_5_p0_5);
957 double max4 = fabs(q0_0_p0_0);
958 if( (max4 < fabs(q0_1_p0_1)) )
960 max4 = fabs(q0_1_p0_1);
962 if( (max4 < fabs(q0_2_p0_2)) )
964 max4 = fabs(q0_2_p0_2);
966 if( (max4 < fabs(q0_3_p0_3)) )
968 max4 = fabs(q0_3_p0_3);
970 if( (max4 < fabs(q0_4_p0_4)) )
972 max4 = fabs(q0_4_p0_4);
974 if( (max4 < fabs(q0_5_p0_5)) )
976 max4 = fabs(q0_5_p0_5);
978 if( (max4 < fabs(q1_0_p0_0)) )
980 max4 = fabs(q1_0_p0_0);
982 if( (max4 < fabs(q1_1_p0_1)) )
984 max4 = fabs(q1_1_p0_1);
986 if( (max4 < fabs(q1_2_p0_2)) )
988 max4 = fabs(q1_2_p0_2);
990 if( (max4 < fabs(q1_3_p0_3)) )
992 max4 = fabs(q1_3_p0_3);
994 if( (max4 < fabs(q1_4_p0_4)) )
996 max4 = fabs(q1_4_p0_4);
998 if( (max4 < fabs(q1_5_p0_5)) )
1000 max4 = fabs(q1_5_p0_5);
1002 double max5 = fabs(q1_0_p0_0);
1003 if( (max5 < fabs(q1_1_p0_1)) )
1005 max5 = fabs(q1_1_p0_1);
1007 if( (max5 < fabs(q1_2_p0_2)) )
1009 max5 = fabs(q1_2_p0_2);
1011 if( (max5 < fabs(q1_3_p0_3)) )
1013 max5 = fabs(q1_3_p0_3);
1015 if( (max5 < fabs(q1_4_p0_4)) )
1017 max5 = fabs(q1_4_p0_4);
1019 if( (max5 < fabs(q1_5_p0_5)) )
1021 max5 = fabs(q1_5_p0_5);
1023 if( (max5 < fabs(q2_0_p0_0)) )
1025 max5 = fabs(q2_0_p0_0);
1027 if( (max5 < fabs(q2_1_p0_1)) )
1029 max5 = fabs(q2_1_p0_1);
1031 if( (max5 < fabs(q2_2_p0_2)) )
1033 max5 = fabs(q2_2_p0_2);
1035 if( (max5 < fabs(q2_3_p0_3)) )
1037 max5 = fabs(q2_3_p0_3);
1039 if( (max5 < fabs(q2_4_p0_4)) )
1041 max5 = fabs(q2_4_p0_4);
1043 if( (max5 < fabs(q2_5_p0_5)) )
1045 max5 = fabs(q2_5_p0_5);
1047 double max6 = fabs(q2_0_p0_0);
1048 if( (max6 < fabs(q2_1_p0_1)) )
1050 max6 = fabs(q2_1_p0_1);
1052 if( (max6 < fabs(q2_2_p0_2)) )
1054 max6 = fabs(q2_2_p0_2);
1056 if( (max6 < fabs(q2_3_p0_3)) )
1058 max6 = fabs(q2_3_p0_3);
1060 if( (max6 < fabs(q2_4_p0_4)) )
1062 max6 = fabs(q2_4_p0_4);
1064 if( (max6 < fabs(q2_5_p0_5)) )
1066 max6 = fabs(q2_5_p0_5);
1068 if( (max6 < fabs(q3_0_p0_0)) )
1070 max6 = fabs(q3_0_p0_0);
1072 if( (max6 < fabs(q3_1_p0_1)) )
1074 max6 = fabs(q3_1_p0_1);
1076 if( (max6 < fabs(q3_2_p0_2)) )
1078 max6 = fabs(q3_2_p0_2);
1080 if( (max6 < fabs(q3_3_p0_3)) )
1082 max6 = fabs(q3_3_p0_3);
1084 if( (max6 < fabs(q3_4_p0_4)) )
1086 max6 = fabs(q3_4_p0_4);
1088 if( (max6 < fabs(q3_5_p0_5)) )
1090 max6 = fabs(q3_5_p0_5);
1092 double lower_bound_1;
1093 double upper_bound_1;
1096 lower_bound_1 = max4;
1097 upper_bound_1 = max4;
1098 if( (max1 < lower_bound_1) )
1100 lower_bound_1 = max1;
1104 if( (max1 > upper_bound_1) )
1106 upper_bound_1 = max1;
1109 if( (max2 < lower_bound_1) )
1111 lower_bound_1 = max2;
1115 if( (max2 > upper_bound_1) )
1117 upper_bound_1 = max2;
1120 if( (max3 < lower_bound_1) )
1122 lower_bound_1 = max3;
1126 if( (max3 > upper_bound_1) )
1128 upper_bound_1 = max3;
1131 if( (max5 < lower_bound_1) )
1133 lower_bound_1 = max5;
1137 if( (max5 > upper_bound_1) )
1139 upper_bound_1 = max5;
1142 if( (max6 < lower_bound_1) )
1144 lower_bound_1 = max6;
1148 if( (max6 > upper_bound_1) )
1150 upper_bound_1 = max6;
1153 if( (lower_bound_1 < 3.31864264949884013629e-50) )
1155 return FPG_UNCERTAIN_VALUE;
1159 if( (upper_bound_1 > 4.83570327845851562508e+24) )
1161 return FPG_UNCERTAIN_VALUE;
1163 eps = (1.66564133587113165316e-11 * (((((max3 * max4) * max2) * max5) * max1) * max6));
1170 if( (Delta < -eps) )
1172 int_tmp_result = -1;
1176 return FPG_UNCERTAIN_VALUE;
1180 Delta_sign = int_tmp_result;
1181 int int_tmp_result_FFWKCAA;
1192 if( (max8 < fabs(p4_5_p0_5)) )
1194 max8 = fabs(p4_5_p0_5);
1196 if( (max8 < fabs(p4_0_p0_0)) )
1198 max8 = fabs(p4_0_p0_0);
1200 if( (max8 < fabs(p4_1_p0_1)) )
1202 max8 = fabs(p4_1_p0_1);
1204 if( (max8 < fabs(p4_2_p0_2)) )
1206 max8 = fabs(p4_2_p0_2);
1208 if( (max8 < fabs(p4_3_p0_3)) )
1210 max8 = fabs(p4_3_p0_3);
1212 if( (max8 < fabs(p4_4_p0_4)) )
1214 max8 = fabs(p4_4_p0_4);
1233 double max10 = max4;
1234 if( (max10 < max3) )
1238 if( (max10 < max5) )
1242 double max11 = max4;
1243 if( (max11 < max5) )
1247 if( (max10 < max11) )
1251 if( (max10 < max6) )
1255 lower_bound_1 = max1;
1256 upper_bound_1 = max1;
1257 if( (max2 < lower_bound_1) )
1259 lower_bound_1 = max2;
1261 if( (max7 < lower_bound_1) )
1263 lower_bound_1 = max7;
1267 if( (max7 > upper_bound_1) )
1269 upper_bound_1 = max7;
1272 if( (max8 < lower_bound_1) )
1274 lower_bound_1 = max8;
1276 if( (max9 < lower_bound_1) )
1278 lower_bound_1 = max9;
1282 if( (max9 > upper_bound_1) )
1284 upper_bound_1 = max9;
1287 if( (max10 < lower_bound_1) )
1289 lower_bound_1 = max10;
1293 if( (max10 > upper_bound_1) )
1295 upper_bound_1 = max10;
1298 if( (max11 < lower_bound_1) )
1300 lower_bound_1 = max11;
1304 if( (max11 > upper_bound_1) )
1306 upper_bound_1 = max11;
1309 if( (lower_bound_1 < 4.87975611107819181771e-38) )
1311 return FPG_UNCERTAIN_VALUE;
1315 if( (upper_bound_1 > 4.83570327845851562508e+24) )
1317 return FPG_UNCERTAIN_VALUE;
1319 eps = (6.92085098542795335117e-10 * (((((((max8 * max11) * max2) * max10) * max1) * max10) * max9) * max7));
1322 int_tmp_result_FFWKCAA = 1;
1328 int_tmp_result_FFWKCAA = -1;
1332 return FPG_UNCERTAIN_VALUE;
1336 return (Delta_sign * int_tmp_result_FFWKCAA);
1340inline 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]);
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]);
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]);
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]);
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]);
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]);
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]);
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]);
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)));
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)));
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)));
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)));
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)));
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)));
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)));
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)));
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)));
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)));
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)));
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)));
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)));
1438 b00 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
1440 b01 = -((((a22 * a33) - (a23 * a32)) + ((a23 * a31) - (a21 * a33))) + ((a21 * a32) - (a22 * a31)));
1442 b02 = ((((a12 * a33) - (a13 * a32)) + ((a13 * a31) - (a11 * a33))) + ((a11 * a32) - (a12 * a31)));
1444 b03 = -((((a12 * a23) - (a13 * a22)) + ((a13 * a21) - (a11 * a23))) + ((a11 * a22) - (a12 * a21)));
1446 b10 = -(((a10 * ((a22 * a33) - (a23 * a32))) - (a20 * ((a12 * a33) - (a13 * a32)))) + (a30 * ((a12 * a23) - (a13 * a22))));
1448 b11 = ((((a22 * a33) - (a23 * a32)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a32) - (a22 * a30)));
1450 b12 = -((((a12 * a33) - (a13 * a32)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a32) - (a12 * a30)));
1452 b13 = ((((a12 * a23) - (a13 * a22)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a22) - (a12 * a20)));
1454 b20 = (((a10 * ((a21 * a33) - (a23 * a31))) - (a20 * ((a11 * a33) - (a13 * a31)))) + (a30 * ((a11 * a23) - (a13 * a21))));
1456 b21 = -((((a21 * a33) - (a23 * a31)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a31) - (a21 * a30)));
1458 b22 = ((((a11 * a33) - (a13 * a31)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a31) - (a11 * a30)));
1460 b23 = -((((a11 * a23) - (a13 * a21)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a21) - (a11 * a20)));
1462 b30 = -(((a10 * ((a21 * a32) - (a22 * a31))) - (a20 * ((a11 * a32) - (a12 * a31)))) + (a30 * ((a11 * a22) - (a12 * a21))));
1464 b31 = ((((a21 * a32) - (a22 * a31)) + ((a22 * a30) - (a20 * a32))) + ((a20 * a31) - (a21 * a30)));
1466 b32 = -((((a11 * a32) - (a12 * a31)) + ((a12 * a30) - (a10 * a32))) + ((a10 * a31) - (a11 * a30)));
1468 b33 = ((((a11 * a22) - (a12 * a21)) + ((a12 * a20) - (a10 * a22))) + ((a10 * a21) - (a11 * a20)));
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);
1480 r = ((Delta * l4) - ((((a40 * DeltaLambda0) + (a41 * DeltaLambda1)) + (a42 * DeltaLambda2)) + (a43 * DeltaLambda3)));
1482 double max1 = fabs(p1_0_p0_0);
1483 if( (max1 < fabs(p1_1_p0_1)) )
1485 max1 = fabs(p1_1_p0_1);
1487 if( (max1 < fabs(p1_2_p0_2)) )
1489 max1 = fabs(p1_2_p0_2);
1491 if( (max1 < fabs(p1_3_p0_3)) )
1493 max1 = fabs(p1_3_p0_3);
1495 if( (max1 < fabs(p1_4_p0_4)) )
1497 max1 = fabs(p1_4_p0_4);
1499 if( (max1 < fabs(p1_5_p0_5)) )
1501 max1 = fabs(p1_5_p0_5);
1503 if( (max1 < fabs(p1_6_p0_6)) )
1505 max1 = fabs(p1_6_p0_6);
1507 double max2 = fabs(p3_0_p0_0);
1508 if( (max2 < fabs(p3_4_p0_4)) )
1510 max2 = fabs(p3_4_p0_4);
1512 if( (max2 < fabs(p3_2_p0_2)) )
1514 max2 = fabs(p3_2_p0_2);
1516 if( (max2 < fabs(p3_1_p0_1)) )
1518 max2 = fabs(p3_1_p0_1);
1520 if( (max2 < fabs(p3_3_p0_3)) )
1522 max2 = fabs(p3_3_p0_3);
1524 if( (max2 < fabs(p3_5_p0_5)) )
1526 max2 = fabs(p3_5_p0_5);
1528 if( (max2 < fabs(p3_6_p0_6)) )
1530 max2 = fabs(p3_6_p0_6);
1532 double max3 = fabs(p2_5_p0_5);
1533 if( (max3 < fabs(p2_2_p0_2)) )
1535 max3 = fabs(p2_2_p0_2);
1537 if( (max3 < fabs(p2_3_p0_3)) )
1539 max3 = fabs(p2_3_p0_3);
1541 if( (max3 < fabs(p2_0_p0_0)) )
1543 max3 = fabs(p2_0_p0_0);
1545 if( (max3 < fabs(p2_1_p0_1)) )
1547 max3 = fabs(p2_1_p0_1);
1549 if( (max3 < fabs(p2_6_p0_6)) )
1551 max3 = fabs(p2_6_p0_6);
1553 if( (max3 < fabs(p2_4_p0_4)) )
1555 max3 = fabs(p2_4_p0_4);
1557 double max4 = fabs(q0_0_p0_0);
1558 if( (max4 < fabs(q0_1_p0_1)) )
1560 max4 = fabs(q0_1_p0_1);
1562 if( (max4 < fabs(q0_2_p0_2)) )
1564 max4 = fabs(q0_2_p0_2);
1566 if( (max4 < fabs(q0_3_p0_3)) )
1568 max4 = fabs(q0_3_p0_3);
1570 if( (max4 < fabs(q0_4_p0_4)) )
1572 max4 = fabs(q0_4_p0_4);
1574 if( (max4 < fabs(q0_5_p0_5)) )
1576 max4 = fabs(q0_5_p0_5);
1578 if( (max4 < fabs(q0_6_p0_6)) )
1580 max4 = fabs(q0_6_p0_6);
1582 if( (max4 < fabs(q1_0_p0_0)) )
1584 max4 = fabs(q1_0_p0_0);
1586 if( (max4 < fabs(q1_1_p0_1)) )
1588 max4 = fabs(q1_1_p0_1);
1590 if( (max4 < fabs(q1_2_p0_2)) )
1592 max4 = fabs(q1_2_p0_2);
1594 if( (max4 < fabs(q1_3_p0_3)) )
1596 max4 = fabs(q1_3_p0_3);
1598 if( (max4 < fabs(q1_4_p0_4)) )
1600 max4 = fabs(q1_4_p0_4);
1602 if( (max4 < fabs(q1_5_p0_5)) )
1604 max4 = fabs(q1_5_p0_5);
1606 if( (max4 < fabs(q1_6_p0_6)) )
1608 max4 = fabs(q1_6_p0_6);
1610 double max5 = fabs(q1_0_p0_0);
1611 if( (max5 < fabs(q1_1_p0_1)) )
1613 max5 = fabs(q1_1_p0_1);
1615 if( (max5 < fabs(q1_2_p0_2)) )
1617 max5 = fabs(q1_2_p0_2);
1619 if( (max5 < fabs(q1_3_p0_3)) )
1621 max5 = fabs(q1_3_p0_3);
1623 if( (max5 < fabs(q1_4_p0_4)) )
1625 max5 = fabs(q1_4_p0_4);
1627 if( (max5 < fabs(q1_5_p0_5)) )
1629 max5 = fabs(q1_5_p0_5);
1631 if( (max5 < fabs(q1_6_p0_6)) )
1633 max5 = fabs(q1_6_p0_6);
1635 if( (max5 < fabs(q2_0_p0_0)) )
1637 max5 = fabs(q2_0_p0_0);
1639 if( (max5 < fabs(q2_1_p0_1)) )
1641 max5 = fabs(q2_1_p0_1);
1643 if( (max5 < fabs(q2_2_p0_2)) )
1645 max5 = fabs(q2_2_p0_2);
1647 if( (max5 < fabs(q2_3_p0_3)) )
1649 max5 = fabs(q2_3_p0_3);
1651 if( (max5 < fabs(q2_4_p0_4)) )
1653 max5 = fabs(q2_4_p0_4);
1655 if( (max5 < fabs(q2_5_p0_5)) )
1657 max5 = fabs(q2_5_p0_5);
1659 if( (max5 < fabs(q2_6_p0_6)) )
1661 max5 = fabs(q2_6_p0_6);
1663 double max6 = fabs(q2_0_p0_0);
1664 if( (max6 < fabs(q2_1_p0_1)) )
1666 max6 = fabs(q2_1_p0_1);
1668 if( (max6 < fabs(q2_2_p0_2)) )
1670 max6 = fabs(q2_2_p0_2);
1672 if( (max6 < fabs(q2_3_p0_3)) )
1674 max6 = fabs(q2_3_p0_3);
1676 if( (max6 < fabs(q2_4_p0_4)) )
1678 max6 = fabs(q2_4_p0_4);
1680 if( (max6 < fabs(q2_5_p0_5)) )
1682 max6 = fabs(q2_5_p0_5);
1684 if( (max6 < fabs(q2_6_p0_6)) )
1686 max6 = fabs(q2_6_p0_6);
1688 if( (max6 < fabs(q3_0_p0_0)) )
1690 max6 = fabs(q3_0_p0_0);
1692 if( (max6 < fabs(q3_1_p0_1)) )
1694 max6 = fabs(q3_1_p0_1);
1696 if( (max6 < fabs(q3_2_p0_2)) )
1698 max6 = fabs(q3_2_p0_2);
1700 if( (max6 < fabs(q3_3_p0_3)) )
1702 max6 = fabs(q3_3_p0_3);
1704 if( (max6 < fabs(q3_4_p0_4)) )
1706 max6 = fabs(q3_4_p0_4);
1708 if( (max6 < fabs(q3_5_p0_5)) )
1710 max6 = fabs(q3_5_p0_5);
1712 if( (max6 < fabs(q3_6_p0_6)) )
1714 max6 = fabs(q3_6_p0_6);
1716 double lower_bound_1;
1717 double upper_bound_1;
1720 lower_bound_1 = max3;
1721 upper_bound_1 = max3;
1722 if( (max4 < lower_bound_1) )
1724 lower_bound_1 = max4;
1728 if( (max4 > upper_bound_1) )
1730 upper_bound_1 = max4;
1733 if( (max1 < lower_bound_1) )
1735 lower_bound_1 = max1;
1739 if( (max1 > upper_bound_1) )
1741 upper_bound_1 = max1;
1744 if( (max2 < lower_bound_1) )
1746 lower_bound_1 = max2;
1750 if( (max2 > upper_bound_1) )
1752 upper_bound_1 = max2;
1755 if( (max5 < lower_bound_1) )
1757 lower_bound_1 = max5;
1761 if( (max5 > upper_bound_1) )
1763 upper_bound_1 = max5;
1766 if( (max6 < lower_bound_1) )
1768 lower_bound_1 = max6;
1772 if( (max6 > upper_bound_1) )
1774 upper_bound_1 = max6;
1777 if( (lower_bound_1 < 3.04548303565602498901e-50) )
1779 return FPG_UNCERTAIN_VALUE;
1783 if( (upper_bound_1 > 4.83570327845851562508e+24) )
1785 return FPG_UNCERTAIN_VALUE;
1787 eps = (2.78873548804336160566e-11 * (((((max1 * max4) * max3) * max5) * max2) * max6));
1794 if( (Delta < -eps) )
1796 int_tmp_result = -1;
1800 return FPG_UNCERTAIN_VALUE;
1804 Delta_sign = int_tmp_result;
1805 int int_tmp_result_FFWKCAA;
1820 if( (max8 < fabs(p4_4_p0_4)) )
1822 max8 = fabs(p4_4_p0_4);
1824 if( (max8 < fabs(p4_1_p0_1)) )
1826 max8 = fabs(p4_1_p0_1);
1828 if( (max8 < fabs(p4_5_p0_5)) )
1830 max8 = fabs(p4_5_p0_5);
1832 if( (max8 < fabs(p4_0_p0_0)) )
1834 max8 = fabs(p4_0_p0_0);
1836 if( (max8 < fabs(p4_2_p0_2)) )
1838 max8 = fabs(p4_2_p0_2);
1840 if( (max8 < fabs(p4_3_p0_3)) )
1842 max8 = fabs(p4_3_p0_3);
1844 if( (max8 < fabs(p4_6_p0_6)) )
1846 max8 = fabs(p4_6_p0_6);
1861 double max10 = max4;
1862 if( (max10 < max1) )
1866 if( (max10 < max5) )
1870 if( (max10 < max6) )
1874 double max11 = max4;
1875 if( (max11 < max5) )
1879 if( (max10 < max11) )
1883 lower_bound_1 = max3;
1884 upper_bound_1 = max3;
1885 if( (max2 < lower_bound_1) )
1887 lower_bound_1 = max2;
1889 if( (max7 < lower_bound_1) )
1891 lower_bound_1 = max7;
1895 if( (max7 > upper_bound_1) )
1897 upper_bound_1 = max7;
1900 if( (max8 < lower_bound_1) )
1902 lower_bound_1 = max8;
1904 if( (max9 < lower_bound_1) )
1906 lower_bound_1 = max9;
1910 if( (max9 > upper_bound_1) )
1912 upper_bound_1 = max9;
1915 if( (max10 < lower_bound_1) )
1917 lower_bound_1 = max10;
1921 if( (max10 > upper_bound_1) )
1923 upper_bound_1 = max10;
1926 if( (max11 < lower_bound_1) )
1928 lower_bound_1 = max11;
1932 if( (max11 > upper_bound_1) )
1934 upper_bound_1 = max11;
1937 if( (lower_bound_1 < 4.48906690519700369396e-38) )
1939 return FPG_UNCERTAIN_VALUE;
1943 if( (upper_bound_1 > 4.83570327845851562508e+24) )
1945 return FPG_UNCERTAIN_VALUE;
1947 eps = (1.34926049830188433875e-09 * (((((((max8 * max11) * max3) * max10) * max2) * max10) * max9) * max7));
1950 int_tmp_result_FFWKCAA = 1;
1956 int_tmp_result_FFWKCAA = -1;
1960 return FPG_UNCERTAIN_VALUE;
1964 return (Delta_sign * int_tmp_result_FFWKCAA);
1968inline 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 double p1_0_p0_0 = (p1[0] - p0[0]);
1970 double p1_1_p0_1 = (p1[1] - p0[1]);
1971 double p1_2_p0_2 = (p1[2] - p0[2]);
1972 double p1_3_p0_3 = (p1[3] - p0[3]);
1973 double p1_4_p0_4 = (p1[4] - p0[4]);
1974 double p1_5_p0_5 = (p1[5] - p0[5]);
1975 double p1_6_p0_6 = (p1[6] - p0[6]);
1976 double p1_7_p0_7 = (p1[7] - p0[7]);
1978 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 double p2_0_p0_0 = (p2[0] - p0[0]);
1980 double p2_1_p0_1 = (p2[1] - p0[1]);
1981 double p2_2_p0_2 = (p2[2] - p0[2]);
1982 double p2_3_p0_3 = (p2[3] - p0[3]);
1983 double p2_4_p0_4 = (p2[4] - p0[4]);
1984 double p2_5_p0_5 = (p2[5] - p0[5]);
1985 double p2_6_p0_6 = (p2[6] - p0[6]);
1986 double p2_7_p0_7 = (p2[7] - p0[7]);
1988 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 double p3_0_p0_0 = (p3[0] - p0[0]);
1990 double p3_1_p0_1 = (p3[1] - p0[1]);
1991 double p3_2_p0_2 = (p3[2] - p0[2]);
1992 double p3_3_p0_3 = (p3[3] - p0[3]);
1993 double p3_4_p0_4 = (p3[4] - p0[4]);
1994 double p3_5_p0_5 = (p3[5] - p0[5]);
1995 double p3_6_p0_6 = (p3[6] - p0[6]);
1996 double p3_7_p0_7 = (p3[7] - p0[7]);
1998 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 double p4_0_p0_0 = (p4[0] - p0[0]);
2000 double p4_1_p0_1 = (p4[1] - p0[1]);
2001 double p4_2_p0_2 = (p4[2] - p0[2]);
2002 double p4_3_p0_3 = (p4[3] - p0[3]);
2003 double p4_4_p0_4 = (p4[4] - p0[4]);
2004 double p4_5_p0_5 = (p4[5] - p0[5]);
2005 double p4_6_p0_6 = (p4[6] - p0[6]);
2006 double p4_7_p0_7 = (p4[7] - p0[7]);
2008 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 double q0_0_p0_0 = (q0[0] - p0[0]);
2010 double q0_1_p0_1 = (q0[1] - p0[1]);
2011 double q0_2_p0_2 = (q0[2] - p0[2]);
2012 double q0_3_p0_3 = (q0[3] - p0[3]);
2013 double q0_4_p0_4 = (q0[4] - p0[4]);
2014 double q0_5_p0_5 = (q0[5] - p0[5]);
2015 double q0_6_p0_6 = (q0[6] - p0[6]);
2016 double q0_7_p0_7 = (q0[7] - p0[7]);
2018 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 double q1_0_p0_0 = (q1[0] - p0[0]);
2020 double q1_1_p0_1 = (q1[1] - p0[1]);
2021 double q1_2_p0_2 = (q1[2] - p0[2]);
2022 double q1_3_p0_3 = (q1[3] - p0[3]);
2023 double q1_4_p0_4 = (q1[4] - p0[4]);
2024 double q1_5_p0_5 = (q1[5] - p0[5]);
2025 double q1_6_p0_6 = (q1[6] - p0[6]);
2026 double q1_7_p0_7 = (q1[7] - p0[7]);
2028 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 double q2_0_p0_0 = (q2[0] - p0[0]);
2030 double q2_1_p0_1 = (q2[1] - p0[1]);
2031 double q2_2_p0_2 = (q2[2] - p0[2]);
2032 double q2_3_p0_3 = (q2[3] - p0[3]);
2033 double q2_4_p0_4 = (q2[4] - p0[4]);
2034 double q2_5_p0_5 = (q2[5] - p0[5]);
2035 double q2_6_p0_6 = (q2[6] - p0[6]);
2036 double q2_7_p0_7 = (q2[7] - p0[7]);
2038 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 double q3_0_p0_0 = (q3[0] - p0[0]);
2040 double q3_1_p0_1 = (q3[1] - p0[1]);
2041 double q3_2_p0_2 = (q3[2] - p0[2]);
2042 double q3_3_p0_3 = (q3[3] - p0[3]);
2043 double q3_4_p0_4 = (q3[4] - p0[4]);
2044 double q3_5_p0_5 = (q3[5] - p0[5]);
2045 double q3_6_p0_6 = (q3[6] - p0[6]);
2046 double q3_7_p0_7 = (q3[7] - p0[7]);
2048 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)));
2050 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)));
2052 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)));
2054 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)));
2056 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)));
2058 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)));
2060 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)));
2062 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)));
2064 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)));
2066 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)));
2068 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)));
2070 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)));
2072 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)));
2074 b00 = (((a11 * ((a22 * a33) - (a23 * a32))) - (a21 * ((a12 * a33) - (a13 * a32)))) + (a31 * ((a12 * a23) - (a13 * a22))));
2076 b01 = -((((a22 * a33) - (a23 * a32)) + ((a23 * a31) - (a21 * a33))) + ((a21 * a32) - (a22 * a31)));
2078 b02 = ((((a12 * a33) - (a13 * a32)) + ((a13 * a31) - (a11 * a33))) + ((a11 * a32) - (a12 * a31)));
2080 b03 = -((((a12 * a23) - (a13 * a22)) + ((a13 * a21) - (a11 * a23))) + ((a11 * a22) - (a12 * a21)));
2082 b10 = -(((a10 * ((a22 * a33) - (a23 * a32))) - (a20 * ((a12 * a33) - (a13 * a32)))) + (a30 * ((a12 * a23) - (a13 * a22))));
2084 b11 = ((((a22 * a33) - (a23 * a32)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a32) - (a22 * a30)));
2086 b12 = -((((a12 * a33) - (a13 * a32)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a32) - (a12 * a30)));
2088 b13 = ((((a12 * a23) - (a13 * a22)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a22) - (a12 * a20)));
2090 b20 = (((a10 * ((a21 * a33) - (a23 * a31))) - (a20 * ((a11 * a33) - (a13 * a31)))) + (a30 * ((a11 * a23) - (a13 * a21))));
2092 b21 = -((((a21 * a33) - (a23 * a31)) + ((a23 * a30) - (a20 * a33))) + ((a20 * a31) - (a21 * a30)));
2094 b22 = ((((a11 * a33) - (a13 * a31)) + ((a13 * a30) - (a10 * a33))) + ((a10 * a31) - (a11 * a30)));
2096 b23 = -((((a11 * a23) - (a13 * a21)) + ((a13 * a20) - (a10 * a23))) + ((a10 * a21) - (a11 * a20)));
2098 b30 = -(((a10 * ((a21 * a32) - (a22 * a31))) - (a20 * ((a11 * a32) - (a12 * a31)))) + (a30 * ((a11 * a22) - (a12 * a21))));
2100 b31 = ((((a21 * a32) - (a22 * a31)) + ((a22 * a30) - (a20 * a32))) + ((a20 * a31) - (a21 * a30)));
2102 b32 = -((((a11 * a32) - (a12 * a31)) + ((a12 * a30) - (a10 * a32))) + ((a10 * a31) - (a11 * a30)));
2104 b33 = ((((a11 * a22) - (a12 * a21)) + ((a12 * a20) - (a10 * a22))) + ((a10 * a21) - (a11 * a20)));
2106 Delta = (((b00 + b10) + b20) + b30);
2107 double DeltaLambda0;
2108 DeltaLambda0 = ((((b01 * l1) + (b02 * l2)) + (b03 * l3)) + b00);
2109 double DeltaLambda1;
2110 DeltaLambda1 = ((((b11 * l1) + (b12 * l2)) + (b13 * l3)) + b10);
2111 double DeltaLambda2;
2112 DeltaLambda2 = ((((b21 * l1) + (b22 * l2)) + (b23 * l3)) + b20);
2113 double DeltaLambda3;
2114 DeltaLambda3 = ((((b31 * l1) + (b32 * l2)) + (b33 * l3)) + b30);
2116 r = ((Delta * l4) - ((((a40 * DeltaLambda0) + (a41 * DeltaLambda1)) + (a42 * DeltaLambda2)) + (a43 * DeltaLambda3)));
2118 double max1 = fabs(p2_5_p0_5);
2119 if( (max1 < fabs(p2_3_p0_3)) )
2121 max1 = fabs(p2_3_p0_3);
2123 if( (max1 < fabs(p2_0_p0_0)) )
2125 max1 = fabs(p2_0_p0_0);
2127 if( (max1 < fabs(p2_1_p0_1)) )
2129 max1 = fabs(p2_1_p0_1);
2131 if( (max1 < fabs(p2_6_p0_6)) )
2133 max1 = fabs(p2_6_p0_6);
2135 if( (max1 < fabs(p2_2_p0_2)) )
2137 max1 = fabs(p2_2_p0_2);
2139 if( (max1 < fabs(p2_4_p0_4)) )
2141 max1 = fabs(p2_4_p0_4);
2143 if( (max1 < fabs(p2_7_p0_7)) )
2145 max1 = fabs(p2_7_p0_7);
2147 double max2 = fabs(p1_4_p0_4);
2148 if( (max2 < fabs(p1_3_p0_3)) )
2150 max2 = fabs(p1_3_p0_3);
2152 if( (max2 < fabs(p1_7_p0_7)) )
2154 max2 = fabs(p1_7_p0_7);
2156 if( (max2 < fabs(p1_0_p0_0)) )
2158 max2 = fabs(p1_0_p0_0);
2160 if( (max2 < fabs(p1_2_p0_2)) )
2162 max2 = fabs(p1_2_p0_2);
2164 if( (max2 < fabs(p1_5_p0_5)) )
2166 max2 = fabs(p1_5_p0_5);
2168 if( (max2 < fabs(p1_1_p0_1)) )
2170 max2 = fabs(p1_1_p0_1);
2172 if( (max2 < fabs(p1_6_p0_6)) )
2174 max2 = fabs(p1_6_p0_6);
2176 double max3 = fabs(p3_3_p0_3);
2177 if( (max3 < fabs(p3_0_p0_0)) )
2179 max3 = fabs(p3_0_p0_0);
2181 if( (max3 < fabs(p3_1_p0_1)) )
2183 max3 = fabs(p3_1_p0_1);
2185 if( (max3 < fabs(p3_2_p0_2)) )
2187 max3 = fabs(p3_2_p0_2);
2189 if( (max3 < fabs(p3_4_p0_4)) )
2191 max3 = fabs(p3_4_p0_4);
2193 if( (max3 < fabs(p3_5_p0_5)) )
2195 max3 = fabs(p3_5_p0_5);
2197 if( (max3 < fabs(p3_6_p0_6)) )
2199 max3 = fabs(p3_6_p0_6);
2201 if( (max3 < fabs(p3_7_p0_7)) )
2203 max3 = fabs(p3_7_p0_7);
2205 double max4 = fabs(q0_0_p0_0);
2206 if( (max4 < fabs(q0_1_p0_1)) )
2208 max4 = fabs(q0_1_p0_1);
2210 if( (max4 < fabs(q0_2_p0_2)) )
2212 max4 = fabs(q0_2_p0_2);
2214 if( (max4 < fabs(q0_3_p0_3)) )
2216 max4 = fabs(q0_3_p0_3);
2218 if( (max4 < fabs(q0_4_p0_4)) )
2220 max4 = fabs(q0_4_p0_4);
2222 if( (max4 < fabs(q0_5_p0_5)) )
2224 max4 = fabs(q0_5_p0_5);
2226 if( (max4 < fabs(q0_6_p0_6)) )
2228 max4 = fabs(q0_6_p0_6);
2230 if( (max4 < fabs(q0_7_p0_7)) )
2232 max4 = fabs(q0_7_p0_7);
2234 if( (max4 < fabs(q1_0_p0_0)) )
2236 max4 = fabs(q1_0_p0_0);
2238 if( (max4 < fabs(q1_1_p0_1)) )
2240 max4 = fabs(q1_1_p0_1);
2242 if( (max4 < fabs(q1_2_p0_2)) )
2244 max4 = fabs(q1_2_p0_2);
2246 if( (max4 < fabs(q1_3_p0_3)) )
2248 max4 = fabs(q1_3_p0_3);
2250 if( (max4 < fabs(q1_4_p0_4)) )
2252 max4 = fabs(q1_4_p0_4);
2254 if( (max4 < fabs(q1_5_p0_5)) )
2256 max4 = fabs(q1_5_p0_5);
2258 if( (max4 < fabs(q1_6_p0_6)) )
2260 max4 = fabs(q1_6_p0_6);
2262 if( (max4 < fabs(q1_7_p0_7)) )
2264 max4 = fabs(q1_7_p0_7);
2266 double max5 = fabs(q1_0_p0_0);
2267 if( (max5 < fabs(q1_1_p0_1)) )
2269 max5 = fabs(q1_1_p0_1);
2271 if( (max5 < fabs(q1_2_p0_2)) )
2273 max5 = fabs(q1_2_p0_2);
2275 if( (max5 < fabs(q1_3_p0_3)) )
2277 max5 = fabs(q1_3_p0_3);
2279 if( (max5 < fabs(q1_4_p0_4)) )
2281 max5 = fabs(q1_4_p0_4);
2283 if( (max5 < fabs(q1_5_p0_5)) )
2285 max5 = fabs(q1_5_p0_5);
2287 if( (max5 < fabs(q1_6_p0_6)) )
2289 max5 = fabs(q1_6_p0_6);
2291 if( (max5 < fabs(q1_7_p0_7)) )
2293 max5 = fabs(q1_7_p0_7);
2295 if( (max5 < fabs(q2_0_p0_0)) )
2297 max5 = fabs(q2_0_p0_0);
2299 if( (max5 < fabs(q2_1_p0_1)) )
2301 max5 = fabs(q2_1_p0_1);
2303 if( (max5 < fabs(q2_2_p0_2)) )
2305 max5 = fabs(q2_2_p0_2);
2307 if( (max5 < fabs(q2_3_p0_3)) )
2309 max5 = fabs(q2_3_p0_3);
2311 if( (max5 < fabs(q2_4_p0_4)) )
2313 max5 = fabs(q2_4_p0_4);
2315 if( (max5 < fabs(q2_5_p0_5)) )
2317 max5 = fabs(q2_5_p0_5);
2319 if( (max5 < fabs(q2_6_p0_6)) )
2321 max5 = fabs(q2_6_p0_6);
2323 if( (max5 < fabs(q2_7_p0_7)) )
2325 max5 = fabs(q2_7_p0_7);
2327 double max6 = fabs(q2_0_p0_0);
2328 if( (max6 < fabs(q2_1_p0_1)) )
2330 max6 = fabs(q2_1_p0_1);
2332 if( (max6 < fabs(q2_2_p0_2)) )
2334 max6 = fabs(q2_2_p0_2);
2336 if( (max6 < fabs(q2_3_p0_3)) )
2338 max6 = fabs(q2_3_p0_3);
2340 if( (max6 < fabs(q2_4_p0_4)) )
2342 max6 = fabs(q2_4_p0_4);
2344 if( (max6 < fabs(q2_5_p0_5)) )
2346 max6 = fabs(q2_5_p0_5);
2348 if( (max6 < fabs(q2_6_p0_6)) )
2350 max6 = fabs(q2_6_p0_6);
2352 if( (max6 < fabs(q2_7_p0_7)) )
2354 max6 = fabs(q2_7_p0_7);
2356 if( (max6 < fabs(q3_0_p0_0)) )
2358 max6 = fabs(q3_0_p0_0);
2360 if( (max6 < fabs(q3_1_p0_1)) )
2362 max6 = fabs(q3_1_p0_1);
2364 if( (max6 < fabs(q3_2_p0_2)) )
2366 max6 = fabs(q3_2_p0_2);
2368 if( (max6 < fabs(q3_3_p0_3)) )
2370 max6 = fabs(q3_3_p0_3);
2372 if( (max6 < fabs(q3_4_p0_4)) )
2374 max6 = fabs(q3_4_p0_4);
2376 if( (max6 < fabs(q3_5_p0_5)) )
2378 max6 = fabs(q3_5_p0_5);
2380 if( (max6 < fabs(q3_6_p0_6)) )
2382 max6 = fabs(q3_6_p0_6);
2384 if( (max6 < fabs(q3_7_p0_7)) )
2386 max6 = fabs(q3_7_p0_7);
2388 double lower_bound_1;
2389 double upper_bound_1;
2392 lower_bound_1 = max4;
2393 upper_bound_1 = max4;
2394 if( (max1 < lower_bound_1) )
2396 lower_bound_1 = max1;
2400 if( (max1 > upper_bound_1) )
2402 upper_bound_1 = max1;
2405 if( (max3 < lower_bound_1) )
2407 lower_bound_1 = max3;
2411 if( (max3 > upper_bound_1) )
2413 upper_bound_1 = max3;
2416 if( (max2 < lower_bound_1) )
2418 lower_bound_1 = max2;
2422 if( (max2 > upper_bound_1) )
2424 upper_bound_1 = max2;
2427 if( (max6 < lower_bound_1) )
2429 lower_bound_1 = max6;
2433 if( (max6 > upper_bound_1) )
2435 upper_bound_1 = max6;
2438 if( (max5 < lower_bound_1) )
2440 lower_bound_1 = max5;
2444 if( (max5 > upper_bound_1) )
2446 upper_bound_1 = max5;
2449 if( (lower_bound_1 < 2.82528483194754087282e-50) )
2451 return FPG_UNCERTAIN_VALUE;
2455 if( (upper_bound_1 > 4.83570327845851562508e+24) )
2457 return FPG_UNCERTAIN_VALUE;
2459 eps = (4.37492894694731040560e-11 * (((((max2 * max4) * max1) * max5) * max3) * max6));
2466 if( (Delta < -eps) )
2468 int_tmp_result = -1;
2472 return FPG_UNCERTAIN_VALUE;
2476 Delta_sign = int_tmp_result;
2477 int int_tmp_result_FFWKCAA;
2484 if( (max8 < fabs(p4_3_p0_3)) )
2486 max8 = fabs(p4_3_p0_3);
2488 if( (max8 < fabs(p4_1_p0_1)) )
2490 max8 = fabs(p4_1_p0_1);
2492 if( (max8 < fabs(p4_4_p0_4)) )
2494 max8 = fabs(p4_4_p0_4);
2496 if( (max8 < fabs(p4_0_p0_0)) )
2498 max8 = fabs(p4_0_p0_0);
2500 if( (max8 < fabs(p4_2_p0_2)) )
2502 max8 = fabs(p4_2_p0_2);
2504 if( (max8 < fabs(p4_5_p0_5)) )
2506 max8 = fabs(p4_5_p0_5);
2508 if( (max8 < fabs(p4_6_p0_6)) )
2510 max8 = fabs(p4_6_p0_6);
2512 if( (max8 < fabs(p4_7_p0_7)) )
2514 max8 = fabs(p4_7_p0_7);
2537 double max10 = max4;
2538 double max11 = max4;
2539 if( (max11 < max5) )
2543 if( (max10 < max11) )
2547 if( (max10 < max2) )
2551 if( (max10 < max6) )
2555 if( (max10 < max5) )
2559 lower_bound_1 = max11;
2560 upper_bound_1 = max11;
2561 if( (max9 < lower_bound_1) )
2563 lower_bound_1 = max9;
2567 if( (max9 > upper_bound_1) )
2569 upper_bound_1 = max9;
2572 if( (max7 < lower_bound_1) )
2574 lower_bound_1 = max7;
2578 if( (max7 > upper_bound_1) )
2580 upper_bound_1 = max7;
2583 if( (max1 < lower_bound_1) )
2585 lower_bound_1 = max1;
2587 if( (max3 < lower_bound_1) )
2589 lower_bound_1 = max3;
2591 if( (max8 < lower_bound_1) )
2593 lower_bound_1 = max8;
2595 if( (max10 < lower_bound_1) )
2597 lower_bound_1 = max10;
2601 if( (max10 > upper_bound_1) )
2603 upper_bound_1 = max10;
2606 if( (lower_bound_1 < 4.17402518597284772324e-38) )
2608 return FPG_UNCERTAIN_VALUE;
2612 if( (upper_bound_1 > 4.83570327845851562508e+24) )
2614 return FPG_UNCERTAIN_VALUE;
2616 eps = (2.41492645607254025015e-09 * (((((((max8 * max11) * max1) * max10) * max3) * max10) * max9) * max7));
2619 int_tmp_result_FFWKCAA = 1;
2625 int_tmp_result_FFWKCAA = -1;
2629 return FPG_UNCERTAIN_VALUE;
2633 return (Delta_sign * int_tmp_result_FFWKCAA);