robust-predicates 2.0.2 → 3.0.1

package/umd/predicates.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];
@@ -139,53 +151,52 @@ const B = vec(4);
 const C1 = vec(8);
 const C2 = vec(12);
 const D = vec(16);
-const u = vec(4);
+const u$2 = vec(4);
 
 function orient2dadapt(ax, ay, bx, by, cx, cy, detsum) {
     let acxtail, acytail, bcxtail, bcytail;
-    let detleft, detright, detlefttail, detrighttail;
-    let B3, u3, s1, t1, s0, t0;
-    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0;
+    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0, s1, s0, t1, t0, u3;
 
     const acx = ax - cx;
     const bcx = bx - cx;
     const acy = ay - cy;
     const bcy = by - cy;
 
-    detleft = acx * bcy;
+    s1 = acx * bcy;
     c = splitter * acx;
     ahi = c - (c - acx);
     alo = acx - ahi;
     c = splitter * bcy;
     bhi = c - (c - bcy);
     blo = bcy - bhi;
-    detlefttail = alo * blo - (detleft - ahi * bhi - alo * bhi - ahi * blo);
-    detright = acy * bcx;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = acy * bcx;
     c = splitter * acy;
     ahi = c - (c - acy);
     alo = acy - ahi;
     c = splitter * bcx;
     bhi = c - (c - bcx);
     blo = bcx - bhi;
-    detrighttail = alo * blo - (detright - ahi * bhi - alo * bhi - ahi * blo);
-
-    _i = detlefttail - detrighttail;
-    bvirt = detlefttail - _i;
-    B[0] = detlefttail - (_i + bvirt) + (bvirt - detrighttail);
-    _j = detleft + _i;
-    bvirt = _j - detleft;
-    _0 = detleft - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - detright;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    B[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    B[1] = _0 - (_i + bvirt) + (bvirt - detright);
-    B3 = _j + _i;
-    bvirt = B3 - _j;
-    B[2] = _j - (B3 - bvirt) + (_i - bvirt);
-    B[3] = B3;
+    B[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    B[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    B[3] = u3;
 
     let det = estimate(4, B);
     let errbound = ccwerrboundB * detsum;
-    if (det >= errbound || -det >= errbound) return det;
+    if (det >= errbound || -det >= errbound) {
+        return det;
+    }
 
     bvirt = ax - acx;
     acxtail = ax - (acx + bvirt) + (bvirt - cx);
@@ -196,13 +207,13 @@ function orient2dadapt(ax, ay, bx, by, cx, cy, detsum) {
     bvirt = by - bcy;
     bcytail = by - (bcy + bvirt) + (bvirt - cy);
 
-    if (acxtail === 0 && acytail === 0 && bcxtail === 0 && bcytail === 0) return det;
+    if (acxtail === 0 && acytail === 0 && bcxtail === 0 && bcytail === 0) {
+        return det;
+    }
 
     errbound = ccwerrboundC * detsum + resulterrbound * Math.abs(det);
     det += (acx * bcytail + bcy * acxtail) - (acy * bcxtail + bcx * acytail);
-    if (det >= errbound || -det >= errbound) {
-        return det;
-    }
+    if (det >= errbound || -det >= errbound) return det;
 
     s1 = acxtail * bcy;
     c = splitter * acxtail;
@@ -222,18 +233,18 @@ function orient2dadapt(ax, ay, bx, by, cx, cy, detsum) {
     t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
     _i = s0 - t0;
     bvirt = s0 - _i;
-    u[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    u$2[0] = s0 - (_i + bvirt) + (bvirt - t0);
     _j = s1 + _i;
     bvirt = _j - s1;
     _0 = s1 - (_j - bvirt) + (_i - bvirt);
     _i = _0 - t1;
     bvirt = _0 - _i;
-    u[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u$2[1] = _0 - (_i + bvirt) + (bvirt - t1);
     u3 = _j + _i;
     bvirt = u3 - _j;
-    u[2] = _j - (u3 - bvirt) + (_i - bvirt);
-    u[3] = u3;
-    const C1length = fast_expansion_sum_zeroelim(4, B, 4, u, C1);
+    u$2[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    u$2[3] = u3;
+    const C1len = sum(4, B, 4, u$2, C1);
 
     s1 = acx * bcytail;
     c = splitter * acx;
@@ -253,18 +264,18 @@ function orient2dadapt(ax, ay, bx, by, cx, cy, detsum) {
     t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
     _i = s0 - t0;
     bvirt = s0 - _i;
-    u[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    u$2[0] = s0 - (_i + bvirt) + (bvirt - t0);
     _j = s1 + _i;
     bvirt = _j - s1;
     _0 = s1 - (_j - bvirt) + (_i - bvirt);
     _i = _0 - t1;
     bvirt = _0 - _i;
-    u[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u$2[1] = _0 - (_i + bvirt) + (bvirt - t1);
     u3 = _j + _i;
     bvirt = u3 - _j;
-    u[2] = _j - (u3 - bvirt) + (_i - bvirt);
-    u[3] = u3;
-    const C2length = fast_expansion_sum_zeroelim(C1length, C1, 4, u, C2);
+    u$2[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    u$2[3] = u3;
+    const C2len = sum(C1len, C1, 4, u$2, C2);
 
     s1 = acxtail * bcytail;
     c = splitter * acxtail;
@@ -284,48 +295,31 @@ function orient2dadapt(ax, ay, bx, by, cx, cy, detsum) {
     t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
     _i = s0 - t0;
     bvirt = s0 - _i;
-    u[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    u$2[0] = s0 - (_i + bvirt) + (bvirt - t0);
     _j = s1 + _i;
     bvirt = _j - s1;
     _0 = s1 - (_j - bvirt) + (_i - bvirt);
     _i = _0 - t1;
     bvirt = _0 - _i;
-    u[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u$2[1] = _0 - (_i + bvirt) + (bvirt - t1);
     u3 = _j + _i;
     bvirt = u3 - _j;
-    u[2] = _j - (u3 - bvirt) + (_i - bvirt);
-    u[3] = u3;
-    const Dlength = fast_expansion_sum_zeroelim(C2length, C2, 4, u, D);
+    u$2[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    u$2[3] = u3;
+    const Dlen = sum(C2len, C2, 4, u$2, D);
 
-    return D[Dlength - 1];
+    return D[Dlen - 1];
 }
 
 function orient2d(ax, ay, bx, by, cx, cy) {
     const detleft = (ay - cy) * (bx - cx);
     const detright = (ax - cx) * (by - cy);
     const det = detleft - detright;
-    let detsum;
-
-    if (detleft > 0) {
-        if (detright <= 0) {
-            return det;
-        }
-        detsum = detleft + detright;
-
-    } else if (detleft < 0) {
-        if (detright >= 0) {
-            return det;
-        }
-        detsum = -detleft - detright;
 
-    } else {
-        return det;
-    }
+    if (detleft === 0 || detright === 0 || (detleft > 0) !== (detright > 0)) return det;
 
-    const errbound = ccwerrboundA * detsum;
-    if (det >= errbound || -det >= errbound) {
-        return det;
-    }
+    const detsum = Math.abs(detleft + detright);
+    if (Math.abs(det) >= ccwerrboundA * detsum) return det;
 
     return -orient2dadapt(ax, ay, bx, by, cx, cy, detsum);
 }
@@ -338,59 +332,210 @@ const o3derrboundA = (7 + 56 * epsilon) * epsilon;
 const o3derrboundB = (3 + 28 * epsilon) * epsilon;
 const o3derrboundC = (26 + 288 * epsilon) * epsilon * epsilon;
 
-const bc = vec(4);
-const ca = vec(4);
-const ab = vec(4);
-const adet = vec(8);
-const bdet = vec(8);
-const cdet = vec(8);
-const abdet = vec(16);
-const fin1 = vec(192);
-const fin2 = vec(192);
+const bc$2 = vec(4);
+const ca$1 = vec(4);
+const ab$2 = vec(4);
 const at_b = vec(4);
 const at_c = vec(4);
 const bt_c = vec(4);
 const bt_a = vec(4);
 const ct_a = vec(4);
 const ct_b = vec(4);
-const bct = vec(8);
-const cat = vec(8);
-const abt = vec(8);
+const bct$1 = vec(8);
+const cat$1 = vec(8);
+const abt$1 = vec(8);
 const u$1 = vec(4);
-const v = vec(12);
-const w = vec(16);
 
-function orient3dadapt(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz, permanent) {
-    let bdxcdy1, cdxbdy1, cdxady1, adxcdy1, adxbdy1, bdxady1;
-    let bdxcdy0, cdxbdy0, cdxady0, adxcdy0, adxbdy0, bdxady0;
-    let bc3, ca3, ab3;
-    let finnow, finother, finswap;
-    let finlength;
+const _8$2 = vec(8);
+const _8b$1 = vec(8);
+const _16$2 = vec(8);
+const _12 = vec(12);
+
+let fin$2 = vec(192);
+let fin2$1 = vec(192);
+
+function finadd$1(finlen, alen, a) {
+    finlen = sum(finlen, fin$2, alen, a, fin2$1);
+    const tmp = fin$2; fin$2 = fin2$1; fin2$1 = tmp;
+    return finlen;
+}
+
+function tailinit(xtail, ytail, ax, ay, bx, by, a, b) {
+    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0, s1, s0, t1, t0, u3, negate;
+    if (xtail === 0) {
+        if (ytail === 0) {
+            a[0] = 0;
+            b[0] = 0;
+            return 1;
+        } else {
+            negate = -ytail;
+            s1 = negate * ax;
+            c = splitter * negate;
+            ahi = c - (c - negate);
+            alo = negate - ahi;
+            c = splitter * ax;
+            bhi = c - (c - ax);
+            blo = ax - bhi;
+            a[0] = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            a[1] = s1;
+            s1 = ytail * bx;
+            c = splitter * ytail;
+            ahi = c - (c - ytail);
+            alo = ytail - ahi;
+            c = splitter * bx;
+            bhi = c - (c - bx);
+            blo = bx - bhi;
+            b[0] = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            b[1] = s1;
+            return 2;
+        }
+    } else {
+        if (ytail === 0) {
+            s1 = xtail * ay;
+            c = splitter * xtail;
+            ahi = c - (c - xtail);
+            alo = xtail - ahi;
+            c = splitter * ay;
+            bhi = c - (c - ay);
+            blo = ay - bhi;
+            a[0] = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            a[1] = s1;
+            negate = -xtail;
+            s1 = negate * by;
+            c = splitter * negate;
+            ahi = c - (c - negate);
+            alo = negate - ahi;
+            c = splitter * by;
+            bhi = c - (c - by);
+            blo = by - bhi;
+            b[0] = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            b[1] = s1;
+            return 2;
+        } else {
+            s1 = xtail * ay;
+            c = splitter * xtail;
+            ahi = c - (c - xtail);
+            alo = xtail - ahi;
+            c = splitter * ay;
+            bhi = c - (c - ay);
+            blo = ay - bhi;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = ytail * ax;
+            c = splitter * ytail;
+            ahi = c - (c - ytail);
+            alo = ytail - ahi;
+            c = splitter * ax;
+            bhi = c - (c - ax);
+            blo = ax - bhi;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 - t0;
+            bvirt = s0 - _i;
+            a[0] = s0 - (_i + bvirt) + (bvirt - t0);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 - t1;
+            bvirt = _0 - _i;
+            a[1] = _0 - (_i + bvirt) + (bvirt - t1);
+            u3 = _j + _i;
+            bvirt = u3 - _j;
+            a[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            a[3] = u3;
+            s1 = ytail * bx;
+            c = splitter * ytail;
+            ahi = c - (c - ytail);
+            alo = ytail - ahi;
+            c = splitter * bx;
+            bhi = c - (c - bx);
+            blo = bx - bhi;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = xtail * by;
+            c = splitter * xtail;
+            ahi = c - (c - xtail);
+            alo = xtail - ahi;
+            c = splitter * by;
+            bhi = c - (c - by);
+            blo = by - bhi;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 - t0;
+            bvirt = s0 - _i;
+            b[0] = s0 - (_i + bvirt) + (bvirt - t0);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 - t1;
+            bvirt = _0 - _i;
+            b[1] = _0 - (_i + bvirt) + (bvirt - t1);
+            u3 = _j + _i;
+            bvirt = u3 - _j;
+            b[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            b[3] = u3;
+            return 4;
+        }
+    }
+}
+
+function tailadd(finlen, a, b, k, z) {
+    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _k, _0, s1, s0, u3;
+    s1 = a * b;
+    c = splitter * a;
+    ahi = c - (c - a);
+    alo = a - ahi;
+    c = splitter * b;
+    bhi = c - (c - b);
+    blo = b - bhi;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    c = splitter * k;
+    bhi = c - (c - k);
+    blo = k - bhi;
+    _i = s0 * k;
+    c = splitter * s0;
+    ahi = c - (c - s0);
+    alo = s0 - ahi;
+    u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
+    _j = s1 * k;
+    c = splitter * s1;
+    ahi = c - (c - s1);
+    alo = s1 - ahi;
+    _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
+    _k = _i + _0;
+    bvirt = _k - _i;
+    u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
+    u3 = _j + _k;
+    u$1[2] = _k - (u3 - _j);
+    u$1[3] = u3;
+    finlen = finadd$1(finlen, 4, u$1);
+    if (z !== 0) {
+        c = splitter * z;
+        bhi = c - (c - z);
+        blo = z - bhi;
+        _i = s0 * z;
+        c = splitter * s0;
+        ahi = c - (c - s0);
+        alo = s0 - ahi;
+        u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
+        _j = s1 * z;
+        c = splitter * s1;
+        ahi = c - (c - s1);
+        alo = s1 - ahi;
+        _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
+        _k = _i + _0;
+        bvirt = _k - _i;
+        u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
+        u3 = _j + _k;
+        u$1[2] = _k - (u3 - _j);
+        u$1[3] = u3;
+        finlen = finadd$1(finlen, 4, u$1);
+    }
+    return finlen;
+}
 
+function orient3dadapt(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz, permanent) {
+    let finlen;
     let adxtail, bdxtail, cdxtail;
     let adytail, bdytail, cdytail;
     let adztail, bdztail, cdztail;
-    let at_blarge, at_clarge;
-    let bt_clarge, bt_alarge;
-    let ct_alarge, ct_blarge;
-    let at_blen, at_clen, bt_clen, bt_alen, ct_alen, ct_blen;
-    let bdxt_cdy1, cdxt_bdy1, cdxt_ady1;
-    let adxt_cdy1, adxt_bdy1, bdxt_ady1;
-    let bdxt_cdy0, cdxt_bdy0, cdxt_ady0;
-    let adxt_cdy0, adxt_bdy0, bdxt_ady0;
-    let bdyt_cdx1, cdyt_bdx1, cdyt_adx1;
-    let adyt_cdx1, adyt_bdx1, bdyt_adx1;
-    let bdyt_cdx0, cdyt_bdx0, cdyt_adx0;
-    let adyt_cdx0, adyt_bdx0, bdyt_adx0;
-    let bdxt_cdyt1, cdxt_bdyt1, cdxt_adyt1;
-    let adxt_cdyt1, adxt_bdyt1, bdxt_adyt1;
-    let bdxt_cdyt0, cdxt_bdyt0, cdxt_adyt0;
-    let adxt_cdyt0, adxt_bdyt0, bdxt_adyt0;
-    let u3;
-    let vlength, wlength;
-    let negate;
-
-    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _k, _0;
+    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0, s1, s0, t1, t0, u3;
 
     const adx = ax - dx;
     const bdx = bx - dx;
@@ -402,103 +547,101 @@ function orient3dadapt(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz, permanent
     const bdz = bz - dz;
     const cdz = cz - dz;
 
-    bdxcdy1 = bdx * cdy;
+    s1 = bdx * cdy;
     c = splitter * bdx;
     ahi = c - (c - bdx);
     alo = bdx - ahi;
     c = splitter * cdy;
     bhi = c - (c - cdy);
     blo = cdy - bhi;
-    bdxcdy0 = alo * blo - (bdxcdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    cdxbdy1 = cdx * bdy;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = cdx * bdy;
     c = splitter * cdx;
     ahi = c - (c - cdx);
     alo = cdx - ahi;
     c = splitter * bdy;
     bhi = c - (c - bdy);
     blo = bdy - bhi;
-    cdxbdy0 = alo * blo - (cdxbdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = bdxcdy0 - cdxbdy0;
-    bvirt = bdxcdy0 - _i;
-    bc[0] = bdxcdy0 - (_i + bvirt) + (bvirt - cdxbdy0);
-    _j = bdxcdy1 + _i;
-    bvirt = _j - bdxcdy1;
-    _0 = bdxcdy1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - cdxbdy1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    bc$2[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 - cdxbdy1);
-    bc3 = _j + _i;
-    bvirt = bc3 - _j;
-    bc[2] = _j - (bc3 - bvirt) + (_i - bvirt);
-    bc[3] = bc3;
-    const alen = scale_expansion_zeroelim(4, bc, adz, adet);
-
-    cdxady1 = cdx * ady;
+    bc$2[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    bc$2[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    bc$2[3] = u3;
+    s1 = cdx * ady;
     c = splitter * cdx;
     ahi = c - (c - cdx);
     alo = cdx - ahi;
     c = splitter * ady;
     bhi = c - (c - ady);
     blo = ady - bhi;
-    cdxady0 = alo * blo - (cdxady1 - ahi * bhi - alo * bhi - ahi * blo);
-    adxcdy1 = adx * cdy;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = adx * cdy;
     c = splitter * adx;
     ahi = c - (c - adx);
     alo = adx - ahi;
     c = splitter * cdy;
     bhi = c - (c - cdy);
     blo = cdy - bhi;
-    adxcdy0 = alo * blo - (adxcdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = cdxady0 - adxcdy0;
-    bvirt = cdxady0 - _i;
-    ca[0] = cdxady0 - (_i + bvirt) + (bvirt - adxcdy0);
-    _j = cdxady1 + _i;
-    bvirt = _j - cdxady1;
-    _0 = cdxady1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - adxcdy1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ca$1[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    ca[1] = _0 - (_i + bvirt) + (bvirt - adxcdy1);
-    ca3 = _j + _i;
-    bvirt = ca3 - _j;
-    ca[2] = _j - (ca3 - bvirt) + (_i - bvirt);
-    ca[3] = ca3;
-    const blen = scale_expansion_zeroelim(4, ca, bdz, bdet);
-
-    adxbdy1 = adx * bdy;
+    ca$1[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    ca$1[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    ca$1[3] = u3;
+    s1 = adx * bdy;
     c = splitter * adx;
     ahi = c - (c - adx);
     alo = adx - ahi;
     c = splitter * bdy;
     bhi = c - (c - bdy);
     blo = bdy - bhi;
-    adxbdy0 = alo * blo - (adxbdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    bdxady1 = bdx * ady;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = bdx * ady;
     c = splitter * bdx;
     ahi = c - (c - bdx);
     alo = bdx - ahi;
     c = splitter * ady;
     bhi = c - (c - ady);
     blo = ady - bhi;
-    bdxady0 = alo * blo - (bdxady1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = adxbdy0 - bdxady0;
-    bvirt = adxbdy0 - _i;
-    ab[0] = adxbdy0 - (_i + bvirt) + (bvirt - bdxady0);
-    _j = adxbdy1 + _i;
-    bvirt = _j - adxbdy1;
-    _0 = adxbdy1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - bdxady1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ab$2[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 - bdxady1);
-    ab3 = _j + _i;
-    bvirt = ab3 - _j;
-    ab[2] = _j - (ab3 - bvirt) + (_i - bvirt);
-    ab[3] = ab3;
-    const clen = scale_expansion_zeroelim(4, ab, cdz, cdet);
+    ab$2[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    ab$2[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    ab$2[3] = u3;
 
-    const ablen = fast_expansion_sum_zeroelim(alen, adet, blen, bdet, abdet);
-    finlength = fast_expansion_sum_zeroelim(ablen, abdet, clen, cdet, fin1);
+    finlen = sum(
+        sum(
+            scale(4, bc$2, adz, _8$2), _8$2,
+            scale(4, ca$1, bdz, _8b$1), _8b$1, _16$2), _16$2,
+        scale(4, ab$2, cdz, _8$2), _8$2, fin$2);
 
-    let det = estimate(finlength, fin1);
+    let det = estimate(finlen, fin$2);
     let errbound = o3derrboundB * permanent;
     if (det >= errbound || -det >= errbound) {
         return det;
@@ -538,737 +681,58 @@ function orient3dadapt(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz, permanent
         return det;
     }
 
-    finnow = fin1;
-    finother = fin2;
-
-    if (adxtail === 0) {
-        if (adytail === 0) {
-            at_b[0] = 0;
-            at_blen = 1;
-            at_c[0] = 0;
-            at_clen = 1;
-        } else {
-            negate = -adytail;
-            at_blarge = negate * bdx;
-            c = splitter * negate;
-            ahi = c - (c - negate);
-            alo = negate - ahi;
-            c = splitter * bdx;
-            bhi = c - (c - bdx);
-            blo = bdx - bhi;
-            at_b[0] = alo * blo - (at_blarge - ahi * bhi - alo * bhi - ahi * blo);
-            at_b[1] = at_blarge;
-            at_blen = 2;
-            at_clarge = adytail * cdx;
-            c = splitter * adytail;
-            ahi = c - (c - adytail);
-            alo = adytail - ahi;
-            c = splitter * cdx;
-            bhi = c - (c - cdx);
-            blo = cdx - bhi;
-            at_c[0] = alo * blo - (at_clarge - ahi * bhi - alo * bhi - ahi * blo);
-            at_c[1] = at_clarge;
-            at_clen = 2;
-        }
-    } else {
-        if (adytail === 0) {
-            at_blarge = adxtail * bdy;
-            c = splitter * adxtail;
-            ahi = c - (c - adxtail);
-            alo = adxtail - ahi;
-            c = splitter * bdy;
-            bhi = c - (c - bdy);
-            blo = bdy - bhi;
-            at_b[0] = alo * blo - (at_blarge - ahi * bhi - alo * bhi - ahi * blo);
-            at_b[1] = at_blarge;
-            at_blen = 2;
-            negate = -adxtail;
-            at_clarge = negate * cdy;
-            c = splitter * negate;
-            ahi = c - (c - negate);
-            alo = negate - ahi;
-            c = splitter * cdy;
-            bhi = c - (c - cdy);
-            blo = cdy - bhi;
-            at_c[0] = alo * blo - (at_clarge - ahi * bhi - alo * bhi - ahi * blo);
-            at_c[1] = at_clarge;
-            at_clen = 2;
-        } else {
-            adxt_bdy1 = adxtail * bdy;
-            c = splitter * adxtail;
-            ahi = c - (c - adxtail);
-            alo = adxtail - ahi;
-            c = splitter * bdy;
-            bhi = c - (c - bdy);
-            blo = bdy - bhi;
-            adxt_bdy0 = alo * blo - (adxt_bdy1 - ahi * bhi - alo * bhi - ahi * blo);
-            adyt_bdx1 = adytail * bdx;
-            c = splitter * adytail;
-            ahi = c - (c - adytail);
-            alo = adytail - ahi;
-            c = splitter * bdx;
-            bhi = c - (c - bdx);
-            blo = bdx - bhi;
-            adyt_bdx0 = alo * blo - (adyt_bdx1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = adxt_bdy0 - adyt_bdx0;
-            bvirt = adxt_bdy0 - _i;
-            at_b[0] = adxt_bdy0 - (_i + bvirt) + (bvirt - adyt_bdx0);
-            _j = adxt_bdy1 + _i;
-            bvirt = _j - adxt_bdy1;
-            _0 = adxt_bdy1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 - adyt_bdx1;
-            bvirt = _0 - _i;
-            at_b[1] = _0 - (_i + bvirt) + (bvirt - adyt_bdx1);
-            at_blarge = _j + _i;
-            bvirt = at_blarge - _j;
-            at_b[2] = _j - (at_blarge - bvirt) + (_i - bvirt);
-            at_b[3] = at_blarge;
-            at_blen = 4;
-            adyt_cdx1 = adytail * cdx;
-            c = splitter * adytail;
-            ahi = c - (c - adytail);
-            alo = adytail - ahi;
-            c = splitter * cdx;
-            bhi = c - (c - cdx);
-            blo = cdx - bhi;
-            adyt_cdx0 = alo * blo - (adyt_cdx1 - ahi * bhi - alo * bhi - ahi * blo);
-            adxt_cdy1 = adxtail * cdy;
-            c = splitter * adxtail;
-            ahi = c - (c - adxtail);
-            alo = adxtail - ahi;
-            c = splitter * cdy;
-            bhi = c - (c - cdy);
-            blo = cdy - bhi;
-            adxt_cdy0 = alo * blo - (adxt_cdy1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = adyt_cdx0 - adxt_cdy0;
-            bvirt = adyt_cdx0 - _i;
-            at_c[0] = adyt_cdx0 - (_i + bvirt) + (bvirt - adxt_cdy0);
-            _j = adyt_cdx1 + _i;
-            bvirt = _j - adyt_cdx1;
-            _0 = adyt_cdx1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 - adxt_cdy1;
-            bvirt = _0 - _i;
-            at_c[1] = _0 - (_i + bvirt) + (bvirt - adxt_cdy1);
-            at_clarge = _j + _i;
-            bvirt = at_clarge - _j;
-            at_c[2] = _j - (at_clarge - bvirt) + (_i - bvirt);
-            at_c[3] = at_clarge;
-            at_clen = 4;
-        }
-    }
-    if (bdxtail === 0) {
-        if (bdytail === 0) {
-            bt_c[0] = 0;
-            bt_clen = 1;
-            bt_a[0] = 0;
-            bt_alen = 1;
-        } else {
-            negate = -bdytail;
-            bt_clarge = negate * cdx;
-            c = splitter * negate;
-            ahi = c - (c - negate);
-            alo = negate - ahi;
-            c = splitter * cdx;
-            bhi = c - (c - cdx);
-            blo = cdx - bhi;
-            bt_c[0] = alo * blo - (bt_clarge - ahi * bhi - alo * bhi - ahi * blo);
-            bt_c[1] = bt_clarge;
-            bt_clen = 2;
-            bt_alarge = bdytail * adx;
-            c = splitter * bdytail;
-            ahi = c - (c - bdytail);
-            alo = bdytail - ahi;
-            c = splitter * adx;
-            bhi = c - (c - adx);
-            blo = adx - bhi;
-            bt_a[0] = alo * blo - (bt_alarge - ahi * bhi - alo * bhi - ahi * blo);
-            bt_a[1] = bt_alarge;
-            bt_alen = 2;
-        }
-    } else {
-        if (bdytail === 0) {
-            bt_clarge = bdxtail * cdy;
-            c = splitter * bdxtail;
-            ahi = c - (c - bdxtail);
-            alo = bdxtail - ahi;
-            c = splitter * cdy;
-            bhi = c - (c - cdy);
-            blo = cdy - bhi;
-            bt_c[0] = alo * blo - (bt_clarge - ahi * bhi - alo * bhi - ahi * blo);
-            bt_c[1] = bt_clarge;
-            bt_clen = 2;
-            negate = -bdxtail;
-            bt_alarge = negate * ady;
-            c = splitter * negate;
-            ahi = c - (c - negate);
-            alo = negate - ahi;
-            c = splitter * ady;
-            bhi = c - (c - ady);
-            blo = ady - bhi;
-            bt_a[0] = alo * blo - (bt_alarge - ahi * bhi - alo * bhi - ahi * blo);
-            bt_a[1] = bt_alarge;
-            bt_alen = 2;
-        } else {
-            bdxt_cdy1 = bdxtail * cdy;
-            c = splitter * bdxtail;
-            ahi = c - (c - bdxtail);
-            alo = bdxtail - ahi;
-            c = splitter * cdy;
-            bhi = c - (c - cdy);
-            blo = cdy - bhi;
-            bdxt_cdy0 = alo * blo - (bdxt_cdy1 - ahi * bhi - alo * bhi - ahi * blo);
-            bdyt_cdx1 = bdytail * cdx;
-            c = splitter * bdytail;
-            ahi = c - (c - bdytail);
-            alo = bdytail - ahi;
-            c = splitter * cdx;
-            bhi = c - (c - cdx);
-            blo = cdx - bhi;
-            bdyt_cdx0 = alo * blo - (bdyt_cdx1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = bdxt_cdy0 - bdyt_cdx0;
-            bvirt = bdxt_cdy0 - _i;
-            bt_c[0] = bdxt_cdy0 - (_i + bvirt) + (bvirt - bdyt_cdx0);
-            _j = bdxt_cdy1 + _i;
-            bvirt = _j - bdxt_cdy1;
-            _0 = bdxt_cdy1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 - bdyt_cdx1;
-            bvirt = _0 - _i;
-            bt_c[1] = _0 - (_i + bvirt) + (bvirt - bdyt_cdx1);
-            bt_clarge = _j + _i;
-            bvirt = bt_clarge - _j;
-            bt_c[2] = _j - (bt_clarge - bvirt) + (_i - bvirt);
-            bt_c[3] = bt_clarge;
-            bt_clen = 4;
-            bdyt_adx1 = bdytail * adx;
-            c = splitter * bdytail;
-            ahi = c - (c - bdytail);
-            alo = bdytail - ahi;
-            c = splitter * adx;
-            bhi = c - (c - adx);
-            blo = adx - bhi;
-            bdyt_adx0 = alo * blo - (bdyt_adx1 - ahi * bhi - alo * bhi - ahi * blo);
-            bdxt_ady1 = bdxtail * ady;
-            c = splitter * bdxtail;
-            ahi = c - (c - bdxtail);
-            alo = bdxtail - ahi;
-            c = splitter * ady;
-            bhi = c - (c - ady);
-            blo = ady - bhi;
-            bdxt_ady0 = alo * blo - (bdxt_ady1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = bdyt_adx0 - bdxt_ady0;
-            bvirt = bdyt_adx0 - _i;
-            bt_a[0] = bdyt_adx0 - (_i + bvirt) + (bvirt - bdxt_ady0);
-            _j = bdyt_adx1 + _i;
-            bvirt = _j - bdyt_adx1;
-            _0 = bdyt_adx1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 - bdxt_ady1;
-            bvirt = _0 - _i;
-            bt_a[1] = _0 - (_i + bvirt) + (bvirt - bdxt_ady1);
-            bt_alarge = _j + _i;
-            bvirt = bt_alarge - _j;
-            bt_a[2] = _j - (bt_alarge - bvirt) + (_i - bvirt);
-            bt_a[3] = bt_alarge;
-            bt_alen = 4;
-        }
-    }
-    if (cdxtail === 0) {
-        if (cdytail === 0) {
-            ct_a[0] = 0;
-            ct_alen = 1;
-            ct_b[0] = 0;
-            ct_blen = 1;
-        } else {
-            negate = -cdytail;
-            ct_alarge = negate * adx;
-            c = splitter * negate;
-            ahi = c - (c - negate);
-            alo = negate - ahi;
-            c = splitter * adx;
-            bhi = c - (c - adx);
-            blo = adx - bhi;
-            ct_a[0] = alo * blo - (ct_alarge - ahi * bhi - alo * bhi - ahi * blo);
-            ct_a[1] = ct_alarge;
-            ct_alen = 2;
-            ct_blarge = cdytail * bdx;
-            c = splitter * cdytail;
-            ahi = c - (c - cdytail);
-            alo = cdytail - ahi;
-            c = splitter * bdx;
-            bhi = c - (c - bdx);
-            blo = bdx - bhi;
-            ct_b[0] = alo * blo - (ct_blarge - ahi * bhi - alo * bhi - ahi * blo);
-            ct_b[1] = ct_blarge;
-            ct_blen = 2;
-        }
-    } else {
-        if (cdytail === 0) {
-            ct_alarge = cdxtail * ady;
-            c = splitter * cdxtail;
-            ahi = c - (c - cdxtail);
-            alo = cdxtail - ahi;
-            c = splitter * ady;
-            bhi = c - (c - ady);
-            blo = ady - bhi;
-            ct_a[0] = alo * blo - (ct_alarge - ahi * bhi - alo * bhi - ahi * blo);
-            ct_a[1] = ct_alarge;
-            ct_alen = 2;
-            negate = -cdxtail;
-            ct_blarge = negate * bdy;
-            c = splitter * negate;
-            ahi = c - (c - negate);
-            alo = negate - ahi;
-            c = splitter * bdy;
-            bhi = c - (c - bdy);
-            blo = bdy - bhi;
-            ct_b[0] = alo * blo - (ct_blarge - ahi * bhi - alo * bhi - ahi * blo);
-            ct_b[1] = ct_blarge;
-            ct_blen = 2;
-        } else {
-            cdxt_ady1 = cdxtail * ady;
-            c = splitter * cdxtail;
-            ahi = c - (c - cdxtail);
-            alo = cdxtail - ahi;
-            c = splitter * ady;
-            bhi = c - (c - ady);
-            blo = ady - bhi;
-            cdxt_ady0 = alo * blo - (cdxt_ady1 - ahi * bhi - alo * bhi - ahi * blo);
-            cdyt_adx1 = cdytail * adx;
-            c = splitter * cdytail;
-            ahi = c - (c - cdytail);
-            alo = cdytail - ahi;
-            c = splitter * adx;
-            bhi = c - (c - adx);
-            blo = adx - bhi;
-            cdyt_adx0 = alo * blo - (cdyt_adx1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = cdxt_ady0 - cdyt_adx0;
-            bvirt = cdxt_ady0 - _i;
-            ct_a[0] = cdxt_ady0 - (_i + bvirt) + (bvirt - cdyt_adx0);
-            _j = cdxt_ady1 + _i;
-            bvirt = _j - cdxt_ady1;
-            _0 = cdxt_ady1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 - cdyt_adx1;
-            bvirt = _0 - _i;
-            ct_a[1] = _0 - (_i + bvirt) + (bvirt - cdyt_adx1);
-            ct_alarge = _j + _i;
-            bvirt = ct_alarge - _j;
-            ct_a[2] = _j - (ct_alarge - bvirt) + (_i - bvirt);
-            ct_a[3] = ct_alarge;
-            ct_alen = 4;
-            cdyt_bdx1 = cdytail * bdx;
-            c = splitter * cdytail;
-            ahi = c - (c - cdytail);
-            alo = cdytail - ahi;
-            c = splitter * bdx;
-            bhi = c - (c - bdx);
-            blo = bdx - bhi;
-            cdyt_bdx0 = alo * blo - (cdyt_bdx1 - ahi * bhi - alo * bhi - ahi * blo);
-            cdxt_bdy1 = cdxtail * bdy;
-            c = splitter * cdxtail;
-            ahi = c - (c - cdxtail);
-            alo = cdxtail - ahi;
-            c = splitter * bdy;
-            bhi = c - (c - bdy);
-            blo = bdy - bhi;
-            cdxt_bdy0 = alo * blo - (cdxt_bdy1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = cdyt_bdx0 - cdxt_bdy0;
-            bvirt = cdyt_bdx0 - _i;
-            ct_b[0] = cdyt_bdx0 - (_i + bvirt) + (bvirt - cdxt_bdy0);
-            _j = cdyt_bdx1 + _i;
-            bvirt = _j - cdyt_bdx1;
-            _0 = cdyt_bdx1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 - cdxt_bdy1;
-            bvirt = _0 - _i;
-            ct_b[1] = _0 - (_i + bvirt) + (bvirt - cdxt_bdy1);
-            ct_blarge = _j + _i;
-            bvirt = ct_blarge - _j;
-            ct_b[2] = _j - (ct_blarge - bvirt) + (_i - bvirt);
-            ct_b[3] = ct_blarge;
-            ct_blen = 4;
-        }
-    }
+    const at_len = tailinit(adxtail, adytail, bdx, bdy, cdx, cdy, at_b, at_c);
+    const bt_len = tailinit(bdxtail, bdytail, cdx, cdy, adx, ady, bt_c, bt_a);
+    const ct_len = tailinit(cdxtail, cdytail, adx, ady, bdx, bdy, ct_a, ct_b);
 
-    const bctlen = fast_expansion_sum_zeroelim(bt_clen, bt_c, ct_blen, ct_b, bct);
-    wlength = scale_expansion_zeroelim(bctlen, bct, adz, w);
-    finlength = fast_expansion_sum_zeroelim(finlength, finnow, wlength, w, finother);
-    finswap = finnow; finnow = finother; finother = finswap;
+    const bctlen = sum(bt_len, bt_c, ct_len, ct_b, bct$1);
+    finlen = finadd$1(finlen, scale(bctlen, bct$1, adz, _16$2), _16$2);
 
-    const catlen = fast_expansion_sum_zeroelim(ct_alen, ct_a, at_clen, at_c, cat);
-    wlength = scale_expansion_zeroelim(catlen, cat, bdz, w);
-    finlength = fast_expansion_sum_zeroelim(finlength, finnow, wlength, w, finother);
-    finswap = finnow; finnow = finother; finother = finswap;
+    const catlen = sum(ct_len, ct_a, at_len, at_c, cat$1);
+    finlen = finadd$1(finlen, scale(catlen, cat$1, bdz, _16$2), _16$2);
 
-    const abtlen = fast_expansion_sum_zeroelim(at_blen, at_b, bt_alen, bt_a, abt);
-    wlength = scale_expansion_zeroelim(abtlen, abt, cdz, w);
-    finlength = fast_expansion_sum_zeroelim(finlength, finnow, wlength, w, finother);
-    finswap = finnow; finnow = finother; finother = finswap;
+    const abtlen = sum(at_len, at_b, bt_len, bt_a, abt$1);
+    finlen = finadd$1(finlen, scale(abtlen, abt$1, cdz, _16$2), _16$2);
 
     if (adztail !== 0) {
-        vlength = scale_expansion_zeroelim(4, bc, adztail, v);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, vlength, v, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
+        finlen = finadd$1(finlen, scale(4, bc$2, adztail, _12), _12);
+        finlen = finadd$1(finlen, scale(bctlen, bct$1, adztail, _16$2), _16$2);
     }
     if (bdztail !== 0) {
-        vlength = scale_expansion_zeroelim(4, ca, bdztail, v);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, vlength, v, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
+        finlen = finadd$1(finlen, scale(4, ca$1, bdztail, _12), _12);
+        finlen = finadd$1(finlen, scale(catlen, cat$1, bdztail, _16$2), _16$2);
     }
     if (cdztail !== 0) {
-        vlength = scale_expansion_zeroelim(4, ab, cdztail, v);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, vlength, v, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
+        finlen = finadd$1(finlen, scale(4, ab$2, cdztail, _12), _12);
+        finlen = finadd$1(finlen, scale(abtlen, abt$1, cdztail, _16$2), _16$2);
     }
 
     if (adxtail !== 0) {
         if (bdytail !== 0) {
-            adxt_bdyt1 = adxtail * bdytail;
-            c = splitter * adxtail;
-            ahi = c - (c - adxtail);
-            alo = adxtail - ahi;
-            c = splitter * bdytail;
-            bhi = c - (c - bdytail);
-            blo = bdytail - bhi;
-            adxt_bdyt0 = alo * blo - (adxt_bdyt1 - ahi * bhi - alo * bhi - ahi * blo);
-            c = splitter * cdz;
-            bhi = c - (c - cdz);
-            blo = cdz - bhi;
-            _i = adxt_bdyt0 * cdz;
-            c = splitter * adxt_bdyt0;
-            ahi = c - (c - adxt_bdyt0);
-            alo = adxt_bdyt0 - ahi;
-            u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-            _j = adxt_bdyt1 * cdz;
-            c = splitter * adxt_bdyt1;
-            ahi = c - (c - adxt_bdyt1);
-            alo = adxt_bdyt1 - ahi;
-            _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-            _k = _i + _0;
-            bvirt = _k - _i;
-            u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-            u3 = _j + _k;
-            u$1[2] = _k - (u3 - _j);
-            u$1[3] = u3;
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
-            if (cdztail !== 0) {
-                c = splitter * cdztail;
-                bhi = c - (c - cdztail);
-                blo = cdztail - bhi;
-                _i = adxt_bdyt0 * cdztail;
-                c = splitter * adxt_bdyt0;
-                ahi = c - (c - adxt_bdyt0);
-                alo = adxt_bdyt0 - ahi;
-                u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-                _j = adxt_bdyt1 * cdztail;
-                c = splitter * adxt_bdyt1;
-                ahi = c - (c - adxt_bdyt1);
-                alo = adxt_bdyt1 - ahi;
-                _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-                _k = _i + _0;
-                bvirt = _k - _i;
-                u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-                u3 = _j + _k;
-                u$1[2] = _k - (u3 - _j);
-                u$1[3] = u3;
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
-            }
+            finlen = tailadd(finlen, adxtail, bdytail, cdz, cdztail);
         }
         if (cdytail !== 0) {
-            negate = -adxtail;
-            adxt_cdyt1 = negate * cdytail;
-            c = splitter * negate;
-            ahi = c - (c - negate);
-            alo = negate - ahi;
-            c = splitter * cdytail;
-            bhi = c - (c - cdytail);
-            blo = cdytail - bhi;
-            adxt_cdyt0 = alo * blo - (adxt_cdyt1 - ahi * bhi - alo * bhi - ahi * blo);
-            c = splitter * bdz;
-            bhi = c - (c - bdz);
-            blo = bdz - bhi;
-            _i = adxt_cdyt0 * bdz;
-            c = splitter * adxt_cdyt0;
-            ahi = c - (c - adxt_cdyt0);
-            alo = adxt_cdyt0 - ahi;
-            u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-            _j = adxt_cdyt1 * bdz;
-            c = splitter * adxt_cdyt1;
-            ahi = c - (c - adxt_cdyt1);
-            alo = adxt_cdyt1 - ahi;
-            _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-            _k = _i + _0;
-            bvirt = _k - _i;
-            u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-            u3 = _j + _k;
-            u$1[2] = _k - (u3 - _j);
-            u$1[3] = u3;
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
-            if (bdztail !== 0) {
-                c = splitter * bdztail;
-                bhi = c - (c - bdztail);
-                blo = bdztail - bhi;
-                _i = adxt_cdyt0 * bdztail;
-                c = splitter * adxt_cdyt0;
-                ahi = c - (c - adxt_cdyt0);
-                alo = adxt_cdyt0 - ahi;
-                u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-                _j = adxt_cdyt1 * bdztail;
-                c = splitter * adxt_cdyt1;
-                ahi = c - (c - adxt_cdyt1);
-                alo = adxt_cdyt1 - ahi;
-                _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-                _k = _i + _0;
-                bvirt = _k - _i;
-                u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-                u3 = _j + _k;
-                u$1[2] = _k - (u3 - _j);
-                u$1[3] = u3;
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
-            }
+            finlen = tailadd(finlen, -adxtail, cdytail, bdz, bdztail);
         }
     }
     if (bdxtail !== 0) {
         if (cdytail !== 0) {
-            bdxt_cdyt1 = bdxtail * cdytail;
-            c = splitter * bdxtail;
-            ahi = c - (c - bdxtail);
-            alo = bdxtail - ahi;
-            c = splitter * cdytail;
-            bhi = c - (c - cdytail);
-            blo = cdytail - bhi;
-            bdxt_cdyt0 = alo * blo - (bdxt_cdyt1 - ahi * bhi - alo * bhi - ahi * blo);
-            c = splitter * adz;
-            bhi = c - (c - adz);
-            blo = adz - bhi;
-            _i = bdxt_cdyt0 * adz;
-            c = splitter * bdxt_cdyt0;
-            ahi = c - (c - bdxt_cdyt0);
-            alo = bdxt_cdyt0 - ahi;
-            u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-            _j = bdxt_cdyt1 * adz;
-            c = splitter * bdxt_cdyt1;
-            ahi = c - (c - bdxt_cdyt1);
-            alo = bdxt_cdyt1 - ahi;
-            _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-            _k = _i + _0;
-            bvirt = _k - _i;
-            u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-            u3 = _j + _k;
-            u$1[2] = _k - (u3 - _j);
-            u$1[3] = u3;
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
-            if (adztail !== 0) {
-                c = splitter * adztail;
-                bhi = c - (c - adztail);
-                blo = adztail - bhi;
-                _i = bdxt_cdyt0 * adztail;
-                c = splitter * bdxt_cdyt0;
-                ahi = c - (c - bdxt_cdyt0);
-                alo = bdxt_cdyt0 - ahi;
-                u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-                _j = bdxt_cdyt1 * adztail;
-                c = splitter * bdxt_cdyt1;
-                ahi = c - (c - bdxt_cdyt1);
-                alo = bdxt_cdyt1 - ahi;
-                _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-                _k = _i + _0;
-                bvirt = _k - _i;
-                u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-                u3 = _j + _k;
-                u$1[2] = _k - (u3 - _j);
-                u$1[3] = u3;
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
-            }
+            finlen = tailadd(finlen, bdxtail, cdytail, adz, adztail);
         }
         if (adytail !== 0) {
-            negate = -bdxtail;
-            bdxt_adyt1 = negate * adytail;
-            c = splitter * negate;
-            ahi = c - (c - negate);
-            alo = negate - ahi;
-            c = splitter * adytail;
-            bhi = c - (c - adytail);
-            blo = adytail - bhi;
-            bdxt_adyt0 = alo * blo - (bdxt_adyt1 - ahi * bhi - alo * bhi - ahi * blo);
-            c = splitter * cdz;
-            bhi = c - (c - cdz);
-            blo = cdz - bhi;
-            _i = bdxt_adyt0 * cdz;
-            c = splitter * bdxt_adyt0;
-            ahi = c - (c - bdxt_adyt0);
-            alo = bdxt_adyt0 - ahi;
-            u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-            _j = bdxt_adyt1 * cdz;
-            c = splitter * bdxt_adyt1;
-            ahi = c - (c - bdxt_adyt1);
-            alo = bdxt_adyt1 - ahi;
-            _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-            _k = _i + _0;
-            bvirt = _k - _i;
-            u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-            u3 = _j + _k;
-            u$1[2] = _k - (u3 - _j);
-            u$1[3] = u3;
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
-            if (cdztail !== 0) {
-                c = splitter * cdztail;
-                bhi = c - (c - cdztail);
-                blo = cdztail - bhi;
-                _i = bdxt_adyt0 * cdztail;
-                c = splitter * bdxt_adyt0;
-                ahi = c - (c - bdxt_adyt0);
-                alo = bdxt_adyt0 - ahi;
-                u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-                _j = bdxt_adyt1 * cdztail;
-                c = splitter * bdxt_adyt1;
-                ahi = c - (c - bdxt_adyt1);
-                alo = bdxt_adyt1 - ahi;
-                _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-                _k = _i + _0;
-                bvirt = _k - _i;
-                u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-                u3 = _j + _k;
-                u$1[2] = _k - (u3 - _j);
-                u$1[3] = u3;
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
-            }
+            finlen = tailadd(finlen, -bdxtail, adytail, cdz, cdztail);
         }
     }
     if (cdxtail !== 0) {
         if (adytail !== 0) {
-            cdxt_adyt1 = cdxtail * adytail;
-            c = splitter * cdxtail;
-            ahi = c - (c - cdxtail);
-            alo = cdxtail - ahi;
-            c = splitter * adytail;
-            bhi = c - (c - adytail);
-            blo = adytail - bhi;
-            cdxt_adyt0 = alo * blo - (cdxt_adyt1 - ahi * bhi - alo * bhi - ahi * blo);
-            c = splitter * bdz;
-            bhi = c - (c - bdz);
-            blo = bdz - bhi;
-            _i = cdxt_adyt0 * bdz;
-            c = splitter * cdxt_adyt0;
-            ahi = c - (c - cdxt_adyt0);
-            alo = cdxt_adyt0 - ahi;
-            u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-            _j = cdxt_adyt1 * bdz;
-            c = splitter * cdxt_adyt1;
-            ahi = c - (c - cdxt_adyt1);
-            alo = cdxt_adyt1 - ahi;
-            _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-            _k = _i + _0;
-            bvirt = _k - _i;
-            u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-            u3 = _j + _k;
-            u$1[2] = _k - (u3 - _j);
-            u$1[3] = u3;
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
-            if (bdztail !== 0) {
-                c = splitter * bdztail;
-                bhi = c - (c - bdztail);
-                blo = bdztail - bhi;
-                _i = cdxt_adyt0 * bdztail;
-                c = splitter * cdxt_adyt0;
-                ahi = c - (c - cdxt_adyt0);
-                alo = cdxt_adyt0 - ahi;
-                u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-                _j = cdxt_adyt1 * bdztail;
-                c = splitter * cdxt_adyt1;
-                ahi = c - (c - cdxt_adyt1);
-                alo = cdxt_adyt1 - ahi;
-                _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-                _k = _i + _0;
-                bvirt = _k - _i;
-                u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-                u3 = _j + _k;
-                u$1[2] = _k - (u3 - _j);
-                u$1[3] = u3;
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
-            }
+            finlen = tailadd(finlen, cdxtail, adytail, bdz, bdztail);
         }
         if (bdytail !== 0) {
-            negate = -cdxtail;
-            cdxt_bdyt1 = negate * bdytail;
-            c = splitter * negate;
-            ahi = c - (c - negate);
-            alo = negate - ahi;
-            c = splitter * bdytail;
-            bhi = c - (c - bdytail);
-            blo = bdytail - bhi;
-            cdxt_bdyt0 = alo * blo - (cdxt_bdyt1 - ahi * bhi - alo * bhi - ahi * blo);
-            c = splitter * adz;
-            bhi = c - (c - adz);
-            blo = adz - bhi;
-            _i = cdxt_bdyt0 * adz;
-            c = splitter * cdxt_bdyt0;
-            ahi = c - (c - cdxt_bdyt0);
-            alo = cdxt_bdyt0 - ahi;
-            u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-            _j = cdxt_bdyt1 * adz;
-            c = splitter * cdxt_bdyt1;
-            ahi = c - (c - cdxt_bdyt1);
-            alo = cdxt_bdyt1 - ahi;
-            _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-            _k = _i + _0;
-            bvirt = _k - _i;
-            u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-            u3 = _j + _k;
-            u$1[2] = _k - (u3 - _j);
-            u$1[3] = u3;
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
-            if (adztail !== 0) {
-                c = splitter * adztail;
-                bhi = c - (c - adztail);
-                blo = adztail - bhi;
-                _i = cdxt_bdyt0 * adztail;
-                c = splitter * cdxt_bdyt0;
-                ahi = c - (c - cdxt_bdyt0);
-                alo = cdxt_bdyt0 - ahi;
-                u$1[0] = alo * blo - (_i - ahi * bhi - alo * bhi - ahi * blo);
-                _j = cdxt_bdyt1 * adztail;
-                c = splitter * cdxt_bdyt1;
-                ahi = c - (c - cdxt_bdyt1);
-                alo = cdxt_bdyt1 - ahi;
-                _0 = alo * blo - (_j - ahi * bhi - alo * bhi - ahi * blo);
-                _k = _i + _0;
-                bvirt = _k - _i;
-                u$1[1] = _i - (_k - bvirt) + (_0 - bvirt);
-                u3 = _j + _k;
-                u$1[2] = _k - (u3 - _j);
-                u$1[3] = u3;
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, 4, u$1, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
-            }
+            finlen = tailadd(finlen, -cdxtail, bdytail, adz, adztail);
         }
     }
 
-    if (adztail !== 0) {
-        wlength = scale_expansion_zeroelim(bctlen, bct, adztail, w);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, wlength, w, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
-    }
-    if (bdztail !== 0) {
-        wlength = scale_expansion_zeroelim(catlen, cat, bdztail, w);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, wlength, w, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
-    }
-    if (cdztail !== 0) {
-        wlength = scale_expansion_zeroelim(abtlen, abt, cdztail, w);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, wlength, w, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
-    }
-
-    return finnow[finlength - 1];
+    return fin$2[finlen - 1];
 }
 
 function orient3d(ax, ay, az, bx, by, bz, cx, cy, cz, dx, dy, dz) {
@@ -1330,98 +794,53 @@ const iccerrboundB = (4 + 48 * epsilon) * epsilon;
 const iccerrboundC = (44 + 576 * epsilon) * epsilon * epsilon;
 
 const bc$1 = vec(4);
-const ca$1 = vec(4);
+const ca = vec(4);
 const ab$1 = vec(4);
-const axbc = vec(8);
-const axxbc = vec(16);
-const aybc = vec(8);
-const ayybc = vec(16);
-const adet$1 = vec(32);
-const bxca = vec(8);
-const bxxca = vec(16);
-const byca = vec(8);
-const byyca = vec(16);
-const bdet$1 = vec(32);
-const cxab = vec(8);
-const cxxab = vec(16);
-const cyab = vec(8);
-const cyyab = vec(16);
-const cdet$1 = vec(32);
-const abdet$1 = vec(64);
-const fin1$1 = vec(1152);
-const fin2$1 = vec(1152);
 const aa = vec(4);
 const bb = vec(4);
 const cc = vec(4);
-const u$2 = vec(4);
-const v$1 = vec(4);
-const temp8 = vec(8);
-const temp16a = vec(16);
-const temp16b = vec(16);
-const temp16c = vec(16);
-const temp32a = vec(32);
-const temp32b = vec(32);
-const temp48 = vec(48);
-const temp64 = vec(64);
-const axtbb = vec(8);
-const axtcc = vec(8);
-const aytbb = vec(8);
-const aytcc = vec(8);
-const bxtaa = vec(8);
-const bxtcc = vec(8);
-const bytaa = vec(8);
-const bytcc = vec(8);
-const cxtaa = vec(8);
-const cxtbb = vec(8);
-const cytaa = vec(8);
-const cytbb = vec(8);
+const u = vec(4);
+const v = vec(4);
 const axtbc = vec(8);
 const aytbc = vec(8);
 const bxtca = vec(8);
 const bytca = vec(8);
 const cxtab = vec(8);
 const cytab = vec(8);
-const axtbct = vec(16);
-const aytbct = vec(16);
-const bxtcat = vec(16);
-const bytcat = vec(16);
-const cxtabt = vec(16);
-const cytabt = vec(16);
-const axtbctt = vec(8);
-const aytbctt = vec(8);
-const bxtcatt = vec(8);
-const bytcatt = vec(8);
-const cxtabtt = vec(8);
-const cytabtt = vec(8);
-const abt$1 = vec(8);
-const bct$1 = vec(8);
-const cat$1 = vec(8);
+const abt = vec(8);
+const bct = vec(8);
+const cat = vec(8);
 const abtt = vec(4);
 const bctt = vec(4);
 const catt = vec(4);
 
-function incircleadapt(ax, ay, bx, by, cx, cy, dx, dy, permanent) {
-    let bdxcdy1, cdxbdy1, cdxady1, adxcdy1, adxbdy1, bdxady1;
-    let bdxcdy0, cdxbdy0, cdxady0, adxcdy0, adxbdy0, bdxady0;
-    let bc3, ca3, ab3;
-    let finnow, finother, finswap, finlength;
+const _8$1 = vec(8);
+const _16$1 = vec(16);
+const _16b = vec(16);
+const _16c = vec(16);
+const _32 = vec(32);
+const _32b = vec(32);
+const _48$1 = vec(48);
+const _64 = vec(64);
+
+let fin$1 = vec(1152);
+let fin2 = vec(1152);
+
+function finadd(finlen, a, alen) {
+    finlen = sum(finlen, fin$1, a, alen, fin2);
+    const tmp = fin$1; fin$1 = fin2; fin2 = tmp;
+    return finlen;
+}
 
+function incircleadapt(ax, ay, bx, by, cx, cy, dx, dy, permanent) {
+    let finlen;
     let adxtail, bdxtail, cdxtail, adytail, bdytail, cdytail;
-    let adxadx1, adyady1, bdxbdx1, bdybdy1, cdxcdx1, cdycdy1;
-    let adxadx0, adyady0, bdxbdx0, bdybdy0, cdxcdx0, cdycdy0;
-    let aa3, bb3, cc3;
-    let ti1, tj1;
-    let ti0, tj0;
-    let u3, v3;
-    let temp8len, temp16alen, temp16blen, temp16clen;
-    let temp32alen, temp32blen, temp48len, temp64len;
     let axtbclen, aytbclen, bxtcalen, bytcalen, cxtablen, cytablen;
     let abtlen, bctlen, catlen;
     let abttlen, bcttlen, cattlen;
-    let abtt3, bctt3, catt3;
-    let negate;
+    let n1, n0;
 
-    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0;
+    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0, s1, s0, t1, t0, u3;
 
     const adx = ax - dx;
     const bdx = bx - dx;
@@ -1430,115 +849,107 @@ function incircleadapt(ax, ay, bx, by, cx, cy, dx, dy, permanent) {
     const bdy = by - dy;
     const cdy = cy - dy;
 
-    bdxcdy1 = bdx * cdy;
+    s1 = bdx * cdy;
     c = splitter * bdx;
     ahi = c - (c - bdx);
     alo = bdx - ahi;
     c = splitter * cdy;
     bhi = c - (c - cdy);
     blo = cdy - bhi;
-    bdxcdy0 = alo * blo - (bdxcdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    cdxbdy1 = cdx * bdy;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = cdx * bdy;
     c = splitter * cdx;
     ahi = c - (c - cdx);
     alo = cdx - ahi;
     c = splitter * bdy;
     bhi = c - (c - bdy);
     blo = bdy - bhi;
-    cdxbdy0 = alo * blo - (cdxbdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = bdxcdy0 - cdxbdy0;
-    bvirt = bdxcdy0 - _i;
-    bc$1[0] = bdxcdy0 - (_i + bvirt) + (bvirt - cdxbdy0);
-    _j = bdxcdy1 + _i;
-    bvirt = _j - bdxcdy1;
-    _0 = bdxcdy1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - cdxbdy1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    bc$1[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[1] = _0 - (_i + bvirt) + (bvirt - cdxbdy1);
-    bc3 = _j + _i;
-    bvirt = bc3 - _j;
-    bc$1[2] = _j - (bc3 - bvirt) + (_i - bvirt);
-    bc$1[3] = bc3;
-    const axbclen = scale_expansion_zeroelim(4, bc$1, adx, axbc);
-    const axxbclen = scale_expansion_zeroelim(axbclen, axbc, adx, axxbc);
-    const aybclen = scale_expansion_zeroelim(4, bc$1, ady, aybc);
-    const ayybclen = scale_expansion_zeroelim(aybclen, aybc, ady, ayybc);
-    const alen = fast_expansion_sum_zeroelim(axxbclen, axxbc, ayybclen, ayybc, adet$1);
-
-    cdxady1 = cdx * ady;
+    bc$1[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    bc$1[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    bc$1[3] = u3;
+    s1 = cdx * ady;
     c = splitter * cdx;
     ahi = c - (c - cdx);
     alo = cdx - ahi;
     c = splitter * ady;
     bhi = c - (c - ady);
     blo = ady - bhi;
-    cdxady0 = alo * blo - (cdxady1 - ahi * bhi - alo * bhi - ahi * blo);
-    adxcdy1 = adx * cdy;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = adx * cdy;
     c = splitter * adx;
     ahi = c - (c - adx);
     alo = adx - ahi;
     c = splitter * cdy;
     bhi = c - (c - cdy);
     blo = cdy - bhi;
-    adxcdy0 = alo * blo - (adxcdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = cdxady0 - adxcdy0;
-    bvirt = cdxady0 - _i;
-    ca$1[0] = cdxady0 - (_i + bvirt) + (bvirt - adxcdy0);
-    _j = cdxady1 + _i;
-    bvirt = _j - cdxady1;
-    _0 = cdxady1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - adxcdy1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ca[0] = s0 - (_i + bvirt) + (bvirt - t0);
+    _j = s1 + _i;
+    bvirt = _j - s1;
+    _0 = s1 - (_j - bvirt) + (_i - bvirt);
+    _i = _0 - t1;
     bvirt = _0 - _i;
-    ca$1[1] = _0 - (_i + bvirt) + (bvirt - adxcdy1);
-    ca3 = _j + _i;
-    bvirt = ca3 - _j;
-    ca$1[2] = _j - (ca3 - bvirt) + (_i - bvirt);
-    ca$1[3] = ca3;
-    const bxcalen = scale_expansion_zeroelim(4, ca$1, bdx, bxca);
-    const bxxcalen = scale_expansion_zeroelim(bxcalen, bxca, bdx, bxxca);
-    const bycalen = scale_expansion_zeroelim(4, ca$1, bdy, byca);
-    const byycalen = scale_expansion_zeroelim(bycalen, byca, bdy, byyca);
-    const blen = fast_expansion_sum_zeroelim(bxxcalen, bxxca, byycalen, byyca, bdet$1);
-
-    adxbdy1 = adx * bdy;
+    ca[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    ca[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    ca[3] = u3;
+    s1 = adx * bdy;
     c = splitter * adx;
     ahi = c - (c - adx);
     alo = adx - ahi;
     c = splitter * bdy;
     bhi = c - (c - bdy);
     blo = bdy - bhi;
-    adxbdy0 = alo * blo - (adxbdy1 - ahi * bhi - alo * bhi - ahi * blo);
-    bdxady1 = bdx * ady;
+    s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+    t1 = bdx * ady;
     c = splitter * bdx;
     ahi = c - (c - bdx);
     alo = bdx - ahi;
     c = splitter * ady;
     bhi = c - (c - ady);
     blo = ady - bhi;
-    bdxady0 = alo * blo - (bdxady1 - ahi * bhi - alo * bhi - ahi * blo);
-    _i = adxbdy0 - bdxady0;
-    bvirt = adxbdy0 - _i;
-    ab$1[0] = adxbdy0 - (_i + bvirt) + (bvirt - bdxady0);
-    _j = adxbdy1 + _i;
-    bvirt = _j - adxbdy1;
-    _0 = adxbdy1 - (_j - bvirt) + (_i - bvirt);
-    _i = _0 - bdxady1;
+    t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+    _i = s0 - t0;
+    bvirt = s0 - _i;
+    ab$1[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[1] = _0 - (_i + bvirt) + (bvirt - bdxady1);
-    ab3 = _j + _i;
-    bvirt = ab3 - _j;
-    ab$1[2] = _j - (ab3 - bvirt) + (_i - bvirt);
-    ab$1[3] = ab3;
-    const cxablen = scale_expansion_zeroelim(4, ab$1, cdx, cxab);
-    const cxxablen = scale_expansion_zeroelim(cxablen, cxab, cdx, cxxab);
-    const cyablen = scale_expansion_zeroelim(4, ab$1, cdy, cyab);
-    const cyyablen = scale_expansion_zeroelim(cyablen, cyab, cdy, cyyab);
-    const clen = fast_expansion_sum_zeroelim(cxxablen, cxxab, cyyablen, cyyab, cdet$1);
-
-    const ablen = fast_expansion_sum_zeroelim(alen, adet$1, blen, bdet$1, abdet$1);
-    finlength = fast_expansion_sum_zeroelim(ablen, abdet$1, clen, cdet$1, fin1$1);
-
-    let det = estimate(finlength, fin1$1);
+    ab$1[1] = _0 - (_i + bvirt) + (bvirt - t1);
+    u3 = _j + _i;
+    bvirt = u3 - _j;
+    ab$1[2] = _j - (u3 - bvirt) + (_i - bvirt);
+    ab$1[3] = u3;
+
+    finlen = sum(
+        sum(
+            sum(
+                scale(scale(4, bc$1, adx, _8$1), _8$1, adx, _16$1), _16$1,
+                scale(scale(4, bc$1, ady, _8$1), _8$1, ady, _16b), _16b, _32), _32,
+            sum(
+                scale(scale(4, ca, bdx, _8$1), _8$1, bdx, _16$1), _16$1,
+                scale(scale(4, ca, bdy, _8$1), _8$1, bdy, _16b), _16b, _32b), _32b, _64), _64,
+        sum(
+            scale(scale(4, ab$1, cdx, _8$1), _8$1, cdx, _16$1), _16$1,
+            scale(scale(4, ab$1, cdy, _8$1), _8$1, cdy, _16b), _16b, _32), _32, fin$1);
+
+    let det = estimate(finlen, fin$1);
     let errbound = iccerrboundB * permanent;
     if (det >= errbound || -det >= errbound) {
         return det;
@@ -1572,624 +983,518 @@ function incircleadapt(ax, ay, bx, by, cx, cy, dx, dy, permanent) {
         return det;
     }
 
-    finnow = fin1$1;
-    finother = fin2$1;
-
     if (bdxtail !== 0 || bdytail !== 0 || cdxtail !== 0 || cdytail !== 0) {
-        adxadx1 = adx * adx;
+        s1 = adx * adx;
         c = splitter * adx;
         ahi = c - (c - adx);
         alo = adx - ahi;
-        adxadx0 = alo * alo - (adxadx1 - ahi * ahi - (ahi + ahi) * alo);
-        adyady1 = ady * ady;
+        s0 = alo * alo - (s1 - ahi * ahi - (ahi + ahi) * alo);
+        t1 = ady * ady;
         c = splitter * ady;
         ahi = c - (c - ady);
         alo = ady - ahi;
-        adyady0 = alo * alo - (adyady1 - ahi * ahi - (ahi + ahi) * alo);
-        _i = adxadx0 + adyady0;
-        bvirt = _i - adxadx0;
-        aa[0] = adxadx0 - (_i - bvirt) + (adyady0 - bvirt);
-        _j = adxadx1 + _i;
-        bvirt = _j - adxadx1;
-        _0 = adxadx1 - (_j - bvirt) + (_i - bvirt);
-        _i = _0 + adyady1;
+        t0 = alo * alo - (t1 - ahi * ahi - (ahi + ahi) * alo);
+        _i = s0 + t0;
+        bvirt = _i - s0;
+        aa[0] = s0 - (_i - bvirt) + (t0 - bvirt);
+        _j = s1 + _i;
+        bvirt = _j - s1;
+        _0 = s1 - (_j - bvirt) + (_i - bvirt);
+        _i = _0 + t1;
         bvirt = _i - _0;
-        aa[1] = _0 - (_i - bvirt) + (adyady1 - bvirt);
-        aa3 = _j + _i;
-        bvirt = aa3 - _j;
-        aa[2] = _j - (aa3 - bvirt) + (_i - bvirt);
-        aa[3] = aa3;
+        aa[1] = _0 - (_i - bvirt) + (t1 - bvirt);
+        u3 = _j + _i;
+        bvirt = u3 - _j;
+        aa[2] = _j - (u3 - bvirt) + (_i - bvirt);
+        aa[3] = u3;
     }
     if (cdxtail !== 0 || cdytail !== 0 || adxtail !== 0 || adytail !== 0) {
-        bdxbdx1 = bdx * bdx;
+        s1 = bdx * bdx;
         c = splitter * bdx;
         ahi = c - (c - bdx);
         alo = bdx - ahi;
-        bdxbdx0 = alo * alo - (bdxbdx1 - ahi * ahi - (ahi + ahi) * alo);
-        bdybdy1 = bdy * bdy;
+        s0 = alo * alo - (s1 - ahi * ahi - (ahi + ahi) * alo);
+        t1 = bdy * bdy;
         c = splitter * bdy;
         ahi = c - (c - bdy);
         alo = bdy - ahi;
-        bdybdy0 = alo * alo - (bdybdy1 - ahi * ahi - (ahi + ahi) * alo);
-        _i = bdxbdx0 + bdybdy0;
-        bvirt = _i - bdxbdx0;
-        bb[0] = bdxbdx0 - (_i - bvirt) + (bdybdy0 - bvirt);
-        _j = bdxbdx1 + _i;
-        bvirt = _j - bdxbdx1;
-        _0 = bdxbdx1 - (_j - bvirt) + (_i - bvirt);
-        _i = _0 + bdybdy1;
+        t0 = alo * alo - (t1 - ahi * ahi - (ahi + ahi) * alo);
+        _i = s0 + t0;
+        bvirt = _i - s0;
+        bb[0] = s0 - (_i - bvirt) + (t0 - bvirt);
+        _j = s1 + _i;
+        bvirt = _j - s1;
+        _0 = s1 - (_j - bvirt) + (_i - bvirt);
+        _i = _0 + t1;
         bvirt = _i - _0;
-        bb[1] = _0 - (_i - bvirt) + (bdybdy1 - bvirt);
-        bb3 = _j + _i;
-        bvirt = bb3 - _j;
-        bb[2] = _j - (bb3 - bvirt) + (_i - bvirt);
-        bb[3] = bb3;
+        bb[1] = _0 - (_i - bvirt) + (t1 - bvirt);
+        u3 = _j + _i;
+        bvirt = u3 - _j;
+        bb[2] = _j - (u3 - bvirt) + (_i - bvirt);
+        bb[3] = u3;
     }
     if (adxtail !== 0 || adytail !== 0 || bdxtail !== 0 || bdytail !== 0) {
-        cdxcdx1 = cdx * cdx;
+        s1 = cdx * cdx;
         c = splitter * cdx;
         ahi = c - (c - cdx);
         alo = cdx - ahi;
-        cdxcdx0 = alo * alo - (cdxcdx1 - ahi * ahi - (ahi + ahi) * alo);
-        cdycdy1 = cdy * cdy;
+        s0 = alo * alo - (s1 - ahi * ahi - (ahi + ahi) * alo);
+        t1 = cdy * cdy;
         c = splitter * cdy;
         ahi = c - (c - cdy);
         alo = cdy - ahi;
-        cdycdy0 = alo * alo - (cdycdy1 - ahi * ahi - (ahi + ahi) * alo);
-        _i = cdxcdx0 + cdycdy0;
-        bvirt = _i - cdxcdx0;
-        cc[0] = cdxcdx0 - (_i - bvirt) + (cdycdy0 - bvirt);
-        _j = cdxcdx1 + _i;
-        bvirt = _j - cdxcdx1;
-        _0 = cdxcdx1 - (_j - bvirt) + (_i - bvirt);
-        _i = _0 + cdycdy1;
+        t0 = alo * alo - (t1 - ahi * ahi - (ahi + ahi) * alo);
+        _i = s0 + t0;
+        bvirt = _i - s0;
+        cc[0] = s0 - (_i - bvirt) + (t0 - bvirt);
+        _j = s1 + _i;
+        bvirt = _j - s1;
+        _0 = s1 - (_j - bvirt) + (_i - bvirt);
+        _i = _0 + t1;
         bvirt = _i - _0;
-        cc[1] = _0 - (_i - bvirt) + (cdycdy1 - bvirt);
-        cc3 = _j + _i;
-        bvirt = cc3 - _j;
-        cc[2] = _j - (cc3 - bvirt) + (_i - bvirt);
-        cc[3] = cc3;
+        cc[1] = _0 - (_i - bvirt) + (t1 - bvirt);
+        u3 = _j + _i;
+        bvirt = u3 - _j;
+        cc[2] = _j - (u3 - bvirt) + (_i - bvirt);
+        cc[3] = u3;
     }
 
     if (adxtail !== 0) {
-        axtbclen = scale_expansion_zeroelim(4, bc$1, adxtail, axtbc);
-        temp16alen = scale_expansion_zeroelim(axtbclen, axtbc, 2 * adx, temp16a);
-
-        const axtcclen = scale_expansion_zeroelim(4, cc, adxtail, axtcc);
-        temp16blen = scale_expansion_zeroelim(axtcclen, axtcc, bdy, temp16b);
-
-        const axtbblen = scale_expansion_zeroelim(4, bb, adxtail, axtbb);
-        temp16clen = scale_expansion_zeroelim(axtbblen, axtbb, -cdy, temp16c);
-
-        temp32alen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32a);
-        temp48len = fast_expansion_sum_zeroelim(temp16clen, temp16c, temp32alen, temp32a, temp48);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
+        axtbclen = scale(4, bc$1, adxtail, axtbc);
+        finlen = finadd(finlen, sum_three(
+            scale(axtbclen, axtbc, 2 * adx, _16$1), _16$1,
+            scale(scale(4, cc, adxtail, _8$1), _8$1, bdy, _16b), _16b,
+            scale(scale(4, bb, adxtail, _8$1), _8$1, -cdy, _16c), _16c, _32, _48$1), _48$1);
     }
     if (adytail !== 0) {
-        aytbclen = scale_expansion_zeroelim(4, bc$1, adytail, aytbc);
-        temp16alen = scale_expansion_zeroelim(aytbclen, aytbc, 2 * ady, temp16a);
-
-        const aytbblen = scale_expansion_zeroelim(4, bb, adytail, aytbb);
-        temp16blen = scale_expansion_zeroelim(aytbblen, aytbb, cdx, temp16b);
-
-        const aytcclen = scale_expansion_zeroelim(4, cc, adytail, aytcc);
-        temp16clen = scale_expansion_zeroelim(aytcclen, aytcc, -bdx, temp16c);
-
-        temp32alen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32a);
-        temp48len = fast_expansion_sum_zeroelim(temp16clen, temp16c, temp32alen, temp32a, temp48);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
+        aytbclen = scale(4, bc$1, adytail, aytbc);
+        finlen = finadd(finlen, sum_three(
+            scale(aytbclen, aytbc, 2 * ady, _16$1), _16$1,
+            scale(scale(4, bb, adytail, _8$1), _8$1, cdx, _16b), _16b,
+            scale(scale(4, cc, adytail, _8$1), _8$1, -bdx, _16c), _16c, _32, _48$1), _48$1);
     }
     if (bdxtail !== 0) {
-        bxtcalen = scale_expansion_zeroelim(4, ca$1, bdxtail, bxtca);
-        temp16alen = scale_expansion_zeroelim(bxtcalen, bxtca, 2 * bdx, temp16a);
-
-        const bxtaalen = scale_expansion_zeroelim(4, aa, bdxtail, bxtaa);
-        temp16blen = scale_expansion_zeroelim(bxtaalen, bxtaa, cdy, temp16b);
-
-        const bxtcclen = scale_expansion_zeroelim(4, cc, bdxtail, bxtcc);
-        temp16clen = scale_expansion_zeroelim(bxtcclen, bxtcc, -ady, temp16c);
-
-        temp32alen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32a);
-        temp48len = fast_expansion_sum_zeroelim(temp16clen, temp16c, temp32alen, temp32a, temp48);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
+        bxtcalen = scale(4, ca, bdxtail, bxtca);
+        finlen = finadd(finlen, sum_three(
+            scale(bxtcalen, bxtca, 2 * bdx, _16$1), _16$1,
+            scale(scale(4, aa, bdxtail, _8$1), _8$1, cdy, _16b), _16b,
+            scale(scale(4, cc, bdxtail, _8$1), _8$1, -ady, _16c), _16c, _32, _48$1), _48$1);
     }
     if (bdytail !== 0) {
-        bytcalen = scale_expansion_zeroelim(4, ca$1, bdytail, bytca);
-        temp16alen = scale_expansion_zeroelim(bytcalen, bytca, 2 * bdy, temp16a);
-
-        const bytcclen = scale_expansion_zeroelim(4, cc, bdytail, bytcc);
-        temp16blen = scale_expansion_zeroelim(bytcclen, bytcc, adx, temp16b);
-
-        const bytaalen = scale_expansion_zeroelim(4, aa, bdytail, bytaa);
-        temp16clen = scale_expansion_zeroelim(bytaalen, bytaa, -cdx, temp16c);
-
-        temp32alen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32a);
-        temp48len = fast_expansion_sum_zeroelim(temp16clen, temp16c, temp32alen, temp32a, temp48);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
+        bytcalen = scale(4, ca, bdytail, bytca);
+        finlen = finadd(finlen, sum_three(
+            scale(bytcalen, bytca, 2 * bdy, _16$1), _16$1,
+            scale(scale(4, cc, bdytail, _8$1), _8$1, adx, _16b), _16b,
+            scale(scale(4, aa, bdytail, _8$1), _8$1, -cdx, _16c), _16c, _32, _48$1), _48$1);
     }
     if (cdxtail !== 0) {
-        cxtablen = scale_expansion_zeroelim(4, ab$1, cdxtail, cxtab);
-        temp16alen = scale_expansion_zeroelim(cxtablen, cxtab, 2 * cdx, temp16a);
-
-        const cxtbblen = scale_expansion_zeroelim(4, bb, cdxtail, cxtbb);
-        temp16blen = scale_expansion_zeroelim(cxtbblen, cxtbb, ady, temp16b);
-
-        const cxtaalen = scale_expansion_zeroelim(4, aa, cdxtail, cxtaa);
-        temp16clen = scale_expansion_zeroelim(cxtaalen, cxtaa, -bdy, temp16c);
-
-        temp32alen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32a);
-        temp48len = fast_expansion_sum_zeroelim(temp16clen, temp16c, temp32alen, temp32a, temp48);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
+        cxtablen = scale(4, ab$1, cdxtail, cxtab);
+        finlen = finadd(finlen, sum_three(
+            scale(cxtablen, cxtab, 2 * cdx, _16$1), _16$1,
+            scale(scale(4, bb, cdxtail, _8$1), _8$1, ady, _16b), _16b,
+            scale(scale(4, aa, cdxtail, _8$1), _8$1, -bdy, _16c), _16c, _32, _48$1), _48$1);
     }
     if (cdytail !== 0) {
-        cytablen = scale_expansion_zeroelim(4, ab$1, cdytail, cytab);
-        temp16alen = scale_expansion_zeroelim(cytablen, cytab, 2 * cdy, temp16a);
-
-        const cytaalen = scale_expansion_zeroelim(4, aa, cdytail, cytaa);
-        temp16blen = scale_expansion_zeroelim(cytaalen, cytaa, bdx, temp16b);
-
-        const cytbblen = scale_expansion_zeroelim(4, bb, cdytail, cytbb);
-        temp16clen = scale_expansion_zeroelim(cytbblen, cytbb, -adx, temp16c);
-
-        temp32alen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32a);
-        temp48len = fast_expansion_sum_zeroelim(temp16clen, temp16c, temp32alen, temp32a, temp48);
-        finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-        finswap = finnow; finnow = finother; finother = finswap;
+        cytablen = scale(4, ab$1, cdytail, cytab);
+        finlen = finadd(finlen, sum_three(
+            scale(cytablen, cytab, 2 * cdy, _16$1), _16$1,
+            scale(scale(4, aa, cdytail, _8$1), _8$1, bdx, _16b), _16b,
+            scale(scale(4, bb, cdytail, _8$1), _8$1, -adx, _16c), _16c, _32, _48$1), _48$1);
     }
 
     if (adxtail !== 0 || adytail !== 0) {
         if (bdxtail !== 0 || bdytail !== 0 || cdxtail !== 0 || cdytail !== 0) {
-            ti1 = bdxtail * cdy;
+            s1 = bdxtail * cdy;
             c = splitter * bdxtail;
             ahi = c - (c - bdxtail);
             alo = bdxtail - ahi;
             c = splitter * cdy;
             bhi = c - (c - cdy);
             blo = cdy - bhi;
-            ti0 = alo * blo - (ti1 - ahi * bhi - alo * bhi - ahi * blo);
-            tj1 = bdx * cdytail;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = bdx * cdytail;
             c = splitter * bdx;
             ahi = c - (c - bdx);
             alo = bdx - ahi;
             c = splitter * cdytail;
             bhi = c - (c - cdytail);
             blo = cdytail - bhi;
-            tj0 = alo * blo - (tj1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = ti0 + tj0;
-            bvirt = _i - ti0;
-            u$2[0] = ti0 - (_i - bvirt) + (tj0 - bvirt);
-            _j = ti1 + _i;
-            bvirt = _j - ti1;
-            _0 = ti1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 + tj1;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 + t0;
+            bvirt = _i - s0;
+            u[0] = s0 - (_i - bvirt) + (t0 - bvirt);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 + t1;
             bvirt = _i - _0;
-            u$2[1] = _0 - (_i - bvirt) + (tj1 - bvirt);
+            u[1] = _0 - (_i - bvirt) + (t1 - bvirt);
             u3 = _j + _i;
             bvirt = u3 - _j;
-            u$2[2] = _j - (u3 - bvirt) + (_i - bvirt);
-            u$2[3] = u3;
-            negate = -bdy;
-            ti1 = cdxtail * negate;
+            u[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            u[3] = u3;
+            s1 = cdxtail * -bdy;
             c = splitter * cdxtail;
             ahi = c - (c - cdxtail);
             alo = cdxtail - ahi;
-            c = splitter * negate;
-            bhi = c - (c - negate);
-            blo = negate - bhi;
-            ti0 = alo * blo - (ti1 - ahi * bhi - alo * bhi - ahi * blo);
-            negate = -bdytail;
-            tj1 = cdx * negate;
+            c = splitter * -bdy;
+            bhi = c - (c - -bdy);
+            blo = -bdy - bhi;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = cdx * -bdytail;
             c = splitter * cdx;
             ahi = c - (c - cdx);
             alo = cdx - ahi;
-            c = splitter * negate;
-            bhi = c - (c - negate);
-            blo = negate - bhi;
-            tj0 = alo * blo - (tj1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = ti0 + tj0;
-            bvirt = _i - ti0;
-            v$1[0] = ti0 - (_i - bvirt) + (tj0 - bvirt);
-            _j = ti1 + _i;
-            bvirt = _j - ti1;
-            _0 = ti1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 + tj1;
+            c = splitter * -bdytail;
+            bhi = c - (c - -bdytail);
+            blo = -bdytail - bhi;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 + t0;
+            bvirt = _i - s0;
+            v[0] = s0 - (_i - bvirt) + (t0 - bvirt);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 + t1;
             bvirt = _i - _0;
-            v$1[1] = _0 - (_i - bvirt) + (tj1 - bvirt);
-            v3 = _j + _i;
-            bvirt = v3 - _j;
-            v$1[2] = _j - (v3 - bvirt) + (_i - bvirt);
-            v$1[3] = v3;
-            bctlen = fast_expansion_sum_zeroelim(4, u$2, 4, v$1, bct$1);
-
-            ti1 = bdxtail * cdytail;
+            v[1] = _0 - (_i - bvirt) + (t1 - bvirt);
+            u3 = _j + _i;
+            bvirt = u3 - _j;
+            v[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            v[3] = u3;
+            bctlen = sum(4, u, 4, v, bct);
+            s1 = bdxtail * cdytail;
             c = splitter * bdxtail;
             ahi = c - (c - bdxtail);
             alo = bdxtail - ahi;
             c = splitter * cdytail;
             bhi = c - (c - cdytail);
             blo = cdytail - bhi;
-            ti0 = alo * blo - (ti1 - ahi * bhi - alo * bhi - ahi * blo);
-            tj1 = cdxtail * bdytail;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = cdxtail * bdytail;
             c = splitter * cdxtail;
             ahi = c - (c - cdxtail);
             alo = cdxtail - ahi;
             c = splitter * bdytail;
             bhi = c - (c - bdytail);
             blo = bdytail - bhi;
-            tj0 = alo * blo - (tj1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = ti0 - tj0;
-            bvirt = ti0 - _i;
-            bctt[0] = ti0 - (_i + bvirt) + (bvirt - tj0);
-            _j = ti1 + _i;
-            bvirt = _j - ti1;
-            _0 = ti1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 - tj1;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 - t0;
+            bvirt = s0 - _i;
+            bctt[0] = s0 - (_i + bvirt) + (bvirt - t0);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 - t1;
             bvirt = _0 - _i;
-            bctt[1] = _0 - (_i + bvirt) + (bvirt - tj1);
-            bctt3 = _j + _i;
-            bvirt = bctt3 - _j;
-            bctt[2] = _j - (bctt3 - bvirt) + (_i - bvirt);
-            bctt[3] = bctt3;
+            bctt[1] = _0 - (_i + bvirt) + (bvirt - t1);
+            u3 = _j + _i;
+            bvirt = u3 - _j;
+            bctt[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            bctt[3] = u3;
             bcttlen = 4;
         } else {
-            bct$1[0] = 0;
+            bct[0] = 0;
             bctlen = 1;
             bctt[0] = 0;
             bcttlen = 1;
         }
-
         if (adxtail !== 0) {
-            temp16alen = scale_expansion_zeroelim(axtbclen, axtbc, adxtail, temp16a);
-            const axtbctlen = scale_expansion_zeroelim(bctlen, bct$1, adxtail, axtbct);
-            temp32alen = scale_expansion_zeroelim(axtbctlen, axtbct, 2 * adx, temp32a);
-            temp48len = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp32alen, temp32a, temp48);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
+            const len = scale(bctlen, bct, adxtail, _16c);
+            finlen = finadd(finlen, sum(
+                scale(axtbclen, axtbc, adxtail, _16$1), _16$1,
+                scale(len, _16c, 2 * adx, _32), _32, _48$1), _48$1);
+
+            const len2 = scale(bcttlen, bctt, adxtail, _8$1);
+            finlen = finadd(finlen, sum_three(
+                scale(len2, _8$1, 2 * adx, _16$1), _16$1,
+                scale(len2, _8$1, adxtail, _16b), _16b,
+                scale(len, _16c, adxtail, _32), _32, _32b, _64), _64);
 
             if (bdytail !== 0) {
-                temp8len = scale_expansion_zeroelim(4, cc, adxtail, temp8);
-                temp16alen = scale_expansion_zeroelim(temp8len, temp8, bdytail, temp16a);
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp16alen, temp16a, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
+                finlen = finadd(finlen, scale(scale(4, cc, adxtail, _8$1), _8$1, bdytail, _16$1), _16$1);
             }
             if (cdytail !== 0) {
-                temp8len = scale_expansion_zeroelim(4, bb, -adxtail, temp8);
-                temp16alen = scale_expansion_zeroelim(temp8len, temp8, cdytail, temp16a);
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp16alen, temp16a, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
+                finlen = finadd(finlen, scale(scale(4, bb, -adxtail, _8$1), _8$1, cdytail, _16$1), _16$1);
             }
-
-            temp32alen = scale_expansion_zeroelim(axtbctlen, axtbct, adxtail, temp32a);
-            const axtbcttlen = scale_expansion_zeroelim(bcttlen, bctt, adxtail, axtbctt);
-            temp16alen = scale_expansion_zeroelim(axtbcttlen, axtbctt, 2 * adx, temp16a);
-            temp16blen = scale_expansion_zeroelim(axtbcttlen, axtbctt, adxtail, temp16b);
-            temp32blen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32b);
-            temp64len = fast_expansion_sum_zeroelim(temp32alen, temp32a, temp32blen, temp32b, temp64);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp64len, temp64, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
         }
-
         if (adytail !== 0) {
-            temp16alen = scale_expansion_zeroelim(aytbclen, aytbc, adytail, temp16a);
-            const aytbctlen = scale_expansion_zeroelim(bctlen, bct$1, adytail, aytbct);
-            temp32alen = scale_expansion_zeroelim(aytbctlen, aytbct, 2 * ady, temp32a);
-            temp48len = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp32alen, temp32a, temp48);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
-
-            temp32alen = scale_expansion_zeroelim(aytbctlen, aytbct, adytail, temp32a);
-            const aytbcttlen = scale_expansion_zeroelim(bcttlen, bctt, adytail, aytbctt);
-            temp16alen = scale_expansion_zeroelim(aytbcttlen, aytbctt, 2 * ady, temp16a);
-            temp16blen = scale_expansion_zeroelim(aytbcttlen, aytbctt, adytail, temp16b);
-            temp32blen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32b);
-            temp64len = fast_expansion_sum_zeroelim(temp32alen, temp32a, temp32blen, temp32b, temp64);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp64len, temp64, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
+            const len = scale(bctlen, bct, adytail, _16c);
+            finlen = finadd(finlen, sum(
+                scale(aytbclen, aytbc, adytail, _16$1), _16$1,
+                scale(len, _16c, 2 * ady, _32), _32, _48$1), _48$1);
+
+            const len2 = scale(bcttlen, bctt, adytail, _8$1);
+            finlen = finadd(finlen, sum_three(
+                scale(len2, _8$1, 2 * ady, _16$1), _16$1,
+                scale(len2, _8$1, adytail, _16b), _16b,
+                scale(len, _16c, adytail, _32), _32, _32b, _64), _64);
         }
     }
     if (bdxtail !== 0 || bdytail !== 0) {
         if (cdxtail !== 0 || cdytail !== 0 || adxtail !== 0 || adytail !== 0) {
-            ti1 = cdxtail * ady;
+            s1 = cdxtail * ady;
             c = splitter * cdxtail;
             ahi = c - (c - cdxtail);
             alo = cdxtail - ahi;
             c = splitter * ady;
             bhi = c - (c - ady);
             blo = ady - bhi;
-            ti0 = alo * blo - (ti1 - ahi * bhi - alo * bhi - ahi * blo);
-            tj1 = cdx * adytail;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = cdx * adytail;
             c = splitter * cdx;
             ahi = c - (c - cdx);
             alo = cdx - ahi;
             c = splitter * adytail;
             bhi = c - (c - adytail);
             blo = adytail - bhi;
-            tj0 = alo * blo - (tj1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = ti0 + tj0;
-            bvirt = _i - ti0;
-            u$2[0] = ti0 - (_i - bvirt) + (tj0 - bvirt);
-            _j = ti1 + _i;
-            bvirt = _j - ti1;
-            _0 = ti1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 + tj1;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 + t0;
+            bvirt = _i - s0;
+            u[0] = s0 - (_i - bvirt) + (t0 - bvirt);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 + t1;
             bvirt = _i - _0;
-            u$2[1] = _0 - (_i - bvirt) + (tj1 - bvirt);
+            u[1] = _0 - (_i - bvirt) + (t1 - bvirt);
             u3 = _j + _i;
             bvirt = u3 - _j;
-            u$2[2] = _j - (u3 - bvirt) + (_i - bvirt);
-            u$2[3] = u3;
-            negate = -cdy;
-            ti1 = adxtail * negate;
+            u[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            u[3] = u3;
+            n1 = -cdy;
+            n0 = -cdytail;
+            s1 = adxtail * n1;
             c = splitter * adxtail;
             ahi = c - (c - adxtail);
             alo = adxtail - ahi;
-            c = splitter * negate;
-            bhi = c - (c - negate);
-            blo = negate - bhi;
-            ti0 = alo * blo - (ti1 - ahi * bhi - alo * bhi - ahi * blo);
-            negate = -cdytail;
-            tj1 = adx * negate;
+            c = splitter * n1;
+            bhi = c - (c - n1);
+            blo = n1 - bhi;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = adx * n0;
             c = splitter * adx;
             ahi = c - (c - adx);
             alo = adx - ahi;
-            c = splitter * negate;
-            bhi = c - (c - negate);
-            blo = negate - bhi;
-            tj0 = alo * blo - (tj1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = ti0 + tj0;
-            bvirt = _i - ti0;
-            v$1[0] = ti0 - (_i - bvirt) + (tj0 - bvirt);
-            _j = ti1 + _i;
-            bvirt = _j - ti1;
-            _0 = ti1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 + tj1;
+            c = splitter * n0;
+            bhi = c - (c - n0);
+            blo = n0 - bhi;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 + t0;
+            bvirt = _i - s0;
+            v[0] = s0 - (_i - bvirt) + (t0 - bvirt);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 + t1;
             bvirt = _i - _0;
-            v$1[1] = _0 - (_i - bvirt) + (tj1 - bvirt);
-            v3 = _j + _i;
-            bvirt = v3 - _j;
-            v$1[2] = _j - (v3 - bvirt) + (_i - bvirt);
-            v$1[3] = v3;
-            catlen = fast_expansion_sum_zeroelim(4, u$2, 4, v$1, cat$1);
-
-            ti1 = cdxtail * adytail;
+            v[1] = _0 - (_i - bvirt) + (t1 - bvirt);
+            u3 = _j + _i;
+            bvirt = u3 - _j;
+            v[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            v[3] = u3;
+            catlen = sum(4, u, 4, v, cat);
+            s1 = cdxtail * adytail;
             c = splitter * cdxtail;
             ahi = c - (c - cdxtail);
             alo = cdxtail - ahi;
             c = splitter * adytail;
             bhi = c - (c - adytail);
             blo = adytail - bhi;
-            ti0 = alo * blo - (ti1 - ahi * bhi - alo * bhi - ahi * blo);
-            tj1 = adxtail * cdytail;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = adxtail * cdytail;
             c = splitter * adxtail;
             ahi = c - (c - adxtail);
             alo = adxtail - ahi;
             c = splitter * cdytail;
             bhi = c - (c - cdytail);
             blo = cdytail - bhi;
-            tj0 = alo * blo - (tj1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = ti0 - tj0;
-            bvirt = ti0 - _i;
-            catt[0] = ti0 - (_i + bvirt) + (bvirt - tj0);
-            _j = ti1 + _i;
-            bvirt = _j - ti1;
-            _0 = ti1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 - tj1;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 - t0;
+            bvirt = s0 - _i;
+            catt[0] = s0 - (_i + bvirt) + (bvirt - t0);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 - t1;
             bvirt = _0 - _i;
-            catt[1] = _0 - (_i + bvirt) + (bvirt - tj1);
-            catt3 = _j + _i;
-            bvirt = catt3 - _j;
-            catt[2] = _j - (catt3 - bvirt) + (_i - bvirt);
-            catt[3] = catt3;
+            catt[1] = _0 - (_i + bvirt) + (bvirt - t1);
+            u3 = _j + _i;
+            bvirt = u3 - _j;
+            catt[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            catt[3] = u3;
             cattlen = 4;
         } else {
-            cat$1[0] = 0;
+            cat[0] = 0;
             catlen = 1;
             catt[0] = 0;
             cattlen = 1;
         }
-
         if (bdxtail !== 0) {
-            temp16alen = scale_expansion_zeroelim(bxtcalen, bxtca, bdxtail, temp16a);
-            const bxtcatlen = scale_expansion_zeroelim(catlen, cat$1, bdxtail, bxtcat);
-            temp32alen = scale_expansion_zeroelim(bxtcatlen, bxtcat, 2 * bdx, temp32a);
-            temp48len = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp32alen, temp32a, temp48);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
+            const len = scale(catlen, cat, bdxtail, _16c);
+            finlen = finadd(finlen, sum(
+                scale(bxtcalen, bxtca, bdxtail, _16$1), _16$1,
+                scale(len, _16c, 2 * bdx, _32), _32, _48$1), _48$1);
+
+            const len2 = scale(cattlen, catt, bdxtail, _8$1);
+            finlen = finadd(finlen, sum_three(
+                scale(len2, _8$1, 2 * bdx, _16$1), _16$1,
+                scale(len2, _8$1, bdxtail, _16b), _16b,
+                scale(len, _16c, bdxtail, _32), _32, _32b, _64), _64);
+
             if (cdytail !== 0) {
-                temp8len = scale_expansion_zeroelim(4, aa, bdxtail, temp8);
-                temp16alen = scale_expansion_zeroelim(temp8len, temp8, cdytail, temp16a);
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp16alen, temp16a, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
+                finlen = finadd(finlen, scale(scale(4, aa, bdxtail, _8$1), _8$1, cdytail, _16$1), _16$1);
             }
             if (adytail !== 0) {
-                temp8len = scale_expansion_zeroelim(4, cc, -bdxtail, temp8);
-                temp16alen = scale_expansion_zeroelim(temp8len, temp8, adytail, temp16a);
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp16alen, temp16a, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
+                finlen = finadd(finlen, scale(scale(4, cc, -bdxtail, _8$1), _8$1, adytail, _16$1), _16$1);
             }
-
-            temp32alen = scale_expansion_zeroelim(bxtcatlen, bxtcat, bdxtail, temp32a);
-            const bxtcattlen = scale_expansion_zeroelim(cattlen, catt, bdxtail, bxtcatt);
-            temp16alen = scale_expansion_zeroelim(bxtcattlen, bxtcatt, 2 * bdx, temp16a);
-            temp16blen = scale_expansion_zeroelim(bxtcattlen, bxtcatt, bdxtail, temp16b);
-            temp32blen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32b);
-            temp64len = fast_expansion_sum_zeroelim(temp32alen, temp32a, temp32blen, temp32b, temp64);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp64len, temp64, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
         }
         if (bdytail !== 0) {
-            temp16alen = scale_expansion_zeroelim(bytcalen, bytca, bdytail, temp16a);
-            const bytcatlen = scale_expansion_zeroelim(catlen, cat$1, bdytail, bytcat);
-            temp32alen = scale_expansion_zeroelim(bytcatlen, bytcat, 2 * bdy, temp32a);
-            temp48len = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp32alen, temp32a, temp48);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
-
-            temp32alen = scale_expansion_zeroelim(bytcatlen, bytcat, bdytail, temp32a);
-            const bytcattlen = scale_expansion_zeroelim(cattlen, catt, bdytail, bytcatt);
-            temp16alen = scale_expansion_zeroelim(bytcattlen, bytcatt, 2 * bdy, temp16a);
-            temp16blen = scale_expansion_zeroelim(bytcattlen, bytcatt, bdytail, temp16b);
-            temp32blen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32b);
-            temp64len = fast_expansion_sum_zeroelim(temp32alen, temp32a, temp32blen, temp32b, temp64);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp64len, temp64, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
+            const len = scale(catlen, cat, bdytail, _16c);
+            finlen = finadd(finlen, sum(
+                scale(bytcalen, bytca, bdytail, _16$1), _16$1,
+                scale(len, _16c, 2 * bdy, _32), _32, _48$1), _48$1);
+
+            const len2 = scale(cattlen, catt, bdytail, _8$1);
+            finlen = finadd(finlen, sum_three(
+                scale(len2, _8$1, 2 * bdy, _16$1), _16$1,
+                scale(len2, _8$1, bdytail, _16b), _16b,
+                scale(len, _16c, bdytail, _32), _32,  _32b, _64), _64);
         }
     }
     if (cdxtail !== 0 || cdytail !== 0) {
         if (adxtail !== 0 || adytail !== 0 || bdxtail !== 0 || bdytail !== 0) {
-            ti1 = adxtail * bdy;
+            s1 = adxtail * bdy;
             c = splitter * adxtail;
             ahi = c - (c - adxtail);
             alo = adxtail - ahi;
             c = splitter * bdy;
             bhi = c - (c - bdy);
             blo = bdy - bhi;
-            ti0 = alo * blo - (ti1 - ahi * bhi - alo * bhi - ahi * blo);
-            tj1 = adx * bdytail;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = adx * bdytail;
             c = splitter * adx;
             ahi = c - (c - adx);
             alo = adx - ahi;
             c = splitter * bdytail;
             bhi = c - (c - bdytail);
             blo = bdytail - bhi;
-            tj0 = alo * blo - (tj1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = ti0 + tj0;
-            bvirt = _i - ti0;
-            u$2[0] = ti0 - (_i - bvirt) + (tj0 - bvirt);
-            _j = ti1 + _i;
-            bvirt = _j - ti1;
-            _0 = ti1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 + tj1;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 + t0;
+            bvirt = _i - s0;
+            u[0] = s0 - (_i - bvirt) + (t0 - bvirt);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 + t1;
             bvirt = _i - _0;
-            u$2[1] = _0 - (_i - bvirt) + (tj1 - bvirt);
+            u[1] = _0 - (_i - bvirt) + (t1 - bvirt);
             u3 = _j + _i;
             bvirt = u3 - _j;
-            u$2[2] = _j - (u3 - bvirt) + (_i - bvirt);
-            u$2[3] = u3;
-            negate = -ady;
-            ti1 = bdxtail * negate;
+            u[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            u[3] = u3;
+            n1 = -ady;
+            n0 = -adytail;
+            s1 = bdxtail * n1;
             c = splitter * bdxtail;
             ahi = c - (c - bdxtail);
             alo = bdxtail - ahi;
-            c = splitter * negate;
-            bhi = c - (c - negate);
-            blo = negate - bhi;
-            ti0 = alo * blo - (ti1 - ahi * bhi - alo * bhi - ahi * blo);
-            negate = -adytail;
-            tj1 = bdx * negate;
+            c = splitter * n1;
+            bhi = c - (c - n1);
+            blo = n1 - bhi;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = bdx * n0;
             c = splitter * bdx;
             ahi = c - (c - bdx);
             alo = bdx - ahi;
-            c = splitter * negate;
-            bhi = c - (c - negate);
-            blo = negate - bhi;
-            tj0 = alo * blo - (tj1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = ti0 + tj0;
-            bvirt = _i - ti0;
-            v$1[0] = ti0 - (_i - bvirt) + (tj0 - bvirt);
-            _j = ti1 + _i;
-            bvirt = _j - ti1;
-            _0 = ti1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 + tj1;
+            c = splitter * n0;
+            bhi = c - (c - n0);
+            blo = n0 - bhi;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 + t0;
+            bvirt = _i - s0;
+            v[0] = s0 - (_i - bvirt) + (t0 - bvirt);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 + t1;
             bvirt = _i - _0;
-            v$1[1] = _0 - (_i - bvirt) + (tj1 - bvirt);
-            v3 = _j + _i;
-            bvirt = v3 - _j;
-            v$1[2] = _j - (v3 - bvirt) + (_i - bvirt);
-            v$1[3] = v3;
-            abtlen = fast_expansion_sum_zeroelim(4, u$2, 4, v$1, abt$1);
-
-            ti1 = adxtail * bdytail;
+            v[1] = _0 - (_i - bvirt) + (t1 - bvirt);
+            u3 = _j + _i;
+            bvirt = u3 - _j;
+            v[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            v[3] = u3;
+            abtlen = sum(4, u, 4, v, abt);
+            s1 = adxtail * bdytail;
             c = splitter * adxtail;
             ahi = c - (c - adxtail);
             alo = adxtail - ahi;
             c = splitter * bdytail;
             bhi = c - (c - bdytail);
             blo = bdytail - bhi;
-            ti0 = alo * blo - (ti1 - ahi * bhi - alo * bhi - ahi * blo);
-            tj1 = bdxtail * adytail;
+            s0 = alo * blo - (s1 - ahi * bhi - alo * bhi - ahi * blo);
+            t1 = bdxtail * adytail;
             c = splitter * bdxtail;
             ahi = c - (c - bdxtail);
             alo = bdxtail - ahi;
             c = splitter * adytail;
             bhi = c - (c - adytail);
             blo = adytail - bhi;
-            tj0 = alo * blo - (tj1 - ahi * bhi - alo * bhi - ahi * blo);
-            _i = ti0 - tj0;
-            bvirt = ti0 - _i;
-            abtt[0] = ti0 - (_i + bvirt) + (bvirt - tj0);
-            _j = ti1 + _i;
-            bvirt = _j - ti1;
-            _0 = ti1 - (_j - bvirt) + (_i - bvirt);
-            _i = _0 - tj1;
+            t0 = alo * blo - (t1 - ahi * bhi - alo * bhi - ahi * blo);
+            _i = s0 - t0;
+            bvirt = s0 - _i;
+            abtt[0] = s0 - (_i + bvirt) + (bvirt - t0);
+            _j = s1 + _i;
+            bvirt = _j - s1;
+            _0 = s1 - (_j - bvirt) + (_i - bvirt);
+            _i = _0 - t1;
             bvirt = _0 - _i;
-            abtt[1] = _0 - (_i + bvirt) + (bvirt - tj1);
-            abtt3 = _j + _i;
-            bvirt = abtt3 - _j;
-            abtt[2] = _j - (abtt3 - bvirt) + (_i - bvirt);
-            abtt[3] = abtt3;
+            abtt[1] = _0 - (_i + bvirt) + (bvirt - t1);
+            u3 = _j + _i;
+            bvirt = u3 - _j;
+            abtt[2] = _j - (u3 - bvirt) + (_i - bvirt);
+            abtt[3] = u3;
             abttlen = 4;
         } else {
-            abt$1[0] = 0;
+            abt[0] = 0;
             abtlen = 1;
             abtt[0] = 0;
             abttlen = 1;
         }
-
         if (cdxtail !== 0) {
-            temp16alen = scale_expansion_zeroelim(cxtablen, cxtab, cdxtail, temp16a);
-            const cxtabtlen = scale_expansion_zeroelim(abtlen, abt$1, cdxtail, cxtabt);
-            temp32alen = scale_expansion_zeroelim(cxtabtlen, cxtabt, 2 * cdx, temp32a);
-            temp48len = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp32alen, temp32a, temp48);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
+            const len = scale(abtlen, abt, cdxtail, _16c);
+            finlen = finadd(finlen, sum(
+                scale(cxtablen, cxtab, cdxtail, _16$1), _16$1,
+                scale(len, _16c, 2 * cdx, _32), _32, _48$1), _48$1);
+
+            const len2 = scale(abttlen, abtt, cdxtail, _8$1);
+            finlen = finadd(finlen, sum_three(
+                scale(len2, _8$1, 2 * cdx, _16$1), _16$1,
+                scale(len2, _8$1, cdxtail, _16b), _16b,
+                scale(len, _16c, cdxtail, _32), _32, _32b, _64), _64);
+
             if (adytail !== 0) {
-                temp8len = scale_expansion_zeroelim(4, bb, cdxtail, temp8);
-                temp16alen = scale_expansion_zeroelim(temp8len, temp8, adytail, temp16a);
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp16alen, temp16a, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
+                finlen = finadd(finlen, scale(scale(4, bb, cdxtail, _8$1), _8$1, adytail, _16$1), _16$1);
             }
             if (bdytail !== 0) {
-                temp8len = scale_expansion_zeroelim(4, aa, -cdxtail, temp8);
-                temp16alen = scale_expansion_zeroelim(temp8len, temp8, bdytail, temp16a);
-                finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp16alen, temp16a, finother);
-                finswap = finnow; finnow = finother; finother = finswap;
+                finlen = finadd(finlen, scale(scale(4, aa, -cdxtail, _8$1), _8$1, bdytail, _16$1), _16$1);
             }
-
-            temp32alen = scale_expansion_zeroelim(cxtabtlen, cxtabt, cdxtail, temp32a);
-            const cxtabttlen = scale_expansion_zeroelim(abttlen, abtt, cdxtail, cxtabtt);
-            temp16alen = scale_expansion_zeroelim(cxtabttlen, cxtabtt, 2 * cdx, temp16a);
-            temp16blen = scale_expansion_zeroelim(cxtabttlen, cxtabtt, cdxtail, temp16b);
-            temp32blen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32b);
-            temp64len = fast_expansion_sum_zeroelim(temp32alen, temp32a, temp32blen, temp32b, temp64);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp64len, temp64, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
         }
         if (cdytail !== 0) {
-            temp16alen = scale_expansion_zeroelim(cytablen, cytab, cdytail, temp16a);
-            const cytabtlen = scale_expansion_zeroelim(abtlen, abt$1, cdytail, cytabt);
-            temp32alen = scale_expansion_zeroelim(cytabtlen, cytabt, 2 * cdy, temp32a);
-            temp48len = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp32alen, temp32a, temp48);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp48len, temp48, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
-
-            temp32alen = scale_expansion_zeroelim(cytabtlen, cytabt, cdytail, temp32a);
-            const cytabttlen = scale_expansion_zeroelim(abttlen, abtt, cdytail, cytabtt);
-            temp16alen = scale_expansion_zeroelim(cytabttlen, cytabtt, 2 * cdy, temp16a);
-            temp16blen = scale_expansion_zeroelim(cytabttlen, cytabtt, cdytail, temp16b);
-            temp32blen = fast_expansion_sum_zeroelim(temp16alen, temp16a, temp16blen, temp16b, temp32b);
-            temp64len = fast_expansion_sum_zeroelim(temp32alen, temp32a, temp32blen, temp32b, temp64);
-            finlength = fast_expansion_sum_zeroelim(finlength, finnow, temp64len, temp64, finother);
-            finswap = finnow; finnow = finother; finother = finswap;
+            const len = scale(abtlen, abt, cdytail, _16c);
+            finlen = finadd(finlen, sum(
+                scale(cytablen, cytab, cdytail, _16$1), _16$1,
+                scale(len, _16c, 2 * cdy, _32), _32, _48$1), _48$1);
+
+            const len2 = scale(abttlen, abtt, cdytail, _8$1);
+            finlen = finadd(finlen, sum_three(
+                scale(len2, _8$1, 2 * cdy, _16$1), _16$1,
+                scale(len2, _8$1, cdytail, _16b), _16b,
+                scale(len, _16c, cdytail, _32), _32, _32b, _64), _64);
         }
     }
 
-    return finnow[finlength - 1];
+    return fin$1[finlen - 1];
 }
 
 function incircle(ax, ay, bx, by, cx, cy, dx, dy) {
@@ -2252,8 +1557,8 @@ const isperrboundA = (16 + 224 * epsilon) * epsilon;
 const isperrboundB = (5 + 72 * epsilon) * epsilon;
 const isperrboundC = (71 + 1408 * epsilon) * epsilon * epsilon;
 
-const ab$2 = vec(4);
-const bc$2 = vec(4);
+const ab = vec(4);
+const bc = vec(4);
 const cd = vec(4);
 const de = vec(4);
 const ea = vec(4);
@@ -2262,9 +1567,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);
@@ -2275,503 +1578,386 @@ 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$2 = vec(1152);
-const bdet$2 = vec(1152);
-const cdet$2 = vec(1152);
+
+const adet = vec(1152);
+const bdet = vec(1152);
+const cdet = vec(1152);
 const ddet = vec(1152);
 const edet = vec(1152);
-const abdet$2 = vec(2304);
+const abdet = vec(2304);
 const cddet = vec(2304);
 const cdedet = vec(3456);
 const deter = vec(5760);
 
+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);
+
+    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 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;
-
-    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0;
-
-    axby1 = ax * by;
+    let bvirt, c, ahi, alo, bhi, blo, _i, _j, _0, s1, s0, t1, t0, u3;
+
+    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$2[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$2[1] = _0 - (_i + bvirt) + (bvirt - bxay1);
-    ab$2[3] = _j + _i;
-    bvirt = ab$2[3] - _j;
-    ab$2[2] = _j - (ab$2[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$2[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$2[1] = _0 - (_i + bvirt) + (bvirt - cxby1);
-    bc$2[3] = _j + _i;
-    bvirt = bc$2[3] - _j;
-    bc$2[2] = _j - (bc$2[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$2, 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$2, 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$2, 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$2, 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$2, 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$2, 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$2);
-
-    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$2);
-
-    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$2);
-
-    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$2, blen, bdet$2, abdet$2);
-    const cdlen = fast_expansion_sum_zeroelim(clen, cdet$2, dlen, ddet, cddet);
-    cdelen = fast_expansion_sum_zeroelim(cdlen, cddet, elen, edet, cdedet);
-    const deterlen = fast_expansion_sum_zeroelim(ablen, abdet$2, 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$1 = vec(48);
 const xdet = vec(96);
 const ydet = vec(96);
 const zdet = vec(96);
-const xydet = vec(192);
-const fin1$2 = 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;
@@ -2786,247 +1972,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$2[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$2[1] = _0 - (_i + bvirt) + (bvirt - bexaey1);
+    ab[1] = _0 - (_i + bvirt) + (bvirt - t1);
     ab3 = _j + _i;
     bvirt = ab3 - _j;
-    ab$2[2] = _j - (ab3 - bvirt) + (_i - bvirt);
-    ab$2[3] = ab3;
-
-    bexcey1 = bex * cey;
+    ab[2] = _j - (ab3 - bvirt) + (_i - bvirt);
+    ab[3] = ab3;
+    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$2[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$2[1] = _0 - (_i + bvirt) + (bvirt - cexbey1);
+    bc[1] = _0 - (_i + bvirt) + (bvirt - t1);
     bc3 = _j + _i;
     bvirt = bc3 - _j;
-    bc$2[2] = _j - (bc3 - bvirt) + (_i - bvirt);
-    bc$2[3] = bc3;
-
-    cexdey1 = cex * dey;
+    bc[2] = _j - (bc3 - bvirt) + (_i - bvirt);
+    bc[3] = bc3;
+    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$2, 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$1);
-    xlen = scale_expansion_zeroelim(temp48len, temp48$1, -aex, xdet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, aey, temp48$1);
-    ylen = scale_expansion_zeroelim(temp48len, temp48$1, -aey, ydet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, aez, temp48$1);
-    zlen = scale_expansion_zeroelim(temp48len, temp48$1, -aez, zdet);
-    xylen = fast_expansion_sum_zeroelim(xlen, xdet, ylen, ydet, xydet);
-    const alen = fast_expansion_sum_zeroelim(xylen, xydet, zlen, zdet, adet$2);
-
-    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$1);
-    xlen = scale_expansion_zeroelim(temp48len, temp48$1, bex, xdet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, bey, temp48$1);
-    ylen = scale_expansion_zeroelim(temp48len, temp48$1, bey, ydet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, bez, temp48$1);
-    zlen = scale_expansion_zeroelim(temp48len, temp48$1, bez, zdet);
-    xylen = fast_expansion_sum_zeroelim(xlen, xdet, ylen, ydet, xydet);
-    const blen = fast_expansion_sum_zeroelim(xylen, xydet, zlen, zdet, bdet$2);
-
-    temp8alen = scale_expansion_zeroelim(4, ab$2, 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$1);
-    xlen = scale_expansion_zeroelim(temp48len, temp48$1, -cex, xdet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, cey, temp48$1);
-    ylen = scale_expansion_zeroelim(temp48len, temp48$1, -cey, ydet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, cez, temp48$1);
-    zlen = scale_expansion_zeroelim(temp48len, temp48$1, -cez, zdet);
-    xylen = fast_expansion_sum_zeroelim(xlen, xdet, ylen, ydet, xydet);
-    const clen = fast_expansion_sum_zeroelim(xylen, xydet, zlen, zdet, cdet$2);
-
-    temp8alen = scale_expansion_zeroelim(4, bc$2, aez, temp8a);
-    temp8blen = scale_expansion_zeroelim(4, ac, -bez, temp8b);
-    temp8clen = scale_expansion_zeroelim(4, ab$2, 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$1);
-    xlen = scale_expansion_zeroelim(temp48len, temp48$1, dex, xdet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, dey, temp48$1);
-    ylen = scale_expansion_zeroelim(temp48len, temp48$1, dey, ydet);
-    temp48len = scale_expansion_zeroelim(temp24len, temp24, dez, temp48$1);
-    zlen = scale_expansion_zeroelim(temp48len, temp48$1, dez, zdet);
-    xylen = fast_expansion_sum_zeroelim(xlen, xdet, ylen, ydet, xydet);
-    const dlen = fast_expansion_sum_zeroelim(xylen, xydet, zlen, zdet, ddet);
-
-    const ablen = fast_expansion_sum_zeroelim(alen, adet$2, blen, bdet$2, abdet$2);
-    const cdlen = fast_expansion_sum_zeroelim(clen, cdet$2, dlen, ddet, cddet);
-    const finlength = fast_expansion_sum_zeroelim(ablen, abdet$2, cdlen, cddet, fin1$2);
-
-    let det = estimate(finlength, fin1$2);
+    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);
+
+    let det = estimate(finlen, fin);
     let errbound = isperrboundB * permanent;
     if (det >= errbound || -det >= errbound) {
         return det;
@@ -3116,7 +2245,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;
@@ -3183,7 +2311,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;
 
@@ -3211,4 +2338,4 @@ exports.orient3dfast = orient3dfast;
 
 Object.defineProperty(exports, '__esModule', { value: true });
 
-}));
+})));