17 int orient2d(
const mpq2 &
a,
const mpq2 &
b,
const mpq2 &
c)
19 mpq_class detleft = (
a[0] -
c[0]) * (
b[1] -
c[1]);
20 mpq_class detright = (
a[1] -
c[1]) * (
b[0] -
c[0]);
21 mpq_class det = detleft - detright;
25 int incircle(
const mpq2 &
a,
const mpq2 &
b,
const mpq2 &
c,
const mpq2 &d)
27 mpq_class adx =
a[0] -
d[0];
28 mpq_class bdx =
b[0] -
d[0];
29 mpq_class cdx =
c[0] -
d[0];
30 mpq_class ady =
a[1] -
d[1];
31 mpq_class bdy =
b[1] -
d[1];
32 mpq_class cdy =
c[1] -
d[1];
34 mpq_class bdxcdy = bdx * cdy;
35 mpq_class cdxbdy = cdx * bdy;
36 mpq_class alift = adx * adx + ady * ady;
38 mpq_class cdxady = cdx * ady;
39 mpq_class adxcdy = adx * cdy;
40 mpq_class blift = bdx * bdx + bdy * bdy;
42 mpq_class adxbdy = adx * bdy;
43 mpq_class bdxady = bdx * ady;
44 mpq_class clift = cdx * cdx + cdy * cdy;
46 mpq_class det = alift * (bdxcdy - cdxbdy) + blift * (cdxady - adxcdy) +
47 clift * (adxbdy - bdxady);
51 int orient3d(
const mpq3 &
a,
const mpq3 &
b,
const mpq3 &
c,
const mpq3 &d)
53 mpq_class adx =
a[0] -
d[0];
54 mpq_class bdx =
b[0] -
d[0];
55 mpq_class cdx =
c[0] -
d[0];
56 mpq_class ady =
a[1] -
d[1];
57 mpq_class bdy =
b[1] -
d[1];
58 mpq_class cdy =
c[1] -
d[1];
59 mpq_class adz =
a[2] -
d[2];
60 mpq_class bdz =
b[2] -
d[2];
61 mpq_class cdz =
c[2] -
d[2];
63 mpq_class bdxcdy = bdx * cdy;
64 mpq_class cdxbdy = cdx * bdy;
66 mpq_class cdxady = cdx * ady;
67 mpq_class adxcdy = adx * cdy;
69 mpq_class adxbdy = adx * bdy;
70 mpq_class bdxady = bdx * ady;
72 mpq_class det = adz * (bdxcdy - cdxbdy) + bdz * (cdxady - adxcdy) + cdz * (adxbdy - bdxady);
83 namespace robust_pred {
96 double orient2dfast(
const double *pa,
const double *pb,
const double *pc);
97 double orient2d(
const double *pa,
const double *pb,
const double *pc);
98 double orient3dfast(
const double *pa,
const double *pb,
const double *pc,
const double *pd);
99 double orient3d(
const double *pa,
const double *pb,
const double *pc,
const double *pd);
100 double incirclefast(
const double *pa,
const double *pb,
const double *pc,
const double *pd);
101 double incircle(
const double *pa,
const double *pb,
const double *pc,
const double *pd);
103 const double *pa,
const double *pb,
const double *pc,
const double *pd,
const double *pe);
105 const double *pa,
const double *pb,
const double *pc,
const double *pd,
const double *pe);
174 #define Absolute(a) ((a) >= 0.0 ? (a) : -(a))
191 #define Fast_Two_Sum_Tail(a, b, x, y) \
195 #define Fast_Two_Sum(a, b, x, y) \
196 x = (double)(a + b); \
197 Fast_Two_Sum_Tail(a, b, x, y)
199 #define Fast_Two_Diff_Tail(a, b, x, y) \
203 #define Fast_Two_Diff(a, b, x, y) \
204 x = (double)(a - b); \
205 Fast_Two_Diff_Tail(a, b, x, y)
207 #define Two_Sum_Tail(a, b, x, y) \
208 bvirt = (double)(x - a); \
210 bround = b - bvirt; \
211 around = a - avirt; \
214 #define Two_Sum(a, b, x, y) \
215 x = (double)(a + b); \
216 Two_Sum_Tail(a, b, x, y)
218 #define Two_Diff_Tail(a, b, x, y) \
219 bvirt = (double)(a - x); \
221 bround = bvirt - b; \
222 around = a - avirt; \
225 #define Two_Diff(a, b, x, y) \
226 x = (double)(a - b); \
227 Two_Diff_Tail(a, b, x, y)
229 #define Split(a, ahi, alo) \
230 c = (double)(splitter * a); \
231 abig = (double)(c - a); \
235 #define Two_Product_Tail(a, b, x, y) \
236 Split(a, ahi, alo); \
237 Split(b, bhi, blo); \
238 err1 = x - (ahi * bhi); \
239 err2 = err1 - (alo * bhi); \
240 err3 = err2 - (ahi * blo); \
241 y = (alo * blo) - err3
243 #define Two_Product(a, b, x, y) \
244 x = (double)(a * b); \
245 Two_Product_Tail(a, b, x, y)
247 #define Two_Product_Presplit(a, b, bhi, blo, x, y) \
248 x = (double)(a * b); \
249 Split(a, ahi, alo); \
250 err1 = x - (ahi * bhi); \
251 err2 = err1 - (alo * bhi); \
252 err3 = err2 - (ahi * blo); \
253 y = (alo * blo) - err3
255 #define Two_Product_2Presplit(a, ahi, alo, b, bhi, blo, x, y) \
256 x = (double)(a * b); \
257 err1 = x - (ahi * bhi); \
258 err2 = err1 - (alo * bhi); \
259 err3 = err2 - (ahi * blo); \
260 y = (alo * blo) - err3
262 #define Square_Tail(a, x, y) \
263 Split(a, ahi, alo); \
264 err1 = x - (ahi * ahi); \
265 err3 = err1 - ((ahi + ahi) * alo); \
266 y = (alo * alo) - err3
268 #define Square(a, x, y) \
269 x = (double)(a * a); \
272 #define Two_One_Sum(a1, a0, b, x2, x1, x0) \
273 Two_Sum(a0, b, _i, x0); \
274 Two_Sum(a1, _i, x2, x1)
276 #define Two_One_Diff(a1, a0, b, x2, x1, x0) \
277 Two_Diff(a0, b, _i, x0); \
278 Two_Sum(a1, _i, x2, x1)
280 #define Two_Two_Sum(a1, a0, b1, b0, x3, x2, x1, x0) \
281 Two_One_Sum(a1, a0, b0, _j, _0, x0); \
282 Two_One_Sum(_j, _0, b1, x3, x2, x1)
284 #define Two_Two_Diff(a1, a0, b1, b0, x3, x2, x1, x0) \
285 Two_One_Diff(a1, a0, b0, _j, _0, x0); \
286 Two_One_Diff(_j, _0, b1, x3, x2, x1)
288 #define Four_One_Sum(a3, a2, a1, a0, b, x4, x3, x2, x1, x0) \
289 Two_One_Sum(a1, a0, b, _j, x1, x0); \
290 Two_One_Sum(a3, a2, _j, x4, x3, x2)
292 #define Four_Two_Sum(a3, a2, a1, a0, b1, b0, x5, x4, x3, x2, x1, x0) \
293 Four_One_Sum(a3, a2, a1, a0, b0, _k, _2, _1, _0, x0); \
294 Four_One_Sum(_k, _2, _1, _0, b1, x5, x4, x3, x2, x1)
296 #define Four_Four_Sum(a3, a2, a1, a0, b4, b3, b1, b0, x7, x6, x5, x4, x3, x2, x1, x0) \
297 Four_Two_Sum(a3, a2, a1, a0, b1, b0, _l, _2, _1, _0, x1, x0); \
298 Four_Two_Sum(_l, _2, _1, _0, b4, b3, x7, x6, x5, x4, x3, x2)
300 #define Eight_One_Sum(a7, a6, a5, a4, a3, a2, a1, a0, b, x8, x7, x6, x5, x4, x3, x2, x1, x0) \
301 Four_One_Sum(a3, a2, a1, a0, b, _j, x3, x2, x1, x0); \
302 Four_One_Sum(a7, a6, a5, a4, _j, x8, x7, x6, x5, x4)
304 #define Eight_Two_Sum( \
305 a7, a6, a5, a4, a3, a2, a1, a0, b1, b0, x9, x8, x7, x6, x5, x4, x3, x2, x1, x0) \
306 Eight_One_Sum(a7, a6, a5, a4, a3, a2, a1, a0, b0, _k, _6, _5, _4, _3, _2, _1, _0, x0); \
307 Eight_One_Sum(_k, _6, _5, _4, _3, _2, _1, _0, b1, x9, x8, x7, x6, x5, x4, x3, x2, x1)
309 #define Eight_Four_Sum(a7, \
333 Eight_Two_Sum(a7, a6, a5, a4, a3, a2, a1, a0, b1, b0, _l, _6, _5, _4, _3, _2, _1, _0, x1, x0); \
334 Eight_Two_Sum(_l, _6, _5, _4, _3, _2, _1, _0, b4, b3, x11, x10, x9, x8, x7, x6, x5, x4, x3, x2)
336 #define Two_One_Product(a1, a0, b, x3, x2, x1, x0) \
337 Split(b, bhi, blo); \
338 Two_Product_Presplit(a0, b, bhi, blo, _i, x0); \
339 Two_Product_Presplit(a1, b, bhi, blo, _j, _0); \
340 Two_Sum(_i, _0, _k, x1); \
341 Fast_Two_Sum(_j, _k, x3, x2)
343 #define Four_One_Product(a3, a2, a1, a0, b, x7, x6, x5, x4, x3, x2, x1, x0) \
344 Split(b, bhi, blo); \
345 Two_Product_Presplit(a0, b, bhi, blo, _i, x0); \
346 Two_Product_Presplit(a1, b, bhi, blo, _j, _0); \
347 Two_Sum(_i, _0, _k, x1); \
348 Fast_Two_Sum(_j, _k, _i, x2); \
349 Two_Product_Presplit(a2, b, bhi, blo, _j, _0); \
350 Two_Sum(_i, _0, _k, x3); \
351 Fast_Two_Sum(_j, _k, _i, x4); \
352 Two_Product_Presplit(a3, b, bhi, blo, _j, _0); \
353 Two_Sum(_i, _0, _k, x5); \
354 Fast_Two_Sum(_j, _k, x7, x6)
356 #define Two_Two_Product(a1, a0, b1, b0, x7, x6, x5, x4, x3, x2, x1, x0) \
357 Split(a0, a0hi, a0lo); \
358 Split(b0, bhi, blo); \
359 Two_Product_2Presplit(a0, a0hi, a0lo, b0, bhi, blo, _i, x0); \
360 Split(a1, a1hi, a1lo); \
361 Two_Product_2Presplit(a1, a1hi, a1lo, b0, bhi, blo, _j, _0); \
362 Two_Sum(_i, _0, _k, _1); \
363 Fast_Two_Sum(_j, _k, _l, _2); \
364 Split(b1, bhi, blo); \
365 Two_Product_2Presplit(a0, a0hi, a0lo, b1, bhi, blo, _i, _0); \
366 Two_Sum(_1, _0, _k, x1); \
367 Two_Sum(_2, _k, _j, _1); \
368 Two_Sum(_l, _j, _m, _2); \
369 Two_Product_2Presplit(a1, a1hi, a1lo, b1, bhi, blo, _j, _0); \
370 Two_Sum(_i, _0, _n, _0); \
371 Two_Sum(_1, _0, _i, x2); \
372 Two_Sum(_2, _i, _k, _1); \
373 Two_Sum(_m, _k, _l, _2); \
374 Two_Sum(_j, _n, _k, _0); \
375 Two_Sum(_1, _0, _j, x3); \
376 Two_Sum(_2, _j, _i, _1); \
377 Two_Sum(_l, _i, _m, _2); \
378 Two_Sum(_1, _k, _i, x4); \
379 Two_Sum(_2, _i, _k, x5); \
380 Two_Sum(_m, _k, x7, x6)
382 #define Two_Square(a1, a0, x5, x4, x3, x2, x1, x0) \
383 Square(a0, _j, x0); \
385 Two_Product(a1, _0, _k, _1); \
386 Two_One_Sum(_k, _1, _j, _l, _2, x1); \
387 Square(a1, _j, _1); \
388 Two_Two_Sum(_j, _1, _l, _2, x5, x4, x3, x2)
419 double check, lastcheck;
437 every_other = !every_other;
439 }
while (!
ELEM(check, 1.0, lastcheck));
466 int elen,
const double *
e,
int flen,
const double *f,
double *h)
472 double avirt, bround, around;
473 int eindex, findex, hindex;
479 if ((fnow > enow) == (fnow > -enow)) {
488 if ((eindex < elen) && (findex < flen)) {
489 if ((fnow > enow) == (fnow > -enow)) {
501 while ((eindex < elen) && (findex < flen)) {
502 if ((fnow > enow) == (fnow > -enow)) {
516 while (eindex < elen) {
524 while (findex < flen) {
532 if ((Q != 0.0) || (hindex == 0)) {
554 double avirt, bround, around;
557 double ahi, alo, bhi, blo;
558 double err1, err2, err3;
566 for (eindex = 1; eindex < elen; eindex++) {
578 if ((Q != 0.0) || (hindex == 0)) {
591 for (eindex = 1; eindex < elen; eindex++) {
615 double orient2dfast(
const double *pa,
const double *pb,
const double *pc)
617 double acx, bcx, acy, bcy;
623 return acx * bcy - acy * bcx;
626 static double orient2dadapt(
const double *pa,
const double *pb,
const double *pc,
double detsum)
628 INEXACT double acx, acy, bcx, bcy;
629 double acxtail, acytail, bcxtail, bcytail;
630 INEXACT double detleft, detright;
631 double detlefttail, detrighttail;
632 double det, errbound;
633 double B[4], C1[8], C2[12],
D[16];
635 int C1length, C2length, Dlength;
642 double avirt, bround, around;
645 double ahi, alo, bhi, blo;
646 double err1, err2, err3;
650 acx = (
double)(pa[0] - pc[0]);
651 bcx = (
double)(pb[0] - pc[0]);
652 acy = (
double)(pa[1] - pc[1]);
653 bcy = (
double)(pb[1] - pc[1]);
658 Two_Two_Diff(detleft, detlefttail, detright, detrighttail,
B3,
B[2],
B[1],
B[0]);
663 if ((det >= errbound) || (-det >= errbound)) {
672 if ((acxtail == 0.0) && (acytail == 0.0) && (bcxtail == 0.0) && (bcytail == 0.0)) {
677 det += (acx * bcytail + bcy * acxtail) - (acy * bcxtail + bcx * acytail);
678 if ((det >= errbound) || (-det >= errbound)) {
700 return (
D[Dlength - 1]);
703 double orient2d(
const double *pa,
const double *pb,
const double *pc)
705 double detleft, detright, det;
706 double detsum, errbound;
708 detleft = (pa[0] - pc[0]) * (pb[1] - pc[1]);
709 detright = (pa[1] - pc[1]) * (pb[0] - pc[0]);
710 det = detleft - detright;
713 if (detright <= 0.0) {
716 detsum = detleft + detright;
718 else if (detleft < 0.0) {
719 if (detright >= 0.0) {
722 detsum = -detleft - detright;
729 if ((det >= errbound) || (-det >= errbound)) {
758 double orient3dfast(
const double *pa,
const double *pb,
const double *pc,
const double *pd)
760 double adx, bdx, cdx;
761 double ady, bdy, cdy;
762 double adz, bdz, cdz;
774 return adx * (bdy * cdz - bdz * cdy) + bdx * (cdy * adz - cdz * ady) +
775 cdx * (ady * bdz - adz * bdy);
784 const double *pa,
const double *pb,
const double *pc,
const double *pd,
double permanent)
786 INEXACT double adx, bdx, cdx, ady, bdy, cdy, adz, bdz, cdz;
787 double det, errbound;
789 INEXACT double bdxcdy1, cdxbdy1, cdxady1, adxcdy1, adxbdy1, bdxady1;
790 double bdxcdy0, cdxbdy0, cdxady0, adxcdy0, adxbdy0, bdxady0;
791 double bc[4], ca[4], ab[4];
793 double adet[8], bdet[8], cdet[8];
794 int alen, blen, clen;
797 double *finnow, *finother, *finswap;
798 double fin1[192], fin2[192];
801 double adxtail, bdxtail, cdxtail;
802 double adytail, bdytail, cdytail;
803 double adztail, bdztail, cdztail;
804 INEXACT double at_blarge, at_clarge;
805 INEXACT double bt_clarge, bt_alarge;
806 INEXACT double ct_alarge, ct_blarge;
807 double at_b[4], at_c[4], bt_c[4], bt_a[4], ct_a[4], ct_b[4];
808 int at_blen, at_clen, bt_clen, bt_alen, ct_alen, ct_blen;
809 INEXACT double bdxt_cdy1, cdxt_bdy1, cdxt_ady1;
810 INEXACT double adxt_cdy1, adxt_bdy1, bdxt_ady1;
811 double bdxt_cdy0, cdxt_bdy0, cdxt_ady0;
812 double adxt_cdy0, adxt_bdy0, bdxt_ady0;
813 INEXACT double bdyt_cdx1, cdyt_bdx1, cdyt_adx1;
814 INEXACT double adyt_cdx1, adyt_bdx1, bdyt_adx1;
815 double bdyt_cdx0, cdyt_bdx0, cdyt_adx0;
816 double adyt_cdx0, adyt_bdx0, bdyt_adx0;
817 double bct[8], cat[8], abt[8];
818 int bctlen, catlen, abtlen;
819 INEXACT double bdxt_cdyt1, cdxt_bdyt1, cdxt_adyt1;
820 INEXACT double adxt_cdyt1, adxt_bdyt1, bdxt_adyt1;
821 double bdxt_cdyt0, cdxt_bdyt0, cdxt_adyt0;
822 double adxt_cdyt0, adxt_bdyt0, bdxt_adyt0;
823 double u[4],
v[12],
w[16];
825 int vlength, wlength;
829 double avirt, bround, around;
832 double ahi, alo, bhi, blo;
833 double err1, err2, err3;
837 adx = (
double)(pa[0] - pd[0]);
838 bdx = (
double)(pb[0] - pd[0]);
839 cdx = (
double)(pc[0] - pd[0]);
840 ady = (
double)(pa[1] - pd[1]);
841 bdy = (
double)(pb[1] - pd[1]);
842 cdy = (
double)(pc[1] - pd[1]);
843 adz = (
double)(pa[2] - pd[2]);
844 bdz = (
double)(pb[2] - pd[2]);
845 cdz = (
double)(pc[2] - pd[2]);
849 Two_Two_Diff(bdxcdy1, bdxcdy0, cdxbdy1, cdxbdy0, bc3, bc[2], bc[1], bc[0]);
855 Two_Two_Diff(cdxady1, cdxady0, adxcdy1, adxcdy0, ca3, ca[2], ca[1], ca[0]);
861 Two_Two_Diff(adxbdy1, adxbdy0, bdxady1, bdxady0, ab3, ab[2], ab[1], ab[0]);
870 if ((det >= errbound) || (-det >= errbound)) {
884 if ((adxtail == 0.0) && (bdxtail == 0.0) && (cdxtail == 0.0) && (adytail == 0.0) &&
885 (bdytail == 0.0) && (cdytail == 0.0) && (adztail == 0.0) && (bdztail == 0.0) &&
891 det += (adz * ((bdx * cdytail + cdy * bdxtail) - (bdy * cdxtail + cdx * bdytail)) +
892 adztail * (bdx * cdy - bdy * cdx)) +
893 (bdz * ((cdx * adytail + ady * cdxtail) - (cdy * adxtail + adx * cdytail)) +
894 bdztail * (cdx * ady - cdy * adx)) +
895 (cdz * ((adx * bdytail + bdy * adxtail) - (ady * bdxtail + bdx * adytail)) +
896 cdztail * (adx * bdy - ady * bdx));
897 if ((det >= errbound) || (-det >= errbound)) {
904 if (adxtail == 0.0) {
905 if (adytail == 0.0) {
922 if (adytail == 0.0) {
935 adxt_bdy1, adxt_bdy0, adyt_bdx1, adyt_bdx0, at_blarge, at_b[2], at_b[1], at_b[0]);
941 adyt_cdx1, adyt_cdx0, adxt_cdy1, adxt_cdy0, at_clarge, at_c[2], at_c[1], at_c[0]);
946 if (bdxtail == 0.0) {
947 if (bdytail == 0.0) {
964 if (bdytail == 0.0) {
977 bdxt_cdy1, bdxt_cdy0, bdyt_cdx1, bdyt_cdx0, bt_clarge, bt_c[2], bt_c[1], bt_c[0]);
983 bdyt_adx1, bdyt_adx0, bdxt_ady1, bdxt_ady0, bt_alarge, bt_a[2], bt_a[1], bt_a[0]);
988 if (cdxtail == 0.0) {
989 if (cdytail == 0.0) {
1001 ct_b[1] = ct_blarge;
1006 if (cdytail == 0.0) {
1008 ct_a[1] = ct_alarge;
1012 ct_b[1] = ct_blarge;
1019 cdxt_ady1, cdxt_ady0, cdyt_adx1, cdyt_adx0, ct_alarge, ct_a[2], ct_a[1], ct_a[0]);
1020 ct_a[3] = ct_alarge;
1025 cdyt_bdx1, cdyt_bdx0, cdxt_bdy1, cdxt_bdy0, ct_blarge, ct_b[2], ct_b[1], ct_b[0]);
1026 ct_b[3] = ct_blarge;
1052 if (adztail != 0.0) {
1059 if (bdztail != 0.0) {
1066 if (cdztail != 0.0) {
1074 if (adxtail != 0.0) {
1075 if (bdytail != 0.0) {
1076 Two_Product(adxtail, bdytail, adxt_bdyt1, adxt_bdyt0);
1083 if (cdztail != 0.0) {
1084 Two_One_Product(adxt_bdyt1, adxt_bdyt0, cdztail, u3, u[2], u[1], u[0]);
1092 if (cdytail != 0.0) {
1094 Two_Product(negate, cdytail, adxt_cdyt1, adxt_cdyt0);
1101 if (bdztail != 0.0) {
1102 Two_One_Product(adxt_cdyt1, adxt_cdyt0, bdztail, u3, u[2], u[1], u[0]);
1111 if (bdxtail != 0.0) {
1112 if (cdytail != 0.0) {
1113 Two_Product(bdxtail, cdytail, bdxt_cdyt1, bdxt_cdyt0);
1120 if (adztail != 0.0) {
1121 Two_One_Product(bdxt_cdyt1, bdxt_cdyt0, adztail, u3, u[2], u[1], u[0]);
1129 if (adytail != 0.0) {
1131 Two_Product(negate, adytail, bdxt_adyt1, bdxt_adyt0);
1138 if (cdztail != 0.0) {
1139 Two_One_Product(bdxt_adyt1, bdxt_adyt0, cdztail, u3, u[2], u[1], u[0]);
1148 if (cdxtail != 0.0) {
1149 if (adytail != 0.0) {
1150 Two_Product(cdxtail, adytail, cdxt_adyt1, cdxt_adyt0);
1157 if (bdztail != 0.0) {
1158 Two_One_Product(cdxt_adyt1, cdxt_adyt0, bdztail, u3, u[2], u[1], u[0]);
1166 if (bdytail != 0.0) {
1168 Two_Product(negate, bdytail, cdxt_bdyt1, cdxt_bdyt0);
1175 if (adztail != 0.0) {
1176 Two_One_Product(cdxt_bdyt1, cdxt_bdyt0, adztail, u3, u[2], u[1], u[0]);
1186 if (adztail != 0.0) {
1193 if (bdztail != 0.0) {
1200 if (cdztail != 0.0) {
1208 return finnow[finlength - 1];
1211 double orient3d(
const double *pa,
const double *pb,
const double *pc,
const double *pd)
1213 double adx, bdx, cdx, ady, bdy, cdy, adz, bdz, cdz;
1214 double bdxcdy, cdxbdy, cdxady, adxcdy, adxbdy, bdxady;
1216 double permanent, errbound;
1218 adx = pa[0] - pd[0];
1219 bdx = pb[0] - pd[0];
1220 cdx = pc[0] - pd[0];
1221 ady = pa[1] - pd[1];
1222 bdy = pb[1] - pd[1];
1223 cdy = pc[1] - pd[1];
1224 adz = pa[2] - pd[2];
1225 bdz = pb[2] - pd[2];
1226 cdz = pc[2] - pd[2];
1237 det = adz * (bdxcdy - cdxbdy) + bdz * (cdxady - adxcdy) + cdz * (adxbdy - bdxady);
1243 if ((det > errbound) || (-det > errbound)) {
1269 double incirclefast(
const double *pa,
const double *pb,
const double *pc,
const double *pd)
1271 double adx, ady, bdx, bdy, cdx, cdy;
1272 double abdet, bcdet, cadet;
1273 double alift, blift, clift;
1275 adx = pa[0] - pd[0];
1276 ady = pa[1] - pd[1];
1277 bdx = pb[0] - pd[0];
1278 bdy = pb[1] - pd[1];
1279 cdx = pc[0] - pd[0];
1280 cdy = pc[1] - pd[1];
1282 abdet = adx * bdy - bdx * ady;
1283 bcdet = bdx * cdy - cdx * bdy;
1284 cadet = cdx * ady - adx * cdy;
1285 alift = adx * adx + ady * ady;
1286 blift = bdx * bdx + bdy * bdy;
1287 clift = cdx * cdx + cdy * cdy;
1289 return alift * bcdet + blift * cadet + clift * abdet;
1298 const double *pa,
const double *pb,
const double *pc,
const double *pd,
double permanent)
1300 INEXACT double adx, bdx, cdx, ady, bdy, cdy;
1301 double det, errbound;
1303 INEXACT double bdxcdy1, cdxbdy1, cdxady1, adxcdy1, adxbdy1, bdxady1;
1304 double bdxcdy0, cdxbdy0, cdxady0, adxcdy0, adxbdy0, bdxady0;
1305 double bc[4], ca[4], ab[4];
1307 double axbc[8], axxbc[16], aybc[8], ayybc[16], adet[32];
1308 int axbclen, axxbclen, aybclen, ayybclen, alen;
1309 double bxca[8], bxxca[16], byca[8], byyca[16], bdet[32];
1310 int bxcalen, bxxcalen, bycalen, byycalen, blen;
1311 double cxab[8], cxxab[16], cyab[8], cyyab[16], cdet[32];
1312 int cxablen, cxxablen, cyablen, cyyablen, clen;
1315 double fin1[1152], fin2[1152];
1316 double *finnow, *finother, *finswap;
1319 double adxtail, bdxtail, cdxtail, adytail, bdytail, cdytail;
1320 INEXACT double adxadx1, adyady1, bdxbdx1, bdybdy1, cdxcdx1, cdycdy1;
1321 double adxadx0, adyady0, bdxbdx0, bdybdy0, cdxcdx0, cdycdy0;
1322 double aa[4], bb[4], cc[4];
1328 double temp8[8], temp16a[16], temp16b[16], temp16c[16];
1329 double temp32a[32], temp32b[32], temp48[48], temp64[64];
1330 int temp8len, temp16alen, temp16blen, temp16clen;
1331 int temp32alen, temp32blen, temp48len, temp64len;
1332 double axtbb[8], axtcc[8], aytbb[8], aytcc[8];
1333 int axtbblen, axtcclen, aytbblen, aytcclen;
1334 double bxtaa[8], bxtcc[8], bytaa[8], bytcc[8];
1335 int bxtaalen, bxtcclen, bytaalen, bytcclen;
1336 double cxtaa[8], cxtbb[8], cytaa[8], cytbb[8];
1337 int cxtaalen, cxtbblen, cytaalen, cytbblen;
1338 double axtbc[8], aytbc[8], bxtca[8], bytca[8], cxtab[8], cytab[8];
1339 int axtbclen, aytbclen, bxtcalen, bytcalen, cxtablen, cytablen;
1340 double axtbct[16], aytbct[16], bxtcat[16], bytcat[16], cxtabt[16], cytabt[16];
1341 int axtbctlen, aytbctlen, bxtcatlen, bytcatlen, cxtabtlen, cytabtlen;
1342 double axtbctt[8], aytbctt[8], bxtcatt[8];
1343 double bytcatt[8], cxtabtt[8], cytabtt[8];
1344 int axtbcttlen, aytbcttlen, bxtcattlen, bytcattlen, cxtabttlen, cytabttlen;
1345 double abt[8], bct[8], cat[8];
1346 int abtlen, bctlen, catlen;
1347 double abtt[4], bctt[4], catt[4];
1348 int abttlen, bcttlen, cattlen;
1349 INEXACT double abtt3, bctt3, catt3;
1353 double avirt, bround, around;
1356 double ahi, alo, bhi, blo;
1357 double err1, err2, err3;
1361 adx = (
double)(pa[0] - pd[0]);
1362 bdx = (
double)(pb[0] - pd[0]);
1363 cdx = (
double)(pc[0] - pd[0]);
1364 ady = (
double)(pa[1] - pd[1]);
1365 bdy = (
double)(pb[1] - pd[1]);
1366 cdy = (
double)(pc[1] - pd[1]);
1370 Two_Two_Diff(bdxcdy1, bdxcdy0, cdxbdy1, cdxbdy0, bc3, bc[2], bc[1], bc[0]);
1380 Two_Two_Diff(cdxady1, cdxady0, adxcdy1, adxcdy0, ca3, ca[2], ca[1], ca[0]);
1390 Two_Two_Diff(adxbdy1, adxbdy0, bdxady1, bdxady0, ab3, ab[2], ab[1], ab[0]);
1403 if ((det >= errbound) || (-det >= errbound)) {
1413 if ((adxtail == 0.0) && (bdxtail == 0.0) && (cdxtail == 0.0) && (adytail == 0.0) &&
1414 (bdytail == 0.0) && (cdytail == 0.0)) {
1419 det += ((adx * adx + ady * ady) *
1420 ((bdx * cdytail + cdy * bdxtail) - (bdy * cdxtail + cdx * bdytail)) +
1421 2.0 * (adx * adxtail + ady * adytail) * (bdx * cdy - bdy * cdx)) +
1422 ((bdx * bdx + bdy * bdy) *
1423 ((cdx * adytail + ady * cdxtail) - (cdy * adxtail + adx * cdytail)) +
1424 2.0 * (bdx * bdxtail + bdy * bdytail) * (cdx * ady - cdy * adx)) +
1425 ((cdx * cdx + cdy * cdy) *
1426 ((adx * bdytail + bdy * adxtail) - (ady * bdxtail + bdx * adytail)) +
1427 2.0 * (cdx * cdxtail + cdy * cdytail) * (adx * bdy - ady * bdx));
1428 if ((det >= errbound) || (-det >= errbound)) {
1435 if ((bdxtail != 0.0) || (bdytail != 0.0) || (cdxtail != 0.0) || (cdytail != 0.0)) {
1436 Square(adx, adxadx1, adxadx0);
1437 Square(ady, adyady1, adyady0);
1438 Two_Two_Sum(adxadx1, adxadx0, adyady1, adyady0, aa3, aa[2], aa[1], aa[0]);
1441 if ((cdxtail != 0.0) || (cdytail != 0.0) || (adxtail != 0.0) || (adytail != 0.0)) {
1442 Square(bdx, bdxbdx1, bdxbdx0);
1443 Square(bdy, bdybdy1, bdybdy0);
1444 Two_Two_Sum(bdxbdx1, bdxbdx0, bdybdy1, bdybdy0, bb3, bb[2], bb[1], bb[0]);
1447 if ((adxtail != 0.0) || (adytail != 0.0) || (bdxtail != 0.0) || (bdytail != 0.0)) {
1448 Square(cdx, cdxcdx1, cdxcdx0);
1449 Square(cdy, cdycdy1, cdycdy0);
1450 Two_Two_Sum(cdxcdx1, cdxcdx0, cdycdy1, cdycdy0, cc3, cc[2], cc[1], cc[0]);
1454 if (adxtail != 0.0) {
1471 if (adytail != 0.0) {
1488 if (bdxtail != 0.0) {
1505 if (bdytail != 0.0) {
1522 if (cdxtail != 0.0) {
1539 if (cdytail != 0.0) {
1557 if ((adxtail != 0.0) || (adytail != 0.0)) {
1558 if ((bdxtail != 0.0) || (bdytail != 0.0) || (cdxtail != 0.0) || (cdytail != 0.0)) {
1561 Two_Two_Sum(ti1, ti0, tj1, tj0, u3, u[2], u[1], u[0]);
1573 Two_Two_Diff(ti1, ti0, tj1, tj0, bctt3, bctt[2], bctt[1], bctt[0]);
1584 if (adxtail != 0.0) {
1593 if (bdytail != 0.0) {
1601 if (cdytail != 0.0) {
1621 if (adytail != 0.0) {
1643 if ((bdxtail != 0.0) || (bdytail != 0.0)) {
1644 if ((cdxtail != 0.0) || (cdytail != 0.0) || (adxtail != 0.0) || (adytail != 0.0)) {
1647 Two_Two_Sum(ti1, ti0, tj1, tj0, u3, u[2], u[1], u[0]);
1659 Two_Two_Diff(ti1, ti0, tj1, tj0, catt3, catt[2], catt[1], catt[0]);
1670 if (bdxtail != 0.0) {
1679 if (cdytail != 0.0) {
1687 if (adytail != 0.0) {
1707 if (bdytail != 0.0) {
1729 if ((cdxtail != 0.0) || (cdytail != 0.0)) {
1730 if ((adxtail != 0.0) || (adytail != 0.0) || (bdxtail != 0.0) || (bdytail != 0.0)) {
1733 Two_Two_Sum(ti1, ti0, tj1, tj0, u3, u[2], u[1], u[0]);
1745 Two_Two_Diff(ti1, ti0, tj1, tj0, abtt3, abtt[2], abtt[1], abtt[0]);
1756 if (cdxtail != 0.0) {
1765 if (adytail != 0.0) {
1773 if (bdytail != 0.0) {
1793 if (cdytail != 0.0) {
1816 return finnow[finlength - 1];
1819 double incircle(
const double *pa,
const double *pb,
const double *pc,
const double *pd)
1821 double adx, bdx, cdx, ady, bdy, cdy;
1822 double bdxcdy, cdxbdy, cdxady, adxcdy, adxbdy, bdxady;
1823 double alift, blift, clift;
1825 double permanent, errbound;
1827 adx = pa[0] - pd[0];
1828 bdx = pb[0] - pd[0];
1829 cdx = pc[0] - pd[0];
1830 ady = pa[1] - pd[1];
1831 bdy = pb[1] - pd[1];
1832 cdy = pc[1] - pd[1];
1836 alift = adx * adx + ady * ady;
1840 blift = bdx * bdx + bdy * bdy;
1844 clift = cdx * cdx + cdy * cdy;
1846 det = alift * (bdxcdy - cdxbdy) + blift * (cdxady - adxcdy) + clift * (adxbdy - bdxady);
1852 if ((det > errbound) || (-det > errbound)) {
1880 const double *pa,
const double *pb,
const double *pc,
const double *pd,
const double *pe)
1882 double aex, bex, cex, dex;
1883 double aey, bey, cey, dey;
1884 double aez, bez, cez, dez;
1885 double alift, blift, clift, dlift;
1886 double ab, bc, cd, da, ac, bd;
1887 double abc, bcd, cda, dab;
1889 aex = pa[0] - pe[0];
1890 bex = pb[0] - pe[0];
1891 cex = pc[0] - pe[0];
1892 dex = pd[0] - pe[0];
1893 aey = pa[1] - pe[1];
1894 bey = pb[1] - pe[1];
1895 cey = pc[1] - pe[1];
1896 dey = pd[1] - pe[1];
1897 aez = pa[2] - pe[2];
1898 bez = pb[2] - pe[2];
1899 cez = pc[2] - pe[2];
1900 dez = pd[2] - pe[2];
1902 ab = aex * bey - bex * aey;
1903 bc = bex * cey - cex * bey;
1904 cd = cex * dey - dex * cey;
1905 da = dex * aey - aex * dey;
1907 ac = aex * cey - cex * aey;
1908 bd = bex * dey - dex * bey;
1910 abc = aez * bc - bez * ac + cez * ab;
1911 bcd = bez * cd - cez * bd + dez * bc;
1912 cda = cez * da + dez * ac + aez * cd;
1913 dab = dez * ab + aez * bd + bez * da;
1915 alift = aex * aex + aey * aey + aez * aez;
1916 blift = bex * bex + bey * bey + bez * bez;
1917 clift = cex * cex + cey * cey + cez * cez;
1918 dlift = dex * dex + dey * dey + dez * dez;
1920 return (dlift * abc - clift * dab) + (blift * cda - alift * bcd);
1924 const double *pa,
const double *pb,
const double *pc,
const double *pd,
const double *pe)
1926 INEXACT double axby1, bxcy1, cxdy1, dxey1, exay1;
1927 INEXACT double bxay1, cxby1, dxcy1, exdy1, axey1;
1928 INEXACT double axcy1, bxdy1, cxey1, dxay1, exby1;
1929 INEXACT double cxay1, dxby1, excy1, axdy1, bxey1;
1930 double axby0, bxcy0, cxdy0, dxey0, exay0;
1931 double bxay0, cxby0, dxcy0, exdy0, axey0;
1932 double axcy0, bxdy0, cxey0, dxay0, exby0;
1933 double cxay0, dxby0, excy0, axdy0, bxey0;
1934 double ab[4], bc[4], cd[4], de[4], ea[4];
1935 double ac[4], bd[4], ce[4], da[4], eb[4];
1936 double temp8a[8], temp8b[8], temp16[16];
1937 int temp8alen, temp8blen, temp16len;
1938 double abc[24], bcd[24], cde[24], dea[24], eab[24];
1939 double abd[24], bce[24], cda[24], deb[24], eac[24];
1940 int abclen, bcdlen, cdelen, dealen, eablen;
1941 int abdlen, bcelen, cdalen, deblen, eaclen;
1942 double temp48a[48], temp48b[48];
1943 int temp48alen, temp48blen;
1944 double abcd[96], bcde[96], cdea[96], deab[96], eabc[96];
1945 int abcdlen, bcdelen, cdealen, deablen, eabclen;
1946 double temp192[192];
1947 double det384x[384], det384y[384], det384z[384];
1948 int xlen, ylen, zlen;
1951 double adet[1152], bdet[1152], cdet[1152], ddet[1152], edet[1152];
1952 int alen, blen, clen, dlen, elen;
1953 double abdet[2304], cddet[2304], cdedet[3456];
1960 double avirt, bround, around;
1963 double ahi, alo, bhi, blo;
1964 double err1, err2, err3;
1970 Two_Two_Diff(axby1, axby0, bxay1, bxay0, ab[3], ab[2], ab[1], ab[0]);
1974 Two_Two_Diff(bxcy1, bxcy0, cxby1, cxby0, bc[3], bc[2], bc[1], bc[0]);
1978 Two_Two_Diff(cxdy1, cxdy0, dxcy1, dxcy0, cd[3], cd[2], cd[1], cd[0]);
1982 Two_Two_Diff(dxey1, dxey0, exdy1, exdy0, de[3], de[2], de[1], de[0]);
1986 Two_Two_Diff(exay1, exay0, axey1, axey0, ea[3], ea[2], ea[1], ea[0]);
1990 Two_Two_Diff(axcy1, axcy0, cxay1, cxay0, ac[3], ac[2], ac[1], ac[0]);
1994 Two_Two_Diff(bxdy1, bxdy0, dxby1, dxby0, bd[3], bd[2], bd[1], bd[0]);
1998 Two_Two_Diff(cxey1, cxey0, excy1, excy0, ce[3], ce[2], ce[1], ce[0]);
2002 Two_Two_Diff(dxay1, dxay0, axdy1, axdy0, da[3], da[2], da[1], da[0]);
2006 Two_Two_Diff(exby1, exby0, bxey1, bxey0, eb[3], eb[2], eb[1], eb[0]);
2070 for (i = 0; i < temp48blen; i++) {
2071 temp48b[i] = -temp48b[i];
2085 for (i = 0; i < temp48blen; i++) {
2086 temp48b[i] = -temp48b[i];
2100 for (i = 0; i < temp48blen; i++) {
2101 temp48b[i] = -temp48b[i];
2115 for (i = 0; i < temp48blen; i++) {
2116 temp48b[i] = -temp48b[i];
2130 for (i = 0; i < temp48blen; i++) {
2131 temp48b[i] = -temp48b[i];
2148 return deter[deterlen - 1];
2158 INEXACT double aex, bex, cex, dex, aey, bey, cey, dey, aez, bez, cez, dez;
2159 double det, errbound;
2161 INEXACT double aexbey1, bexaey1, bexcey1, cexbey1;
2162 INEXACT double cexdey1, dexcey1, dexaey1, aexdey1;
2163 INEXACT double aexcey1, cexaey1, bexdey1, dexbey1;
2164 double aexbey0, bexaey0, bexcey0, cexbey0;
2165 double cexdey0, dexcey0, dexaey0, aexdey0;
2166 double aexcey0, cexaey0, bexdey0, dexbey0;
2167 double ab[4], bc[4], cd[4], da[4], ac[4], bd[4];
2168 INEXACT double ab3, bc3, cd3, da3, ac3, bd3;
2169 double abeps, bceps, cdeps, daeps, aceps, bdeps;
2170 double temp8a[8], temp8b[8], temp8c[8], temp16[16], temp24[24], temp48[48];
2171 int temp8alen, temp8blen, temp8clen, temp16len, temp24len, temp48len;
2172 double xdet[96], ydet[96], zdet[96], xydet[192];
2173 int xlen, ylen, zlen, xylen;
2174 double adet[288], bdet[288], cdet[288], ddet[288];
2175 int alen, blen, clen, dlen;
2176 double abdet[576], cddet[576];
2181 double aextail, bextail, cextail, dextail;
2182 double aeytail, beytail, ceytail, deytail;
2183 double aeztail, beztail, ceztail, deztail;
2186 double avirt, bround, around;
2189 double ahi, alo, bhi, blo;
2190 double err1, err2, err3;
2194 aex = (
double)(pa[0] - pe[0]);
2195 bex = (
double)(pb[0] - pe[0]);
2196 cex = (
double)(pc[0] - pe[0]);
2197 dex = (
double)(pd[0] - pe[0]);
2198 aey = (
double)(pa[1] - pe[1]);
2199 bey = (
double)(pb[1] - pe[1]);
2200 cey = (
double)(pc[1] - pe[1]);
2201 dey = (
double)(pd[1] - pe[1]);
2202 aez = (
double)(pa[2] - pe[2]);
2203 bez = (
double)(pb[2] - pe[2]);
2204 cez = (
double)(pc[2] - pe[2]);
2205 dez = (
double)(pd[2] - pe[2]);
2209 Two_Two_Diff(aexbey1, aexbey0, bexaey1, bexaey0, ab3, ab[2], ab[1], ab[0]);
2214 Two_Two_Diff(bexcey1, bexcey0, cexbey1, cexbey0, bc3, bc[2], bc[1], bc[0]);
2219 Two_Two_Diff(cexdey1, cexdey0, dexcey1, dexcey0, cd3, cd[2], cd[1], cd[0]);
2224 Two_Two_Diff(dexaey1, dexaey0, aexdey1, aexdey0, da3, da[2], da[1], da[0]);
2229 Two_Two_Diff(aexcey1, aexcey0, cexaey1, cexaey0, ac3, ac[2], ac[1], ac[0]);
2234 Two_Two_Diff(bexdey1, bexdey0, dexbey1, dexbey0, bd3, bd[2], bd[1], bd[0]);
2299 if ((det >= errbound) || (-det >= errbound)) {
2315 if ((aextail == 0.0) && (aeytail == 0.0) && (aeztail == 0.0) && (bextail == 0.0) &&
2316 (beytail == 0.0) && (beztail == 0.0) && (cextail == 0.0) && (ceytail == 0.0) &&
2317 (ceztail == 0.0) && (dextail == 0.0) && (deytail == 0.0) && (deztail == 0.0)) {
2322 abeps = (aex * beytail + bey * aextail) - (aey * bextail + bex * aeytail);
2323 bceps = (bex * ceytail + cey * bextail) - (bey * cextail + cex * beytail);
2324 cdeps = (cex * deytail + dey * cextail) - (cey * dextail + dex * ceytail);
2325 daeps = (dex * aeytail + aey * dextail) - (dey * aextail + aex * deytail);
2326 aceps = (aex * ceytail + cey * aextail) - (aey * cextail + cex * aeytail);
2327 bdeps = (bex * deytail + dey * bextail) - (bey * dextail + dex * beytail);
2329 (((bex * bex + bey * bey + bez * bez) * ((cez * daeps + dez * aceps + aez * cdeps) +
2330 (ceztail * da3 + deztail * ac3 + aeztail * cd3)) +
2331 (dex * dex + dey * dey + dez * dez) * ((aez * bceps - bez * aceps + cez * abeps) +
2332 (aeztail * bc3 - beztail * ac3 + ceztail * ab3))) -
2333 ((aex * aex + aey * aey + aez * aez) * ((bez * cdeps - cez * bdeps + dez * bceps) +
2334 (beztail * cd3 - ceztail * bd3 + deztail * bc3)) +
2335 (cex * cex + cey * cey + cez * cez) * ((dez * abeps + aez * bdeps + bez * daeps) +
2336 (deztail * ab3 + aeztail * bd3 + beztail * da3)))) +
2338 (((bex * bextail + bey * beytail + bez * beztail) * (cez * da3 + dez * ac3 + aez * cd3) +
2339 (dex * dextail + dey * deytail + dez * deztail) *
2340 (aez * bc3 - bez * ac3 + cez * ab3)) -
2341 ((aex * aextail + aey * aeytail + aez * aeztail) * (bez * cd3 - cez * bd3 + dez * bc3) +
2342 (cex * cextail + cey * ceytail + cez * ceztail) *
2343 (dez * ab3 + aez * bd3 + bez * da3)));
2344 if ((det >= errbound) || (-det >= errbound)) {
2352 const double *pa,
const double *pb,
const double *pc,
const double *pd,
const double *pe)
2354 double aex, bex, cex, dex;
2355 double aey, bey, cey, dey;
2356 double aez, bez, cez, dez;
2357 double aexbey, bexaey, bexcey, cexbey, cexdey, dexcey, dexaey, aexdey;
2358 double aexcey, cexaey, bexdey, dexbey;
2359 double alift, blift, clift, dlift;
2360 double ab, bc, cd, da, ac, bd;
2361 double abc, bcd, cda, dab;
2362 double aezplus, bezplus, cezplus, dezplus;
2363 double aexbeyplus, bexaeyplus, bexceyplus, cexbeyplus;
2364 double cexdeyplus, dexceyplus, dexaeyplus, aexdeyplus;
2365 double aexceyplus, cexaeyplus, bexdeyplus, dexbeyplus;
2367 double permanent, errbound;
2369 aex = pa[0] - pe[0];
2370 bex = pb[0] - pe[0];
2371 cex = pc[0] - pe[0];
2372 dex = pd[0] - pe[0];
2373 aey = pa[1] - pe[1];
2374 bey = pb[1] - pe[1];
2375 cey = pc[1] - pe[1];
2376 dey = pd[1] - pe[1];
2377 aez = pa[2] - pe[2];
2378 bez = pb[2] - pe[2];
2379 cez = pc[2] - pe[2];
2380 dez = pd[2] - pe[2];
2384 ab = aexbey - bexaey;
2387 bc = bexcey - cexbey;
2390 cd = cexdey - dexcey;
2393 da = dexaey - aexdey;
2397 ac = aexcey - cexaey;
2400 bd = bexdey - dexbey;
2402 abc = aez * bc - bez * ac + cez * ab;
2403 bcd = bez * cd - cez * bd + dez * bc;
2404 cda = cez * da + dez * ac + aez * cd;
2405 dab = dez * ab + aez * bd + bez * da;
2407 alift = aex * aex + aey * aey + aez * aez;
2408 blift = bex * bex + bey * bey + bez * bez;
2409 clift = cex * cex + cey * cey + cez * cez;
2410 dlift = dex * dex + dey * dey + dez * dez;
2412 det = (dlift * abc - clift * dab) + (blift * cda - alift * bcd);
2430 permanent = ((cexdeyplus + dexceyplus) * bezplus + (dexbeyplus + bexdeyplus) * cezplus +
2431 (bexceyplus + cexbeyplus) * dezplus) *
2433 ((dexaeyplus + aexdeyplus) * cezplus + (aexceyplus + cexaeyplus) * dezplus +
2434 (cexdeyplus + dexceyplus) * aezplus) *
2436 ((aexbeyplus + bexaeyplus) * dezplus + (bexdeyplus + dexbeyplus) * aezplus +
2437 (dexaeyplus + aexdeyplus) * bezplus) *
2439 ((bexceyplus + cexbeyplus) * aezplus + (cexaeyplus + aexceyplus) * bezplus +
2440 (aexbeyplus + bexaeyplus) * cezplus) *
2443 if ((det > errbound) || (-det > errbound)) {
2454 return (
x > 0) ? 1 : ((
x < 0) ? -1 : 0);
Math vector functions needed specifically for mesh intersect and boolean.
typedef double(DMatrix)[4][4]
ATTR_WARN_UNUSED_RESULT const BMVert const BMEdge * e
ATTR_WARN_UNUSED_RESULT const BMVert * v
SIMD_FORCE_INLINE const btScalar & w() const
Return the w value.
static T sum(const btAlignedObjectArray< T > &items)
#define Two_Product(a, b, x, y)
#define Two_Product_Presplit(a, b, bhi, blo, x, y)
#define Two_Diff_Tail(a, b, x, y)
#define Two_Sum(a, b, x, y)
#define Two_One_Product(a1, a0, b, x3, x2, x1, x0)
#define Two_Two_Diff(a1, a0, b1, b0, x3, x2, x1, x0)
#define Two_Two_Sum(a1, a0, b1, b0, x3, x2, x1, x0)
#define Fast_Two_Sum(a, b, x, y)
#define Split(a, ahi, alo)
static double B3(double u)
static int fast_expansion_sum_zeroelim(int elen, const double *e, int flen, const double *f, double *h)
static double iccerrboundA
double orient2dfast(const double *pa, const double *pb, const double *pc)
static double resulterrbound
static double ccwerrboundB
static double orient2dadapt(const double *pa, const double *pb, const double *pc, double detsum)
static double o3derrboundC
double incirclefast(const double *pa, const double *pb, const double *pc, const double *pd)
static double isperrboundB
static double o3derrboundB
double insphere(const double *pa, const double *pb, const double *pc, const double *pd, const double *pe)
static double isperrboundC
static double ccwerrboundC
static double iccerrboundB
static double estimate(int elen, const double *e)
double orient3dfast(const double *pa, const double *pb, const double *pc, const double *pd)
static double isperrboundA
static RobustInitCaller init_caller
static double incircleadapt(const double *pa, const double *pb, const double *pc, const double *pd, double permanent)
static double insphereadapt(const double *pa, const double *pb, const double *pc, const double *pd, const double *pe, double permanent)
static double o3derrboundA
static double orient3dadapt(const double *pa, const double *pb, const double *pc, const double *pd, double permanent)
double orient3d(const double *pa, const double *pb, const double *pc, const double *pd)
static double iccerrboundC
double incircle(const double *pa, const double *pb, const double *pc, const double *pd)
static int scale_expansion_zeroelim(int elen, const double *e, double b, double *h)
double inspherefast(const double *pa, const double *pb, const double *pc, const double *pd, const double *pe)
double orient2d(const double *pa, const double *pb, const double *pc)
static double ccwerrboundA
static double insphereexact(const double *pa, const double *pb, const double *pc, const double *pd, const double *pe)
int orient3d_fast(const double3 &a, const double3 &b, const double3 &c, const double3 &d)
int orient2d(const double2 &a, const double2 &b, const double2 &c)
int incircle(const double2 &a, const double2 &b, const double2 &c, const double2 &d)
int orient3d(const double3 &a, const double3 &b, const double3 &c, const double3 &d)
int incircle_fast(const double2 &a, const double2 &b, const double2 &c, const double2 &d)
int orient2d_fast(const double2 &a, const double2 &b, const double2 &c)
int insphere_fast(const double3 &a, const double3 &b, const double3 &c, const double3 &d, const double3 &e)
int insphere(const double3 &a, const double3 &b, const double3 &c, const double3 &d, const double3 &e)
static const pxr::TfToken b("b", pxr::TfToken::Immortal)