reze-engine 0.10.2 → 0.11.0

package/src/shaders/default.ts

@@ -0,0 +1,169 @@
+// Blender 3.6 Principled BSDF defaults + Filmic "Medium High Contrast" tone mapping.
+// Metallic=0, Specular=0.5 (F0=0.04), Roughness=0.5.
+// Tone mapping via LUT sampled from Blender's OCIO pipeline (exposure -0.3 baked in).
+
+export const DEFAULT_SHADER_WGSL = /* wgsl */ `
+
+const PI: f32 = 3.141592653589793;
+const F0_DIELECTRIC: f32 = 0.04;
+const ROUGHNESS: f32 = 0.5;
+
+struct CameraUniforms {
+  view: mat4x4f,
+  projection: mat4x4f,
+  viewPos: vec3f,
+  _padding: f32,
+};
+
+struct Light {
+  direction: vec4f,
+  color: vec4f,
+};
+
+struct LightUniforms {
+  ambientColor: vec4f,
+  lights: array<Light, 4>,
+};
+
+// Per-material uniforms. Add fields here only when a shader actually reads them;
+// preset-specific shaders (face.ts, future hair.ts) share this struct so the
+// engine can use one material bind-group layout.
+struct MaterialUniforms {
+  diffuseColor: vec3f,  // tint; multiplies sampled albedo (unused by current fs, reserved)
+  alpha: f32,            // 0 → discard; <1 → transparent draw call
+};
+
+struct VertexOutput {
+  @builtin(position) position: vec4f,
+  @location(0) normal: vec3f,
+  @location(1) uv: vec2f,
+  @location(2) worldPos: vec3f,
+};
+
+struct LightVP { viewProj: mat4x4f, };
+
+@group(0) @binding(0) var<uniform> camera: CameraUniforms;
+@group(0) @binding(1) var<uniform> light: LightUniforms;
+@group(0) @binding(2) var diffuseSampler: sampler;
+@group(0) @binding(3) var shadowMap: texture_depth_2d;
+@group(0) @binding(4) var shadowSampler: sampler_comparison;
+@group(0) @binding(5) var<uniform> lightVP: LightVP;
+@group(1) @binding(0) var<storage, read> skinMats: array<mat4x4f>;
+@group(2) @binding(0) var diffuseTexture: texture_2d<f32>;
+@group(2) @binding(1) var<uniform> material: MaterialUniforms;
+
+// ─── GGX specular helpers ───────────────────────────────────────────
+
+fn ggx_d(ndoth: f32, a2: f32) -> f32 {
+  let denom = ndoth * ndoth * (a2 - 1.0) + 1.0;
+  return a2 / (PI * denom * denom);
+}
+
+fn smith_g1(ndotx: f32, a2: f32) -> f32 {
+  return 2.0 * ndotx / (ndotx + sqrt(a2 + (1.0 - a2) * ndotx * ndotx));
+}
+
+fn fresnel_schlick(cosTheta: f32, f0: f32) -> f32 {
+  return f0 + (1.0 - f0) * pow(1.0 - cosTheta, 5.0);
+}
+
+// ─── Filmic tone mapping (LUT extracted from Blender 3.6 OCIO) ─────
+// View transform = Filmic, Look = Medium High Contrast, Exposure = -0.3.
+// 14 samples at integer log2 stops from -10 to +3 (inclusive).
+// Extracted via scripts/extract_filmic_lut.py → probe image through scene
+// color management. Input: linear scene-referred. Output: sRGB display.
+
+fn filmic(x: f32) -> f32 {
+  var lut = array<f32, 14>(
+    0.0067, 0.0141, 0.0272, 0.0499, 0.0885, 0.1512, 0.2462,
+    0.3753, 0.5273, 0.6776, 0.8031, 0.8929, 0.9495, 0.9814
+  );
+  let t = clamp(log2(max(x, 1e-10)) + 10.0, 0.0, 13.0);
+  let i = u32(t);
+  let j = min(i + 1u, 13u);
+  return mix(lut[i], lut[j], t - f32(i));
+}
+
+fn tonemap(hdr: vec3f) -> vec3f {
+  return vec3f(filmic(hdr.x), filmic(hdr.y), filmic(hdr.z));
+}
+
+// ─── Shadow sampling (3×3 PCF) ──────────────────────────────────────
+
+fn sampleShadow(worldPos: vec3f, n: vec3f) -> f32 {
+  let biasedPos = worldPos + n * 0.08;
+  let lclip = lightVP.viewProj * vec4f(biasedPos, 1.0);
+  let ndc = lclip.xyz / max(lclip.w, 1e-6);
+  let suv = vec2f(ndc.x * 0.5 + 0.5, 0.5 - ndc.y * 0.5);
+  let cmpZ = ndc.z - 0.001;
+  let ts = 1.0 / 4096.0;
+  var vis = 0.0;
+  for (var y = -1; y <= 1; y++) {
+    for (var x = -1; x <= 1; x++) {
+      vis += textureSampleCompare(shadowMap, shadowSampler, suv + vec2f(f32(x), f32(y)) * ts, cmpZ);
+    }
+  }
+  return vis / 9.0;
+}
+
+// ─── Vertex / Fragment ──────────────────────────────────────────────
+
+@vertex fn vs(
+  @location(0) position: vec3f,
+  @location(1) normal: vec3f,
+  @location(2) uv: vec2f,
+  @location(3) joints0: vec4<u32>,
+  @location(4) weights0: vec4<f32>
+) -> VertexOutput {
+  var output: VertexOutput;
+  let pos4 = vec4f(position, 1.0);
+  let weightSum = weights0.x + weights0.y + weights0.z + weights0.w;
+  let invWeightSum = select(1.0, 1.0 / weightSum, weightSum > 0.0001);
+  let nw = select(vec4f(1.0, 0.0, 0.0, 0.0), weights0 * invWeightSum, weightSum > 0.0001);
+  var skinnedPos = vec4f(0.0);
+  var skinnedNrm = vec3f(0.0);
+  for (var i = 0u; i < 4u; i++) {
+    let m = skinMats[joints0[i]];
+    let w = nw[i];
+    skinnedPos += (m * pos4) * w;
+    skinnedNrm += (mat3x3f(m[0].xyz, m[1].xyz, m[2].xyz) * normal) * w;
+  }
+  output.position = camera.projection * camera.view * vec4f(skinnedPos.xyz, 1.0);
+  output.normal = normalize(skinnedNrm);
+  output.uv = uv;
+  output.worldPos = skinnedPos.xyz;
+  return output;
+}
+
+@fragment fn fs(input: VertexOutput) -> @location(0) vec4f {
+  let alpha = material.alpha;
+  if (alpha < 0.001) { discard; }
+
+  let n = normalize(input.normal);
+  let v = normalize(camera.viewPos - input.worldPos);
+  let albedo = textureSample(diffuseTexture, diffuseSampler, input.uv).rgb;
+
+  let l = -light.lights[0].direction.xyz;
+  let sunColor = light.lights[0].color.xyz * light.lights[0].color.w;
+  let h = normalize(l + v);
+
+  let ndotl = max(dot(n, l), 0.0);
+  let ndotv = max(dot(n, v), 0.001);
+  let ndoth = max(dot(n, h), 0.0);
+  let vdoth = max(dot(v, h), 0.0);
+
+  let a2 = ROUGHNESS * ROUGHNESS;
+  let D = ggx_d(ndoth, a2);
+  let G = smith_g1(ndotl, a2) * smith_g1(ndotv, a2);
+  let F = fresnel_schlick(vdoth, F0_DIELECTRIC);
+  let spec = (D * G * F) / max(4.0 * ndotl * ndotv, 0.001);
+
+  let shadow = sampleShadow(input.worldPos, n);
+  let kd = (1.0 - F) * albedo / PI;
+  let direct = (kd + spec) * sunColor * ndotl * shadow;
+  let ambient = albedo * light.ambientColor.xyz;
+
+  return vec4f(ambient + direct, alpha);
+}
+
+`

package/src/shaders/dfg_lut.ts

@@ -0,0 +1,127 @@
+// One-shot bake pass that precomputes EEVEE's BRDF split-sum DFG LUT.
+// Direct WGSL port of Blender 3.6 source/blender/draw/engines/eevee/shaders/
+//   bsdf_lut_frag.glsl + bsdf_sampling_lib.glsl (VNDF GGX branch).
+//
+// Output texture: 64×64 rg16float, written once at engine init.
+//   R = (1 - Fc) × BRDF integrated over GGX VNDF — scales f0
+//   G = Fc × BRDF integrated over GGX VNDF — scales f90
+// Runtime sampler: see brdf_lut_baked() in nodes.ts. Plug into
+// F_brdf_single_scatter / F_brdf_multi_scatter verbatim.
+
+export const DFG_LUT_SIZE = 64
+export const DFG_LUT_SAMPLE_COUNT = 32
+
+export const DFG_LUT_WGSL = /* wgsl */ `
+const LUT_SIZE: f32 = ${DFG_LUT_SIZE}.0;
+const SAMPLE_COUNT: u32 = ${DFG_LUT_SAMPLE_COUNT}u;
+const M_2PI: f32 = 6.283185307179586;
+
+@vertex fn vs(@builtin(vertex_index) vid: u32) -> @builtin(position) vec4f {
+  // Full-screen triangle covering [-1,1]² in NDC.
+  let x = f32((vid << 1u) & 2u) * 2.0 - 1.0;
+  let y = f32(vid & 2u) * 2.0 - 1.0;
+  return vec4f(x, y, 0.0, 1.0);
+}
+
+// common_math_geom_lib.glsl:165 — make_orthonormal_basis.
+fn orthonormal_basis(N: vec3f) -> mat2x3f {
+  let up = select(vec3f(1.0, 0.0, 0.0), vec3f(0.0, 0.0, 1.0), abs(N.z) < 0.99999);
+  let T = normalize(cross(up, N));
+  let B = cross(N, T);
+  return mat2x3f(T, B);
+}
+
+// bsdf_sampling_lib.glsl:27 — Heitz 2018 VNDF sampling in tangent space.
+fn sample_ggx_vndf(rand: vec3f, alpha: f32, Vt: vec3f) -> vec3f {
+  let Vh = normalize(vec3f(alpha * Vt.xy, Vt.z));
+  let tb = orthonormal_basis(Vh);
+  let Th = tb[0];
+  let Bh = tb[1];
+  let r = sqrt(rand.x);
+  let x = r * rand.y;
+  var y = r * rand.z;
+  let s = 0.5 * (1.0 + Vh.z);
+  y = (1.0 - s) * sqrt(1.0 - x * x) + s * y;
+  let z = sqrt(saturate(1.0 - x * x - y * y));
+  let Hh = x * Th + y * Bh + z * Vh;
+  return normalize(vec3f(alpha * Hh.xy, saturate(Hh.z)));
+}
+
+// bsdf_common_lib.glsl:105 — G1 Smith GGX (Brian Karis opti form).
+fn G1_Smith_GGX_opti(NX: f32, a2: f32) -> f32 {
+  return NX + sqrt(NX * (NX - NX * a2) + a2);
+}
+
+// bsdf_common_lib.glsl:50 — exact dielectric Fresnel (monochromatic).
+fn F_eta(eta: f32, cos_theta: f32) -> f32 {
+  let c = abs(cos_theta);
+  var g = eta * eta - 1.0 + c * c;
+  if (g > 0.0) {
+    g = sqrt(g);
+    let A = (g - c) / (g + c);
+    let B = (c * (g + c) - 1.0) / (c * (g - c) + 1.0);
+    return 0.5 * A * A * (1.0 + B * B);
+  }
+  return 1.0; // total internal reflection
+}
+
+fn f0_from_ior(eta: f32) -> f32 {
+  let A = (eta - 1.0) / (eta + 1.0);
+  return A * A;
+}
+
+// F_color_blend(eta, fresnel, vec3(0)).r — blend factor only.
+fn F_color_blend_zero(eta: f32, fresnel: f32) -> f32 {
+  let f0 = f0_from_ior(eta);
+  return saturate((fresnel - f0) / (1.0 - f0));
+}
+
+@fragment fn fs(@builtin(position) frag: vec4f) -> @location(0) vec2f {
+  let y_uv = floor(frag.y) / (LUT_SIZE - 1.0);
+  let x_uv = floor(frag.x) / (LUT_SIZE - 1.0);
+
+  let NV = clamp(1.0 - y_uv * y_uv, 1e-4, 0.9999);
+  let a = x_uv * x_uv;
+  let a2 = clamp(a * a, 1e-4, 0.9999);
+
+  let V = vec3f(sqrt(1.0 - NV * NV), 0.0, NV);
+
+  // principled specular=1.0 — max value, matches bsdf_lut_frag.glsl:41.
+  let eta = (2.0 / (1.0 - sqrt(0.08 * 1.0))) - 1.0;
+
+  var brdf_accum = 0.0;
+  var fresnel_accum = 0.0;
+  let sc_f = f32(SAMPLE_COUNT);
+  for (var j: u32 = 0u; j < SAMPLE_COUNT; j = j + 1u) {
+    for (var i: u32 = 0u; i < SAMPLE_COUNT; i = i + 1u) {
+      let ix = (f32(i) + 0.5) / sc_f;
+      let iy = (f32(j) + 0.5) / sc_f;
+      // Xi.x = radial, Xi.yz = (cos, sin) of azimuth — bsdf_lut_frag.glsl:22.
+      let Xi = vec3f(ix, cos(iy * M_2PI), sin(iy * M_2PI));
+
+      let H = sample_ggx_vndf(Xi, a, V);
+      let L = -reflect(V, H);
+      let NL = L.z;
+      if (NL > 0.0) {
+        let NH = max(H.z, 0.0);
+        let VH = max(dot(V, H), 0.0);
+
+        // G_smith (divided form): 4·NV·NL / (G1_v·G1_l). See bsdf_common_lib.glsl:105.
+        let G1v = G1_Smith_GGX_opti(NV, a2);
+        let G1l = G1_Smith_GGX_opti(NL, a2);
+        let G_smith = 4.0 * NV * NL / (G1v * G1l);
+
+        let brdf = (G_smith * VH) / (NH * NV);
+
+        let fresnel = F_eta(eta, VH);
+        let Fc = F_color_blend_zero(eta, fresnel);
+
+        brdf_accum = brdf_accum + (1.0 - Fc) * brdf;
+        fresnel_accum = fresnel_accum + Fc * brdf;
+      }
+    }
+  }
+  let n2 = sc_f * sc_f;
+  return vec2f(brdf_accum / n2, fresnel_accum / n2);
+}
+`

package/src/shaders/eye.ts

@@ -0,0 +1,143 @@
+// Eye preset — default Principled BSDF (F0=0.04, Roughness=0.5) + Emission socket set to albedo × 1.5.
+// Matches the published preset's instruction: "keep eyes in the default nodegraph, add emission 1.5".
+// Blender's Principled BSDF Emission socket is added on top of the shaded output (pre-tonemap, feeds bloom).
+
+export const EYE_SHADER_WGSL = /* wgsl */ `
+
+const PI: f32 = 3.141592653589793;
+const F0_DIELECTRIC: f32 = 0.04;
+const ROUGHNESS: f32 = 0.5;
+const EYE_EMISSION_STRENGTH: f32 = 1.5;
+
+struct CameraUniforms {
+  view: mat4x4f,
+  projection: mat4x4f,
+  viewPos: vec3f,
+  _padding: f32,
+};
+
+struct Light {
+  direction: vec4f,
+  color: vec4f,
+};
+
+struct LightUniforms {
+  ambientColor: vec4f,
+  lights: array<Light, 4>,
+};
+
+struct MaterialUniforms {
+  diffuseColor: vec3f,
+  alpha: f32,
+};
+
+struct VertexOutput {
+  @builtin(position) position: vec4f,
+  @location(0) normal: vec3f,
+  @location(1) uv: vec2f,
+  @location(2) worldPos: vec3f,
+};
+
+struct LightVP { viewProj: mat4x4f, };
+
+@group(0) @binding(0) var<uniform> camera: CameraUniforms;
+@group(0) @binding(1) var<uniform> light: LightUniforms;
+@group(0) @binding(2) var diffuseSampler: sampler;
+@group(0) @binding(3) var shadowMap: texture_depth_2d;
+@group(0) @binding(4) var shadowSampler: sampler_comparison;
+@group(0) @binding(5) var<uniform> lightVP: LightVP;
+@group(1) @binding(0) var<storage, read> skinMats: array<mat4x4f>;
+@group(2) @binding(0) var diffuseTexture: texture_2d<f32>;
+@group(2) @binding(1) var<uniform> material: MaterialUniforms;
+
+fn ggx_d(ndoth: f32, a2: f32) -> f32 {
+  let denom = ndoth * ndoth * (a2 - 1.0) + 1.0;
+  return a2 / (PI * denom * denom);
+}
+
+fn smith_g1(ndotx: f32, a2: f32) -> f32 {
+  return 2.0 * ndotx / (ndotx + sqrt(a2 + (1.0 - a2) * ndotx * ndotx));
+}
+
+fn fresnel_schlick(cosTheta: f32, f0: f32) -> f32 {
+  return f0 + (1.0 - f0) * pow(1.0 - cosTheta, 5.0);
+}
+
+fn sampleShadow(worldPos: vec3f, n: vec3f) -> f32 {
+  let biasedPos = worldPos + n * 0.08;
+  let lclip = lightVP.viewProj * vec4f(biasedPos, 1.0);
+  let ndc = lclip.xyz / max(lclip.w, 1e-6);
+  let suv = vec2f(ndc.x * 0.5 + 0.5, 0.5 - ndc.y * 0.5);
+  let cmpZ = ndc.z - 0.001;
+  let ts = 1.0 / 4096.0;
+  var vis = 0.0;
+  for (var y = -1; y <= 1; y++) {
+    for (var x = -1; x <= 1; x++) {
+      vis += textureSampleCompare(shadowMap, shadowSampler, suv + vec2f(f32(x), f32(y)) * ts, cmpZ);
+    }
+  }
+  return vis / 9.0;
+}
+
+@vertex fn vs(
+  @location(0) position: vec3f,
+  @location(1) normal: vec3f,
+  @location(2) uv: vec2f,
+  @location(3) joints0: vec4<u32>,
+  @location(4) weights0: vec4<f32>
+) -> VertexOutput {
+  var output: VertexOutput;
+  let pos4 = vec4f(position, 1.0);
+  let weightSum = weights0.x + weights0.y + weights0.z + weights0.w;
+  let invWeightSum = select(1.0, 1.0 / weightSum, weightSum > 0.0001);
+  let nw = select(vec4f(1.0, 0.0, 0.0, 0.0), weights0 * invWeightSum, weightSum > 0.0001);
+  var skinnedPos = vec4f(0.0);
+  var skinnedNrm = vec3f(0.0);
+  for (var i = 0u; i < 4u; i++) {
+    let m = skinMats[joints0[i]];
+    let w = nw[i];
+    skinnedPos += (m * pos4) * w;
+    skinnedNrm += (mat3x3f(m[0].xyz, m[1].xyz, m[2].xyz) * normal) * w;
+  }
+  output.position = camera.projection * camera.view * vec4f(skinnedPos.xyz, 1.0);
+  output.normal = normalize(skinnedNrm);
+  output.uv = uv;
+  output.worldPos = skinnedPos.xyz;
+  return output;
+}
+
+@fragment fn fs(input: VertexOutput) -> @location(0) vec4f {
+  let alpha = material.alpha;
+  if (alpha < 0.001) { discard; }
+
+  let n = normalize(input.normal);
+  let v = normalize(camera.viewPos - input.worldPos);
+  let albedo = textureSample(diffuseTexture, diffuseSampler, input.uv).rgb;
+
+  let l = -light.lights[0].direction.xyz;
+  let sunColor = light.lights[0].color.xyz * light.lights[0].color.w;
+  let h = normalize(l + v);
+
+  let ndotl = max(dot(n, l), 0.0);
+  let ndotv = max(dot(n, v), 0.001);
+  let ndoth = max(dot(n, h), 0.0);
+  let vdoth = max(dot(v, h), 0.0);
+
+  let a2 = ROUGHNESS * ROUGHNESS;
+  let D = ggx_d(ndoth, a2);
+  let G = smith_g1(ndotl, a2) * smith_g1(ndotv, a2);
+  let F = fresnel_schlick(vdoth, F0_DIELECTRIC);
+  let spec = (D * G * F) / max(4.0 * ndotl * ndotv, 0.001);
+
+  let shadow = sampleShadow(input.worldPos, n);
+  let kd = (1.0 - F) * albedo / PI;
+  let direct = (kd + spec) * sunColor * ndotl * shadow;
+  let ambient = albedo * light.ambientColor.xyz;
+
+  // Principled Emission socket: emissive = emission_color × strength, added on top of shading.
+  let emission = albedo * EYE_EMISSION_STRENGTH;
+
+  return vec4f(ambient + direct + emission, alpha);
+}
+
+`

package/src/shaders/face.ts

@@ -0,0 +1,213 @@
+// M_Face — WGSL trace of 仿深空之眼渲染预设v1.0_by_小绿毛猫_material_graph_dump.json "M_Face"; VALTORGB stops from m_graphs (dump omits curve keys).
+
+import { NODES_WGSL } from "./nodes"
+
+export const FACE_SHADER_WGSL = /* wgsl */ `
+
+${NODES_WGSL}
+
+struct CameraUniforms {
+  view: mat4x4f,
+  projection: mat4x4f,
+  viewPos: vec3f,
+  _padding: f32,
+};
+
+struct Light {
+  direction: vec4f,
+  color: vec4f,
+};
+
+struct LightUniforms {
+  ambientColor: vec4f,
+  lights: array<Light, 4>,
+};
+
+struct MaterialUniforms {
+  diffuseColor: vec3f,
+  alpha: f32,
+};
+
+struct VertexOutput {
+  @builtin(position) position: vec4f,
+  @location(0) normal: vec3f,
+  @location(1) uv: vec2f,
+  @location(2) worldPos: vec3f,
+};
+
+struct LightVP { viewProj: mat4x4f, };
+
+@group(0) @binding(0) var<uniform> camera: CameraUniforms;
+@group(0) @binding(1) var<uniform> light: LightUniforms;
+@group(0) @binding(2) var diffuseSampler: sampler;
+@group(0) @binding(3) var shadowMap: texture_depth_2d;
+@group(0) @binding(4) var shadowSampler: sampler_comparison;
+@group(0) @binding(5) var<uniform> lightVP: LightVP;
+@group(1) @binding(0) var<storage, read> skinMats: array<mat4x4f>;
+@group(2) @binding(0) var diffuseTexture: texture_2d<f32>;
+@group(2) @binding(1) var<uniform> material: MaterialUniforms;
+
+// 3x3 PCF shadow sampling, 4096 map, normal-bias 0.08, depth-bias 0.001
+fn sampleShadow(worldPos: vec3f, n: vec3f) -> f32 {
+  let biasedPos = worldPos + n * 0.08;
+  let lclip = lightVP.viewProj * vec4f(biasedPos, 1.0);
+  let ndc = lclip.xyz / max(lclip.w, 1e-6);
+  let suv = vec2f(ndc.x * 0.5 + 0.5, 0.5 - ndc.y * 0.5);
+  let cmpZ = ndc.z - 0.001;
+  let ts = 1.0 / 4096.0;
+  var vis = 0.0;
+  for (var y = -1; y <= 1; y++) {
+    for (var x = -1; x <= 1; x++) {
+      vis += textureSampleCompare(shadowMap, shadowSampler, suv + vec2f(f32(x), f32(y)) * ts, cmpZ);
+    }
+  }
+  return vis / 9.0;
+}
+
+const PI_F: f32 = 3.141592653589793;
+const FACE_SPECULAR: f32 = 0.5;
+const FACE_ROUGHNESS: f32 = 0.3;
+// Dump M_Face unlinked defaults (math op enum not serialized — warm clamp chain still from m_graphs)
+const FACE_RIM2_POW: f32 = 0.6300000548362732;
+const FACE_RIM2_BG: vec3f = vec3f(1.0, 0.4684903025627136, 0.3698573112487793);
+const FACE_WARM_AO_MUL: f32 = 0.30000001192092896; // 运算.004 MULTIPLY after invert (was 0.5 in older trace)
+const FACE_BRIGHT_TEX_THRESH: f32 = 0.9300000071525574; // 运算.005 GREATER_THAN Value_001
+const FACE_MIX_NPR: f32 = 0.5; // 混合着色器.001 Fac
+
+@vertex fn vs(
+  @location(0) position: vec3f,
+  @location(1) normal: vec3f,
+  @location(2) uv: vec2f,
+  @location(3) joints0: vec4<u32>,
+  @location(4) weights0: vec4<f32>
+) -> VertexOutput {
+  var output: VertexOutput;
+  let pos4 = vec4f(position, 1.0);
+  let weightSum = weights0.x + weights0.y + weights0.z + weights0.w;
+  let invWeightSum = select(1.0, 1.0 / weightSum, weightSum > 0.0001);
+  let nw = select(vec4f(1.0, 0.0, 0.0, 0.0), weights0 * invWeightSum, weightSum > 0.0001);
+  var skinnedPos = vec4f(0.0);
+  var skinnedNrm = vec3f(0.0);
+  for (var i = 0u; i < 4u; i++) {
+    let m = skinMats[joints0[i]];
+    let w = nw[i];
+    skinnedPos += (m * pos4) * w;
+    skinnedNrm += (mat3x3f(m[0].xyz, m[1].xyz, m[2].xyz) * normal) * w;
+  }
+  output.position = camera.projection * camera.view * vec4f(skinnedPos.xyz, 1.0);
+  output.normal = normalize(skinnedNrm);
+  output.uv = uv;
+  output.worldPos = skinnedPos.xyz;
+  return output;
+}
+
+// Fragment: M_Face NPR + Principled hybrid
+// TEX → HueSat shadow/lit → toon gate → BrightContrast → AO chain → emission stack
+// Fresnel rims, warm AO emission, bright-texture gate, noise-bumped Principled
+// Final = mix(Principled, NPR, 0.5)
+@fragment fn fs(input: VertexOutput) -> @location(0) vec4f {
+  let alpha = material.alpha;
+  if (alpha < 0.001) { discard; }
+
+  let n = normalize(input.normal);
+  let v = normalize(camera.viewPos - input.worldPos);
+  let l = -light.lights[0].direction.xyz;
+  let intensity = light.lights[0].color.w;
+  let sun = light.lights[0].color.xyz * intensity;
+
+  let tex_color = textureSample(diffuseTexture, diffuseSampler, input.uv).rgb;
+  let shadow = sampleShadow(input.worldPos, n);
+
+  // ═══ SOURCES ═══
+  // DiffuseBSDF(white) → ShaderToRGB (energy-matched); shadow on direct only
+  let ndotl_raw = shader_to_rgb_diffuse(n, l, sun, light.ambientColor.xyz, shadow);
+  // ramp.008 CONSTANT — edge AA avoids binary fac shimmer / white specks on terminator (fwidth + smoothstep)
+  let toon = ramp_constant_edge_aa(ndotl_raw, 0.2966, vec4f(0,0,0,1), vec4f(1,1,1,1)).r;
+  let ao = ao_fake(n, v);
+
+  // ═══ TOON COLOR ═══
+  let shadow_tint = hue_sat(0.46000000834465027, 2.0, 0.3499999940395355, 1.0, tex_color); // HueSat.002
+  let lit_tint = hue_sat(0.46000000834465027, 1.600000023841858, 1.5, 1.0, tex_color); // HueSat.001
+  let toon_color = mix_blend(toon, shadow_tint, lit_tint); // Mix.004
+  let bc = bright_contrast(toon_color, 0.1, 0.2);
+
+  // ═══ AO CHAIN ═══
+  // ramp CONSTANT [0→white, 0.5995→black]
+  let ao_ramp = ramp_constant(ao, 0.0, vec4f(1,1,1,1), 0.5995, vec4f(0,0,0,1)).r;
+  // Mix.003(Factor=ao_ramp, A=bc, B=reddish tint)
+  let ao_mixed = mix_blend(ao_ramp, bc, vec3f(0.8302, 0.3346, 0.2795));
+
+  // ═══ EMISSION 3 ═══
+  let emission3 = ao_mixed * 2.5; // Emission.003(Strength=2.5)
+
+  // ═══ WARM EMISSION ═══
+  let ao_inv = invert_f(1.0, ao_ramp);
+  let warm_str = ao_inv * FACE_WARM_AO_MUL; // 反转 → 运算.004 MULTIPLY Value_001
+  let warm_input = clamp(toon * 0.5 + 0.5, 0.0, 1.0); // 运算.001→运算.006→Clamp
+  // ramp.003 CARDINAL [0.2409→warm dark, 0.4663→warm light]
+  let warm_color = ramp_cardinal(warm_input, 0.2409,
+    vec4f(0.2426, 0.068, 0.0588, 1.0), 0.4663,
+    vec4f(0.6677, 0.5024, 0.5126, 1.0)).rgb;
+  let warm_emission = warm_color * warm_str; // Emission.001
+
+  // ═══ RIM 1 ═══
+  // Fresnel(IOR=2.0) × LayerWeight.001(Facing, Blend=0.24)
+  let rim1_str = fresnel(2.0, n, v) * layer_weight_facing(0.24, n, v);
+  let rim1 = vec3f(0.984157919883728, 0.6110184788703918, 0.5736401677131653) * rim1_str;
+
+  // ═══ RIM 2 ═══
+  // Fresnel.001(IOR=1.45) × LayerWeight.002(Fresnel output, Blend=0.61)
+  let rim2_raw = fresnel(1.45, n, v) * layer_weight_fresnel(0.61, n, v);
+  let rim2_fac = math_power(rim2_raw, FACE_RIM2_POW);
+  // MixShader.002: Shader=Emission.003, Shader_001=背景
+  let rim2_mixed = mix(emission3, FACE_RIM2_BG, rim2_fac);
+
+  // 转接点.005(tex) → 运算.005 GREATER_THAN Value_001
+  let tex_gate = math_greater_than(tex_color.r, FACE_BRIGHT_TEX_THRESH);
+  let bright_emit = vec3f(tex_gate) * 3.0; // Emission.002(Strength=3.0)
+
+  // ═══ NPR STACK (AddShader chain) ═══
+  let add2 = rim2_mixed + bright_emit; // AddShader.002
+  let add0 = rim1 + add2; // AddShader
+  let npr_stack = add0 + warm_emission; // AddShader.001
+
+  // ═══ PRINCIPLED BSDF ═══
+  // Noise-based bump normal
+  let gen = mapping_point(input.worldPos, vec3f(0.0), vec3f(0.0), vec3f(1.0, 1.0, 1.5));
+  let noise_val = tex_noise(gen, 1.0, 2.0, 0.5, 0.0);
+  let noise_ramp = ramp_linear(noise_val, 0.0, vec4f(0,0,0,1), 1.0, vec4f(1,1,1,1)).r;
+  let bumped_n = bump_lh(0.324644535779953, noise_ramp, n, input.worldPos); // 凹凸 Strength; LH bump
+
+  // Mix.001(Factor=noise_ramp, A=bc, B=dark red)
+  let principled_base = mix_blend(noise_ramp, bc, vec3f(0.6832, 0.1947, 0.1373));
+  // Emission input from reroute.011 (bc), Strength=0.2
+  let p_emission = bc * 0.2;
+  // AO.002 → ramp.005 LINEAR [0.003→black, 1.0→gray] for subsurface approx
+  let ao2 = ao_fake(n, v);
+  let sss = ramp_linear(ao2, 0.003, vec4f(0,0,0,1), 1.0, vec4f(0.0786, 0.0786, 0.0786, 1.0)).r;
+
+  // 原理化BSDF (EEVEE port): metallic=0, specular=0.5, roughness=0.3, specular_tint=0.
+  let NL = max(dot(bumped_n, l), 0.0);
+  let NV = max(dot(bumped_n, v), 1e-4);
+
+  let f0 = vec3f(0.08 * FACE_SPECULAR);
+  let f90 = mix(f0, vec3f(1.0), sqrt(FACE_SPECULAR));
+  let split_sum = brdf_lut_baked(NV, FACE_ROUGHNESS);
+  let reflection_color = F_brdf_multi_scatter(f0, f90, split_sum);
+
+  let spec_direct = bsdf_ggx(bumped_n, l, v, FACE_ROUGHNESS) * sun * shadow * ltc_brdf_scale(NV, FACE_ROUGHNESS);
+  let spec_indirect = light.ambientColor.xyz;
+  let spec_radiance = (spec_direct + spec_indirect) * reflection_color;
+
+  // Indirect diffuse = base_color × L_w per Blender closure_eval_surface_lib.glsl line 302;
+  // probe_evaluate_world_diff returns radiance (SH-projected, not cosine-convolved).
+  let diffuse_radiance = principled_base * (sun * NL * shadow / PI_F + light.ambientColor.xyz);
+  let principled = diffuse_radiance + spec_radiance + p_emission + vec3f(sss);
+
+  // 混合着色器.001: Shader=相加着色器.001, Shader_001=原理化BSDF — Fac blends toward second
+  let final_color = mix(npr_stack, principled, FACE_MIX_NPR);
+
+  return vec4f(final_color, alpha);
+}
+
+`