reze-engine 0.14.0 → 0.15.0

package/dist/physics.js

@@ -1,527 +0,0 @@
-import { Quat, Vec3, Mat4 } from "./math";
-// Physics-local scratch pool for per-frame sync (syncFromBones, applyAmmoRigidbodiesToBones).
-// Each method uses only these slots and completes synchronously before the next is called.
-const _physMat = [
-    new Float32Array(16), new Float32Array(16), new Float32Array(16),
-];
-const _physQuat = new Quat(0, 0, 0, 1);
-import { loadAmmo } from "./ammo-loader";
-export var RigidbodyShape;
-(function (RigidbodyShape) {
-    RigidbodyShape[RigidbodyShape["Sphere"] = 0] = "Sphere";
-    RigidbodyShape[RigidbodyShape["Box"] = 1] = "Box";
-    RigidbodyShape[RigidbodyShape["Capsule"] = 2] = "Capsule";
-})(RigidbodyShape || (RigidbodyShape = {}));
-export var RigidbodyType;
-(function (RigidbodyType) {
-    RigidbodyType[RigidbodyType["Static"] = 0] = "Static";
-    RigidbodyType[RigidbodyType["Dynamic"] = 1] = "Dynamic";
-    RigidbodyType[RigidbodyType["Kinematic"] = 2] = "Kinematic";
-})(RigidbodyType || (RigidbodyType = {}));
-export class Physics {
-    constructor(rigidbodies, joints = [], options) {
-        this.gravity = new Vec3(0, -98, 0); // Gravity acceleration (cm/s²), MMD-style default
-        this.constraintSolverPattern = null;
-        this.ammoInitialized = false;
-        this.ammoPromise = null;
-        this.ammo = null;
-        // eslint-disable-next-line @typescript-eslint/no-explicit-any
-        this.dynamicsWorld = null; // btDiscreteDynamicsWorld
-        // eslint-disable-next-line @typescript-eslint/no-explicit-any
-        this.ammoRigidbodies = []; // btRigidBody instances
-        // eslint-disable-next-line @typescript-eslint/no-explicit-any
-        this.ammoConstraints = []; // btTypedConstraint instances
-        this.rigidbodiesInitialized = false; // bodyOffsetMatrixInverse computed and bodies positioned
-        this.jointsCreated = false; // Joints delayed until after rigidbodies are positioned
-        this.firstFrame = true; // Needed to reposition bodies before creating joints
-        // eslint-disable-next-line @typescript-eslint/no-explicit-any
-        this.zeroVector = null; // Cached zero vector for velocity clearing
-        this.rigidbodies = rigidbodies;
-        this.joints = joints;
-        const keywords = options?.constraintSolverKeywords ?? [];
-        if (keywords.length > 0) {
-            this.constraintSolverPattern = new RegExp(keywords.map(k => k.replace(/[.*+?^${}()|[\]\\]/g, "\\$&")).join("|"), "i");
-        }
-        this.initAmmo();
-    }
-    async initAmmo() {
-        if (this.ammoInitialized || this.ammoPromise)
-            return;
-        this.ammoPromise = loadAmmo();
-        try {
-            this.ammo = await this.ammoPromise;
-            this.createAmmoWorld();
-            this.ammoInitialized = true;
-        }
-        catch (error) {
-            console.error("[Physics] Failed to initialize Ammo:", error);
-            this.ammoPromise = null;
-        }
-    }
-    setGravity(gravity) {
-        this.gravity = gravity;
-        if (this.dynamicsWorld && this.ammo) {
-            const Ammo = this.ammo;
-            const gravityVec = new Ammo.btVector3(gravity.x, gravity.y, gravity.z);
-            this.dynamicsWorld.setGravity(gravityVec);
-            Ammo.destroy(gravityVec);
-        }
-    }
-    getGravity() {
-        return this.gravity;
-    }
-    getRigidbodies() {
-        return this.rigidbodies;
-    }
-    getJoints() {
-        return this.joints;
-    }
-    getRigidbodyTransforms() {
-        const transforms = [];
-        if (!this.ammo || !this.ammoRigidbodies.length) {
-            for (let i = 0; i < this.rigidbodies.length; i++) {
-                transforms.push({
-                    position: new Vec3(this.rigidbodies[i].shapePosition.x, this.rigidbodies[i].shapePosition.y, this.rigidbodies[i].shapePosition.z),
-                    rotation: Quat.fromEuler(this.rigidbodies[i].shapeRotation.x, this.rigidbodies[i].shapeRotation.y, this.rigidbodies[i].shapeRotation.z),
-                });
-            }
-            return transforms;
-        }
-        for (let i = 0; i < this.ammoRigidbodies.length; i++) {
-            const ammoBody = this.ammoRigidbodies[i];
-            if (!ammoBody) {
-                const rb = this.rigidbodies[i];
-                transforms.push({
-                    position: new Vec3(rb.shapePosition.x, rb.shapePosition.y, rb.shapePosition.z),
-                    rotation: Quat.fromEuler(rb.shapeRotation.x, rb.shapeRotation.y, rb.shapeRotation.z),
-                });
-                continue;
-            }
-            const transform = ammoBody.getWorldTransform();
-            const origin = transform.getOrigin();
-            const rotQuat = transform.getRotation();
-            transforms.push({
-                position: new Vec3(origin.x(), origin.y(), origin.z()),
-                rotation: new Quat(rotQuat.x(), rotQuat.y(), rotQuat.z(), rotQuat.w()),
-            });
-        }
-        return transforms;
-    }
-    createAmmoWorld() {
-        if (!this.ammo)
-            return;
-        const Ammo = this.ammo;
-        const collisionConfiguration = new Ammo.btDefaultCollisionConfiguration();
-        const dispatcher = new Ammo.btCollisionDispatcher(collisionConfiguration);
-        const overlappingPairCache = new Ammo.btDbvtBroadphase();
-        const solver = new Ammo.btSequentialImpulseConstraintSolver();
-        this.dynamicsWorld = new Ammo.btDiscreteDynamicsWorld(dispatcher, overlappingPairCache, solver, collisionConfiguration);
-        const gravityVec = new Ammo.btVector3(this.gravity.x, this.gravity.y, this.gravity.z);
-        this.dynamicsWorld.setGravity(gravityVec);
-        Ammo.destroy(gravityVec);
-        this.createAmmoRigidbodies();
-    }
-    createAmmoRigidbodies() {
-        if (!this.ammo || !this.dynamicsWorld)
-            return;
-        const Ammo = this.ammo;
-        this.ammoRigidbodies = [];
-        for (let i = 0; i < this.rigidbodies.length; i++) {
-            const rb = this.rigidbodies[i];
-            // eslint-disable-next-line @typescript-eslint/no-explicit-any
-            let shape = null;
-            const size = rb.size;
-            switch (rb.shape) {
-                case RigidbodyShape.Sphere:
-                    const radius = size.x;
-                    shape = new Ammo.btSphereShape(radius);
-                    break;
-                case RigidbodyShape.Box:
-                    const sizeVector = new Ammo.btVector3(size.x, size.y, size.z);
-                    shape = new Ammo.btBoxShape(sizeVector);
-                    Ammo.destroy(sizeVector);
-                    break;
-                case RigidbodyShape.Capsule:
-                    const capsuleRadius = size.x;
-                    const capsuleHalfHeight = size.y;
-                    shape = new Ammo.btCapsuleShape(capsuleRadius, capsuleHalfHeight);
-                    break;
-                default:
-                    const defaultHalfExtents = new Ammo.btVector3(size.x / 2, size.y / 2, size.z / 2);
-                    shape = new Ammo.btBoxShape(defaultHalfExtents);
-                    Ammo.destroy(defaultHalfExtents);
-                    break;
-            }
-            // Bodies must start at correct position to avoid explosions when joints are created
-            const transform = new Ammo.btTransform();
-            transform.setIdentity();
-            const shapePos = new Ammo.btVector3(rb.shapePosition.x, rb.shapePosition.y, rb.shapePosition.z);
-            transform.setOrigin(shapePos);
-            Ammo.destroy(shapePos);
-            const shapeRotQuat = Quat.fromEuler(rb.shapeRotation.x, rb.shapeRotation.y, rb.shapeRotation.z);
-            const quat = new Ammo.btQuaternion(shapeRotQuat.x, shapeRotQuat.y, shapeRotQuat.z, shapeRotQuat.w);
-            transform.setRotation(quat);
-            Ammo.destroy(quat);
-            // All types use the same motionState constructor
-            const motionState = new Ammo.btDefaultMotionState(transform);
-            const mass = rb.type === RigidbodyType.Dynamic ? rb.mass : 0;
-            const isDynamic = rb.type === RigidbodyType.Dynamic;
-            const localInertia = new Ammo.btVector3(0, 0, 0);
-            if (isDynamic && mass > 0) {
-                shape.calculateLocalInertia(mass, localInertia);
-            }
-            const rbInfo = new Ammo.btRigidBodyConstructionInfo(mass, motionState, shape, localInertia);
-            rbInfo.set_m_restitution(rb.restitution);
-            rbInfo.set_m_friction(rb.friction);
-            rbInfo.set_m_linearDamping(rb.linearDamping);
-            rbInfo.set_m_angularDamping(rb.angularDamping);
-            const body = new Ammo.btRigidBody(rbInfo);
-            body.setSleepingThresholds(0.0, 0.0);
-            // Static (FollowBone) should be kinematic, not static - must follow bones
-            if (rb.type === RigidbodyType.Static || rb.type === RigidbodyType.Kinematic) {
-                body.setCollisionFlags(body.getCollisionFlags() | 2); // CF_KINEMATIC_OBJECT
-                body.setActivationState(4); // DISABLE_DEACTIVATION
-            }
-            const collisionGroup = 1 << rb.group;
-            const collisionMask = rb.collisionMask;
-            const isZeroVolume = (rb.shape === RigidbodyShape.Sphere && rb.size.x === 0) ||
-                (rb.shape === RigidbodyShape.Box && (rb.size.x === 0 || rb.size.y === 0 || rb.size.z === 0)) ||
-                (rb.shape === RigidbodyShape.Capsule && (rb.size.x === 0 || rb.size.y === 0));
-            if (collisionMask === 0 || isZeroVolume) {
-                body.setCollisionFlags(body.getCollisionFlags() | 4); // CF_NO_CONTACT_RESPONSE
-            }
-            this.dynamicsWorld.addRigidBody(body, collisionGroup, collisionMask);
-            this.ammoRigidbodies.push(body);
-            Ammo.destroy(rbInfo);
-            Ammo.destroy(localInertia);
-        }
-    }
-    createAmmoJoints() {
-        if (!this.ammo || !this.dynamicsWorld || this.ammoRigidbodies.length === 0)
-            return;
-        const Ammo = this.ammo;
-        this.ammoConstraints = [];
-        for (const joint of this.joints) {
-            const rbIndexA = joint.rigidbodyIndexA;
-            const rbIndexB = joint.rigidbodyIndexB;
-            if (rbIndexA < 0 ||
-                rbIndexA >= this.ammoRigidbodies.length ||
-                rbIndexB < 0 ||
-                rbIndexB >= this.ammoRigidbodies.length) {
-                console.warn(`[Physics] Invalid joint indices: ${rbIndexA}, ${rbIndexB}`);
-                continue;
-            }
-            const bodyA = this.ammoRigidbodies[rbIndexA];
-            const bodyB = this.ammoRigidbodies[rbIndexB];
-            if (!bodyA || !bodyB) {
-                console.warn(`[Physics] Body not found for joint ${joint.name}: bodyA=${!!bodyA}, bodyB=${!!bodyB}`);
-                continue;
-            }
-            // Compute joint frames using actual current body positions (after repositioning)
-            const bodyATransform = bodyA.getWorldTransform();
-            const bodyBTransform = bodyB.getWorldTransform();
-            const bodyAOrigin = bodyATransform.getOrigin();
-            const bodyARotQuat = bodyATransform.getRotation();
-            const bodyAPos = new Vec3(bodyAOrigin.x(), bodyAOrigin.y(), bodyAOrigin.z());
-            const bodyARot = new Quat(bodyARotQuat.x(), bodyARotQuat.y(), bodyARotQuat.z(), bodyARotQuat.w());
-            const bodyAMat = Mat4.fromPositionRotation(bodyAPos, bodyARot);
-            const bodyBOrigin = bodyBTransform.getOrigin();
-            const bodyBRotQuat = bodyBTransform.getRotation();
-            const bodyBPos = new Vec3(bodyBOrigin.x(), bodyBOrigin.y(), bodyBOrigin.z());
-            const bodyBRot = new Quat(bodyBRotQuat.x(), bodyBRotQuat.y(), bodyBRotQuat.z(), bodyBRotQuat.w());
-            const bodyBMat = Mat4.fromPositionRotation(bodyBPos, bodyBRot);
-            const scalingFactor = 1.0;
-            const jointRotQuat = Quat.fromEuler(joint.rotation.x, joint.rotation.y, joint.rotation.z);
-            const jointPos = new Vec3(joint.position.x * scalingFactor, joint.position.y * scalingFactor, joint.position.z * scalingFactor);
-            const jointTransform = Mat4.fromPositionRotation(jointPos, jointRotQuat);
-            // Transform joint world position to body A's local space
-            const frameInAMat = bodyAMat.inverse().multiply(jointTransform);
-            const framePosA = frameInAMat.getPosition();
-            const frameRotA = frameInAMat.toQuat();
-            // Transform joint world position to body B's local space
-            const frameInBMat = bodyBMat.inverse().multiply(jointTransform);
-            const framePosB = frameInBMat.getPosition();
-            const frameRotB = frameInBMat.toQuat();
-            const frameInA = new Ammo.btTransform();
-            frameInA.setIdentity();
-            const pivotInA = new Ammo.btVector3(framePosA.x, framePosA.y, framePosA.z);
-            frameInA.setOrigin(pivotInA);
-            const quatA = new Ammo.btQuaternion(frameRotA.x, frameRotA.y, frameRotA.z, frameRotA.w);
-            frameInA.setRotation(quatA);
-            const frameInB = new Ammo.btTransform();
-            frameInB.setIdentity();
-            const pivotInB = new Ammo.btVector3(framePosB.x, framePosB.y, framePosB.z);
-            frameInB.setOrigin(pivotInB);
-            const quatB = new Ammo.btQuaternion(frameRotB.x, frameRotB.y, frameRotB.z, frameRotB.w);
-            frameInB.setRotation(quatB);
-            const useLinearReferenceFrameA = true;
-            const constraint = new Ammo.btGeneric6DofSpringConstraint(bodyA, bodyB, frameInA, frameInB, useLinearReferenceFrameA);
-            // Per-joint Bullet 2.75 constraint solver: disable m_useOffsetForConstraintFrame for
-            // joints whose name matches constraintSolverKeywords.
-            if (this.constraintSolverPattern && this.constraintSolverPattern.test(joint.name)) {
-                let jointPtr;
-                if (typeof Ammo.getPointer === "function") {
-                    jointPtr = Ammo.getPointer(constraint);
-                }
-                else {
-                    const constraintWithPtr = constraint;
-                    jointPtr = constraintWithPtr.ptr;
-                }
-                if (jointPtr !== undefined && Ammo.HEAP8) {
-                    const heap8 = Ammo.HEAP8;
-                    if (heap8[jointPtr + 1300] === (useLinearReferenceFrameA ? 1 : 0) && heap8[jointPtr + 1301] === 1) {
-                        heap8[jointPtr + 1301] = 0;
-                    }
-                }
-            }
-            for (let i = 0; i < 6; ++i) {
-                constraint.setParam(2, 0.475, i); // BT_CONSTRAINT_STOP_ERP
-            }
-            const lowerLinear = new Ammo.btVector3(joint.positionMin.x, joint.positionMin.y, joint.positionMin.z);
-            const upperLinear = new Ammo.btVector3(joint.positionMax.x, joint.positionMax.y, joint.positionMax.z);
-            constraint.setLinearLowerLimit(lowerLinear);
-            constraint.setLinearUpperLimit(upperLinear);
-            const lowerAngular = new Ammo.btVector3(this.normalizeAngle(joint.rotationMin.x), this.normalizeAngle(joint.rotationMin.y), this.normalizeAngle(joint.rotationMin.z));
-            const upperAngular = new Ammo.btVector3(this.normalizeAngle(joint.rotationMax.x), this.normalizeAngle(joint.rotationMax.y), this.normalizeAngle(joint.rotationMax.z));
-            constraint.setAngularLowerLimit(lowerAngular);
-            constraint.setAngularUpperLimit(upperAngular);
-            // Linear springs: only enable if stiffness is non-zero
-            if (joint.springPosition.x !== 0) {
-                constraint.setStiffness(0, joint.springPosition.x);
-                constraint.enableSpring(0, true);
-            }
-            else {
-                constraint.enableSpring(0, false);
-            }
-            if (joint.springPosition.y !== 0) {
-                constraint.setStiffness(1, joint.springPosition.y);
-                constraint.enableSpring(1, true);
-            }
-            else {
-                constraint.enableSpring(1, false);
-            }
-            if (joint.springPosition.z !== 0) {
-                constraint.setStiffness(2, joint.springPosition.z);
-                constraint.enableSpring(2, true);
-            }
-            else {
-                constraint.enableSpring(2, false);
-            }
-            // Angular springs: always enable
-            constraint.setStiffness(3, joint.springRotation.x);
-            constraint.enableSpring(3, true);
-            constraint.setStiffness(4, joint.springRotation.y);
-            constraint.enableSpring(4, true);
-            constraint.setStiffness(5, joint.springRotation.z);
-            constraint.enableSpring(5, true);
-            this.dynamicsWorld.addConstraint(constraint, false);
-            this.ammoConstraints.push(constraint);
-            Ammo.destroy(pivotInA);
-            Ammo.destroy(pivotInB);
-            Ammo.destroy(quatA);
-            Ammo.destroy(quatB);
-            Ammo.destroy(lowerLinear);
-            Ammo.destroy(upperLinear);
-            Ammo.destroy(lowerAngular);
-            Ammo.destroy(upperAngular);
-        }
-    }
-    // Normalize angle to [-π, π] range
-    normalizeAngle(angle) {
-        const pi = Math.PI;
-        const twoPi = 2 * pi;
-        angle = angle % twoPi;
-        if (angle < -pi) {
-            angle += twoPi;
-        }
-        else if (angle > pi) {
-            angle -= twoPi;
-        }
-        return angle;
-    }
-    // Re-snap all rigidbodies to current bone poses and zero velocities / forces.
-    // Use when simulation has diverged (explosion, NaN, extreme external teleport).
-    reset(boneWorldMatrices) {
-        if (!this.ammoInitialized || !this.ammo || !this.dynamicsWorld)
-            return;
-        if (!this.rigidbodiesInitialized)
-            return;
-        this.positionBodiesFromBones(boneWorldMatrices, boneWorldMatrices.length);
-        if (this.dynamicsWorld.clearForces) {
-            this.dynamicsWorld.clearForces();
-        }
-        if (this.dynamicsWorld.stepSimulation) {
-            this.dynamicsWorld.stepSimulation(0, 0, 0);
-        }
-    }
-    // Syncs bones to rigidbodies, simulates dynamics, solves constraints
-    // Modifies boneWorldMatrices in-place for dynamic rigidbodies that drive bones
-    step(dt, boneWorldMatrices, boneInverseBindMatrices) {
-        // Wait for Ammo to initialize
-        if (!this.ammoInitialized || !this.ammo || !this.dynamicsWorld) {
-            return;
-        }
-        const boneCount = boneWorldMatrices.length;
-        if (this.firstFrame) {
-            if (!this.rigidbodiesInitialized) {
-                this.computeBodyOffsets(boneInverseBindMatrices, boneCount);
-                this.rigidbodiesInitialized = true;
-            }
-            // Position bodies based on current bone poses (not bind pose) before creating joints
-            this.positionBodiesFromBones(boneWorldMatrices, boneCount);
-            if (!this.jointsCreated) {
-                this.createAmmoJoints();
-                this.jointsCreated = true;
-            }
-            if (this.dynamicsWorld.stepSimulation) {
-                this.dynamicsWorld.stepSimulation(0, 0, 0);
-            }
-            this.firstFrame = false;
-        }
-        // Step order: 1) Sync Static/Kinematic from bones, 2) Step physics, 3) Apply dynamic to bones
-        this.syncFromBones(boneWorldMatrices, boneCount);
-        this.stepAmmoPhysics(dt);
-        this.applyAmmoRigidbodiesToBones(boneWorldMatrices, boneCount);
-    }
-    // Compute bodyOffsetMatrixInverse for all rigidbodies (called once during initialization)
-    computeBodyOffsets(boneInverseBindMatrices, boneCount) {
-        if (!this.ammo || !this.dynamicsWorld)
-            return;
-        for (let i = 0; i < this.rigidbodies.length; i++) {
-            const rb = this.rigidbodies[i];
-            if (rb.boneIndex >= 0 && rb.boneIndex < boneCount) {
-                const boneIdx = rb.boneIndex;
-                const invBindIdx = boneIdx * 16;
-                const invBindMat = new Mat4(boneInverseBindMatrices.subarray(invBindIdx, invBindIdx + 16));
-                // Compute shape transform in bone-local space: shapeLocal = boneInverseBind × shapeWorldBind
-                const shapeRotQuat = Quat.fromEuler(rb.shapeRotation.x, rb.shapeRotation.y, rb.shapeRotation.z);
-                const shapeWorldBind = Mat4.fromPositionRotation(rb.shapePosition, shapeRotQuat);
-                // shapeLocal = boneInverseBind × shapeWorldBind (not shapeWorldBind × boneInverseBind)
-                const bodyOffsetMatrix = invBindMat.multiply(shapeWorldBind);
-                rb.bodyOffsetMatrixInverse = bodyOffsetMatrix.inverse();
-                rb.bodyOffsetMatrix = bodyOffsetMatrix; // Cache non-inverse to avoid expensive inverse() calls
-            }
-            else {
-                rb.bodyOffsetMatrixInverse = Mat4.identity();
-                rb.bodyOffsetMatrix = Mat4.identity(); // Cache non-inverse
-            }
-        }
-    }
-    // Position bodies based on current bone transforms (called on first frame only)
-    positionBodiesFromBones(boneWorldMatrices, boneCount) {
-        if (!this.ammo || !this.dynamicsWorld)
-            return;
-        const Ammo = this.ammo;
-        for (let i = 0; i < this.rigidbodies.length; i++) {
-            const rb = this.rigidbodies[i];
-            const ammoBody = this.ammoRigidbodies[i];
-            if (!ammoBody || rb.boneIndex < 0 || rb.boneIndex >= boneCount)
-                continue;
-            const boneIdx = rb.boneIndex;
-            const boneWorldMat = boneWorldMatrices[boneIdx];
-            // nodeWorld = boneWorld × shapeLocal (not shapeLocal × boneWorld)
-            const bodyOffsetMatrix = rb.bodyOffsetMatrix || rb.bodyOffsetMatrixInverse.inverse();
-            const nodeWorldMatrix = boneWorldMat.multiply(bodyOffsetMatrix);
-            const worldPos = nodeWorldMatrix.getPosition();
-            const worldRot = nodeWorldMatrix.toQuat();
-            const transform = new Ammo.btTransform();
-            const pos = new Ammo.btVector3(worldPos.x, worldPos.y, worldPos.z);
-            const quat = new Ammo.btQuaternion(worldRot.x, worldRot.y, worldRot.z, worldRot.w);
-            transform.setOrigin(pos);
-            transform.setRotation(quat);
-            if (rb.type === RigidbodyType.Static || rb.type === RigidbodyType.Kinematic) {
-                ammoBody.setCollisionFlags(ammoBody.getCollisionFlags() | 2); // CF_KINEMATIC_OBJECT
-                ammoBody.setActivationState(4); // DISABLE_DEACTIVATION
-            }
-            ammoBody.setWorldTransform(transform);
-            ammoBody.getMotionState().setWorldTransform(transform);
-            if (!this.zeroVector) {
-                this.zeroVector = new Ammo.btVector3(0, 0, 0);
-            }
-            ammoBody.setLinearVelocity(this.zeroVector);
-            ammoBody.setAngularVelocity(this.zeroVector);
-            Ammo.destroy(pos);
-            Ammo.destroy(quat);
-            Ammo.destroy(transform);
-        }
-    }
-    // Sync Static (FollowBone) and Kinematic rigidbodies to follow bone transforms
-    syncFromBones(boneWorldMatrices, boneCount) {
-        if (!this.ammo || !this.dynamicsWorld)
-            return;
-        const Ammo = this.ammo;
-        for (let i = 0; i < this.rigidbodies.length; i++) {
-            const rb = this.rigidbodies[i];
-            const ammoBody = this.ammoRigidbodies[i];
-            if (!ammoBody)
-                continue;
-            // Sync both Static (FollowBone) and Kinematic bodies - they both follow bones
-            if ((rb.type === RigidbodyType.Static || rb.type === RigidbodyType.Kinematic) &&
-                rb.boneIndex >= 0 &&
-                rb.boneIndex < boneCount) {
-                const boneIdx = rb.boneIndex;
-                const boneWorldMat = boneWorldMatrices[boneIdx];
-                // Lazy-cache bodyOffsetMatrix on first hit (cold path).
-                if (!rb.bodyOffsetMatrix)
-                    rb.bodyOffsetMatrix = rb.bodyOffsetMatrixInverse.inverse();
-                // nodeWorld = boneWorld × bodyOffsetMatrix → _physMat[0]
-                Mat4.multiplyArrays(boneWorldMat.values, 0, rb.bodyOffsetMatrix.values, 0, _physMat[0], 0);
-                const nodeVals = _physMat[0];
-                const wx = nodeVals[12], wy = nodeVals[13], wz = nodeVals[14];
-                Mat4.toQuatFromArrayInto(nodeVals, 0, _physQuat);
-                const transform = new Ammo.btTransform();
-                const pos = new Ammo.btVector3(wx, wy, wz);
-                const quat = new Ammo.btQuaternion(_physQuat.x, _physQuat.y, _physQuat.z, _physQuat.w);
-                transform.setOrigin(pos);
-                transform.setRotation(quat);
-                ammoBody.setWorldTransform(transform);
-                ammoBody.getMotionState().setWorldTransform(transform);
-                if (!this.zeroVector) {
-                    this.zeroVector = new Ammo.btVector3(0, 0, 0);
-                }
-                ammoBody.setLinearVelocity(this.zeroVector);
-                ammoBody.setAngularVelocity(this.zeroVector);
-                Ammo.destroy(pos);
-                Ammo.destroy(quat);
-                Ammo.destroy(transform);
-            }
-        }
-    }
-    // Step Ammo physics simulation
-    stepAmmoPhysics(dt) {
-        if (!this.ammo || !this.dynamicsWorld)
-            return;
-        const fixedTimeStep = 1 / 75;
-        const maxSubSteps = 10;
-        this.dynamicsWorld.stepSimulation(dt, maxSubSteps, fixedTimeStep);
-    }
-    // Apply dynamic rigidbody world transforms to bone world matrices in-place
-    applyAmmoRigidbodiesToBones(boneWorldMatrices, boneCount) {
-        if (!this.ammo || !this.dynamicsWorld)
-            return;
-        for (let i = 0; i < this.rigidbodies.length; i++) {
-            const rb = this.rigidbodies[i];
-            const ammoBody = this.ammoRigidbodies[i];
-            if (!ammoBody)
-                continue;
-            // Only dynamic rigidbodies drive bones (Static/Kinematic follow bones)
-            if (rb.type === RigidbodyType.Dynamic && rb.boneIndex >= 0 && rb.boneIndex < boneCount) {
-                const boneIdx = rb.boneIndex;
-                const transform = ammoBody.getWorldTransform();
-                const origin = transform.getOrigin();
-                const rotation = transform.getRotation();
-                // nodeWorldMatrix → _physMat[0] (from ammo position/rotation directly)
-                Mat4.fromPositionRotationInto(origin.x(), origin.y(), origin.z(), rotation.x(), rotation.y(), rotation.z(), rotation.w(), _physMat[0]);
-                // boneWorld = nodeWorld × bodyOffsetMatrixInverse → _physMat[1]
-                const boneVals = _physMat[1];
-                Mat4.multiplyArrays(_physMat[0], 0, rb.bodyOffsetMatrixInverse.values, 0, boneVals, 0);
-                if (!isNaN(boneVals[0]) && !isNaN(boneVals[15]) && Math.abs(boneVals[0]) < 1e6 && Math.abs(boneVals[15]) < 1e6) {
-                    boneWorldMatrices[boneIdx].values.set(boneVals);
-                }
-                else {
-                    console.warn(`[Physics] Invalid bone world matrix for rigidbody ${i} (${rb.name}), skipping update`);
-                }
-            }
-        }
-    }
-}

package/dist/shaders/body.d.ts

@@ -1,2 +0,0 @@
-export declare const BODY_SHADER_WGSL = "\n\n\n\n// Baked 64\u00D764 rgba8unorm combined BRDF LUT \u2014 created once at engine init by dfg_lut.ts.\n//   .rg = split-sum DFG (Karis: tint = f0\u00B7x + f90\u00B7y)  \u2192 F_brdf_*_scatter\n//   .ba = Heitz 2016 LTC magnitude (ltc_mag_ggx)       \u2192 ltc_brdf_scale_from_lut\n// Paired with group(0) binding(2) diffuseSampler (linear filter). Sample once per\n// fragment via brdf_lut_sample() \u2014 callers feed .rg and the whole vec4 into the\n// helpers below, halving LUT taps on the default Principled path.\n@group(0) @binding(9) var brdfLut: texture_2d<f32>;\n\n// \u2500\u2500\u2500 RGB \u2194 HSV \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n\nfn rgb_to_hsv(rgb: vec3f) -> vec3f {\n  let c_max = max(rgb.r, max(rgb.g, rgb.b));\n  let c_min = min(rgb.r, min(rgb.g, rgb.b));\n  let delta = c_max - c_min;\n\n  var h = 0.0;\n  if (delta > 1e-6) {\n    if (c_max == rgb.r) {\n      h = (rgb.g - rgb.b) / delta;\n      if (h < 0.0) { h += 6.0; }\n    } else if (c_max == rgb.g) {\n      h = 2.0 + (rgb.b - rgb.r) / delta;\n    } else {\n      h = 4.0 + (rgb.r - rgb.g) / delta;\n    }\n    h /= 6.0;\n  }\n  let s = select(0.0, delta / c_max, c_max > 1e-6);\n  return vec3f(h, s, c_max);\n}\n\nfn hsv_to_rgb(hsv: vec3f) -> vec3f {\n  let h = hsv.x;\n  let s = hsv.y;\n  let v = hsv.z;\n  if (s < 1e-6) { return vec3f(v); }\n\n  let hh = fract(h) * 6.0;\n  let sector = u32(hh);\n  let f = hh - f32(sector);\n  let p = v * (1.0 - s);\n  let q = v * (1.0 - s * f);\n  let t = v * (1.0 - s * (1.0 - f));\n\n  switch (sector) {\n    case 0u: { return vec3f(v, t, p); }\n    case 1u: { return vec3f(q, v, p); }\n    case 2u: { return vec3f(p, v, t); }\n    case 3u: { return vec3f(p, q, v); }\n    case 4u: { return vec3f(t, p, v); }\n    default: { return vec3f(v, p, q); }\n  }\n}\n\n// \u2500\u2500\u2500 HUE_SAT node \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n\nfn hue_sat(hue: f32, saturation: f32, value: f32, fac: f32, color: vec3f) -> vec3f {\n  var hsv = rgb_to_hsv(color);\n  hsv.x = fract(hsv.x + hue - 0.5);\n  hsv.y = clamp(hsv.y * saturation, 0.0, 1.0);\n  hsv.z *= value;\n  return mix(color, hsv_to_rgb(hsv), fac);\n}\n\n// hue_sat specialization for hue=0.5 (identity hue shift \u2014 fract(h + 0.5 - 0.5) = h).\n// Branchless equivalent that skips the rgb_to_hsv \u2192 hsv_to_rgb roundtrip: WebKit's\n// Metal backend serializes the 3-way if chain in rgb_to_hsv and the 6-way switch in\n// hsv_to_rgb, where this form compiles to linear SIMD ops + a single select.\nfn hue_sat_id(saturation: f32, value: f32, fac: f32, color: vec3f) -> vec3f {\n  let m = max(max(color.r, color.g), color.b);\n  let n = min(min(color.r, color.g), color.b);\n  // Unclamped (sat*old_s \u2264 1): reproj = mix(vec3f(m), color, saturation).\n  // Clamped (saturated to 1):   reproj = (color - n) * m / (m - n).\n  let range = max(m - n, 1e-6);\n  let unclamped = mix(vec3f(m), color, saturation);\n  let clamped = (color - vec3f(n)) * m / range;\n  let needs_clamp = (m - n) * saturation >= m;\n  let reproj = select(unclamped, clamped, needs_clamp);\n  return mix(color, reproj * value, fac);\n}\n\n// \u2500\u2500\u2500 BRIGHTCONTRAST node \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n\nfn bright_contrast(color: vec3f, bright: f32, contrast: f32) -> vec3f {\n  let a = 1.0 + contrast;\n  let b = bright - contrast * 0.5;\n  return max(vec3f(0.0), color * a + vec3f(b));\n}\n\n// \u2500\u2500\u2500 INVERT node \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n\nfn invert(fac: f32, color: vec3f) -> vec3f {\n  return mix(color, vec3f(1.0) - color, fac);\n}\n\nfn invert_f(fac: f32, val: f32) -> f32 {\n  return mix(val, 1.0 - val, fac);\n}\n\n// \u2500\u2500\u2500 Color ramp (VALTORGB) \u2014 2-stop variants \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// All 7 presets use exclusively 2-stop ramps.\n\nfn ramp_constant(f: f32, p0: f32, c0: vec4f, p1: f32, c1: vec4f) -> vec4f {\n  return select(c0, c1, f >= p1);\n}\n\n// CONSTANT ramp with screen-space edge AA \u2014 kills sparkle where fwidth(f) straddles a hard step (NPR terminator)\nfn ramp_constant_edge_aa(f: f32, edge: f32, c0: vec4f, c1: vec4f) -> vec4f {\n  let w = max(fwidth(f) * 1.75, 6e-6);\n  let t = smoothstep(edge - w, edge + w, f);\n  return mix(c0, c1, t);\n}\n\nfn ramp_linear(f: f32, p0: f32, c0: vec4f, p1: f32, c1: vec4f) -> vec4f {\n  let t = saturate((f - p0) / max(p1 - p0, 1e-6));\n  return mix(c0, c1, t);\n}\n\nfn ramp_cardinal(f: f32, p0: f32, c0: vec4f, p1: f32, c1: vec4f) -> vec4f {\n  // cardinal spline with 2 stops degrades to smoothstep\n  let t = saturate((f - p0) / max(p1 - p0, 1e-6));\n  let ss = t * t * (3.0 - 2.0 * t);\n  return mix(c0, c1, ss);\n}\n\n// \u2500\u2500\u2500 MATH node operations \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n\nfn math_add(a: f32, b: f32) -> f32 { return a + b; }\nfn math_multiply(a: f32, b: f32) -> f32 { return a * b; }\nfn math_power(a: f32, b: f32) -> f32 { return pow(max(a, 0.0), b); }\nfn math_greater_than(a: f32, b: f32) -> f32 { return select(0.0, 1.0, a > b); }\n\n// Blender's implicit Color \u2192 Float socket conversion uses BT.601 grayscale\n// (rgb_to_grayscale in blenkernel/intern/node.cc). When a material graph plugs a\n// Color output into a Math node's Value input, this is the scalar it actually sees.\nfn color_to_value(c: vec3f) -> f32 {\n  return 0.299 * c.r + 0.587 * c.g + 0.114 * c.b;\n}\n\n// \u2500\u2500\u2500 MIX node (blend_type variants) \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n\nfn mix_blend(fac: f32, a: vec3f, b: vec3f) -> vec3f {\n  return mix(a, b, fac);\n}\n\nfn mix_overlay(fac: f32, a: vec3f, b: vec3f) -> vec3f {\n  let lo = 2.0 * a * b;\n  let hi = vec3f(1.0) - 2.0 * (vec3f(1.0) - a) * (vec3f(1.0) - b);\n  let overlay = select(hi, lo, a < vec3f(0.5));\n  return mix(a, overlay, fac);\n}\n\nfn mix_multiply(fac: f32, a: vec3f, b: vec3f) -> vec3f {\n  return mix(a, a * b, fac);\n}\n\nfn mix_lighten(fac: f32, a: vec3f, b: vec3f) -> vec3f {\n  return mix(a, max(a, b), fac);\n}\n\n// Blender Mix (Color) blend LINEAR_LIGHT: result = mix(A, A + 2*B - 1, Fac)\nfn mix_linear_light(fac: f32, a: vec3f, b: vec3f) -> vec3f {\n  return mix(a, a + 2.0 * b - vec3f(1.0), fac);\n}\n\n// Luminance for Shader\u2192RGB scalar gates (linear RGB, Rec.709 weights)\nfn luminance_rec709_linear(c: vec3f) -> f32 {\n  return dot(max(c, vec3f(0.0)), vec3f(0.2126, 0.7152, 0.0722));\n}\n\n// \u2500\u2500\u2500 FRESNEL node \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// Schlick approximation matching Blender's Fresnel node\n\nfn fresnel(ior: f32, n: vec3f, v: vec3f) -> f32 {\n  let r = (ior - 1.0) / (ior + 1.0);\n  let f0 = r * r;\n  let cos_theta = clamp(dot(n, v), 0.0, 1.0);\n  let m = 1.0 - cos_theta;\n  let m2 = m * m;\n  let m5 = m2 * m2 * m;\n  return f0 + (1.0 - f0) * m5;\n}\n\n// \u2500\u2500\u2500 LAYER_WEIGHT node \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n\nfn layer_weight_fresnel(blend: f32, n: vec3f, v: vec3f) -> f32 {\n  let eta = max(1.0 - blend, 1e-4);\n  let r = (1.0 - eta) / (1.0 + eta);\n  let f0 = r * r;\n  let cos_theta = clamp(abs(dot(n, v)), 0.0, 1.0);\n  let m = 1.0 - cos_theta;\n  let m2 = m * m;\n  let m5 = m2 * m2 * m;\n  return f0 + (1.0 - f0) * m5;\n}\n\nfn layer_weight_facing(blend: f32, n: vec3f, v: vec3f) -> f32 {\n  var facing = abs(dot(n, v));\n  let b = clamp(blend, 0.0, 0.99999);\n  if (b != 0.5) {\n    let exponent = select(2.0 * b, 0.5 / (1.0 - b), b >= 0.5);\n    facing = pow(facing, exponent);\n  }\n  return 1.0 - facing;\n}\n\n// \u2500\u2500\u2500 SHADER_TO_RGB (white DiffuseBSDF) \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// Eevee captures lit diffuse: (albedo/\u03C0)*sun*N\u00B7L*shadow + ambient (linear). Albedo=1.\n// Matches default.ts direct term scale so VALTORGB thresholds from Blender JSON stay valid.\n\nfn shader_to_rgb_diffuse(n: vec3f, l: vec3f, sun_rgb: vec3f, ambient_rgb: vec3f, shadow: f32) -> f32 {\n  const PI_S: f32 = 3.141592653589793;\n  let ndotl = max(dot(n, l), 0.0);\n  let rgb = sun_rgb * (ndotl * shadow / PI_S) + ambient_rgb;\n  return luminance_rec709_linear(rgb);\n}\n\n// \u2500\u2500\u2500 BUMP node \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// Screen-space bump from a scalar height field. Needs dFdx/dFdy which\n// WGSL provides as dpdx/dpdy.\n\nfn bump(strength: f32, height: f32, normal: vec3f, world_pos: vec3f) -> vec3f {\n  let dhdx = dpdx(height);\n  let dhdy = dpdy(height);\n  let dpdx_pos = dpdx(world_pos);\n  let dpdy_pos = dpdy(world_pos);\n  let perturbed = normalize(normal) - strength * (dhdx * normalize(cross(dpdy_pos, normal)) + dhdy * normalize(cross(normal, dpdx_pos)));\n  return normalize(perturbed);\n}\n\n// LH engine + WebGPU fragment Y: flip dhdy contribution so height peaks read as outward bumps vs Blender reference\nfn bump_lh(strength: f32, height: f32, normal: vec3f, world_pos: vec3f) -> vec3f {\n  let dhdx = dpdx(height);\n  let dhdy = dpdy(height);\n  let dpdx_pos = dpdx(world_pos);\n  let dpdy_pos = dpdy(world_pos);\n  let perturbed = normalize(normal) - strength * (dhdx * normalize(cross(dpdy_pos, normal)) - dhdy * normalize(cross(normal, dpdx_pos)));\n  return normalize(perturbed);\n}\n\n// \u2500\u2500\u2500 NOISE texture (Perlin-style) \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// Simplified gradient noise matching Blender's default noise output.\n\n// PCG-style integer hash. Replaces the classic 'fract(sin(q) * LARGE)' trick because\n// WebKit's Metal backend compiles 'sin' to a full transcendental op (slow), while\n// Safari's Apple-GPU scalar ALU handles int muls/xors near free. Inputs arrive as\n// integer-valued floats (floor(p) + unit offsets) from _noise3, so vec3i cast is exact.\nfn _hash33(p: vec3f) -> vec3f {\n  var h = vec3u(vec3i(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  // Mask to 24 bits \u2014 above that f32 loses precision on the u32\u2192f32 convert.\n  let hm = h & vec3u(16777215u);\n  return vec3f(hm) * (2.0 / 16777216.0) - 1.0;\n}\n\nfn _noise3(p: vec3f) -> f32 {\n  let i = floor(p);\n  let f = fract(p);\n  let u = f * f * (3.0 - 2.0 * f);\n\n  return mix(\n    mix(\n      mix(dot(_hash33(i + vec3f(0,0,0)), f - vec3f(0,0,0)),\n          dot(_hash33(i + vec3f(1,0,0)), f - vec3f(1,0,0)), u.x),\n      mix(dot(_hash33(i + vec3f(0,1,0)), f - vec3f(0,1,0)),\n          dot(_hash33(i + vec3f(1,1,0)), f - vec3f(1,1,0)), u.x), u.y),\n    mix(\n      mix(dot(_hash33(i + vec3f(0,0,1)), f - vec3f(0,0,1)),\n          dot(_hash33(i + vec3f(1,0,1)), f - vec3f(1,0,1)), u.x),\n      mix(dot(_hash33(i + vec3f(0,1,1)), f - vec3f(0,1,1)),\n          dot(_hash33(i + vec3f(1,1,1)), f - vec3f(1,1,1)), u.x), u.y),\n    u.z);\n}\n\nfn tex_noise(p: vec3f, scale: f32, detail: f32, roughness: f32, distortion: f32) -> f32 {\n  var q = p;\n  if (abs(distortion) > 1e-6) {\n    let w = _noise3(p * scale * 1.37 + vec3f(2.31, 5.17, 8.09));\n    q = p + (w * 2.0 - 1.0) * distortion;\n  }\n  let coords = q * scale;\n  var value = 0.0;\n  var amplitude = 1.0;\n  var frequency = 1.0;\n  var total_amp = 0.0;\n  let octaves = i32(clamp(detail, 0.0, 15.0)) + 1;\n  for (var i = 0; i < octaves; i++) {\n    value += amplitude * _noise3(coords * frequency);\n    total_amp += amplitude;\n    amplitude *= roughness;\n    frequency *= 2.0;\n  }\n  return value / max(total_amp, 1e-6) * 0.5 + 0.5;\n}\n\n// tex_noise specialization: detail=2.0 (3 octaves), roughness=0.5, distortion=0.\n// WebKit can't unroll tex_noise's for-loop because 'octaves' is a runtime value;\n// this variant is fully unrolled with constants folded (total_amp = 1.75).\nfn tex_noise_d2(p: vec3f, scale: f32) -> f32 {\n  let c = p * scale;\n  let v = _noise3(c) + 0.5 * _noise3(c * 2.0) + 0.25 * _noise3(c * 4.0);\n  return v * (1.0 / 1.75) * 0.5 + 0.5;\n}\n\n// \u2500\u2500\u2500 TEX_GRADIENT (linear) \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// Used by Stockings preset. Maps the input vector's X to a 0\u20131 gradient.\n\nfn tex_gradient_linear(uv: vec3f) -> f32 {\n  return clamp(uv.x, 0.0, 1.0);\n}\n\n// \u2500\u2500\u2500 TEX_VORONOI (distance only) \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// Used by Metal preset. Simplified F1 cell noise.\n\nfn tex_voronoi_f1(p: vec3f, scale: f32) -> f32 {\n  let coords = p * scale;\n  let i = floor(coords);\n  let f = fract(coords);\n  var min_dist = 1e10;\n  for (var z = -1; z <= 1; z++) {\n    for (var y = -1; y <= 1; y++) {\n      for (var x = -1; x <= 1; x++) {\n        let neighbor = vec3f(f32(x), f32(y), f32(z));\n        let point = _hash33(i + neighbor) * 0.5 + 0.5;\n        let diff = neighbor + point - f;\n        min_dist = min(min_dist, dot(diff, diff));\n      }\n    }\n  }\n  return sqrt(min_dist);\n}\n\n// \u2500\u2500\u2500 SEPXYZ node \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n\nfn separate_xyz(v: vec3f) -> vec3f { return v; }\n\n// \u2500\u2500\u2500 VECT_MATH (cross product) \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n\nfn vect_math_cross(a: vec3f, b: vec3f) -> vec3f { return cross(a, b); }\n\n// \u2500\u2500\u2500 MAPPING node \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// Point-type mapping: scale, rotate (euler XYZ), translate.\n\nfn mapping_point(v: vec3f, loc: vec3f, rot: vec3f, scl: vec3f) -> vec3f {\n  var p = v * scl;\n  // simplified: skip rotation when all angles are zero (common case)\n  if (abs(rot.x) + abs(rot.y) + abs(rot.z) > 1e-6) {\n    let cx = cos(rot.x); let sx = sin(rot.x);\n    let cy = cos(rot.y); let sy = sin(rot.y);\n    let cz = cos(rot.z); let sz = sin(rot.z);\n    let rx = vec3f(p.x, cx*p.y - sx*p.z, sx*p.y + cx*p.z);\n    let ry = vec3f(cy*rx.x + sy*rx.z, rx.y, -sy*rx.x + cy*rx.z);\n    p = vec3f(cz*ry.x - sz*ry.y, sz*ry.x + cz*ry.y, ry.z);\n  }\n  return p + loc;\n}\n\n// \u2500\u2500\u2500 NORMAL_MAP node (tangent-space) \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// Applies a tangent-space normal map. Requires TBN from vertex stage.\n\nfn normal_map(strength: f32, map_color: vec3f, normal: vec3f, tangent: vec3f, bitangent: vec3f) -> vec3f {\n  let ts = map_color * 2.0 - 1.0;\n  let perturbed = normalize(tangent * ts.x + bitangent * ts.y + normal * ts.z);\n  return normalize(mix(normal, perturbed, strength));\n}\n\n// \u2500\u2500\u2500 EEVEE Principled BSDF primitives \u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n// Ports from Blender 3.6 source/blender/draw/engines/eevee/shaders/\n//   bsdf_common_lib.glsl + gpu_shader_material_principled.glsl.\n// Usage pattern (see material shaders): direct spec = bsdf_ggx \u00D7 sun \u00D7 shadow\n// (NL baked in, no F yet); ambient spec = probe_radiance; tint both with\n// reflection_color = F_brdf_multi_scatter(f0, f90, split_sum) AFTER summing.\n\nconst EEVEE_PI: f32 = 3.141592653589793;\n\n// Fused analytic GGX specular (direct lights). Returns BRDF \u00D7 NL.\n// 4\u00B7NL\u00B7NV is cancelled via G1_Smith reciprocal form \u2014 see bsdf_common_lib.glsl:115.\n// Caller passes NL, NV (already computed for diffuse + brdf_lut_sample) so WebKit\n// can reuse them instead of recomputing dot products across the function boundary.\nfn bsdf_ggx(N: vec3f, L: vec3f, V: vec3f, NL_in: f32, NV_in: f32, roughness: f32) -> f32 {\n  let a = max(roughness, 1e-4);\n  let a2 = a * a;\n  let H = normalize(L + V);\n  let NH = max(dot(N, H), 1e-8);\n  let NL = max(NL_in, 1e-8);\n  let NV = max(NV_in, 1e-8);\n  // G1_Smith_GGX_opti reciprocal form \u2014 denominator piece only.\n  let G1L = NL + sqrt(NL * (NL - NL * a2) + a2);\n  let G1V = NV + sqrt(NV * (NV - NV * a2) + a2);\n  let G = G1L * G1V;\n  // D_ggx_opti = pi * denom\u00B2 \u2014 reciprocal D \u00D7 a\u00B2.\n  let tmp = (NH * a2 - NH) * NH + 1.0;\n  let D_opti = EEVEE_PI * tmp * tmp;\n  return NL * a2 / (D_opti * G);\n}\n\n// Split-sum DFG LUT \u2014 Karis 2013 curve fit stand-in for the 64\u00D764 baked LUT.\n// Returns (lut.x, lut.y) in Blender convention: tint = f0\u00B7lut.x + f90\u00B7lut.y.\nfn brdf_lut_approx(NV: f32, roughness: f32) -> vec2f {\n  let c0 = vec4f(-1.0, -0.0275, -0.572, 0.022);\n  let c1 = vec4f(1.0, 0.0425, 1.04, -0.04);\n  let r = roughness * c0 + c1;\n  let a004 = min(r.x * r.x, exp2(-9.28 * NV)) * r.x + r.y;\n  return vec2f(-1.04, 1.04) * a004 + r.zw;\n}\n\n// Baked combined BRDF LUT \u2014 exact port of Blender bsdf_lut_frag.glsl packed with\n// ltc_mag_ggx from eevee_lut.c. Single sample returns DFG (.rg) and LTC mag (.ba).\n// Addressed as Blender's common_utiltex_lib.glsl:lut_coords:\n//   coords = (roughness, sqrt(1 - NV)), then half-texel bias for filtering.\n// Requires group(0) binding(9) brdfLut + binding(2) diffuseSampler in the host shader.\nfn brdf_lut_sample(NV: f32, roughness: f32) -> vec4f {\n  let LUT_SIZE: f32 = 64.0;\n  var uv = vec2f(saturate(roughness), sqrt(saturate(1.0 - NV)));\n  uv = uv * ((LUT_SIZE - 1.0) / LUT_SIZE) + 0.5 / LUT_SIZE;\n  return textureSampleLevel(brdfLut, diffuseSampler, uv, 0.0);\n}\n\nfn F_brdf_single_scatter(f0: vec3f, f90: vec3f, lut: vec2f) -> vec3f {\n  return lut.y * f90 + lut.x * f0;\n}\n\n// Fdez-Ag\u00FCera 2019 multi-scatter compensation (EEVEE do_multiscatter=1).\nfn F_brdf_multi_scatter(f0: vec3f, f90: vec3f, lut: vec2f) -> vec3f {\n  let FssEss = lut.y * f90 + lut.x * f0;\n  let Ess = lut.x + lut.y;\n  let Ems = 1.0 - Ess;\n  let Favg = f0 + (1.0 - f0) / 21.0;\n  let Fms = FssEss * Favg / (1.0 - (1.0 - Ess) * Favg);\n  return FssEss + Fms * Ems;\n}\n\n// EEVEE direct-specular energy compensation factor \u2014 closure_eval_glossy_lib.glsl:79-81:\n//   ltc_brdf_scale = (ltc.x + ltc.y) / (split_sum.x + split_sum.y)\n// Blender evaluates direct lights via LTC (Heitz 2016) but indirect via split-sum;\n// direct radiance is rescaled so total-energy matches the split-sum LUT.\n// Takes a pre-sampled vec4f from brdf_lut_sample() to share the fetch with\n// F_brdf_multi_scatter on the same fragment.\nfn ltc_brdf_scale_from_lut(lut: vec4f) -> f32 {\n  return (lut.z + lut.w) / max(lut.x + lut.y, 1e-6);\n}\n\n// Luminance-normalized hue extraction \u2014 Blender tint_from_color (isolates hue+sat).\nfn tint_from_color(color: vec3f) -> vec3f {\n  let lum = dot(color, vec3f(0.3, 0.6, 0.1));\n  return select(vec3f(1.0), color / lum, lum > 0.0);\n}\n\n\n\nstruct CameraUniforms {\n  view: mat4x4f,\n  projection: mat4x4f,\n  viewPos: vec3f,\n  _padding: f32,\n};\n\nstruct Light {\n  direction: vec4f,\n  color: vec4f,\n};\n\nstruct LightUniforms {\n  ambientColor: vec4f,\n  lights: array<Light, 4>,\n};\n\nstruct MaterialUniforms {\n  diffuseColor: vec3f,\n  alpha: f32,\n};\n\nstruct VertexOutput {\n  @builtin(position) position: vec4f,\n  @location(0) normal: vec3f,\n  @location(1) uv: vec2f,\n  @location(2) worldPos: vec3f,\n};\n\nstruct LightVP { viewProj: mat4x4f, };\n\n@group(0) @binding(0) var<uniform> camera: CameraUniforms;\n@group(0) @binding(1) var<uniform> light: LightUniforms;\n@group(0) @binding(2) var diffuseSampler: sampler;\n@group(0) @binding(3) var shadowMap: texture_depth_2d;\n@group(0) @binding(4) var shadowSampler: sampler_comparison;\n@group(0) @binding(5) var<uniform> lightVP: LightVP;\n@group(1) @binding(0) var<storage, read> skinMats: array<mat4x4f>;\n@group(2) @binding(0) var diffuseTexture: texture_2d<f32>;\n@group(2) @binding(1) var<uniform> material: MaterialUniforms;\n\nfn sampleShadow(worldPos: vec3f, n: vec3f) -> f32 {\n  // Back-facing to key light: direct contribution is zero anyway, skip 9 texture samples.\n  if (dot(n, -light.lights[0].direction.xyz) <= 0.0) { return 0.0; }\n  let biasedPos = worldPos + n * 0.08;\n  let lclip = lightVP.viewProj * vec4f(biasedPos, 1.0);\n  let ndc = lclip.xyz / max(lclip.w, 1e-6);\n  let suv = vec2f(ndc.x * 0.5 + 0.5, 0.5 - ndc.y * 0.5);\n  let cmpZ = ndc.z - 0.001;\n  let ts = 1.0 / 2048.0;\n  // 3x3 PCF unrolled \u2014 Safari's Metal backend doesn't unroll nested shadow loops reliably.\n  let s00 = textureSampleCompareLevel(shadowMap, shadowSampler, suv + vec2f(-ts, -ts), cmpZ);\n  let s10 = textureSampleCompareLevel(shadowMap, shadowSampler, suv + vec2f(0.0, -ts), cmpZ);\n  let s20 = textureSampleCompareLevel(shadowMap, shadowSampler, suv + vec2f( ts, -ts), cmpZ);\n  let s01 = textureSampleCompareLevel(shadowMap, shadowSampler, suv + vec2f(-ts, 0.0), cmpZ);\n  let s11 = textureSampleCompareLevel(shadowMap, shadowSampler, suv, cmpZ);\n  let s21 = textureSampleCompareLevel(shadowMap, shadowSampler, suv + vec2f( ts, 0.0), cmpZ);\n  let s02 = textureSampleCompareLevel(shadowMap, shadowSampler, suv + vec2f(-ts,  ts), cmpZ);\n  let s12 = textureSampleCompareLevel(shadowMap, shadowSampler, suv + vec2f(0.0,  ts), cmpZ);\n  let s22 = textureSampleCompareLevel(shadowMap, shadowSampler, suv + vec2f( ts,  ts), cmpZ);\n  return (s00 + s10 + s20 + s01 + s11 + s21 + s02 + s12 + s22) * (1.0 / 9.0);\n}\n\nconst PI_B: f32 = 3.141592653589793;\nconst BODY_ROUGHNESS: f32 = 0.3;\nconst BODY_RIM2_LAYER_BLEND: f32 = 0.20000000298023224;\nconst BODY_RIM2_POW: f32 = 1.4300000667572021;\nconst BODY_RIM2_BG: vec3f = vec3f(1.0, 0.4303792119026184, 0.3315804898738861);\nconst BODY_WARM_STR: f32 = 0.30000001192092896;\nconst BODY_SPECULAR: f32 = 0.5;\nconst BODY_MIX_NPR: f32 = 0.5;\n// EEVEE Light Clamp equivalent \u2014 caps firefly specular from noise-bumped NDF aliasing.\nconst BODY_SPEC_CLAMP: f32 = 10.0;\n\n// smoothstep-based ramp: t*t*(3-2*t) between two color stops\nfn ramp_ease(f: f32, p0: f32, c0: vec4f, p1: f32, c1: vec4f) -> vec4f {\n  let t = saturate((f - p0) / max(p1 - p0, 1e-6));\n  let ss = t * t * (3.0 - 2.0 * t);\n  return mix(c0, c1, ss);\n}\n\n@vertex fn vs(\n  @location(0) position: vec3f,\n  @location(1) normal: vec3f,\n  @location(2) uv: vec2f,\n  @location(3) joints0: vec4<u32>,\n  @location(4) weights0: vec4<f32>\n) -> VertexOutput {\n  var output: VertexOutput;\n  let pos4 = vec4f(position, 1.0);\n  let weightSum = weights0.x + weights0.y + weights0.z + weights0.w;\n  let invWeightSum = select(1.0, 1.0 / weightSum, weightSum > 0.0001);\n  let nw = select(vec4f(1.0, 0.0, 0.0, 0.0), weights0 * invWeightSum, weightSum > 0.0001);\n  var skinnedPos = vec4f(0.0);\n  var skinnedNrm = vec3f(0.0);\n  for (var i = 0u; i < 4u; i++) {\n    let m = skinMats[joints0[i]];\n    let w = nw[i];\n    skinnedPos += (m * pos4) * w;\n    skinnedNrm += (mat3x3f(m[0].xyz, m[1].xyz, m[2].xyz) * normal) * w;\n  }\n  output.position = camera.projection * camera.view * vec4f(skinnedPos.xyz, 1.0);\n  // Skip VS normalize \u2014 interpolation denormalizes anyway, and FS always does normalize(input.normal).\n  output.normal = skinnedNrm;\n  output.uv = uv;\n  output.worldPos = skinnedPos.xyz;\n  return output;\n}\n\nstruct FSOut {\n  @location(0) color: vec4f,\n  @location(1) mask: f32,\n};\n\n@fragment fn fs(input: VertexOutput) -> FSOut {\n  let alpha = material.alpha;\n  if (alpha < 0.001) { discard; }\n\n  let n = normalize(input.normal);\n  let v = normalize(camera.viewPos - input.worldPos);\n  let l = -light.lights[0].direction.xyz;\n  let sun = light.lights[0].color.xyz * light.lights[0].color.w;\n\n  let tex_color = textureSample(diffuseTexture, diffuseSampler, input.uv).rgb;\n  let shadow = sampleShadow(input.worldPos, n);\n\n  let ndotl_raw = shader_to_rgb_diffuse(n, l, sun, light.ambientColor.xyz, shadow);\n  let toon = ramp_constant(ndotl_raw, 0.0, vec4f(0,0,0,1), 0.2966, vec4f(1,1,1,1)).r;\n\n  let shadow_tint = hue_sat_id(2.0, 0.3499999940395355, 1.0, tex_color);\n  let lit_tint = hue_sat_id(1.5, 1.0, 1.0, tex_color);\n  let toon_color = mix_blend(toon, shadow_tint, lit_tint);\n  let bc = bright_contrast(toon_color, 0.1, 0.2);\n\n  let emission3 = bc * 4.0;\n\n  let warm_input = clamp(toon + 0.5, 0.0, 1.0);\n  let warm_color = ramp_cardinal(warm_input, 0.2409,\n    vec4f(0.2426, 0.068, 0.0588, 1.0), 0.4663,\n    vec4f(0.6677, 0.5024, 0.5126, 1.0)).rgb;\n  let warm_emission = warm_color * BODY_WARM_STR;\n\n  let rim1_str = fresnel(2.0, n, v) * layer_weight_facing(0.24000005424022675, n, v);\n  let rim1 = vec3f(0.984157919883728, 0.6110184788703918, 0.5736401677131653) * rim1_str;\n\n  let facing_raw = layer_weight_facing(BODY_RIM2_LAYER_BLEND, n, v);\n  let facing_pow = math_power(facing_raw, BODY_RIM2_POW);\n  let rim2_fac = ramp_ease(facing_pow, 0.0, vec4f(0,0,0,1), 0.5052, vec4f(1,1,1,1)).r;\n  let rim2_mixed = mix(emission3, BODY_RIM2_BG, rim2_fac);\n\n  let npr_stack = rim1 + rim2_mixed + warm_emission;\n\n  // Noise bump \u2014 Mapping loc=rot=0 folds to a plain scale multiply.\n  let noise_val = tex_noise_d2(input.worldPos * vec3f(1.0, 1.0, 1.5), 1.0);\n  let noise_ramp = ramp_linear(noise_val, 0.0, vec4f(0,0,0,1), 1.0, vec4f(1,1,1,1)).r;\n  let bumped_n = bump_lh(0.324644535779953, noise_ramp, n, input.worldPos);\n\n  let principled_base = mix_blend(noise_ramp, bc, vec3f(0.6831911206245422, 0.19474034011363983, 0.13732507824897766));\n  let p_emission = bc * 0.2;\n\n  // Principled BSDF (EEVEE port): metallic=0, specular=0.5, roughness=0.3, specular_tint=0.\n  let NL = max(dot(bumped_n, l), 0.0);\n  let NV = max(dot(bumped_n, v), 1e-4);\n\n  // f0/f90 per gpu_shader_material_principled.glsl \u2014 specular_tint=0 \u2192 dielectric_f0_color=white.\n  let f0 = vec3f(0.08 * BODY_SPECULAR);\n  let f90 = mix(f0, vec3f(1.0), sqrt(BODY_SPECULAR));\n  let brdf_lut = brdf_lut_sample(NV, BODY_ROUGHNESS);\n  let reflection_color = F_brdf_multi_scatter(f0, f90, brdf_lut.xy);\n\n  // Direct glossy \u2014 bsdf_ggx already includes NL; no F applied here (tinted after accum).\n  // ltc_brdf_scale: EEVEE direct path uses LTC; split-sum LUT path is rescaled to match.\n  let spec_direct_raw = bsdf_ggx(bumped_n, l, v, NL, NV, BODY_ROUGHNESS) * sun * shadow * ltc_brdf_scale_from_lut(brdf_lut);\n  let spec_direct = min(spec_direct_raw, vec3f(BODY_SPEC_CLAMP));\n  // Indirect glossy \u2014 flat world probe (solid color). Phase 2 adds cubemap.\n  let spec_indirect = light.ambientColor.xyz;\n  let spec_radiance = (spec_direct + spec_indirect) * reflection_color;\n\n  // Indirect diffuse = base_color \u00D7 L_w per Blender closure_eval_surface_lib.glsl line 302;\n  // probe_evaluate_world_diff returns radiance (SH-projected, not cosine-convolved).\n  // No (1-F) factor per EEVEE \u2014 it doesn't energy-conserve spec<->diffuse.\n  let diffuse_radiance = principled_base * (sun * NL * shadow / PI_B + light.ambientColor.xyz);\n  let principled = diffuse_radiance + spec_radiance + p_emission;\n\n  let final_color = mix(npr_stack, principled, BODY_MIX_NPR);\n\n  var out: FSOut;\n  out.color = vec4f(final_color, alpha);\n  out.mask = 1.0;\n  return out;\n}\n\n";
-//# sourceMappingURL=body.d.ts.map

package/dist/shaders/body.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"body.d.ts","sourceRoot":"","sources":["../../src/shaders/body.ts"],"names":[],"mappings":"AAIA,eAAO,MAAM,gBAAgB,gu7BAoM5B,CAAA"}