@knowvah/dot-engine 2.1.0 → 2.3.0

package/dist/common/libm-hypot.d.ts

@@ -0,0 +1,3 @@
+/** libm `hypot(a, b)` of the native build host. */
+export declare function libmHypot(a: number, b: number): number;
+//# sourceMappingURL=libm-hypot.d.ts.map

package/dist/common/libm-hypot.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"libm-hypot.d.ts","sourceRoot":"","sources":["../../src/common/libm-hypot.ts"],"names":[],"mappings":"AAyBA,mDAAmD;AACnD,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAQtD"}

package/dist/common/poly-sizing.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"poly-sizing.d.ts","sourceRoot":"","sources":["../../src/common/poly-sizing.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,kBAAkB,CAAC;AAM9C,4FAA4F;AAC5F,eAAO,MAAM,GAAG,IAAI,CAAC;AAUrB;;;GAGG;AACH,MAAM,WAAW,cAAc;IAC7B,oDAAoD;IACpD,UAAU,EAAE,KAAK,CAAC;IAClB,2EAA2E;IAC3E,KAAK,EAAE,MAAM,CAAC;IACd,2DAA2D;IAC3D,WAAW,EAAE,MAAM,CAAC;IACpB,8DAA8D;IAC9D,WAAW,EAAE,MAAM,CAAC;IACpB,yBAAyB;IACzB,UAAU,EAAE,MAAM,CAAC;IACnB,mBAAmB;IACnB,IAAI,EAAE,MAAM,CAAC;IACb,qDAAqD;IACrD,OAAO,EAAE,OAAO,CAAC;IACjB,8DAA8D;IAC9D,OAAO,EAAE,OAAO,CAAC;IACjB,qDAAqD;IACrD,OAAO,EAAE,MAAM,CAAC;IAChB,qDAAqD;IACrD,QAAQ,EAAE,MAAM,CAAC;IACjB,kFAAkF;IAClF,WAAW,EAAE,MAAM,CAAC;IACpB,oEAAoE;IACpE,MAAM,EAAE,MAAM,GAAG,SAAS,CAAC;IAC3B,wCAAwC;IACxC,SAAS,EAAE,MAAM,CAAC;IAClB,2CAA2C;IAC3C,QAAQ,EAAE,MAAM,GAAG,SAAS,CAAC;IAC7B,8DAA8D;IAC9D,SAAS,EAAE,MAAM,CAAC;IAClB,wDAAwD;IACxD,IAAI,EAAE,OAAO,CAAC;IACd,qFAAqF;IACrF,QAAQ,CAAC,EAAE,MAAM,CAAC;IAClB;;;;;OAKG;IACH,WAAW,CAAC,EAAE,MAAM,CAAC;IACrB;;;;;OAKG;IACH,SAAS,CAAC,EAAE,OAAO,CAAC;CACrB;AAED,kEAAkE;AAClE,MAAM,WAAW,QAAQ;IACvB,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;CACZ;AAED,0EAA0E;AAC1E,MAAM,WAAW,cAAe,SAAQ,QAAQ;IAC9C,4EAA4E;IAC5E,QAAQ,EAAE,MAAM,CAAC;IACjB,uDAAuD;IACvD,QAAQ,EAAE,MAAM,CAAC;IACjB;;;;;OAKG;IACH,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;IACd;;;;;;OAMG;IACH,UAAU,EAAE,OAAO,CAAC;IACpB;;;;;OAKG;IACH,KAAK,EAAE,KAAK,CAAC;IACb,+DAA+D;IAC/D,MAAM,EAAE,MAAM,CAAC;CAChB;AAUD;;;;GAIG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,cAAc,GAAG;IAAE,KAAK,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAQnF;AAaD;;;;GAIG;AACH,wBAAgB,aAAa,CAC3B,UAAU,EAAE,KAAK,EACjB,MAAM,EAAE,MAAM,GAAG,SAAS,EAC1B,OAAO,EAAE,OAAO,GACf,KAAK,CAaP;AASD;;;;GAIG;AACH;;;;;;GAMG;AACH,wBAAgB,QAAQ,CAAC,EAAE,EAAE,KAAK,GAAG,KAAK,CAQzC;AAED;;;GAGG;AACH,wBAAgB,YAAY,CAAC,EAAE,EAAE,KAAK,GAAG,KAAK,CAE7C;AAED,wBAAgB,cAAc,CAC5B,EAAE,EAAE,KAAK,EACT,KAAK,EAAE,MAAM,EACb,SAAS,EAAE,MAAM,EACjB,MAAM,EAAE,MAAM,GACb,KAAK,CAgBP;AAqBD,2DAA2D;AAC3D,MAAM,WAAW,QAAQ;IACvB,KAAK,EAAE,MAAM,CAAC;IACd,WAAW,EAAE,MAAM,CAAC;IACpB,UAAU,EAAE,MAAM,CAAC;IACnB,IAAI,EAAE,MAAM,CAAC;CACd;AAiBD;;;;GAIG;AACH,wBAAgB,eAAe,CAC7B,EAAE,EAAE,KAAK,EACT,CAAC,EAAE,QAAQ,EACX,KAAK,EAAE,OAAO,GACb;IAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,MAAM,CAAA;CAAE,CA6BhD;AAuDD;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,KAAK,EAAE,MAAM,GAAG,KAAK,EAAE,CAQzE;AAED;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,KAAK,EAAE,CAU3F;AAED;;;;GAIG;AACH,wBAAgB,kBAAkB,CAChC,KAAK,EAAE,KAAK,EAAE,EACd,KAAK,EAAE,MAAM,EACb,CAAC,EAAE;IAAE,WAAW,EAAE,MAAM,CAAC;IAAC,QAAQ,EAAE,MAAM,CAAC;IAAC,UAAU,EAAE,OAAO,CAAA;CAAE,EACjE,EAAE,EAAE,KAAK,GACR;IAAE,EAAE,EAAE,KAAK,CAAC;IAAC,OAAO,EAAE,KAAK,CAAA;CAAE,CAkB/B;AA6DD,gGAAgG;AAChG,wBAAgB,UAAU,CAAC,QAAQ,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,OAAO,GAAG,QAAQ,CAGvF;AAiCD;;;;GAIG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,cAAc,GAAG,cAAc,CA+B1D"}
+{"version":3,"file":"poly-sizing.d.ts","sourceRoot":"","sources":["../../src/common/poly-sizing.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,kBAAkB,CAAC;AAO9C,4FAA4F;AAC5F,eAAO,MAAM,GAAG,IAAI,CAAC;AAUrB;;;GAGG;AACH,MAAM,WAAW,cAAc;IAC7B,oDAAoD;IACpD,UAAU,EAAE,KAAK,CAAC;IAClB,2EAA2E;IAC3E,KAAK,EAAE,MAAM,CAAC;IACd,2DAA2D;IAC3D,WAAW,EAAE,MAAM,CAAC;IACpB,8DAA8D;IAC9D,WAAW,EAAE,MAAM,CAAC;IACpB,yBAAyB;IACzB,UAAU,EAAE,MAAM,CAAC;IACnB,mBAAmB;IACnB,IAAI,EAAE,MAAM,CAAC;IACb,qDAAqD;IACrD,OAAO,EAAE,OAAO,CAAC;IACjB,8DAA8D;IAC9D,OAAO,EAAE,OAAO,CAAC;IACjB,qDAAqD;IACrD,OAAO,EAAE,MAAM,CAAC;IAChB,qDAAqD;IACrD,QAAQ,EAAE,MAAM,CAAC;IACjB,kFAAkF;IAClF,WAAW,EAAE,MAAM,CAAC;IACpB,oEAAoE;IACpE,MAAM,EAAE,MAAM,GAAG,SAAS,CAAC;IAC3B,wCAAwC;IACxC,SAAS,EAAE,MAAM,CAAC;IAClB,2CAA2C;IAC3C,QAAQ,EAAE,MAAM,GAAG,SAAS,CAAC;IAC7B,8DAA8D;IAC9D,SAAS,EAAE,MAAM,CAAC;IAClB,wDAAwD;IACxD,IAAI,EAAE,OAAO,CAAC;IACd,qFAAqF;IACrF,QAAQ,CAAC,EAAE,MAAM,CAAC;IAClB;;;;;OAKG;IACH,WAAW,CAAC,EAAE,MAAM,CAAC;IACrB;;;;;OAKG;IACH,SAAS,CAAC,EAAE,OAAO,CAAC;CACrB;AAED,kEAAkE;AAClE,MAAM,WAAW,QAAQ;IACvB,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;CACZ;AAED,0EAA0E;AAC1E,MAAM,WAAW,cAAe,SAAQ,QAAQ;IAC9C,4EAA4E;IAC5E,QAAQ,EAAE,MAAM,CAAC;IACjB,uDAAuD;IACvD,QAAQ,EAAE,MAAM,CAAC;IACjB;;;;;OAKG;IACH,KAAK,EAAE,MAAM,CAAC;IACd,KAAK,EAAE,MAAM,CAAC;IACd;;;;;;OAMG;IACH,UAAU,EAAE,OAAO,CAAC;IACpB;;;;;OAKG;IACH,KAAK,EAAE,KAAK,CAAC;IACb,+DAA+D;IAC/D,MAAM,EAAE,MAAM,CAAC;CAChB;AAUD;;;;GAIG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,cAAc,GAAG;IAAE,KAAK,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAQnF;AAaD;;;;GAIG;AACH,wBAAgB,aAAa,CAC3B,UAAU,EAAE,KAAK,EACjB,MAAM,EAAE,MAAM,GAAG,SAAS,EAC1B,OAAO,EAAE,OAAO,GACf,KAAK,CAaP;AASD;;;;GAIG;AACH;;;;;;GAMG;AACH,wBAAgB,QAAQ,CAAC,EAAE,EAAE,KAAK,GAAG,KAAK,CAQzC;AAED;;;GAGG;AACH,wBAAgB,YAAY,CAAC,EAAE,EAAE,KAAK,GAAG,KAAK,CAE7C;AAED,wBAAgB,cAAc,CAC5B,EAAE,EAAE,KAAK,EACT,KAAK,EAAE,MAAM,EACb,SAAS,EAAE,MAAM,EACjB,MAAM,EAAE,MAAM,GACb,KAAK,CAgBP;AAqBD,2DAA2D;AAC3D,MAAM,WAAW,QAAQ;IACvB,KAAK,EAAE,MAAM,CAAC;IACd,WAAW,EAAE,MAAM,CAAC;IACpB,UAAU,EAAE,MAAM,CAAC;IACnB,IAAI,EAAE,MAAM,CAAC;CACd;AAgBD;;;;GAIG;AACH,wBAAgB,eAAe,CAC7B,EAAE,EAAE,KAAK,EACT,CAAC,EAAE,QAAQ,EACX,KAAK,EAAE,OAAO,GACb;IAAE,KAAK,EAAE,KAAK,EAAE,CAAC;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,MAAM,CAAA;CAAE,CA6BhD;AAuDD;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,KAAK,EAAE,MAAM,GAAG,KAAK,EAAE,CAQzE;AAED;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,KAAK,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,GAAG,KAAK,EAAE,CAU3F;AAED;;;;GAIG;AACH,wBAAgB,kBAAkB,CAChC,KAAK,EAAE,KAAK,EAAE,EACd,KAAK,EAAE,MAAM,EACb,CAAC,EAAE;IAAE,WAAW,EAAE,MAAM,CAAC;IAAC,QAAQ,EAAE,MAAM,CAAC;IAAC,UAAU,EAAE,OAAO,CAAA;CAAE,EACjE,EAAE,EAAE,KAAK,GACR;IAAE,EAAE,EAAE,KAAK,CAAC;IAAC,OAAO,EAAE,KAAK,CAAA;CAAE,CAkB/B;AA6DD,gGAAgG;AAChG,wBAAgB,UAAU,CAAC,QAAQ,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,OAAO,GAAG,QAAQ,CAGvF;AAiCD;;;;GAIG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,cAAc,GAAG,cAAc,CA+B1D"}

package/dist/index.js

@@ -16909,6 +16909,24 @@ function RADIANS(deg) {
   return deg / 180 * Math.PI;
 }
 
+// src/common/libm-hypot.ts
+var BIG = 2 ** 500;
+var SMALL = 2 ** -500;
+var SCALE_UP = 2 ** 600;
+var SCALE_DOWN = 2 ** -600;
+function libmHypot(a, b) {
+  if (a === Infinity || a === -Infinity || b === Infinity || b === -Infinity) return Infinity;
+  let x = Math.abs(a);
+  let y = Math.abs(b);
+  if (x < y) [x, y] = [y, x];
+  if (x > BIG) return core(x * SCALE_DOWN, y * SCALE_DOWN) * SCALE_UP;
+  if (x < SMALL) return core(x * SCALE_UP, y * SCALE_UP) * SCALE_DOWN;
+  return core(x, y);
+}
+function core(x, y) {
+  return Math.sqrt(fma(y, y, x * x));
+}
+
 // src/common/poly-sizing.ts
 var GAP = 4;
 var XPAD = 4 * GAP;
@@ -16996,16 +17014,16 @@ function ellipsePeripheryBB(bb, peripheries, penwidth, hasOutline) {
   return { bb: b, outline };
 }
 function transformUnitVertex(R2, g, c, bb) {
-  const D = { x: R2.x * (c.skewdist + R2.y * c.gdistortion) + R2.y * c.gskew, y: R2.y };
+  const D = { x: fma(R2.x, fma(R2.y, c.gdistortion, c.skewdist), R2.y * c.gskew), y: R2.y };
   const alpha = RADIANS(g.orientation) + Math.atan2(D.y, D.x);
-  const r = Math.sqrt(D.x * D.x + D.y * D.y);
+  const r = libmHypot(D.x, D.y);
   return { x: r * Math.cos(alpha) * bb.x, y: r * Math.sin(alpha) * bb.y };
 }
 function polygonVertices(bb, g, isBox3) {
   const sectorangle = 2 * Math.PI / g.sides;
   const sidelength = Math.sin(sectorangle / 2);
   const c = {
-    skewdist: Math.hypot(Math.abs(g.distortion) + Math.abs(g.skew), 1),
+    skewdist: libmHypot(Math.abs(g.distortion) + Math.abs(g.skew), 1),
     gdistortion: g.distortion * SQRT2 / Math.cos(sectorangle / 2),
     gskew: g.skew / 2
   };
@@ -17430,13 +17448,13 @@ function evalBezier(sp, t, left2, right2) {
   return pt;
 }
 var HW = 2;
-var SMALL = 1e-10;
+var SMALL2 = 1e-10;
 function ptToLine2(a, b, p) {
   const dx = b.x - a.x;
   const dy = b.y - a.y;
   let a2 = (p.y - a.y) * dx - (p.x - a.x) * dy;
   a2 *= a2;
-  if (a2 < SMALL) return 0;
+  if (a2 < SMALL2) return 0;
   return a2 / (dx * dx + dy * dy);
 }
 function checkControlPoints(cp) {
@@ -19979,7 +19997,7 @@ function feasibleTree(ctx) {
 
 // src/util/bsd-qsort.ts
 var THRESH = 7;
-var BIG = 40;
+var BIG2 = 40;
 function swap(arr, i, j) {
   const t = arr[i];
   arr[i] = arr[j];
@@ -20005,7 +20023,7 @@ function choosePivot(arr, lo, n, cmp) {
   let pm = lo + (n >> 1);
   let pl = lo;
   let pn = lo + (n - 1);
-  if (n > BIG) {
+  if (n > BIG2) {
     const d = n >> 3;
     pl = med3(arr, pl, pl + d, pl + 2 * d, cmp);
     pm = med3(arr, pm - d, pm, pm + d, cmp);
@@ -35806,7 +35824,7 @@ function addLoop(sg, cp, dp, sp) {
 var MARGIN = 36;
 var delta = 1;
 var mu = 500;
-var BIG2 = 16384;
+var BIG3 = 16384;
 function chanSz(w) {
   return (w - 3) / 2;
 }
@@ -35826,7 +35844,7 @@ function updateWt(ep, sz) {
   ep.cnt++;
   if (ep.cnt > sz) {
     ep.cnt = 0;
-    ep.weight += BIG2;
+    ep.weight += BIG3;
   }
 }
 function updateWts(g, cp, ep) {
@@ -35884,12 +35902,12 @@ function createSEdges(cp, g) {
   let vwt = delta * (bb.UR.y - bb.LL.y);
   let wt = (hwt + vwt) / 2 + mu;
   if (isSmall(bb.UR.y - bb.LL.y) && !(cp.flags & MZ_SMALLV)) {
-    hwt = BIG2;
-    wt = BIG2;
+    hwt = BIG3;
+    wt = BIG3;
   }
   if (isSmall(bb.UR.x - bb.LL.x) && !(cp.flags & MZ_SMALLH)) {
-    vwt = BIG2;
-    wt = BIG2;
+    vwt = BIG3;
+    wt = BIG3;
   }
   const L2 = cp.sides[M_LEFT];
   const T2 = cp.sides[M_TOP];
@@ -41494,7 +41512,7 @@ function raySeg2(v, w, a, b) {
   }
   return wind(v, a, w) * wind(v, a, b) >= 0;
 }
-var SMALL2 = 1e-10;
+var SMALL3 = 1e-10;
 function raySegIntersect(v, w, a, b) {
   if (!raySeg2(v, w, a, b)) return null;
   const mv = { x: w.x - v.x, y: w.y - v.y };
@@ -41502,7 +41520,7 @@ function raySegIntersect(v, w, a, b) {
   const ln = { x: -lv.y, y: lv.x };
   const lc = -(ln.x * a.x + ln.y * a.y);
   const dt = ln.x * mv.x + ln.y * mv.y;
-  if (Math.abs(dt) < SMALL2) return null;
+  if (Math.abs(dt) < SMALL3) return null;
   const s = (ln.x * v.x + ln.y * v.y + lc) / dt;
   return { x: v.x - s * mv.x, y: v.y - s * mv.y };
 }
@@ -48218,7 +48236,10 @@ var fdpParms = {
   C: 0,
   Tfact: 1,
   K: -1,
-  T0: -1
+  T0: -1,
+  /** Distance beyond which repulsion is 0; HUGE_VAL unless `dot -Lm` (no library
+   *  equivalent). @see lib/fdpgen/fdp.h:fdpParms_s.Mlimit */
+  Mlimit: Infinity
 };
 var parms = {
   useGrid: 0,
@@ -48546,8 +48567,8 @@ function placeByNeighbors(g, np, ctr) {
     np.info.pos[0] = p.x;
     np.info.pos[1] = p.y;
   } else if (cnt === 1) {
-    np.info.pos[0] = 0.98 * p.x + 0.1 * ctr.x;
-    np.info.pos[1] = 0.9 * p.y + 0.1 * ctr.y;
+    np.info.pos[0] = fma(0.98, p.x, 0.1 * ctr.x);
+    np.info.pos[1] = fma(0.9, p.y, 0.1 * ctr.y);
   } else {
     const angle = Math.PI * 2 * drand48();
     const radius = 0.9 * drand48();
@@ -48563,8 +48584,8 @@ function averagePlacedNeighbors(g, np, p) {
     const op = ep.head === np ? ep.tail : ep.head;
     if (!hasPos3(op)) continue;
     if (cnt) {
-      p.x = (p.x * cnt + op.info.pos[0]) / (cnt + 1);
-      p.y = (p.y * cnt + op.info.pos[1]) / (cnt + 1);
+      p.x = fma(p.x, cnt, op.info.pos[0]) / (cnt + 1);
+      p.y = fma(p.y, cnt, op.info.pos[1]) / (cnt + 1);
     } else {
       p.x = op.info.pos[0];
       p.y = op.info.pos[1];
@@ -48626,8 +48647,8 @@ function sizeBoundary(width, height) {
 function placeWithPorts(g, pp, ctr) {
   for (const port of pp) {
     const np = port.n;
-    np.info.pos[0] = parms.Wd * Math.cos(port.alpha) + ctr.x;
-    np.info.pos[1] = parms.Ht * Math.sin(port.alpha) + ctr.y;
+    np.info.pos[0] = fma(parms.Wd, Math.cos(port.alpha), ctr.x);
+    np.info.pos[1] = fma(parms.Ht, Math.sin(port.alpha), ctr.y);
     dndata(np).pinned = P_SET;
   }
   for (const np of g.nodes.values()) {
@@ -48689,19 +48710,20 @@ function coincidentDelta() {
   } while (xdelta === 0 && ydelta === 0);
   return { xdelta, ydelta };
 }
-function doRep(p, q, xdelta, ydelta, dist22) {
+function doRep(p, q, xdelta, ydelta, dist4) {
   let force;
-  if (dist22 === 0) {
+  if (!(dist4 > 0)) {
     const d = coincidentDelta();
     xdelta = d.xdelta;
     ydelta = d.ydelta;
-    dist22 = fma(xdelta, xdelta, ydelta * ydelta);
+    dist4 = libmHypot(xdelta, ydelta);
   }
-  if (parms.useNew) {
-    const dist4 = Math.sqrt(dist22);
-    force = parms.K * parms.K / (dist4 * dist22);
+  if (dist4 > fdpParms.Mlimit) {
+    force = 0;
+  } else if (parms.useNew) {
+    force = parms.K * parms.K / (dist4 * dist4 * dist4);
   } else {
-    force = parms.K * parms.K / dist22;
+    force = parms.K * parms.K / (dist4 * dist4);
   }
   if (isPort(p) && isPort(q)) force *= 10;
   const dq = disp(q);
@@ -48714,7 +48736,7 @@ function doRep(p, q, xdelta, ydelta, dist22) {
 function applyRep(p, q) {
   const xdelta = q.info.pos[0] - p.info.pos[0];
   const ydelta = q.info.pos[1] - p.info.pos[1];
-  doRep(p, q, xdelta, ydelta, fma(xdelta, xdelta, ydelta * ydelta));
+  doRep(p, q, xdelta, ydelta, libmHypot(xdelta, ydelta));
 }
 function doNeighbor(grid2, i, j, nodes) {
   const cellp = grid2.find(i, j);
@@ -48723,8 +48745,8 @@ function doNeighbor(grid2, i, j, nodes) {
     for (const q of cellp.nodes) {
       const xdelta = q.info.pos[0] - p.info.pos[0];
       const ydelta = q.info.pos[1] - p.info.pos[1];
-      const dist22 = fma(xdelta, xdelta, ydelta * ydelta);
-      if (dist22 < parms.Cell * parms.Cell) doRep(p, q, xdelta, ydelta, dist22);
+      const dist4 = libmHypot(xdelta, ydelta);
+      if (dist4 < parms.Cell) doRep(p, q, xdelta, ydelta, dist4);
     }
   }
 }
@@ -48896,7 +48918,7 @@ function ht2(n) {
   return xMarg.doAdd ? n.info.height / 2 + xMarg.y : n.info.height * xMarg.y / 2;
 }
 function rad(n) {
-  return Math.hypot(wd2(n), ht2(n));
+  return libmHypot(wd2(n), ht2(n));
 }
 function xinitParams(n, xpms) {
   xParams.K = xpms.K;
@@ -48931,37 +48953,38 @@ function cntOverlaps(g) {
 function applyRep2(p, q, rs) {
   let xdelta = q.info.pos[0] - p.info.pos[0];
   let ydelta = q.info.pos[1] - p.info.pos[1];
-  let dist22 = xdelta * xdelta + ydelta * ydelta;
-  if (dist22 === 0) {
+  let dist4 = libmHypot(xdelta, ydelta);
+  if (!(dist4 > 0)) {
     const d = coincidentDelta();
     xdelta = d.xdelta;
     ydelta = d.ydelta;
-    dist22 = xdelta * xdelta + ydelta * ydelta;
+    dist4 = libmHypot(xdelta, ydelta);
   }
   const ov = overlap4(p, q) ? 1 : 0;
-  const force = ov ? rs.ov / dist22 : rs.nonov / dist22;
+  let force = ov ? rs.ov / (dist4 * dist4) : rs.nonov / (dist4 * dist4);
+  if (dist4 > fdpParms.Mlimit) force = 0;
   const dq = disp(q);
   const dp = disp(p);
-  dq[0] += xdelta * force;
-  dq[1] += ydelta * force;
-  dp[0] -= xdelta * force;
-  dp[1] -= ydelta * force;
+  dq[0] = fma(xdelta, force, dq[0]);
+  dq[1] = fma(ydelta, force, dq[1]);
+  dp[0] = fms(xdelta, force, dp[0]);
+  dp[1] = fms(ydelta, force, dp[1]);
   return ov;
 }
 function applyAttr2(p, q) {
   if (overlap4(p, q)) return;
   const xdelta = q.info.pos[0] - p.info.pos[0];
   const ydelta = q.info.pos[1] - p.info.pos[1];
-  const dist4 = Math.hypot(xdelta, ydelta);
+  const dist4 = libmHypot(xdelta, ydelta);
   const din = rad(p) + rad(q);
   const dout = dist4 - din;
   const force = dout * dout / ((xParams.K + din) * dist4);
   const dq = disp(q);
   const dp = disp(p);
-  dq[0] -= xdelta * force;
-  dq[1] -= ydelta * force;
-  dp[0] += xdelta * force;
-  dp[1] += ydelta * force;
+  dq[0] = fms(xdelta, force, dq[0]);
+  dq[1] = fms(ydelta, force, dq[1]);
+  dp[0] = fma(xdelta, force, dp[0]);
+  dp[1] = fma(ydelta, force, dp[1]);
 }
 function accumForces(g, rs) {
   let overlaps2 = 0;
@@ -48983,7 +49006,7 @@ function moveNodes(g, temp) {
     if (dndata(n).pinned === P_PIN) continue;
     const dx = disp(n)[0];
     const dy = disp(n)[1];
-    const len2 = dx * dx + dy * dy;
+    const len2 = fma(dx, dx, dy * dy);
     if (len2 < temp2) {
       n.info.pos[0] += dx;
       n.info.pos[1] += dy;
@@ -49237,7 +49260,7 @@ function deriveClusterNodes(g, dg, infop) {
   }
 }
 function deriveRealNodes(g, dg) {
-  for (const n of g.nodes.values()) {
+  for (const n of nodesInSeq(g)) {
     if (getDnode(n) !== null) continue;
     const par = getParent(n);
     if (par !== null && par !== gdata(g).parent) {
@@ -49279,7 +49302,7 @@ function deriveEdge(dg, byPair, tl, hd) {
   return de;
 }
 function deriveEdges(g, dg, byPair) {
-  for (const n of g.nodes.values()) {
+  for (const n of nodesInSeq(g)) {
     const tl = getDnode(n);
     for (const e of n.outEdges(g)) {
       const hd = getDnode(e.head);