MCPcopy Create free account
hub / github.com/PyMesh/PyMesh / incircleslow

Function incircleslow

tools/Predicates/predicates.c:2492–2650  ·  view source on GitHub ↗
(pa, pb, pc, pd)

Source from the content-addressed store, hash-verified

2490}
2491
2492REAL incircleslow(pa, pb, pc, pd)
2493REAL *pa;
2494REAL *pb;
2495REAL *pc;
2496REAL *pd;
2497{
2498 INEXACT REAL adx, bdx, cdx, ady, bdy, cdy;
2499 REAL adxtail, bdxtail, cdxtail;
2500 REAL adytail, bdytail, cdytail;
2501 REAL negate, negatetail;
2502 INEXACT REAL axby7, bxcy7, axcy7, bxay7, cxby7, cxay7;
2503 REAL axby[8], bxcy[8], axcy[8], bxay[8], cxby[8], cxay[8];
2504 REAL temp16[16];
2505 int temp16len;
2506 REAL detx[32], detxx[64], detxt[32], detxxt[64], detxtxt[64];
2507 int xlen, xxlen, xtlen, xxtlen, xtxtlen;
2508 REAL x1[128], x2[192];
2509 int x1len, x2len;
2510 REAL dety[32], detyy[64], detyt[32], detyyt[64], detytyt[64];
2511 int ylen, yylen, ytlen, yytlen, ytytlen;
2512 REAL y1[128], y2[192];
2513 int y1len, y2len;
2514 REAL adet[384], bdet[384], cdet[384], abdet[768], deter[1152];
2515 int alen, blen, clen, ablen, deterlen;
2516 int i;
2517
2518 INEXACT REAL bvirt;
2519 REAL avirt, bround, around;
2520 INEXACT REAL c;
2521 INEXACT REAL abig;
2522 REAL a0hi, a0lo, a1hi, a1lo, bhi, blo;
2523 REAL err1, err2, err3;
2524 INEXACT REAL _i, _j, _k, _l, _m, _n;
2525 REAL _0, _1, _2;
2526
2527 Two_Diff(pa[0], pd[0], adx, adxtail);
2528 Two_Diff(pa[1], pd[1], ady, adytail);
2529 Two_Diff(pb[0], pd[0], bdx, bdxtail);
2530 Two_Diff(pb[1], pd[1], bdy, bdytail);
2531 Two_Diff(pc[0], pd[0], cdx, cdxtail);
2532 Two_Diff(pc[1], pd[1], cdy, cdytail);
2533
2534 Two_Two_Product(adx, adxtail, bdy, bdytail,
2535 axby7, axby[6], axby[5], axby[4],
2536 axby[3], axby[2], axby[1], axby[0]);
2537 axby[7] = axby7;
2538 negate = -ady;
2539 negatetail = -adytail;
2540 Two_Two_Product(bdx, bdxtail, negate, negatetail,
2541 bxay7, bxay[6], bxay[5], bxay[4],
2542 bxay[3], bxay[2], bxay[1], bxay[0]);
2543 bxay[7] = bxay7;
2544 Two_Two_Product(bdx, bdxtail, cdy, cdytail,
2545 bxcy7, bxcy[6], bxcy[5], bxcy[4],
2546 bxcy[3], bxcy[2], bxcy[1], bxcy[0]);
2547 bxcy[7] = bxcy7;
2548 negate = -bdy;
2549 negatetail = -bdytail;

Callers

nothing calls this directly

Calls 2

scale_expansion_zeroelimFunction · 0.85

Tested by

no test coverage detected