reze-engine 0.50.10 → 0.50.12

package/dist/physics/autofit.js

@@ -0,0 +1,501 @@
+import { Vec3, Quat, Mat4 } from "../math";
+import { RigidbodyShape, RigidbodyType } from "./types";
+const DEFAULTS = {
+    minProtrusion: 0.05,
+    maxProtrusionRatio: 1.0,
+    minVertices: 24,
+    quantile: 0.02,
+    minGain: 0.1,
+    minDirectionality: 0.5,
+    lobeCos: 0.5,
+    minConcentration: 1.3,
+    maxShapesPerBody: 3,
+    radiusQuantile: 0.5,
+    lobeDepthFraction: 0.45,
+};
+const VERTEX_STRIDE = 8;
+/**
+ * Fit supplementary colliders for every bone-following body whose own mesh
+ * escapes it.
+ *
+ * `vertices` is the model's interleaved vertex buffer (stride 8, position
+ * first) and `skinning` its joints/weights — both in PMX bind pose, the same
+ * space as `shapePosition` / `shapeRotation`, so nothing has to be skinned.
+ *
+ * Returns one entry per body it chose to supplement; the caller appends
+ * `.body` to the rigid body list before constructing RezePhysics.
+ */
+export function fitSupplementaryColliders(rigidbodies, vertices, skinning, options = {}) {
+    const opt = { ...DEFAULTS, ...options };
+    const vertexCount = vertices.length / VERTEX_STRIDE;
+    if (vertexCount === 0)
+        return [];
+    // Bucket vertices by their dominant bone once. Every body then reads only the
+    // mesh it stands in for, instead of the pass being O(bodies × vertices).
+    const byBone = new Map();
+    for (let i = 0; i < vertexCount; i++) {
+        let bestW = -1;
+        let bestB = -1;
+        for (let k = 0; k < 4; k++) {
+            const w = skinning.weights[i * 4 + k];
+            if (w > bestW) {
+                bestW = w;
+                bestB = skinning.joints[i * 4 + k];
+            }
+        }
+        if (bestB < 0)
+            continue;
+        const list = byBone.get(bestB);
+        if (list === undefined)
+            byBone.set(bestB, [i]);
+        else
+            list.push(i);
+    }
+    const out = [];
+    for (let b = 0; b < rigidbodies.length; b++) {
+        const rb = rigidbodies[b];
+        // Only bone-following proxies. A dynamic body IS the cloth — it has no
+        // skinned surface of its own to escape from.
+        if (rb.type === RigidbodyType.Dynamic && rb.mass > 0)
+            continue;
+        if (rb.boneIndex < 0)
+            continue;
+        const owned = byBone.get(rb.boneIndex);
+        if (owned === undefined || owned.length < opt.minVertices)
+            continue;
+        // Fit repeatedly, each pass measuring against what the earlier ones cover.
+        const covered = [];
+        for (let pass = 0; pass < opt.maxShapesPerBody; pass++) {
+            const fitted = fitOne(rb, b, owned, vertices, opt, covered);
+            if (fitted === null)
+                break;
+            covered.push(fitted.local);
+            out.push(fitted.result);
+        }
+    }
+    return out;
+}
+// Scratch, reused across bodies.
+const _local = [];
+function fitOne(rb, index, owned, vertices, opt, covered) {
+    const q = Quat.fromEuler(rb.shapeRotation.x, rb.shapeRotation.y, rb.shapeRotation.z);
+    const radius = shapeRadius(rb);
+    if (radius <= 0)
+        return null;
+    const maxProtrusion = radius * opt.maxProtrusionRatio;
+    // Collect protruding vertices in the SOURCE BODY's local frame. Working in
+    // that frame rather than the bone's means the fitted box inherits the
+    // proxy's own orientation, which is already aligned with the limb — no
+    // covariance fit, and no second convention to keep straight.
+    //
+    // Each one is stored with the direction it escapes in (the gradient of the
+    // signed distance), because which way the mesh escapes is what decides
+    // whether this is a feature worth fitting at all.
+    _local.length = 0;
+    let deepest = 0;
+    let seedX = 0, seedY = 1, seedZ = 0;
+    let sumProtrusion = 0;
+    for (let n = 0; n < owned.length; n++) {
+        const o = owned[n] * VERTEX_STRIDE;
+        _v.x = vertices[o] - rb.shapePosition.x;
+        _v.y = vertices[o + 1] - rb.shapePosition.y;
+        _v.z = vertices[o + 2] - rb.shapePosition.z;
+        // Into the body's local frame: rotate by the conjugate.
+        Quat.rotateVecInvInto(q, _v, _r);
+        const px = _r.x, py = _r.y, pz = _r.z;
+        const sd = coveredSignedDistance(rb, covered, px, py, pz);
+        if (sd <= opt.minProtrusion || sd > maxProtrusion)
+            continue;
+        coveredOutwardDir(rb, covered, px, py, pz, _r);
+        // Seed the lobe on the DEEPEST vertex, not the mean direction. A pelvis
+        // escapes its capsule at the back by 0.53 and at the front by 0.14, so the
+        // mean cancels below any usable threshold and the butt — the whole point —
+        // is never fitted. The deepest point is always ON the feature.
+        if (sd > deepest) {
+            deepest = sd;
+            seedX = _r.x;
+            seedY = _r.y;
+            seedZ = _r.z;
+        }
+        sumProtrusion += sd;
+        _local.push(px, py, pz, _r.x, _r.y, _r.z, sd);
+    }
+    let found = _local.length / 7;
+    if (found < opt.minVertices)
+        return null;
+    // Averaged over the PROTRUDING vertices, deliberately — the question being
+    // asked is "is the deepest lobe deeper than protrusion typically is on this
+    // body", which is exactly the test that separates a butt from a hairstyle.
+    //
+    // Averaging over all the body's vertices instead was tried and is worse: it
+    // lifts every ratio, and measured on 托特 it scored the head's hair (2.4–4.2)
+    // ABOVE the pelvis (2.85), so no threshold could separate them and the pass
+    // fitted 33 bodies including three around the hair.
+    //
+    // The known blind spot is a mesh that escapes in exactly one place at uniform
+    // depth: the lobe IS the whole protruding set, so the ratio is 1 and nothing
+    // is fitted. Real skin always has a gradient, and for a pass that adds bodies
+    // to someone else's rig, declining when the evidence is ambiguous is the
+    // right way to be wrong.
+    const meanProtrusion = sumProtrusion / found;
+    // Settle the lobe axis: take everything within lobeCos of the seed, recentre
+    // on that set's depth-weighted mean direction, repeat. Two passes is enough —
+    // the seed already sits on the feature, this only centres it.
+    let mdx = seedX, mdy = seedY, mdz = seedZ;
+    for (let pass = 0; pass < 2; pass++) {
+        let ax = 0, ay = 0, az = 0;
+        for (let i = 0; i < found; i++) {
+            const s = i * 7;
+            if (_local[s + 3] * mdx + _local[s + 4] * mdy + _local[s + 5] * mdz < opt.lobeCos)
+                continue;
+            const wgt = _local[s + 6];
+            ax += _local[s + 3] * wgt;
+            ay += _local[s + 4] * wgt;
+            az += _local[s + 5] * wgt;
+        }
+        const len = Math.sqrt(ax * ax + ay * ay + az * az);
+        if (len < 1e-9)
+            break;
+        mdx = ax / len;
+        mdy = ay / len;
+        mdz = az / len;
+    }
+    // Keep only the lobe. Fitting the whole protruding shell would wrap the limb
+    // and put box corners out in open air on the side that was already covered.
+    // Depth floor for the lobe's core, relative to the deepest point IN the lobe
+    // (not the body's overall deepest, which may sit in a different direction).
+    let lobeDeepest = 0;
+    for (let i = 0; i < found; i++) {
+        const s = i * 7;
+        if (_local[s + 3] * mdx + _local[s + 4] * mdy + _local[s + 5] * mdz < opt.lobeCos)
+            continue;
+        if (_local[s + 6] > lobeDeepest)
+            lobeDeepest = _local[s + 6];
+    }
+    const depthFloor = lobeDeepest * opt.lobeDepthFraction;
+    let w = 0;
+    let lobeSum = 0;
+    let dirX = 0, dirY = 0, dirZ = 0;
+    for (let i = 0; i < found; i++) {
+        const s = i * 7;
+        if (_local[s + 3] * mdx + _local[s + 4] * mdy + _local[s + 5] * mdz < opt.lobeCos)
+            continue;
+        if (_local[s + 6] < depthFloor)
+            continue;
+        lobeSum += _local[s + 6];
+        dirX += _local[s + 3];
+        dirY += _local[s + 4];
+        dirZ += _local[s + 5];
+        _local[w * 3 + 0] = _local[s + 0];
+        _local[w * 3 + 1] = _local[s + 1];
+        _local[w * 3 + 2] = _local[s + 2];
+        w++;
+    }
+    found = w;
+    if (found < opt.minVertices)
+        return null;
+    // Concentration: is this lobe actually a FEATURE, or is the whole mesh
+    // uniformly outside its proxy? Hair around a skull escapes by a similar depth
+    // in every direction, so its deepest lobe is no deeper than typical and the
+    // ratio sits near 1. A butt is markedly deeper than the pelvis average. This
+    // is what stops the pass from boxing a hairstyle.
+    //
+    // Only the FIRST pass is judged on it. The question it answers — is this
+    // proxy standing in for this mesh at all — is settled once per body. Later
+    // passes are filling in what the first capsule could not reach, and by then
+    // the remaining protrusion is evenly spread by construction, so re-applying
+    // the test would reject every one of them and cap coverage at a single
+    // capsule. They still have to clear minGain and the vertex count.
+    const concentration = meanProtrusion > 1e-6 ? lobeSum / found / meanProtrusion : 0;
+    if (covered.length === 0 && concentration < opt.minConcentration)
+        return null;
+    const directionality = Math.sqrt(dirX * dirX + dirY * dirY + dirZ * dirZ) / found;
+    if (directionality < opt.minDirectionality)
+        return null;
+    // Quantile-clipped extent per axis.
+    const lo = [0, 0, 0];
+    const hi = [0, 0, 0];
+    const axis = new Array(found);
+    for (let a = 0; a < 3; a++) {
+        for (let i = 0; i < found; i++)
+            axis[i] = _local[i * 3 + a];
+        axis.sort((x, y) => x - y);
+        const k = Math.min(found - 1, Math.floor(found * opt.quantile));
+        lo[a] = axis[k];
+        hi[a] = axis[found - 1 - k];
+    }
+    const hx = (hi[0] - lo[0]) * 0.5;
+    const hy = (hi[1] - lo[1]) * 0.5;
+    const hz = (hi[2] - lo[2]) * 0.5;
+    if (hx <= 0 || hy <= 0 || hz <= 0)
+        return null;
+    const cx = (hi[0] + lo[0]) * 0.5;
+    const cy = (hi[1] + lo[1]) * 0.5;
+    const cz = (hi[2] + lo[2]) * 0.5;
+    // Fit a CAPSULE, not a box, and keep the source body's own orientation.
+    //
+    // A box fitted to a curved lobe puts its eight corners out in open air — on
+    // 托特's butt the worst corner measured 0.81 units from the nearest skin
+    // vertex, and cloth resting on that corner would float further off the body
+    // than the clipping it was added to fix. A capsule has no corners.
+    //
+    // Reusing the source rotation verbatim is what makes this cheap AND correct.
+    // PMX capsules are Y-aligned in their own frame, and a proxy's axis already
+    // runs along its limb — 下半身 is a transverse pill across the hips, which is
+    // exactly the way a butt lobe runs. Inheriting the rotation means no
+    // covariance fit and, more to the point, no quaternion→Euler conversion whose
+    // convention could disagree with the loader's.
+    // Radius from the lobe's own perpendicular spread, at a quantile — NOT from
+    // the enclosing half-extent.
+    //
+    // This is what decides whether the result looks right. A capsule inflates
+    // equally in every direction perpendicular to its axis, so sizing the radius
+    // to ENCLOSE the lobe makes a wide shallow patch into a fat capsule: it only
+    // needed to reach the butt's depth, but it also grew to the butt's width, and
+    // the surplus is surface standing off the skin everywhere else. Measured on
+    // 托特's pelvis that was 0.68 units of standoff, which reads as a skirt
+    // hovering in idle — trading a clipping artifact for a worse one.
+    //
+    // Taking a quantile of the actual distances puts the capsule's surface
+    // through the lobe rather than around it. The deepest points end up slightly
+    // outside it, which the next pass picks up.
+    const perp = new Array(found);
+    for (let i = 0; i < found; i++) {
+        const dx = _local[i * 3] - cx;
+        const dz = _local[i * 3 + 2] - cz;
+        perp[i] = Math.sqrt(dx * dx + dz * dz);
+    }
+    perp.sort((a, b) => a - b);
+    const fitRadius = perp[Math.min(found - 1, Math.floor(found * opt.radiusQuantile))];
+    if (fitRadius <= 0)
+        return null;
+    // size.y is the FULL cylinder length between cap centres, so the capsule
+    // spans h + 2r along the axis. Solve for the lobe's own extent.
+    const cyl = Math.max(0, hy * 2 - fitRadius * 2);
+    // Gain and overshoot at the capsule's six axis extremes: how far it reaches
+    // past the shape it supplements, and how far it sits from the surface it was
+    // fitted to. Overshoot is the number to watch if a result looks puffy.
+    let gain = 0;
+    let overshoot = 0;
+    const half = cyl * 0.5;
+    const probes = [
+        [cx + fitRadius, cy, cz], [cx - fitRadius, cy, cz],
+        [cx, cy + half + fitRadius, cz], [cx, cy - half - fitRadius, cz],
+        [cx, cy, cz + fitRadius], [cx, cy, cz - fitRadius],
+    ];
+    for (const [px, py, pz] of probes) {
+        const sd = coveredSignedDistance(rb, covered, px, py, pz);
+        if (sd > gain)
+            gain = sd;
+        // Only the part of the capsule that reaches OUTSIDE the original shape is
+        // new collision surface. A probe pointing inward has no skin vertex near it
+        // by construction — the lobe is a patch on the outside — and counting it
+        // would report overshoot for geometry that is buried in the character.
+        if (sd <= 0)
+            continue;
+        let nearest = Infinity;
+        for (let i = 0; i < found; i++) {
+            const dx = _local[i * 3] - px;
+            const dy = _local[i * 3 + 1] - py;
+            const dz = _local[i * 3 + 2] - pz;
+            const d2 = dx * dx + dy * dy + dz * dz;
+            if (d2 < nearest)
+                nearest = d2;
+        }
+        const d = Math.sqrt(nearest);
+        if (d > overshoot)
+            overshoot = d;
+    }
+    if (gain < opt.minGain)
+        return null;
+    // Back out to bind-pose world. Orientation is copied verbatim from the source
+    // rather than converted — the capsule is expressed in that frame, and
+    // re-deriving Euler angles would only introduce a convention to get wrong.
+    _v.x = cx;
+    _v.y = cy;
+    _v.z = cz;
+    Quat.rotateVecInto(q, _v, _r);
+    const wx = rb.shapePosition.x + _r.x;
+    const wy = rb.shapePosition.y + _r.y;
+    const wz = rb.shapePosition.z + _r.z;
+    const result = {
+        sourceIndex: index,
+        sourceName: rb.name,
+        vertexCount: owned.length,
+        protrudingCount: found,
+        maxProtrusion: deepest,
+        gain,
+        directionality,
+        concentration,
+        overshoot,
+        body: {
+            name: `${rb.name}__fit`,
+            englishName: `${rb.englishName}__fit`,
+            boneIndex: rb.boneIndex,
+            // Same group and mask as the body it supplements: it must be visible to
+            // exactly the cloth that was already meant to collide with this limb, and
+            // invisible to whatever the rigger excluded.
+            group: rb.group,
+            collisionMask: rb.collisionMask,
+            shape: RigidbodyShape.Capsule,
+            // PMX capsule semantics: size.x is the radius, size.y the FULL cylinder
+            // length between cap centres (Bullet halves it internally), size.z unused.
+            size: new Vec3(fitRadius, cyl, fitRadius),
+            shapePosition: new Vec3(wx, wy, wz),
+            shapeRotation: new Vec3(rb.shapeRotation.x, rb.shapeRotation.y, rb.shapeRotation.z),
+            // Static: it follows its bone and never integrates. Mass 0 keeps it out
+            // of the dynamic path entirely.
+            mass: 0,
+            linearDamping: 0,
+            angularDamping: 0,
+            restitution: rb.restitution,
+            friction: rb.friction,
+            type: RigidbodyType.Static,
+            bodyOffsetMatrixInverse: Mat4.identity(),
+        },
+    };
+    return { result, local: { cx, cy, cz, radius: fitRadius, cyl } };
+}
+/** Smallest characteristic dimension — the "how wide is this proxy" scale the
+ *  protrusion cap is judged against. */
+function shapeRadius(rb) {
+    switch (rb.shape) {
+        case RigidbodyShape.Sphere:
+            return rb.size.x;
+        case RigidbodyShape.Capsule:
+            return rb.size.x;
+        case RigidbodyShape.Box:
+            return Math.min(rb.size.x, rb.size.y, rb.size.z);
+        default:
+            return 0;
+    }
+}
+/** Signed distance to a Y-axis capsule centred at (cx,cy,cz), in the same
+ *  local frame. `cyl` is the full length between cap centres. */
+function capsuleSignedDistance(c, px, py, pz) {
+    const half = c.cyl * 0.5;
+    const dy = py - c.cy;
+    const t = dy > half ? half : dy < -half ? -half : dy;
+    const ax = px - c.cx;
+    const az = pz - c.cz;
+    const by = dy - t;
+    return Math.sqrt(ax * ax + by * by + az * az) - c.radius;
+}
+/** Distance to the UNION of the source shape and everything already fitted to
+ *  it. Union is a min over signed distances, so a later pass sees only what is
+ *  still exposed. */
+function coveredSignedDistance(rb, covered, px, py, pz) {
+    let best = localSignedDistance(rb, px, py, pz);
+    for (let i = 0; i < covered.length; i++) {
+        const d = capsuleSignedDistance(covered[i], px, py, pz);
+        if (d < best)
+            best = d;
+    }
+    return best;
+}
+/** Outward direction from whichever member of the union is nearest — the
+ *  gradient of coveredSignedDistance. */
+function coveredOutwardDir(rb, covered, px, py, pz, out) {
+    let best = localSignedDistance(rb, px, py, pz);
+    let nearest = -1;
+    for (let i = 0; i < covered.length; i++) {
+        const d = capsuleSignedDistance(covered[i], px, py, pz);
+        if (d < best) {
+            best = d;
+            nearest = i;
+        }
+    }
+    if (nearest < 0) {
+        localOutwardDir(rb, px, py, pz, out);
+        return;
+    }
+    // Same construction as localOutwardDir, for the winning capsule.
+    const c = covered[nearest];
+    const half = c.cyl * 0.5;
+    const dy = py - c.cy;
+    const t = dy > half ? half : dy < -half ? -half : dy;
+    const dx = px - c.cx;
+    const dz = pz - c.cz;
+    const by = dy - t;
+    const len = Math.sqrt(dx * dx + by * by + dz * dz);
+    if (len > 1e-9) {
+        out.x = dx / len;
+        out.y = by / len;
+        out.z = dz / len;
+    }
+    else {
+        out.x = 0;
+        out.y = 1;
+        out.z = 0;
+    }
+}
+/** Direction a body-local point escapes the shape in — the gradient of the
+ *  signed distance, i.e. the unit vector from its closest surface point. */
+function localOutwardDir(rb, px, py, pz, out) {
+    let dx, dy, dz;
+    switch (rb.shape) {
+        case RigidbodyShape.Capsule: {
+            const half = rb.size.y * 0.5;
+            const t = py > half ? half : py < -half ? -half : py;
+            dx = px;
+            dy = py - t;
+            dz = pz;
+            break;
+        }
+        case RigidbodyShape.Box: {
+            dx = px - Math.max(-rb.size.x, Math.min(rb.size.x, px));
+            dy = py - Math.max(-rb.size.y, Math.min(rb.size.y, py));
+            dz = pz - Math.max(-rb.size.z, Math.min(rb.size.z, pz));
+            break;
+        }
+        default: {
+            dx = px;
+            dy = py;
+            dz = pz;
+            break;
+        }
+    }
+    const len = Math.sqrt(dx * dx + dy * dy + dz * dz);
+    if (len > 1e-9) {
+        out.x = dx / len;
+        out.y = dy / len;
+        out.z = dz / len;
+    }
+    else {
+        out.x = 0;
+        out.y = 1;
+        out.z = 0;
+    }
+}
+/** Signed distance from a body-local point to the shape's surface; > 0 outside. */
+function localSignedDistance(rb, px, py, pz) {
+    switch (rb.shape) {
+        case RigidbodyShape.Sphere:
+            return Math.sqrt(px * px + py * py + pz * pz) - rb.size.x;
+        case RigidbodyShape.Capsule: {
+            // PMX capsules are Y-aligned and size.y is the FULL cylinder length
+            // between cap centres, so the segment runs ±size.y/2.
+            const half = rb.size.y * 0.5;
+            const t = py > half ? half : py < -half ? -half : py;
+            const dy = py - t;
+            return Math.sqrt(px * px + dy * dy + pz * pz) - rb.size.x;
+        }
+        case RigidbodyShape.Box: {
+            // size is half-extents for boxes.
+            const qx = Math.abs(px) - rb.size.x;
+            const qy = Math.abs(py) - rb.size.y;
+            const qz = Math.abs(pz) - rb.size.z;
+            const ox = Math.max(qx, 0), oy = Math.max(qy, 0), oz = Math.max(qz, 0);
+            const outside = Math.sqrt(ox * ox + oy * oy + oz * oz);
+            return outside + Math.min(Math.max(qx, qy, qz), 0);
+        }
+        default:
+            return -Infinity;
+    }
+}
+// Rotation goes through Quat.rotateVecInto / rotateVecInvInto rather than
+// open-coded component formulas — this pass runs once at load, so there is
+// nothing to win by hand-rolling it and a sign error here would misplace every
+// fitted body.
+const _v = new Vec3(0, 0, 0);
+const _r = new Vec3(0, 0, 0);

package/dist/shaders/anchor-table.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"anchor-table.d.ts","sourceRoot":"","sources":["../../src/shaders/anchor-table.ts"],"names":[],"mappings":"AAsBA,UAAU,aAAa;IACrB,IAAI,EAAE,MAAM,CAAA;IACZ,KAAK,EAAE,OAAO,CAAA;CACf;AAED,4EAA4E;AAC5E,eAAO,MAAM,kBAAkB,EAAE,WAAqD,CAAA;AAEtF,MAAM,WAAW,WAAW;IAC1B,+EAA+E;IAC/E,OAAO,EAAE,aAAa,EAAE,CAAA;IACxB,8EAA8E;IAC9E,KAAK,EAAE,MAAM,EAAE,EAAE,CAAA;IACjB,6EAA6E;IAC7E,OAAO,EAAE;QAAE,MAAM,EAAE,MAAM,CAAC;QAAC,IAAI,EAAE,MAAM,CAAA;KAAE,EAAE,CAAA;CAC5C;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,gBAAgB,CAAC,QAAQ,EAAE,aAAa,EAAE,EAAE,EAAE,GAAG,EAAE,MAAM,GAAG,WAAW,CAiCtF;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,cAAc,CAAC,UAAU,EAAE,MAAM,EAAE,GAAG,MAAM,CAkB3D;AAED;;;;;;;GAOG;AACH,wBAAgB,eAAe,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,MAAM,CAkBvD"}
+{"version":3,"file":"anchor-table.d.ts","sourceRoot":"","sources":["../../src/shaders/anchor-table.ts"],"names":[],"mappings":"AAsBA,UAAU,aAAa;IACrB,IAAI,EAAE,MAAM,CAAA;IACZ,KAAK,EAAE,OAAO,CAAA;CACf;AAED,4EAA4E;AAC5E,eAAO,MAAM,kBAAkB,EAAE,WAAqD,CAAA;AAEtF,MAAM,WAAW,WAAW;IAC1B,+EAA+E;IAC/E,OAAO,EAAE,aAAa,EAAE,CAAA;IACxB,8EAA8E;IAC9E,KAAK,EAAE,MAAM,EAAE,EAAE,CAAA;IACjB,6EAA6E;IAC7E,OAAO,EAAE;QAAE,MAAM,EAAE,MAAM,CAAC;QAAC,IAAI,EAAE,MAAM,CAAA;KAAE,EAAE,CAAA;CAC5C;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,gBAAgB,CAAC,QAAQ,EAAE,aAAa,EAAE,EAAE,EAAE,GAAG,EAAE,MAAM,GAAG,WAAW,CAiCtF;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,cAAc,CAAC,UAAU,EAAE,MAAM,EAAE,GAAG,MAAM,CAkB3D;AAED;;;;;;;GAOG;AACH,wBAAgB,eAAe,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,MAAM,CAmCvD"}

package/dist/shaders/anchor-table.js

@@ -108,6 +108,23 @@ ${cases}
  * folds away.
  */
 export function anchorAliasWgsl(alias) {
+    // NO @anchor AT ALL. Every slot is unmapped, and rzAnchor() must report
+    // .valid false for all of them.
+    //
+    // This used to fall through to the identity branch below — `[].every()` is
+    // true — so an effect that declared no bone and called rzAnchor(c, 0) got
+    // slot 0 of the SCENE's table: whichever bone the first effect that did
+    // declare one happened to name. Alone in a scene it read a zeroed buffer and
+    // looked correct; add a second effect and it silently anchored to that
+    // effect's wrist. Cross-effect bleed, and invisible from the file you were
+    // reading.
+    if (alias.length === 0) {
+        return `
+/** No @anchor in this effect: every slot is unmapped, so rzAnchor() is invalid
+ *  rather than reading whichever bone another effect declared first. */
+fn _rzSlot(local: i32) -> i32 { return -1; }
+`;
+    }
     const identity = alias.every((g, i) => g === i);
     if (identity) {
         return `

package/dist/shaders/materials/common.d.ts

@@ -19,6 +19,6 @@ export declare function commonFsOutWgsl(): string;
  * clumps of a few millimetres rather than as single-pixel snow, which at any
  * distance is just a fade.
  */
-export declare const DISSOLVE_WGSL = "\n/** Clump size, in object-space units \u2014 MMD's are roughly centimetres. */\nconst RZ_DISSOLVE_GRAIN: f32 = 0.42;\n/** Grit cell, as a fraction of a clump. See rz_dissolve_threshold: this is what\n *  keeps a material SMALLER than a clump from having one threshold for all of\n *  it. Seven millimetres or so on a normal model. */\nconst RZ_DISSOLVE_GRIT: f32 = 0.16;\n/** How wide the glowing front is, in threshold units. Thin: the front is a\n *  rim, and a wide one lights whole limbs at once. */\nconst RZ_DISSOLVE_EDGE: f32 = 0.11;\n/** What the front burns with. Emission, added after lighting and after the\n *  graph \u2014 a colour that ramps or takes a shadow reads as paint, not as heat.\n *  Above 1 on purpose: it is meant to reach the bloom pyramid. Blue-violet,\n *  the same light the motes leaving her are drawn and lit with. */\nconst RZ_BURN_COLOR: vec3f = vec3f(0.45, 0.62, 2.40);\n\nfn rz_dissolve_hash(p: vec3f) -> f32 {\n  let q = fract(p * 0.3183099 + vec3f(0.1, 0.2, 0.3));\n  let r = q * 17.0 * (q + 34.0);\n  return fract(r.x * r.y * r.z);\n}\n\n/** Value noise over that hash \u2014 smooth within a flake, uncorrelated between. */\nfn rz_dissolve_field(p: vec3f) -> f32 {\n  let i = floor(p);\n  let f = fract(p);\n  let u = f * f * (3.0 - 2.0 * f);\n  let c000 = rz_dissolve_hash(i);\n  let c100 = rz_dissolve_hash(i + vec3f(1.0, 0.0, 0.0));\n  let c010 = rz_dissolve_hash(i + vec3f(0.0, 1.0, 0.0));\n  let c110 = rz_dissolve_hash(i + vec3f(1.0, 1.0, 0.0));\n  let c001 = rz_dissolve_hash(i + vec3f(0.0, 0.0, 1.0));\n  let c101 = rz_dissolve_hash(i + vec3f(1.0, 0.0, 1.0));\n  let c011 = rz_dissolve_hash(i + vec3f(0.0, 1.0, 1.0));\n  let c111 = rz_dissolve_hash(i + vec3f(1.0, 1.0, 1.0));\n  let x00 = mix(c000, c100, u.x);\n  let x10 = mix(c010, c110, u.x);\n  let x01 = mix(c001, c101, u.x);\n  let x11 = mix(c011, c111, u.x);\n  return mix(mix(x00, x10, u.y), mix(x01, x11, u.y), u.z);\n}\n\n/**\n * The threshold this fragment is measured against.\n *\n * NOISE ALONE, so the whole body comes apart at once \u2014 a hand, a shoulder and\n * an ankle all thinning together rather than a line travelling up her. It was\n * tilted by height at first, on the argument that a sweep reads as something\n * happening TO her; what it actually reads as is a wipe, and a wipe has a\n * direction, an edge and a speed the rest of the effect knows nothing about.\n * Everywhere at once is also the honest match for what leaves: the motes are\n * shed from every limb in the same breath.\n *\n * TWO SCALES, and the second one is not decoration \u2014 it is the fix for a real\n * failure. Clumps alone meant every fragment of a material SMALLER than one\n * cell shared a single threshold: an eye highlight, a brow, a button, a hair\n * clip. All of it sat inside the glowing front at once, so the piece flashed\n * white as a whole and then popped out of existence \u2014 the two things a\n * disintegration must never do. The grit is a plain per-cell hash at a fraction\n * of the clump, so a part a few millimetres across still comes apart in pieces.\n *\n * Weighted toward the clumps, because they are the LOOK: flakes a centimetre\n * across, in an order fixed to the surface, so the same piece of sleeve goes at\n * the same moment every time she leaves. The grit only breaks ties.\n */\nfn rz_dissolve_threshold(restPos: vec3f) -> f32 {\n  let clump = rz_dissolve_field(restPos / RZ_DISSOLVE_GRAIN);\n  let grit = rz_dissolve_hash(floor(restPos / (RZ_DISSOLVE_GRAIN * RZ_DISSOLVE_GRIT)));\n  return clamp(clump * 0.7 + grit * 0.3, 0.0, 1.0);\n}\n";
+export declare const DISSOLVE_WGSL = "\n/** Clump size, in object-space units \u2014 MMD's are roughly centimetres. */\nconst RZ_DISSOLVE_GRAIN: f32 = 0.42;\n/** How much the clumps pull a fragment off its own grit value \u2014 the whole of\n *  the LOOK, and none of the rate. See rz_dissolve_threshold. */\nconst RZ_DISSOLVE_CLUMP: f32 = 0.6;\n/** Grit cell, as a fraction of a clump. See rz_dissolve_threshold: this is what\n *  keeps a material SMALLER than a clump from having one threshold for all of\n *  it, and it has to be smaller than the SMALLEST material a model carries \u2014\n *  an eye highlight, a button, a buckle \u2014 not merely small. Two and a half\n *  millimetres on a normal model, which is under all of them. */\nconst RZ_DISSOLVE_GRIT: f32 = 0.06;\n/** How wide the glowing front is, in threshold units. Thin: the front is a\n *  rim, and a wide one lights whole limbs at once. */\nconst RZ_DISSOLVE_EDGE: f32 = 0.11;\n/** What the front burns with. Emission, added after lighting and after the\n *  graph \u2014 a colour that ramps or takes a shadow reads as paint, not as heat.\n *  Above 1 on purpose: it is meant to reach the bloom pyramid. Blue-violet,\n *  the same light the motes leaving her are drawn and lit with. */\nconst RZ_BURN_COLOR: vec3f = vec3f(0.45, 0.62, 2.40);\n\nfn rz_dissolve_hash(p: vec3f) -> f32 {\n  // INTEGER, for exactly the reason _hash33 in nodes.ts is, and its comment\n  // says it: above 24 bits an f32 loses the low ones. The float hash this\n  // replaces multiplied three terms into the hundreds of millions and took\n  // fract() of the product \u2014 which past 2^24 is not noise, it is a quantised\n  // staircase, and wherever that staircase landed on zero the threshold was\n  // zero. A surface whose threshold is zero never crosses it: it stays behind\n  // while the rest of her leaves, and it sits inside the burning front the\n  // whole time she is gone. Those were the white pieces that would not go.\n  //\n  // Floored first, unlike _hash33: vec3i() truncates toward zero, so -0.5 and\n  // +0.5 are the same cell, and half this model is at negative x.\n  var h = vec3u(vec3i(floor(p)) + vec3i(32768));\n  h = h * vec3u(1664525u, 1013904223u, 2654435761u);\n  h = (h.yzx ^ h) * vec3u(2246822519u, 3266489917u, 668265263u);\n  h = h ^ (h >> vec3u(16u));\n  return f32((h.x ^ h.y ^ h.z) & 16777215u) * (1.0 / 16777216.0);\n}\n\n/** Value noise over that hash \u2014 smooth within a flake, uncorrelated between. */\nfn rz_dissolve_field(p: vec3f) -> f32 {\n  let i = floor(p);\n  let f = fract(p);\n  let u = f * f * (3.0 - 2.0 * f);\n  let c000 = rz_dissolve_hash(i);\n  let c100 = rz_dissolve_hash(i + vec3f(1.0, 0.0, 0.0));\n  let c010 = rz_dissolve_hash(i + vec3f(0.0, 1.0, 0.0));\n  let c110 = rz_dissolve_hash(i + vec3f(1.0, 1.0, 0.0));\n  let c001 = rz_dissolve_hash(i + vec3f(0.0, 0.0, 1.0));\n  let c101 = rz_dissolve_hash(i + vec3f(1.0, 0.0, 1.0));\n  let c011 = rz_dissolve_hash(i + vec3f(0.0, 1.0, 1.0));\n  let c111 = rz_dissolve_hash(i + vec3f(1.0, 1.0, 1.0));\n  let x00 = mix(c000, c100, u.x);\n  let x10 = mix(c010, c110, u.x);\n  let x01 = mix(c001, c101, u.x);\n  let x11 = mix(c011, c111, u.x);\n  return mix(mix(x00, x10, u.y), mix(x01, x11, u.y), u.z);\n}\n\n/**\n * The threshold this fragment is measured against.\n *\n * NOISE ALONE, so the whole body comes apart at once \u2014 a hand, a shoulder and\n * an ankle all thinning together rather than a line travelling up her. It was\n * tilted by height at first, on the argument that a sweep reads as something\n * happening TO her; what it actually reads as is a wipe, and a wipe has a\n * direction, an edge and a speed the rest of the effect knows nothing about.\n * Everywhere at once is also the honest match for what leaves: the motes are\n * shed from every limb in the same breath.\n *\n * TWO SCALES. Clumps alone meant every fragment of a material SMALLER than one\n * cell shared a single threshold: an eye highlight, a brow, a button, a hair\n * clip. All of it sat inside the glowing front at once, so the piece flashed\n * white as a whole and then popped \u2014 the two things a disintegration must never\n * do. The grit is a per-cell hash at a fraction of the clump, so a part a few\n * millimetres across still comes apart in pieces.\n *\n * THE GRIT IS THE BASE, though, and the clumps only push a fragment off it.\n * That ordering is what makes the dissolve smooth, and it is not obvious: the\n * two used to be WEIGHTED AND ADDED, and a sum is not uniformly distributed\n * however uniform each half is. It piles up in the middle, so most of the\n * surface had a threshold near a half \u2014 all of it went as the dissolve swept\n * through the middle, and what remained was the thin tail near zero. That tail\n * was the second stage everyone could see: a body that mostly vanished, then a\n * scatter of glowing scraps, then nothing.\n *\n * WRAPPED, not added. Offsetting a uniform value and taking fract() leaves it\n * EXACTLY uniform, where scaling a sum back into range only flattens the middle\n * \u2014 measured over four hundred thousand samples, the sum removed 12% of her per\n * step at the halfway mark and 0.1% at each end (a 110x swing), the scaled sum\n * still swung 89x, and the wrap is flat: 5.0% per step from the first flake to\n * the last, a ratio of 1.0.\n *\n * That flatness IS the smoothness. An uneven rate is what produced the two\n * stages you could see \u2014 a body that mostly vanished at once, then a scatter of\n * scraps that took as long again to follow.\n *\n * The clumping survives the wrap: neighbouring fragments share a clump value,\n * so the offset moves them together. What wraps past 1 comes back at 0 and goes\n * at the other end instead, which costs a little speckle in a region and buys\n * an even rate everywhere.\n */\nfn rz_dissolve_threshold(restPos: vec3f) -> f32 {\n  let grit = rz_dissolve_hash(floor(restPos / (RZ_DISSOLVE_GRAIN * RZ_DISSOLVE_GRIT)));\n  let clump = rz_dissolve_field(restPos / RZ_DISSOLVE_GRAIN) - 0.5;\n  return fract(grit + clump * RZ_DISSOLVE_CLUMP);\n}\n";
 export declare const COMMON_MATERIAL_PRELUDE_WGSL: string;
 //# sourceMappingURL=common.d.ts.map

package/dist/shaders/materials/common.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"common.d.ts","sourceRoot":"","sources":["../../../src/shaders/materials/common.ts"],"names":[],"mappings":"AA2BA,eAAO,MAAM,oBAAoB,QAuHhC,CAAC;AAgIF,wBAAgB,eAAe,IAAI,MAAM,CAExC;AAGD;;;;;;;;;;;;;;;;;;GAkBG;AACH,eAAO,MAAM,aAAa,mgHAsEzB,CAAA;AAQD,eAAO,MAAM,4BAA4B,QAC0D,CAAA"}
+{"version":3,"file":"common.d.ts","sourceRoot":"","sources":["../../../src/shaders/materials/common.ts"],"names":[],"mappings":"AA2BA,eAAO,MAAM,oBAAoB,QAuHhC,CAAC;AAgIF,wBAAgB,eAAe,IAAI,MAAM,CAExC;AAGD;;;;;;;;;;;;;;;;;;GAkBG;AACH,eAAO,MAAM,aAAa,k8LA4GzB,CAAA;AAQD,eAAO,MAAM,4BAA4B,QAC0D,CAAA"}

package/dist/shaders/materials/common.js

@@ -288,10 +288,15 @@ export function commonFsOutWgsl() {
 export const DISSOLVE_WGSL = /* wgsl */ `
 /** Clump size, in object-space units — MMD's are roughly centimetres. */
 const RZ_DISSOLVE_GRAIN: f32 = 0.42;
+/** How much the clumps pull a fragment off its own grit value — the whole of
+ *  the LOOK, and none of the rate. See rz_dissolve_threshold. */
+const RZ_DISSOLVE_CLUMP: f32 = 0.6;
 /** Grit cell, as a fraction of a clump. See rz_dissolve_threshold: this is what
  *  keeps a material SMALLER than a clump from having one threshold for all of
- *  it. Seven millimetres or so on a normal model. */
-const RZ_DISSOLVE_GRIT: f32 = 0.16;
+ *  it, and it has to be smaller than the SMALLEST material a model carries —
+ *  an eye highlight, a button, a buckle — not merely small. Two and a half
+ *  millimetres on a normal model, which is under all of them. */
+const RZ_DISSOLVE_GRIT: f32 = 0.06;
 /** How wide the glowing front is, in threshold units. Thin: the front is a
  *  rim, and a wide one lights whole limbs at once. */
 const RZ_DISSOLVE_EDGE: f32 = 0.11;
@@ -302,9 +307,22 @@ const RZ_DISSOLVE_EDGE: f32 = 0.11;
 const RZ_BURN_COLOR: vec3f = vec3f(0.45, 0.62, 2.40);
 
 fn rz_dissolve_hash(p: vec3f) -> f32 {
-  let q = fract(p * 0.3183099 + vec3f(0.1, 0.2, 0.3));
-  let r = q * 17.0 * (q + 34.0);
-  return fract(r.x * r.y * r.z);
+  // INTEGER, for exactly the reason _hash33 in nodes.ts is, and its comment
+  // says it: above 24 bits an f32 loses the low ones. The float hash this
+  // replaces multiplied three terms into the hundreds of millions and took
+  // fract() of the product — which past 2^24 is not noise, it is a quantised
+  // staircase, and wherever that staircase landed on zero the threshold was
+  // zero. A surface whose threshold is zero never crosses it: it stays behind
+  // while the rest of her leaves, and it sits inside the burning front the
+  // whole time she is gone. Those were the white pieces that would not go.
+  //
+  // Floored first, unlike _hash33: vec3i() truncates toward zero, so -0.5 and
+  // +0.5 are the same cell, and half this model is at negative x.
+  var h = vec3u(vec3i(floor(p)) + vec3i(32768));
+  h = h * vec3u(1664525u, 1013904223u, 2654435761u);
+  h = (h.yzx ^ h) * vec3u(2246822519u, 3266489917u, 668265263u);
+  h = h ^ (h >> vec3u(16u));
+  return f32((h.x ^ h.y ^ h.z) & 16777215u) * (1.0 / 16777216.0);
 }
 
 /** Value noise over that hash — smooth within a flake, uncorrelated between. */
@@ -338,22 +356,42 @@ fn rz_dissolve_field(p: vec3f) -> f32 {
  * Everywhere at once is also the honest match for what leaves: the motes are
  * shed from every limb in the same breath.
  *
- * TWO SCALES, and the second one is not decoration — it is the fix for a real
- * failure. Clumps alone meant every fragment of a material SMALLER than one
+ * TWO SCALES. Clumps alone meant every fragment of a material SMALLER than one
  * cell shared a single threshold: an eye highlight, a brow, a button, a hair
  * clip. All of it sat inside the glowing front at once, so the piece flashed
- * white as a whole and then popped out of existence — the two things a
- * disintegration must never do. The grit is a plain per-cell hash at a fraction
- * of the clump, so a part a few millimetres across still comes apart in pieces.
+ * white as a whole and then popped — the two things a disintegration must never
+ * do. The grit is a per-cell hash at a fraction of the clump, so a part a few
+ * millimetres across still comes apart in pieces.
+ *
+ * THE GRIT IS THE BASE, though, and the clumps only push a fragment off it.
+ * That ordering is what makes the dissolve smooth, and it is not obvious: the
+ * two used to be WEIGHTED AND ADDED, and a sum is not uniformly distributed
+ * however uniform each half is. It piles up in the middle, so most of the
+ * surface had a threshold near a half — all of it went as the dissolve swept
+ * through the middle, and what remained was the thin tail near zero. That tail
+ * was the second stage everyone could see: a body that mostly vanished, then a
+ * scatter of glowing scraps, then nothing.
+ *
+ * WRAPPED, not added. Offsetting a uniform value and taking fract() leaves it
+ * EXACTLY uniform, where scaling a sum back into range only flattens the middle
+ * — measured over four hundred thousand samples, the sum removed 12% of her per
+ * step at the halfway mark and 0.1% at each end (a 110x swing), the scaled sum
+ * still swung 89x, and the wrap is flat: 5.0% per step from the first flake to
+ * the last, a ratio of 1.0.
+ *
+ * That flatness IS the smoothness. An uneven rate is what produced the two
+ * stages you could see — a body that mostly vanished at once, then a scatter of
+ * scraps that took as long again to follow.
  *
- * Weighted toward the clumps, because they are the LOOK: flakes a centimetre
- * across, in an order fixed to the surface, so the same piece of sleeve goes at
- * the same moment every time she leaves. The grit only breaks ties.
+ * The clumping survives the wrap: neighbouring fragments share a clump value,
+ * so the offset moves them together. What wraps past 1 comes back at 0 and goes
+ * at the other end instead, which costs a little speckle in a region and buys
+ * an even rate everywhere.
  */
 fn rz_dissolve_threshold(restPos: vec3f) -> f32 {
-  let clump = rz_dissolve_field(restPos / RZ_DISSOLVE_GRAIN);
   let grit = rz_dissolve_hash(floor(restPos / (RZ_DISSOLVE_GRAIN * RZ_DISSOLVE_GRIT)));
-  return clamp(clump * 0.7 + grit * 0.3, 0.0, 1.0);
+  let clump = rz_dissolve_field(restPos / RZ_DISSOLVE_GRAIN) - 0.5;
+  return fract(grit + clump * RZ_DISSOLVE_CLUMP);
 }
 `;
 // ─── Convenience: full shared prelude ───────────────────────────────