@knowvah/dot-engine 2.1.0 → 2.2.0

package/dist/layout/fdp/derive.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"derive.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/derive.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAC7C,OAAO,EAAE,IAAI,EAAE,MAAM,qBAAqB,CAAC;AAwB3C,oEAAoE;AACpE,MAAM,WAAW,UAAU;IACzB,oDAAoD;IACpD,KAAK,EAAE,KAAK,CAAC;IACb,kEAAkE;IAClE,SAAS,EAAE,OAAO,CAAC;IACnB,GAAG,EAAE,MAAM,CAAC;IACZ,6CAA6C;IAC7C,IAAI,EAAE,OAAO,kBAAkB,EAAE,QAAQ,CAAC;CAC3C;AAMD;;;GAGG;AACH,wBAAgB,YAAY,CAAC,EAAE,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM,GAAG,IAAI,CAM1D;AAsLD;;;;GAIG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,EAAE,UAAU,GAAG,KAAK,GAAG,IAAI,CAiBrE"}
+{"version":3,"file":"derive.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/derive.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAC7C,OAAO,EAAE,IAAI,EAAE,MAAM,qBAAqB,CAAC;AAyB3C,oEAAoE;AACpE,MAAM,WAAW,UAAU;IACzB,oDAAoD;IACpD,KAAK,EAAE,KAAK,CAAC;IACb,kEAAkE;IAClE,SAAS,EAAE,OAAO,CAAC;IACnB,GAAG,EAAE,MAAM,CAAC;IACZ,6CAA6C;IAC7C,IAAI,EAAE,OAAO,kBAAkB,EAAE,QAAQ,CAAC;CAC3C;AAMD;;;GAGG;AACH,wBAAgB,YAAY,CAAC,EAAE,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM,GAAG,IAAI,CAM1D;AAwLD;;;;GAIG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,KAAK,EAAE,KAAK,EAAE,UAAU,GAAG,KAAK,GAAG,IAAI,CAiBrE"}

package/dist/layout/fdp/fdp-model.d.ts

@@ -2,11 +2,11 @@
  * fdp engine data model — gdata, dndata, bport_t, xparams, and the
  * accessor macros of the C private interface.
  *
- * Spec read at the 15.0.0 tag (post-tag fdpgen commits add an Mlimit
- * branch and float reorderings the golden refs do not have).
+ * Matches Graphviz after 15.0.0 (fdpParms.Mlimit lives in
+ * tlayout-parms.ts).
  *
- * @see lib/fdpgen/fdp.h (15.0.0)
- * @see lib/fdpgen/xlayout.h (15.0.0)
+ * @see lib/fdpgen/fdp.h
+ * @see lib/fdpgen/xlayout.h
  */
 import type { Graph } from '../../model/graph.js';
 import type { Node } from '../../model/node.js';

package/dist/layout/fdp/tlayout-init.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"tlayout-init.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/tlayout-init.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;GAWG;AAEH,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAGlD,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAEjD,OAAO,EACL,KAAK,KAAK,EAMX,MAAM,gBAAgB,CAAC;AA4KxB;;;GAGG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,GAAG,IAAI,GAAG,KAAK,CA0BjE"}
+{"version":3,"file":"tlayout-init.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/tlayout-init.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;GAWG;AAEH,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAGlD,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,qBAAqB,CAAC;AAGjD,OAAO,EACL,KAAK,KAAK,EAMX,MAAM,gBAAgB,CAAC;AA6KxB;;;GAGG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,GAAG,IAAI,GAAG,KAAK,CA0BjE"}

package/dist/layout/fdp/tlayout-parms.d.ts

@@ -7,7 +7,7 @@
  *
  * @see lib/fdpgen/tlayout.c (15.0.0) — parms_t, fdp_initParams,
  *      init_params, cool, reset_params
- * @see lib/common/globals.c (15.0.0) — fdpParms defaults
+ * @see lib/common/globals.c — fdpParms defaults
  */
 import type { Graph } from '../../model/graph.js';
 import { type XParams } from './fdp-model.js';
@@ -36,6 +36,9 @@ export declare const fdpParms: {
     Tfact: number;
     K: number;
     T0: number;
+    /** Distance beyond which repulsion is 0; HUGE_VAL unless `dot -Lm` (no library
+     *  equivalent). @see lib/fdpgen/fdp.h:fdpParms_s.Mlimit */
+    Mlimit: number;
 };
 /** @see lib/fdpgen/tlayout.c:parms_t */
 export interface TParms {

package/dist/layout/fdp/tlayout-parms.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"tlayout-parms.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/tlayout-parms.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;GAUG;AAEH,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAElD,OAAO,EAAE,KAAK,OAAO,EAAc,MAAM,gBAAgB,CAAC;AAM1D,0CAA0C;AAC1C,eAAO,MAAM,SAAS,MAAM,CAAC;AAU7B,0CAA0C;AAC1C,eAAO,MAAM,SAAS,IAAI,CAAC;AAC3B,6CAA6C;AAC7C,eAAO,MAAM,YAAY,IAAI,CAAC;AAC9B,4CAA4C;AAC5C,eAAO,MAAM,WAAW,IAAI,CAAC;AAQ7B;;;;;;;GAOG;AACH,eAAO,MAAM,QAAQ;;;;;;;;;CASpB,CAAC;AAMF,wCAAwC;AACxC,MAAM,WAAW,MAAM;IACrB,OAAO,EAAE,MAAM,CAAC;IAChB,MAAM,EAAE,MAAM,CAAC;IACf,IAAI,EAAE,MAAM,CAAC;IACb,QAAQ,EAAE,MAAM,CAAC;IACjB,QAAQ,EAAE,MAAM,CAAC;IACjB,QAAQ,EAAE,MAAM,CAAC;IACjB,CAAC,EAAE,MAAM,CAAC;IACV,KAAK,EAAE,MAAM,CAAC;IACd,CAAC,EAAE,MAAM,CAAC;IACV,EAAE,EAAE,MAAM,CAAC;IACX,KAAK,EAAE,MAAM,CAAC;IACd,IAAI,EAAE,MAAM,CAAC;IACb,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,KAAK,EAAE,MAAM,CAAC;IACd,OAAO,EAAE,MAAM,CAAC;CACjB;AAED,+EAA+E;AAC/E,eAAO,MAAM,KAAK,EAAE,MAInB,CAAC;AAEF,8DAA8D;AAC9D,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEtC;AAED,6CAA6C;AAC7C,wBAAgB,WAAW,IAAI,IAAI,CAElC;AAMD;;;;;;;;;GASG;AACH,wBAAgB,UAAU,CACxB,CAAC,EAAE,KAAK,EACR,IAAI,EAAE,MAAM,EACZ,OAAO,EAAE;IAAE,KAAK,EAAE,MAAM,CAAA;CAAE,GACzB,MAAM,CAaR;AAED;;;GAGG;AACH,wBAAgB,aAAa,CAC3B,CAAC,EAAE,MAAM,EACT,IAAI,EAAE,MAAM,GACX;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,MAAM,CAAA;CAAE,CAOhC;AAMD;;;GAGG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,KAAK,GAAG,IAAI,CAsB5C;AAeD;;;;;GAKG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,OAAO,GAAG,OAAO,CA2B3D"}
+{"version":3,"file":"tlayout-parms.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/tlayout-parms.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;GAUG;AAEH,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAElD,OAAO,EAAE,KAAK,OAAO,EAAc,MAAM,gBAAgB,CAAC;AAM1D,0CAA0C;AAC1C,eAAO,MAAM,SAAS,MAAM,CAAC;AAU7B,0CAA0C;AAC1C,eAAO,MAAM,SAAS,IAAI,CAAC;AAC3B,6CAA6C;AAC7C,eAAO,MAAM,YAAY,IAAI,CAAC;AAC9B,4CAA4C;AAC5C,eAAO,MAAM,WAAW,IAAI,CAAC;AAQ7B;;;;;;;GAOG;AACH,eAAO,MAAM,QAAQ;;;;;;;;;IASnB;+DAC2D;;CAE5D,CAAC;AAMF,wCAAwC;AACxC,MAAM,WAAW,MAAM;IACrB,OAAO,EAAE,MAAM,CAAC;IAChB,MAAM,EAAE,MAAM,CAAC;IACf,IAAI,EAAE,MAAM,CAAC;IACb,QAAQ,EAAE,MAAM,CAAC;IACjB,QAAQ,EAAE,MAAM,CAAC;IACjB,QAAQ,EAAE,MAAM,CAAC;IACjB,CAAC,EAAE,MAAM,CAAC;IACV,KAAK,EAAE,MAAM,CAAC;IACd,CAAC,EAAE,MAAM,CAAC;IACV,EAAE,EAAE,MAAM,CAAC;IACX,KAAK,EAAE,MAAM,CAAC;IACd,IAAI,EAAE,MAAM,CAAC;IACb,EAAE,EAAE,MAAM,CAAC;IACX,EAAE,EAAE,MAAM,CAAC;IACX,KAAK,EAAE,MAAM,CAAC;IACd,OAAO,EAAE,MAAM,CAAC;CACjB;AAED,+EAA+E;AAC/E,eAAO,MAAM,KAAK,EAAE,MAInB,CAAC;AAEF,8DAA8D;AAC9D,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEtC;AAED,6CAA6C;AAC7C,wBAAgB,WAAW,IAAI,IAAI,CAElC;AAMD;;;;;;;;;GASG;AACH,wBAAgB,UAAU,CACxB,CAAC,EAAE,KAAK,EACR,IAAI,EAAE,MAAM,EACZ,OAAO,EAAE;IAAE,KAAK,EAAE,MAAM,CAAA;CAAE,GACzB,MAAM,CAaR;AAED;;;GAGG;AACH,wBAAgB,aAAa,CAC3B,CAAC,EAAE,MAAM,EACT,IAAI,EAAE,MAAM,GACX;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,MAAM,CAAA;CAAE,CAOhC;AAMD;;;GAGG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,KAAK,GAAG,IAAI,CAsB5C;AAeD;;;;;GAKG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,OAAO,GAAG,OAAO,CA2B3D"}

package/dist/layout/fdp/tlayout.d.ts

@@ -2,9 +2,8 @@
  * fdp initial layout — grid-accelerated Fruchterman–Reingold on point
  * nodes with boundary ports.
  *
- * Spec read at the 15.0.0 tag: the post-tag tree reorders doRep floats
- * (hypot vs sqrt(x²·y²)) and adds an Mlimit cutoff; the golden refs are
- * 15.0.0 output, so neither is ported.
+ * Repulsion follows Graphviz after 15.0.0: doRep takes the hypot distance
+ * (the host libm's, see libm-hypot.ts) and honours the Mlimit cutoff.
  *
  * All math is double precision (unlike neato's float32 stress kernel).
  * Force-accumulation order is load-bearing: node iteration is the
@@ -14,7 +13,7 @@
  * Parameter state lives in tlayout-parms.ts; initial placement in
  * tlayout-init.ts.
  *
- * @see lib/fdpgen/tlayout.c (15.0.0)
+ * @see lib/fdpgen/tlayout.c
  */
 import type { Graph } from '../../model/graph.js';
 import { type XParams } from './fdp-model.js';

package/dist/layout/fdp/tlayout.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"tlayout.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/tlayout.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;;GAiBG;AAEH,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAKlD,OAAO,EAEL,KAAK,OAAO,EAMb,MAAM,gBAAgB,CAAC;AAKxB,OAAO,EAAE,aAAa,EAAE,QAAQ,EAAE,UAAU,EAAE,WAAW,EAAE,YAAY,EAAE,SAAS,EAAE,MAAM,oBAAoB,CAAC;AAM/G;;;;;;;;;;;GAWG;AACH,wBAAgB,eAAe,IAAI;IAAE,MAAM,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAQpE;AA8ND;;;GAGG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,OAAO,GAAG,IAAI,CAyBxD"}
+{"version":3,"file":"tlayout.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/tlayout.ts"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAMlD,OAAO,EAEL,KAAK,OAAO,EAMb,MAAM,gBAAgB,CAAC;AAKxB,OAAO,EAAE,aAAa,EAAE,QAAQ,EAAE,UAAU,EAAE,WAAW,EAAE,YAAY,EAAE,SAAS,EAAE,MAAM,oBAAoB,CAAC;AAM/G;;;;;;;;;;;GAWG;AACH,wBAAgB,eAAe,IAAI;IAAE,MAAM,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAQpE;AA+ND;;;GAGG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,OAAO,GAAG,IAAI,CAyBxD"}

package/dist/layout/fdp/xlayout.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"xlayout.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/xlayout.ts"],"names":[],"mappings":"AAoBA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAKlD,OAAO,EACL,KAAK,OAAO,EAKb,MAAM,gBAAgB,CAAC;AAuQxB;;;;;GAKG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,OAAO,GAAG,IAAI,CAMxD"}
+{"version":3,"file":"xlayout.d.ts","sourceRoot":"","sources":["../../../src/layout/fdp/xlayout.ts"],"names":[],"mappings":"AAoBA,OAAO,KAAK,EAAE,KAAK,EAAE,MAAM,sBAAsB,CAAC;AAKlD,OAAO,EACL,KAAK,OAAO,EAKb,MAAM,gBAAgB,CAAC;AA2QxB;;;;;GAKG;AACH,wBAAgB,UAAU,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,OAAO,GAAG,IAAI,CAMxD"}

package/dist/render.js

@@ -10513,6 +10513,24 @@ function RADIANS(deg) {
   return deg / 180 * Math.PI;
 }
 
+// src/common/libm-hypot.ts
+var BIG = 2 ** 500;
+var SMALL2 = 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 < SMALL2) 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;
@@ -10600,16 +10618,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
   };
@@ -21484,7 +21502,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];
@@ -21510,7 +21528,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);
@@ -31834,7 +31852,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;
 }
@@ -31854,7 +31872,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) {
@@ -31912,12 +31930,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];
@@ -37522,7 +37540,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 };
@@ -37530,7 +37548,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 };
 }
@@ -44246,7 +44264,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,
@@ -44574,8 +44595,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();
@@ -44591,8 +44612,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];
@@ -44654,8 +44675,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()) {
@@ -44717,19 +44738,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);
@@ -44742,7 +44764,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);
@@ -44751,8 +44773,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);
     }
   }
 }
@@ -44924,7 +44946,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;
@@ -44959,37 +44981,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;
@@ -45011,7 +45034,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;
@@ -45265,7 +45288,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) {
@@ -45307,7 +45330,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);