robust-predicates 2.0.2 → 3.0.1

package/umd/insphere.js

@@ -1,14 +1,15 @@
 (function (global, factory) {
 typeof exports === 'object' && typeof module !== 'undefined' ? factory(exports) :
 typeof define === 'function' && define.amd ? define(['exports'], factory) :
-(global = global || self, factory(global.predicates = {}));
-}(this, function (exports) { 'use strict';
+(global = typeof globalThis !== 'undefined' ? globalThis : global || self, factory(global.predicates = {}));
+}(this, (function (exports) { 'use strict';
 
 const epsilon = 1.1102230246251565e-16;
 const splitter = 134217729;
 const resulterrbound = (3 + 8 * epsilon) * epsilon;
 
-function fast_expansion_sum_zeroelim(elen, e, flen, f, h) {
+// fast_expansion_sum_zeroelim routine from oritinal code
+function sum(elen, e, flen, f, h) {
     let Q, Qnew, hh, bvirt;
     let enow = e[0];
     let fnow = f[0];
@@ -22,7 +23,7 @@ function fast_expansion_sum_zeroelim(elen, e, flen, f, h) {
         fnow = f[++findex];
     }
     let hindex = 0;
-    if ((eindex < elen) && (findex < flen)) {
+    if (eindex < elen && findex < flen) {
         if ((fnow > enow) === (fnow > -enow)) {
             Qnew = enow + Q;
             hh = Q - (Qnew - enow);
@@ -36,7 +37,7 @@ function fast_expansion_sum_zeroelim(elen, e, flen, f, h) {
         if (hh !== 0) {
             h[hindex++] = hh;
         }
-        while ((eindex < elen) && (findex < flen)) {
+        while (eindex < elen && findex < flen) {
             if ((fnow > enow) === (fnow > -enow)) {
                 Qnew = Q + enow;
                 bvirt = Qnew - Q;
@@ -80,24 +81,30 @@ function fast_expansion_sum_zeroelim(elen, e, flen, f, h) {
     return hindex;
 }
 
-function scale_expansion_zeroelim(elen, e, b, h) {
+function sum_three(alen, a, blen, b, clen, c, tmp, out) {
+    return sum(sum(alen, a, blen, b, tmp), tmp, clen, c, out);
+}
+
+// scale_expansion_zeroelim routine from oritinal code
+function scale(elen, e, b, h) {
     let Q, sum, hh, product1, product0;
     let bvirt, c, ahi, alo, bhi, blo;
 
     c = splitter * b;
     bhi = c - (c - b);
     blo = b - bhi;
-    Q = e[0] * b;
-    c = splitter * e[0];
-    ahi = c - (c - e[0]);
-    alo = e[0] - ahi;
+    let enow = e[0];
+    Q = enow * b;
+    c = splitter * enow;
+    ahi = c - (c - enow);
+    alo = enow - ahi;
     hh = alo * blo - (Q - ahi * bhi - alo * bhi - ahi * blo);
     let hindex = 0;
     if (hh !== 0) {
         h[hindex++] = hh;
     }
-    for (let eindex = 1; eindex < elen; eindex++) {
-        const enow = e[eindex];
+    for (let i = 1; i < elen; i++) {
+        enow = e[i];
         product1 = enow * b;
         c = splitter * enow;
         ahi = c - (c - enow);
@@ -115,12 +122,17 @@ function scale_expansion_zeroelim(elen, e, b, h) {
             h[hindex++] = hh;
         }
     }
-    if (Q !== 0.0 || hindex === 0) {
+    if (Q !== 0 || hindex === 0) {
         h[hindex++] = Q;
     }
     return hindex;
 }
 
+function negate(elen, e) {
+    for (let i = 0; i < elen; i++) e[i] = -e[i];
+    return elen;
+}
+
 function estimate(elen, e) {
     let Q = e[0];
     for (let i = 1; i < elen; i++) Q += e[i];
@@ -145,9 +157,7 @@ const bd = vec(4);
 const ce = vec(4);
 const da = vec(4);
 const eb = vec(4);
-const temp8a = vec(8);
-const temp8b = vec(8);
-const temp16 = vec(16);
+
 const abc = vec(24);
 const bcd = vec(24);
 const cde = vec(24);
@@ -158,18 +168,7 @@ const bce = vec(24);
 const cda = vec(24);
 const deb = vec(24);
 const eac = vec(24);
-const temp48a = vec(48);
-const temp48b = vec(48);
-const abcd = vec(96);
-const bcde = vec(96);
-const cdea = vec(96);
-const deab = vec(96);
-const eabc = vec(96);
-const temp192 = vec(192);
-const det384x = vec(384);
-const det384y = vec(384);
-const det384z = vec(384);
-const detxy = vec(768);
+
 const adet = vec(1152);
 const bdet = vec(1152);
 const cdet = vec(1152);
@@ -180,481 +179,375 @@ const cddet = vec(2304);
 const cdedet = vec(3456);
 const deter = vec(5760);
 
-function insphereexact(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz, ex, ey, ez) {
-    let axby1, bxcy1, cxdy1, dxey1, exay1;
-    let bxay1, cxby1, dxcy1, exdy1, axey1;
-    let axcy1, bxdy1, cxey1, dxay1, exby1;
-    let cxay1, dxby1, excy1, axdy1, bxey1;
-    let axby0, bxcy0, cxdy0, dxey0, exay0;
-    let bxay0, cxby0, dxcy0, exdy0, axey0;
-    let axcy0, bxdy0, cxey0, dxay0, exby0;
-    let cxay0, dxby0, excy0, axdy0, bxey0;
-    let temp8alen, temp8blen, temp16len;
-    let temp48alen, temp48blen;
-    let xlen, ylen, zlen;
-    let xylen;
-    let i;
+const _8 = vec(8);
+const _8b = vec(8);
+const _8c = vec(8);
+const _16 = vec(16);
+const _24 = vec(24);
+const _48 = vec(48);
+const _48b = vec(48);
+const _96 = vec(96);
+const _192 = vec(192);
+const _384x = vec(384);
+const _384y = vec(384);
+const _384z = vec(384);
+const _768 = vec(768);
+
+function sum_three_scale(a, b, c, az, bz, cz, out) {
+    return sum_three(
+        scale(4, a, az, _8), _8,
+        scale(4, b, bz, _8b), _8b,
+        scale(4, c, cz, _8c), _8c, _16, out);
+}
+
+function liftexact(alen, a, blen, b, clen, c, dlen, d, x, y, z, out) {
+    const len = sum(
+        sum(alen, a, blen, b, _48), _48,
+        negate(sum(clen, c, dlen, d, _48b), _48b), _48b, _96);
 
-    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0;
+    return sum_three(
+        scale(scale(len, _96, x, _192), _192, x, _384x), _384x,
+        scale(scale(len, _96, y, _192), _192, y, _384y), _384y,
+        scale(scale(len, _96, z, _192), _192, z, _384z), _384z, _768, out);
+}
+
+function insphereexact(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz, ex, ey, ez) {
+    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0, s1, s0, t1, t0, u3;
 
-    axby1 = ax * by;
+    s1 = ax * by;
     c = splitter * ax;
     ahi = c - (c - ax);
     alo = ax - ahi;
     c = splitter * by;
     bhi = c - (c - by);
     blo = by - bhi;
-    axby0 = alo * blo - (axby1 - ahi * bhi - alo * bhi - ahi * blo);
-    bxay1 = bx * ay;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = bx * ay;
     c = splitter * bx;
     ahi = c - (c - bx);
     alo = bx - ahi;
     c = splitter * ay;
     bhi = c - (c - ay);
     blo = ay - bhi;
-    bxay0 = alo * blo - (bxay1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = axby0 - bxay0;
-    bvirt = axby0 - _i;
-    ab[0] = axby0 - (_i + bvirt) + (bvirt - bxay0);
-    _j = axby1 + _i;
-    bvirt = _j - axby1;
-    _0 = axby1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - bxay1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ab[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    ab[1] = _0 - (_i + bvirt) + (bvirt - bxay1);
-    ab[3] = _j + _i;
-    bvirt = ab[3] - _j;
-    ab[2] = _j - (ab[3] - bvirt) + (_i - bvirt);
-
-    bxcy1 = bx * cy;
+    ab[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    ab[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    ab[3] = u3;
+    s1 = bx * cy;
     c = splitter * bx;
     ahi = c - (c - bx);
     alo = bx - ahi;
     c = splitter * cy;
     bhi = c - (c - cy);
     blo = cy - bhi;
-    bxcy0 = alo * blo - (bxcy1 - ahi * bhi - alo * bhi - ahi * blo);
-    cxby1 = cx * by;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = cx * by;
     c = splitter * cx;
     ahi = c - (c - cx);
     alo = cx - ahi;
     c = splitter * by;
     bhi = c - (c - by);
     blo = by - bhi;
-    cxby0 = alo * blo - (cxby1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = bxcy0 - cxby0;
-    bvirt = bxcy0 - _i;
-    bc[0] = bxcy0 - (_i + bvirt) + (bvirt - cxby0);
-    _j = bxcy1 + _i;
-    bvirt = _j - bxcy1;
-    _0 = bxcy1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - cxby1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    bc[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    bc[1] = _0 - (_i + bvirt) + (bvirt - cxby1);
-    bc[3] = _j + _i;
-    bvirt = bc[3] - _j;
-    bc[2] = _j - (bc[3] - bvirt) + (_i - bvirt);
-
-    cxdy1 = cx * dy;
+    bc[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    bc[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    bc[3] = u3;
+    s1 = cx * dy;
     c = splitter * cx;
     ahi = c - (c - cx);
     alo = cx - ahi;
     c = splitter * dy;
     bhi = c - (c - dy);
     blo = dy - bhi;
-    cxdy0 = alo * blo - (cxdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    dxcy1 = dx * cy;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = dx * cy;
     c = splitter * dx;
     ahi = c - (c - dx);
     alo = dx - ahi;
     c = splitter * cy;
     bhi = c - (c - cy);
     blo = cy - bhi;
-    dxcy0 = alo * blo - (dxcy1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = cxdy0 - dxcy0;
-    bvirt = cxdy0 - _i;
-    cd[0] = cxdy0 - (_i + bvirt) + (bvirt - dxcy0);
-    _j = cxdy1 + _i;
-    bvirt = _j - cxdy1;
-    _0 = cxdy1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - dxcy1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    cd[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    cd[1] = _0 - (_i + bvirt) + (bvirt - dxcy1);
-    cd[3] = _j + _i;
-    bvirt = cd[3] - _j;
-    cd[2] = _j - (cd[3] - bvirt) + (_i - bvirt);
-
-    dxey1 = dx * ey;
+    cd[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    cd[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    cd[3] = u3;
+    s1 = dx * ey;
     c = splitter * dx;
     ahi = c - (c - dx);
     alo = dx - ahi;
     c = splitter * ey;
     bhi = c - (c - ey);
     blo = ey - bhi;
-    dxey0 = alo * blo - (dxey1 - ahi * bhi - alo * bhi - ahi * blo);
-    exdy1 = ex * dy;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = ex * dy;
     c = splitter * ex;
     ahi = c - (c - ex);
     alo = ex - ahi;
     c = splitter * dy;
     bhi = c - (c - dy);
     blo = dy - bhi;
-    exdy0 = alo * blo - (exdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = dxey0 - exdy0;
-    bvirt = dxey0 - _i;
-    de[0] = dxey0 - (_i + bvirt) + (bvirt - exdy0);
-    _j = dxey1 + _i;
-    bvirt = _j - dxey1;
-    _0 = dxey1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - exdy1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    de[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    de[1] = _0 - (_i + bvirt) + (bvirt - exdy1);
-    de[3] = _j + _i;
-    bvirt = de[3] - _j;
-    de[2] = _j - (de[3] - bvirt) + (_i - bvirt);
-
-    exay1 = ex * ay;
+    de[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    de[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    de[3] = u3;
+    s1 = ex * ay;
     c = splitter * ex;
     ahi = c - (c - ex);
     alo = ex - ahi;
     c = splitter * ay;
     bhi = c - (c - ay);
     blo = ay - bhi;
-    exay0 = alo * blo - (exay1 - ahi * bhi - alo * bhi - ahi * blo);
-    axey1 = ax * ey;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = ax * ey;
     c = splitter * ax;
     ahi = c - (c - ax);
     alo = ax - ahi;
     c = splitter * ey;
     bhi = c - (c - ey);
     blo = ey - bhi;
-    axey0 = alo * blo - (axey1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = exay0 - axey0;
-    bvirt = exay0 - _i;
-    ea[0] = exay0 - (_i + bvirt) + (bvirt - axey0);
-    _j = exay1 + _i;
-    bvirt = _j - exay1;
-    _0 = exay1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - axey1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ea[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    ea[1] = _0 - (_i + bvirt) + (bvirt - axey1);
-    ea[3] = _j + _i;
-    bvirt = ea[3] - _j;
-    ea[2] = _j - (ea[3] - bvirt) + (_i - bvirt);
-
-    axcy1 = ax * cy;
+    ea[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    ea[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    ea[3] = u3;
+    s1 = ax * cy;
     c = splitter * ax;
     ahi = c - (c - ax);
     alo = ax - ahi;
     c = splitter * cy;
     bhi = c - (c - cy);
     blo = cy - bhi;
-    axcy0 = alo * blo - (axcy1 - ahi * bhi - alo * bhi - ahi * blo);
-    cxay1 = cx * ay;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = cx * ay;
     c = splitter * cx;
     ahi = c - (c - cx);
     alo = cx - ahi;
     c = splitter * ay;
     bhi = c - (c - ay);
     blo = ay - bhi;
-    cxay0 = alo * blo - (cxay1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = axcy0 - cxay0;
-    bvirt = axcy0 - _i;
-    ac[0] = axcy0 - (_i + bvirt) + (bvirt - cxay0);
-    _j = axcy1 + _i;
-    bvirt = _j - axcy1;
-    _0 = axcy1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - cxay1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ac[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    ac[1] = _0 - (_i + bvirt) + (bvirt - cxay1);
-    ac[3] = _j + _i;
-    bvirt = ac[3] - _j;
-    ac[2] = _j - (ac[3] - bvirt) + (_i - bvirt);
-
-    bxdy1 = bx * dy;
+    ac[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    ac[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    ac[3] = u3;
+    s1 = bx * dy;
     c = splitter * bx;
     ahi = c - (c - bx);
     alo = bx - ahi;
     c = splitter * dy;
     bhi = c - (c - dy);
     blo = dy - bhi;
-    bxdy0 = alo * blo - (bxdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    dxby1 = dx * by;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = dx * by;
     c = splitter * dx;
     ahi = c - (c - dx);
     alo = dx - ahi;
     c = splitter * by;
     bhi = c - (c - by);
     blo = by - bhi;
-    dxby0 = alo * blo - (dxby1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = bxdy0 - dxby0;
-    bvirt = bxdy0 - _i;
-    bd[0] = bxdy0 - (_i + bvirt) + (bvirt - dxby0);
-    _j = bxdy1 + _i;
-    bvirt = _j - bxdy1;
-    _0 = bxdy1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - dxby1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    bd[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    bd[1] = _0 - (_i + bvirt) + (bvirt - dxby1);
-    bd[3] = _j + _i;
-    bvirt = bd[3] - _j;
-    bd[2] = _j - (bd[3] - bvirt) + (_i - bvirt);
-
-    cxey1 = cx * ey;
+    bd[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    bd[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    bd[3] = u3;
+    s1 = cx * ey;
     c = splitter * cx;
     ahi = c - (c - cx);
     alo = cx - ahi;
     c = splitter * ey;
     bhi = c - (c - ey);
     blo = ey - bhi;
-    cxey0 = alo * blo - (cxey1 - ahi * bhi - alo * bhi - ahi * blo);
-    excy1 = ex * cy;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = ex * cy;
     c = splitter * ex;
     ahi = c - (c - ex);
     alo = ex - ahi;
     c = splitter * cy;
     bhi = c - (c - cy);
     blo = cy - bhi;
-    excy0 = alo * blo - (excy1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = cxey0 - excy0;
-    bvirt = cxey0 - _i;
-    ce[0] = cxey0 - (_i + bvirt) + (bvirt - excy0);
-    _j = cxey1 + _i;
-    bvirt = _j - cxey1;
-    _0 = cxey1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - excy1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ce[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    ce[1] = _0 - (_i + bvirt) + (bvirt - excy1);
-    ce[3] = _j + _i;
-    bvirt = ce[3] - _j;
-    ce[2] = _j - (ce[3] - bvirt) + (_i - bvirt);
-
-    dxay1 = dx * ay;
+    ce[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    ce[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    ce[3] = u3;
+    s1 = dx * ay;
     c = splitter * dx;
     ahi = c - (c - dx);
     alo = dx - ahi;
     c = splitter * ay;
     bhi = c - (c - ay);
     blo = ay - bhi;
-    dxay0 = alo * blo - (dxay1 - ahi * bhi - alo * bhi - ahi * blo);
-    axdy1 = ax * dy;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = ax * dy;
     c = splitter * ax;
     ahi = c - (c - ax);
     alo = ax - ahi;
     c = splitter * dy;
     bhi = c - (c - dy);
     blo = dy - bhi;
-    axdy0 = alo * blo - (axdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = dxay0 - axdy0;
-    bvirt = dxay0 - _i;
-    da[0] = dxay0 - (_i + bvirt) + (bvirt - axdy0);
-    _j = dxay1 + _i;
-    bvirt = _j - dxay1;
-    _0 = dxay1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - axdy1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    da[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    da[1] = _0 - (_i + bvirt) + (bvirt - axdy1);
-    da[3] = _j + _i;
-    bvirt = da[3] - _j;
-    da[2] = _j - (da[3] - bvirt) + (_i - bvirt);
-
-    exby1 = ex * by;
+    da[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    da[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    da[3] = u3;
+    s1 = ex * by;
     c = splitter * ex;
     ahi = c - (c - ex);
     alo = ex - ahi;
     c = splitter * by;
     bhi = c - (c - by);
     blo = by - bhi;
-    exby0 = alo * blo - (exby1 - ahi * bhi - alo * bhi - ahi * blo);
-    bxey1 = bx * ey;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = bx * ey;
     c = splitter * bx;
     ahi = c - (c - bx);
     alo = bx - ahi;
     c = splitter * ey;
     bhi = c - (c - ey);
     blo = ey - bhi;
-    bxey0 = alo * blo - (bxey1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = exby0 - bxey0;
-    bvirt = exby0 - _i;
-    eb[0] = exby0 - (_i + bvirt) + (bvirt - bxey0);
-    _j = exby1 + _i;
-    bvirt = _j - exby1;
-    _0 = exby1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - bxey1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    eb[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    eb[1] = _0 - (_i + bvirt) + (bvirt - bxey1);
-    eb[3] = _j + _i;
-    bvirt = eb[3] - _j;
-    eb[2] = _j - (eb[3] - bvirt) + (_i - bvirt);
-
-    temp8alen = scale_expansion_zeroelim(4, bc, az, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, ac, -bz, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, ab, cz, temp8a);
-    const abclen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, abc);
-
-    temp8alen = scale_expansion_zeroelim(4, cd, bz, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, bd, -cz, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, bc, dz, temp8a);
-    const bcdlen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, bcd);
-
-    temp8alen = scale_expansion_zeroelim(4, de, cz, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, ce, -dz, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, cd, ez, temp8a);
-    let cdelen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, cde);
-
-    temp8alen = scale_expansion_zeroelim(4, ea, dz, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, da, -ez, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, de, az, temp8a);
-    const dealen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, dea);
-
-    temp8alen = scale_expansion_zeroelim(4, ab, ez, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, eb, -az, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, ea, bz, temp8a);
-    const eablen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, eab);
-
-    temp8alen = scale_expansion_zeroelim(4, bd, az, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, da, bz, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, ab, dz, temp8a);
-    const abdlen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, abd);
-
-    temp8alen = scale_expansion_zeroelim(4, ce, bz, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, eb, cz, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, bc, ez, temp8a);
-    const bcelen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, bce);
-
-    temp8alen = scale_expansion_zeroelim(4, da, cz, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, ac, dz, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, cd, az, temp8a);
-    const cdalen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, cda);
-
-    temp8alen = scale_expansion_zeroelim(4, eb, dz, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, bd, ez, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, de, bz, temp8a);
-    const deblen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, deb);
-
-    temp8alen = scale_expansion_zeroelim(4, ac, ez, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, ce, az, temp8b);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp8alen = scale_expansion_zeroelim(4, ea, cz, temp8a);
-    const eaclen = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp16len, temp16, eac);
-
-    temp48alen = fast_expansion_sum_zeroelim(cdelen, cde, bcelen, bce, temp48a);
-    temp48blen = fast_expansion_sum_zeroelim(deblen, deb, bcdlen, bcd, temp48b);
-    for (i = 0; i < temp48blen; i++) {
-        temp48b[i] = -temp48b[i];
-    }
-    const bcdelen = fast_expansion_sum_zeroelim(temp48alen, temp48a, temp48blen, temp48b, bcde);
-    xlen = scale_expansion_zeroelim(bcdelen, bcde, ax, temp192);
-    xlen = scale_expansion_zeroelim(xlen, temp192, ax, det384x);
-    ylen = scale_expansion_zeroelim(bcdelen, bcde, ay, temp192);
-    ylen = scale_expansion_zeroelim(ylen, temp192, ay, det384y);
-    zlen = scale_expansion_zeroelim(bcdelen, bcde, az, temp192);
-    zlen = scale_expansion_zeroelim(zlen, temp192, az, det384z);
-    xylen = fast_expansion_sum_zeroelim(xlen, det384x, ylen, det384y, detxy);
-    const alen = fast_expansion_sum_zeroelim(xylen, detxy, zlen, det384z, adet);
-
-    temp48alen = fast_expansion_sum_zeroelim(dealen, dea, cdalen, cda, temp48a);
-    temp48blen = fast_expansion_sum_zeroelim(eaclen, eac, cdelen, cde, temp48b);
-    for (i = 0; i < temp48blen; i++) {
-        temp48b[i] = -temp48b[i];
-    }
-    const cdealen = fast_expansion_sum_zeroelim(temp48alen, temp48a, temp48blen, temp48b, cdea);
-    xlen = scale_expansion_zeroelim(cdealen, cdea, bx, temp192);
-    xlen = scale_expansion_zeroelim(xlen, temp192, bx, det384x);
-    ylen = scale_expansion_zeroelim(cdealen, cdea, by, temp192);
-    ylen = scale_expansion_zeroelim(ylen, temp192, by, det384y);
-    zlen = scale_expansion_zeroelim(cdealen, cdea, bz, temp192);
-    zlen = scale_expansion_zeroelim(zlen, temp192, bz, det384z);
-    xylen = fast_expansion_sum_zeroelim(xlen, det384x, ylen, det384y, detxy);
-    const blen = fast_expansion_sum_zeroelim(xylen, detxy, zlen, det384z, bdet);
-
-    temp48alen = fast_expansion_sum_zeroelim(eablen, eab, deblen, deb, temp48a);
-    temp48blen = fast_expansion_sum_zeroelim(abdlen, abd, dealen, dea, temp48b);
-    for (i = 0; i < temp48blen; i++) {
-        temp48b[i] = -temp48b[i];
-    }
-    const deablen = fast_expansion_sum_zeroelim(temp48alen, temp48a, temp48blen, temp48b, deab);
-    xlen = scale_expansion_zeroelim(deablen, deab, cx, temp192);
-    xlen = scale_expansion_zeroelim(xlen, temp192, cx, det384x);
-    ylen = scale_expansion_zeroelim(deablen, deab, cy, temp192);
-    ylen = scale_expansion_zeroelim(ylen, temp192, cy, det384y);
-    zlen = scale_expansion_zeroelim(deablen, deab, cz, temp192);
-    zlen = scale_expansion_zeroelim(zlen, temp192, cz, det384z);
-    xylen = fast_expansion_sum_zeroelim(xlen, det384x, ylen, det384y, detxy);
-    const clen = fast_expansion_sum_zeroelim(xylen, detxy, zlen, det384z, cdet);
-
-    temp48alen = fast_expansion_sum_zeroelim(abclen, abc, eaclen, eac, temp48a);
-    temp48blen = fast_expansion_sum_zeroelim(bcelen, bce, eablen, eab, temp48b);
-    for (i = 0; i < temp48blen; i++) {
-        temp48b[i] = -temp48b[i];
-    }
-    const eabclen = fast_expansion_sum_zeroelim(temp48alen, temp48a, temp48blen, temp48b, eabc);
-    xlen = scale_expansion_zeroelim(eabclen, eabc, dx, temp192);
-    xlen = scale_expansion_zeroelim(xlen, temp192, dx, det384x);
-    ylen = scale_expansion_zeroelim(eabclen, eabc, dy, temp192);
-    ylen = scale_expansion_zeroelim(ylen, temp192, dy, det384y);
-    zlen = scale_expansion_zeroelim(eabclen, eabc, dz, temp192);
-    zlen = scale_expansion_zeroelim(zlen, temp192, dz, det384z);
-    xylen = fast_expansion_sum_zeroelim(xlen, det384x, ylen, det384y, detxy);
-    const dlen = fast_expansion_sum_zeroelim(xylen, detxy, zlen, det384z, ddet);
-
-    temp48alen = fast_expansion_sum_zeroelim(bcdlen, bcd, abdlen, abd, temp48a);
-    temp48blen = fast_expansion_sum_zeroelim(cdalen, cda, abclen, abc, temp48b);
-    for (i = 0; i < temp48blen; i++) {
-        temp48b[i] = -temp48b[i];
-    }
-    const abcdlen = fast_expansion_sum_zeroelim(temp48alen, temp48a, temp48blen, temp48b, abcd);
-    xlen = scale_expansion_zeroelim(abcdlen, abcd, ex, temp192);
-    xlen = scale_expansion_zeroelim(xlen, temp192, ex, det384x);
-    ylen = scale_expansion_zeroelim(abcdlen, abcd, ey, temp192);
-    ylen = scale_expansion_zeroelim(ylen, temp192, ey, det384y);
-    zlen = scale_expansion_zeroelim(abcdlen, abcd, ez, temp192);
-    zlen = scale_expansion_zeroelim(zlen, temp192, ez, det384z);
-    xylen = fast_expansion_sum_zeroelim(xlen, det384x, ylen, det384y, detxy);
-    const elen = fast_expansion_sum_zeroelim(xylen, detxy, zlen, det384z, edet);
-
-    const ablen = fast_expansion_sum_zeroelim(alen, adet, blen, bdet, abdet);
-    const cdlen = fast_expansion_sum_zeroelim(clen, cdet, dlen, ddet, cddet);
-    cdelen = fast_expansion_sum_zeroelim(cdlen, cddet, elen, edet, cdedet);
-    const deterlen = fast_expansion_sum_zeroelim(ablen, abdet, cdelen, cdedet, deter);
+    eb[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    eb[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    eb[3] = u3;
+
+    const abclen = sum_three_scale(ab, bc, ac, cz, az, -bz, abc);
+    const bcdlen = sum_three_scale(bc, cd, bd, dz, bz, -cz, bcd);
+    const cdelen = sum_three_scale(cd, de, ce, ez, cz, -dz, cde);
+    const dealen = sum_three_scale(de, ea, da, az, dz, -ez, dea);
+    const eablen = sum_three_scale(ea, ab, eb, bz, ez, -az, eab);
+    const abdlen = sum_three_scale(ab, bd, da, dz, az, bz, abd);
+    const bcelen = sum_three_scale(bc, ce, eb, ez, bz, cz, bce);
+    const cdalen = sum_three_scale(cd, da, ac, az, cz, dz, cda);
+    const deblen = sum_three_scale(de, eb, bd, bz, dz, ez, deb);
+    const eaclen = sum_three_scale(ea, ac, ce, cz, ez, az, eac);
+
+    const deterlen = sum_three(
+        liftexact(cdelen, cde, bcelen, bce, deblen, deb, bcdlen, bcd, ax, ay, az, adet), adet,
+        liftexact(dealen, dea, cdalen, cda, eaclen, eac, cdelen, cde, bx, by, bz, bdet), bdet,
+        sum_three(
+            liftexact(eablen, eab, deblen, deb, abdlen, abd, dealen, dea, cx, cy, cz, cdet), cdet,
+            liftexact(abclen, abc, eaclen, eac, bcelen, bce, eablen, eab, dx, dy, dz, ddet), ddet,
+            liftexact(bcdlen, bcd, abdlen, abd, cdalen, cda, abclen, abc, ex, ey, ez, edet), edet, cddet, cdedet), cdedet, abdet, deter);
 
     return deter[deterlen - 1];
 }
 
-const temp8c = vec(8);
-const temp24 = vec(24);
-const temp48 = vec(48);
 const xdet = vec(96);
 const ydet = vec(96);
 const zdet = vec(96);
-const xydet = vec(192);
-const fin1 = vec(1152);
+const fin = vec(1152);
+
+function liftadapt(a, b, c, az, bz, cz, x, y, z, out) {
+    const len = sum_three_scale(a, b, c, az, bz, cz, _24);
+    return sum_three(
+        scale(scale(len, _24, x, _48), _48, x, xdet), xdet,
+        scale(scale(len, _24, y, _48), _48, y, ydet), ydet,
+        scale(scale(len, _24, z, _48), _48, z, zdet), zdet, _192, out);
+}
 
 function insphereadapt(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz, ex, ey, ez, permanent) {
-    let aexbey1, bexaey1, bexcey1, cexbey1;
-    let cexdey1, dexcey1, dexaey1, aexdey1;
-    let aexcey1, cexaey1, bexdey1, dexbey1;
-    let aexbey0, bexaey0, bexcey0, cexbey0;
-    let cexdey0, dexcey0, dexaey0, aexdey0;
-    let aexcey0, cexaey0, bexdey0, dexbey0;
     let ab3, bc3, cd3, da3, ac3, bd3;
-    let temp8alen, temp8blen, temp8clen, temp16len, temp24len, temp48len;
-    let xlen, ylen, zlen, xylen;
 
     let aextail, bextail, cextail, dextail;
     let aeytail, beytail, ceytail, deytail;
     let aeztail, beztail, ceztail, deztail;
 
-    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0;
+    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0, s1, s0, t1, t0;
 
     const aex = ax - ex;
     const bex = bx - ex;
@@ -669,247 +562,190 @@ function insphereadapt(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz, ex, ey, e
     const cez = cz - ez;
     const dez = dz - ez;
 
-    aexbey1 = aex * bey;
+    s1 = aex * bey;
     c = splitter * aex;
     ahi = c - (c - aex);
     alo = aex - ahi;
     c = splitter * bey;
     bhi = c - (c - bey);
     blo = bey - bhi;
-    aexbey0 = alo * blo - (aexbey1 - ahi * bhi - alo * bhi - ahi * blo);
-    bexaey1 = bex * aey;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = bex * aey;
     c = splitter * bex;
     ahi = c - (c - bex);
     alo = bex - ahi;
     c = splitter * aey;
     bhi = c - (c - aey);
     blo = aey - bhi;
-    bexaey0 = alo * blo - (bexaey1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = aexbey0 - bexaey0;
-    bvirt = aexbey0 - _i;
-    ab[0] = aexbey0 - (_i + bvirt) + (bvirt - bexaey0);
-    _j = aexbey1 + _i;
-    bvirt = _j - aexbey1;
-    _0 = aexbey1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - bexaey1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ab[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    ab[1] = _0 - (_i + bvirt) + (bvirt - bexaey1);
+    ab[1] = _0 - (_i + bvirt) + (bvirt - t1);
     ab3 = _j + _i;
     bvirt = ab3 - _j;
     ab[2] = _j - (ab3 - bvirt) + (_i - bvirt);
     ab[3] = ab3;
-
-    bexcey1 = bex * cey;
+    s1 = bex * cey;
     c = splitter * bex;
     ahi = c - (c - bex);
     alo = bex - ahi;
     c = splitter * cey;
     bhi = c - (c - cey);
     blo = cey - bhi;
-    bexcey0 = alo * blo - (bexcey1 - ahi * bhi - alo * bhi - ahi * blo);
-    cexbey1 = cex * bey;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = cex * bey;
     c = splitter * cex;
     ahi = c - (c - cex);
     alo = cex - ahi;
     c = splitter * bey;
     bhi = c - (c - bey);
     blo = bey - bhi;
-    cexbey0 = alo * blo - (cexbey1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = bexcey0 - cexbey0;
-    bvirt = bexcey0 - _i;
-    bc[0] = bexcey0 - (_i + bvirt) + (bvirt - cexbey0);
-    _j = bexcey1 + _i;
-    bvirt = _j - bexcey1;
-    _0 = bexcey1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - cexbey1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    bc[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    bc[1] = _0 - (_i + bvirt) + (bvirt - cexbey1);
+    bc[1] = _0 - (_i + bvirt) + (bvirt - t1);
     bc3 = _j + _i;
     bvirt = bc3 - _j;
     bc[2] = _j - (bc3 - bvirt) + (_i - bvirt);
     bc[3] = bc3;
-
-    cexdey1 = cex * dey;
+    s1 = cex * dey;
     c = splitter * cex;
     ahi = c - (c - cex);
     alo = cex - ahi;
     c = splitter * dey;
     bhi = c - (c - dey);
     blo = dey - bhi;
-    cexdey0 = alo * blo - (cexdey1 - ahi * bhi - alo * bhi - ahi * blo);
-    dexcey1 = dex * cey;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = dex * cey;
     c = splitter * dex;
     ahi = c - (c - dex);
     alo = dex - ahi;
     c = splitter * cey;
     bhi = c - (c - cey);
     blo = cey - bhi;
-    dexcey0 = alo * blo - (dexcey1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = cexdey0 - dexcey0;
-    bvirt = cexdey0 - _i;
-    cd[0] = cexdey0 - (_i + bvirt) + (bvirt - dexcey0);
-    _j = cexdey1 + _i;
-    bvirt = _j - cexdey1;
-    _0 = cexdey1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - dexcey1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    cd[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    cd[1] = _0 - (_i + bvirt) + (bvirt - dexcey1);
+    cd[1] = _0 - (_i + bvirt) + (bvirt - t1);
     cd3 = _j + _i;
     bvirt = cd3 - _j;
     cd[2] = _j - (cd3 - bvirt) + (_i - bvirt);
     cd[3] = cd3;
-
-    dexaey1 = dex * aey;
+    s1 = dex * aey;
     c = splitter * dex;
     ahi = c - (c - dex);
     alo = dex - ahi;
     c = splitter * aey;
     bhi = c - (c - aey);
     blo = aey - bhi;
-    dexaey0 = alo * blo - (dexaey1 - ahi * bhi - alo * bhi - ahi * blo);
-    aexdey1 = aex * dey;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = aex * dey;
     c = splitter * aex;
     ahi = c - (c - aex);
     alo = aex - ahi;
     c = splitter * dey;
     bhi = c - (c - dey);
     blo = dey - bhi;
-    aexdey0 = alo * blo - (aexdey1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = dexaey0 - aexdey0;
-    bvirt = dexaey0 - _i;
-    da[0] = dexaey0 - (_i + bvirt) + (bvirt - aexdey0);
-    _j = dexaey1 + _i;
-    bvirt = _j - dexaey1;
-    _0 = dexaey1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - aexdey1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    da[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    da[1] = _0 - (_i + bvirt) + (bvirt - aexdey1);
+    da[1] = _0 - (_i + bvirt) + (bvirt - t1);
     da3 = _j + _i;
     bvirt = da3 - _j;
     da[2] = _j - (da3 - bvirt) + (_i - bvirt);
     da[3] = da3;
-
-    aexcey1 = aex * cey;
+    s1 = aex * cey;
     c = splitter * aex;
     ahi = c - (c - aex);
     alo = aex - ahi;
     c = splitter * cey;
     bhi = c - (c - cey);
     blo = cey - bhi;
-    aexcey0 = alo * blo - (aexcey1 - ahi * bhi - alo * bhi - ahi * blo);
-    cexaey1 = cex * aey;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = cex * aey;
     c = splitter * cex;
     ahi = c - (c - cex);
     alo = cex - ahi;
     c = splitter * aey;
     bhi = c - (c - aey);
     blo = aey - bhi;
-    cexaey0 = alo * blo - (cexaey1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = aexcey0 - cexaey0;
-    bvirt = aexcey0 - _i;
-    ac[0] = aexcey0 - (_i + bvirt) + (bvirt - cexaey0);
-    _j = aexcey1 + _i;
-    bvirt = _j - aexcey1;
-    _0 = aexcey1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - cexaey1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ac[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    ac[1] = _0 - (_i + bvirt) + (bvirt - cexaey1);
+    ac[1] = _0 - (_i + bvirt) + (bvirt - t1);
     ac3 = _j + _i;
     bvirt = ac3 - _j;
     ac[2] = _j - (ac3 - bvirt) + (_i - bvirt);
     ac[3] = ac3;
-
-    bexdey1 = bex * dey;
+    s1 = bex * dey;
     c = splitter * bex;
     ahi = c - (c - bex);
     alo = bex - ahi;
     c = splitter * dey;
     bhi = c - (c - dey);
     blo = dey - bhi;
-    bexdey0 = alo * blo - (bexdey1 - ahi * bhi - alo * bhi - ahi * blo);
-    dexbey1 = dex * bey;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = dex * bey;
     c = splitter * dex;
     ahi = c - (c - dex);
     alo = dex - ahi;
     c = splitter * bey;
     bhi = c - (c - bey);
     blo = bey - bhi;
-    dexbey0 = alo * blo - (dexbey1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = bexdey0 - dexbey0;
-    bvirt = bexdey0 - _i;
-    bd[0] = bexdey0 - (_i + bvirt) + (bvirt - dexbey0);
-    _j = bexdey1 + _i;
-    bvirt = _j - bexdey1;
-    _0 = bexdey1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - dexbey1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    bd[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    bd[1] = _0 - (_i + bvirt) + (bvirt - dexbey1);
+    bd[1] = _0 - (_i + bvirt) + (bvirt - t1);
     bd3 = _j + _i;
     bvirt = bd3 - _j;
     bd[2] = _j - (bd3 - bvirt) + (_i - bvirt);
     bd[3] = bd3;
 
-    temp8alen = scale_expansion_zeroelim(4, cd, bez, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, bd, -cez, temp8b);
-    temp8clen = scale_expansion_zeroelim(4, bc, dez, temp8c);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp24len = fast_expansion_sum_zeroelim(temp8clen, temp8c, temp16len, temp16, temp24);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, aex, temp48);
-    xlen = scale_expansion_zeroelim(temp48len, temp48, -aex, xdet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, aey, temp48);
-    ylen = scale_expansion_zeroelim(temp48len, temp48, -aey, ydet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, aez, temp48);
-    zlen = scale_expansion_zeroelim(temp48len, temp48, -aez, zdet);
-    xylen = fast_expansion_sum_zeroelim(xlen, xdet, ylen, ydet, xydet);
-    const alen = fast_expansion_sum_zeroelim(xylen, xydet, zlen, zdet, adet);
-
-    temp8alen = scale_expansion_zeroelim(4, da, cez, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, ac, dez, temp8b);
-    temp8clen = scale_expansion_zeroelim(4, cd, aez, temp8c);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp24len = fast_expansion_sum_zeroelim(temp8clen, temp8c, temp16len, temp16, temp24);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, bex, temp48);
-    xlen = scale_expansion_zeroelim(temp48len, temp48, bex, xdet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, bey, temp48);
-    ylen = scale_expansion_zeroelim(temp48len, temp48, bey, ydet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, bez, temp48);
-    zlen = scale_expansion_zeroelim(temp48len, temp48, bez, zdet);
-    xylen = fast_expansion_sum_zeroelim(xlen, xdet, ylen, ydet, xydet);
-    const blen = fast_expansion_sum_zeroelim(xylen, xydet, zlen, zdet, bdet);
-
-    temp8alen = scale_expansion_zeroelim(4, ab, dez, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, bd, aez, temp8b);
-    temp8clen = scale_expansion_zeroelim(4, da, bez, temp8c);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp24len = fast_expansion_sum_zeroelim(temp8clen, temp8c, temp16len, temp16, temp24);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, cex, temp48);
-    xlen = scale_expansion_zeroelim(temp48len, temp48, -cex, xdet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, cey, temp48);
-    ylen = scale_expansion_zeroelim(temp48len, temp48, -cey, ydet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, cez, temp48);
-    zlen = scale_expansion_zeroelim(temp48len, temp48, -cez, zdet);
-    xylen = fast_expansion_sum_zeroelim(xlen, xdet, ylen, ydet, xydet);
-    const clen = fast_expansion_sum_zeroelim(xylen, xydet, zlen, zdet, cdet);
-
-    temp8alen = scale_expansion_zeroelim(4, bc, aez, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, ac, -bez, temp8b);
-    temp8clen = scale_expansion_zeroelim(4, ab, cez, temp8c);
-    temp16len = fast_expansion_sum_zeroelim(temp8alen, temp8a, temp8blen, temp8b, temp16);
-    temp24len = fast_expansion_sum_zeroelim(temp8clen, temp8c, temp16len, temp16, temp24);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, dex, temp48);
-    xlen = scale_expansion_zeroelim(temp48len, temp48, dex, xdet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, dey, temp48);
-    ylen = scale_expansion_zeroelim(temp48len, temp48, dey, ydet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, dez, temp48);
-    zlen = scale_expansion_zeroelim(temp48len, temp48, dez, zdet);
-    xylen = fast_expansion_sum_zeroelim(xlen, xdet, ylen, ydet, xydet);
-    const dlen = fast_expansion_sum_zeroelim(xylen, xydet, zlen, zdet, ddet);
+    const finlen = sum(
+        sum(
+            negate(liftadapt(bc, cd, bd, dez, bez, -cez, aex, aey, aez, adet), adet), adet,
+            liftadapt(cd, da, ac, aez, cez, dez, bex, bey, bez, bdet), bdet, abdet), abdet,
+        sum(
+            negate(liftadapt(da, ab, bd, bez, dez, aez, cex, cey, cez, cdet), cdet), cdet,
+            liftadapt(ab, bc, ac, cez, aez, -bez, dex, dey, dez, ddet), ddet, cddet), cddet, fin);
 
-    const ablen = fast_expansion_sum_zeroelim(alen, adet, blen, bdet, abdet);
-    const cdlen = fast_expansion_sum_zeroelim(clen, cdet, dlen, ddet, cddet);
-    const finlength = fast_expansion_sum_zeroelim(ablen, abdet, cdlen, cddet, fin1);
-
-    let det = estimate(finlength, fin1);
+    let det = estimate(finlen, fin);
     let errbound = isperrboundB * permanent;
     if (det >= errbound || -det >= errbound) {
         return det;
@@ -999,7 +835,6 @@ function insphere(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz, ex, ey, ez) {
     const dexaey = dex * aey;
     const aexdey = aex * dey;
     const da = dexaey - aexdey;
-
     const aexcey = aex * cey;
     const cexaey = cex * aey;
     const ac = aexcey - cexaey;
@@ -1066,7 +901,6 @@ function inspherefast(pax, pay, paz, pbx, pby, pbz, pcx, pcy, pcz, pdx, pdy, pdz
     const bc = bex * cey - cex * bey;
     const cd = cex * dey - dex * cey;
     const da = dex * aey - aex * dey;
-
     const ac = aex * cey - cex * aey;
     const bd = bex * dey - dex * bey;
 
@@ -1088,4 +922,4 @@ exports.inspherefast = inspherefast;
 
 Object.defineProperty(exports, '__esModule', { value: true });
 
-}));
+})));