@woosh/meep-engine 3.43.0 → 3.45.0

package/src/engine/physics/cloth/solver/cloth_worst_tether_ratio.js

@@ -0,0 +1,56 @@
+import { v3_length } from "../../../../core/geom/vec3/v3_length.js";
+import { max2_exact } from "../../../../core/math/max2_exact.js";
+
+/**
+ * How far the cloth has stretched away from its root, as the worst ratio over
+ * nodes of straight-line distance from node 0 to the fabric between them:
+ * `max ‖x_node − x_0‖ / d_geo(node)` over nodes `1 .. node_count − 1` — §7 P2's
+ * measure. At most 1 for a cloth that has not over-stretched anywhere, since no
+ * straight line is longer than the fabric it spans; `1.05` is a node five
+ * percent further out than the cloth reaches.
+ *
+ * Everything is measured from node 0, the root of the seeder's tree, which is
+ * where the rigs this was written for are held. A cloth held in several places
+ * wants each node measured from its own hold instead, and this is not that
+ * measure.
+ *
+ * A node at a NaN position reads NaN and so does the answer, whatever the other
+ * nodes say, so a broken state cannot pass for a well-behaved one.
+ *
+ * Takes the arrays rather than a `ClothInstance` so the solver side does not
+ * depend on the ECS one.
+ *
+ * @param {Float32Array|Float64Array} position 3 per particle
+ * @param {Int32Array|Uint32Array|number[]} particle_of_node particle index per node
+ * @param {Float32Array|Float64Array|number[]} geodesic fabric distance from node 0 per node, metres
+ * @param {number} node_count
+ * @returns {number} dimensionless; 0 for a cloth of fewer than two nodes
+ *
+ * @author Alex Goldring
+ * @copyright Company Named Limited (c) 2026
+ */
+export function cloth_worst_tether_ratio(position, particle_of_node, geodesic, node_count) {
+    const root = particle_of_node[0] * 3;
+
+    const root_x = position[root];
+    const root_y = position[root + 1];
+    const root_z = position[root + 2];
+
+    let worst = 0;
+
+    for (let node = 1; node < node_count; node++) {
+        const p = particle_of_node[node] * 3;
+
+        const distance = v3_length(
+            position[p] - root_x,
+            position[p + 1] - root_y,
+            position[p + 2] - root_z
+        );
+
+        const ratio = distance / geodesic[node];
+
+        worst = max2_exact(worst, ratio);
+    }
+
+    return worst;
+}

package/src/engine/physics/cloth/solver/constraint/cloth_contact_row.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"cloth_contact_row.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/physics/cloth/solver/constraint/cloth_contact_row.js"],"names":[],"mappings":"AAuHA;;;;;;;;GAQG;AACH,sDANW,YAAY,+BAEZ,MAAM,SACN,YAAY,GAAC,YAAY,QA+LnC;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,mEANW,MAAM,SACN,MAAM,eACN,MAAM,eACN,MAAM,qBACN,MAAM,QA8ChB;AAED;;;;;;;;;;;;;;;GAeG;AACH,oEALW,YAAY,GAAC,YAAY,QACzB,MAAM,eACN,MAAM,qBACN,MAAM,QAqGhB;AAxdD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2DG;AAEH;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,kCAFU,MAAM,CAEqB;AAErC;;;;;GAKG;AACH,qCAFU,MAAM,CAE2B"}
+{"version":3,"file":"cloth_contact_row.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/physics/cloth/solver/constraint/cloth_contact_row.js"],"names":[],"mappings":"AAwHA;;;;;;;;GAQG;AACH,sDANW,YAAY,+BAEZ,MAAM,SACN,YAAY,GAAC,YAAY,QA+LnC;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,mEANW,MAAM,SACN,MAAM,eACN,MAAM,eACN,MAAM,qBACN,MAAM,QA8ChB;AAED;;;;;;;;;;;;;;;GAeG;AACH,oEALW,YAAY,GAAC,YAAY,QACzB,MAAM,eACN,MAAM,qBACN,MAAM,QAqGhB;AAxdD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2DG;AAEH;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,kCAFU,MAAM,CAEqB;AAErC;;;;;GAKG;AACH,qCAFU,MAAM,CAE2B"}

package/src/engine/physics/cloth/solver/constraint/cloth_contact_row.js

@@ -1,3 +1,4 @@
+import { clamp } from "../../../../../core/math/clamp.js";
 import {
     CCR_FEATURE,
     CCR_FRICTION,
@@ -373,8 +374,8 @@ export function cloth_contact_warm_start(state, alpha, gamma, k_min_ratio, k_max
             lambda[slot] *= decay;
             friction_lambda[slot] *= decay;
 
-            penalty[slot] = clamp_band(penalty[slot] * gamma, lo, hi);
-            friction_penalty[slot] = clamp_band(friction_penalty[slot] * gamma, lo, hi);
+            penalty[slot] = clamp(penalty[slot] * gamma, lo, hi);
+            friction_penalty[slot] = clamp(friction_penalty[slot] * gamma, lo, hi);
         }
     }
 }
@@ -496,21 +497,6 @@ export function cloth_contact_dual_update(state, table, beta, k_max_ratio, inver
     }
 }
 
-/**
- * @param {number} value
- * @param {number} lo
- * @param {number} hi
- * @returns {number}
- * @private
- */
-function clamp_band(value, lo, hi) {
-    if (value < lo) {
-        return lo;
-    }
-
-    return value > hi ? hi : value;
-}
-
 /**
  * @type {Float64Array}
  * @private

package/src/engine/physics/cloth/solver/constraint/cloth_leash_row.d.ts

@@ -11,7 +11,7 @@
  * | | `s` | `r` |
  * |---|---|---|
  * | slack | the particle's **animated position** | `slack · d_geo,i` |
- * | tether (LRA) | the **anchor's** position | `d_geo,i` |
+ * | tether (LRA) | the nearest **held particle's** position | `d_geo,i` |
  *
  * So they are one module, parameterised by mode, exactly as
  * `physics/constraint/solve_constraints.js` does for its 6-DOF joints at the

package/src/engine/physics/cloth/solver/constraint/cloth_leash_row.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"cloth_leash_row.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/physics/cloth/solver/constraint/cloth_leash_row.js"],"names":[],"mappings":"AAcA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+DG;AAEH;;;;;;;;;;;GAWG;AACH,+CANW,YAAY,GAAC,MAAM,EAAE,mBACrB,MAAM,0BAEN,MAAM,GACJ,MAAM,CA8BlB;AAED;;;;;;GAMG;AACH,oDAJW,YAAY,+BAEZ,MAAM,QA8GhB;AAED;;;;;;;;;;;;;;GAcG;AACH,iEANW,MAAM,SACN,MAAM,eACN,MAAM,eACN,MAAM,qBACN,MAAM,QA0ChB;AAED;;;;;;;;GAQG;AACH,iEAJW,MAAM,eACN,MAAM,qBACN,MAAM,QA2ChB"}
+{"version":3,"file":"cloth_leash_row.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/physics/cloth/solver/constraint/cloth_leash_row.js"],"names":[],"mappings":"AAcA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+DG;AAEH;;;;;;;;;;;GAWG;AACH,+CANW,YAAY,GAAC,MAAM,EAAE,mBACrB,MAAM,0BAEN,MAAM,GACJ,MAAM,CA8BlB;AAED;;;;;;GAMG;AACH,oDAJW,YAAY,+BAEZ,MAAM,QAqHhB;AAED;;;;;;;;;;;;;;GAcG;AACH,iEANW,MAAM,SACN,MAAM,eACN,MAAM,eACN,MAAM,qBACN,MAAM,QA0ChB;AAED;;;;;;;;GAQG;AACH,iEAJW,MAAM,eACN,MAAM,qBACN,MAAM,QA2ChB"}