@pibbl/core 0.0.1 → 0.0.3-alpha.109.1

package/dist/features/textures/lib/chemical-field.d.ts

@@ -0,0 +1,56 @@
+/** Private periodic finite-difference state; presentation adds a duplicate seam row/column. */
+export interface ChemicalConfig {
+    readonly a: number;
+    readonly b: number;
+    readonly epsilon: number;
+    readonly diffusionU: number;
+    readonly diffusionV: number;
+    readonly seed: number;
+    readonly initialNoise: number;
+    readonly warmupSteps: number;
+    readonly speed: number;
+    readonly paused: boolean;
+}
+interface Model {
+    readonly equilibrium: (c: ChemicalConfig) => readonly [number, number];
+    readonly reaction: (u: number, v: number, c: ChemicalConfig) => readonly [number, number];
+    readonly excite: boolean;
+}
+export declare const schnakenbergModel: Model;
+export declare const brusselatorModel: Model;
+export declare const barkleyModel: Model;
+export declare class ChemicalField {
+    private readonly options;
+    u: Float32Array;
+    v: Float32Array;
+    private nextU;
+    private nextV;
+    private maxU;
+    private maxV;
+    readonly width: number;
+    readonly height: number;
+    readonly baseline: number;
+    constructor(options: {
+        width: number;
+        height: number;
+        config: ChemicalConfig;
+        model: Model;
+    });
+    advance(count: number): void;
+    paint(options: {
+        u: number;
+        v: number;
+        radius: number;
+        strength: number;
+    }): void;
+    stir(options: {
+        u: number;
+        v: number;
+        du: number;
+        dv: number;
+        radius: number;
+    }): void;
+    intensity(index: number): number;
+    dispose(): void;
+}
+export {};

package/dist/features/textures/lib/chemical-field.js

@@ -0,0 +1,154 @@
+const schnakenbergModel = {
+  equilibrium: (c) => [c.a + c.b, c.b / (c.a + c.b) ** 2],
+  reaction: (u, v, c) => {
+    const q = u * u * v;
+    return [c.a - u + q, c.b - q];
+  },
+  excite: false
+};
+const brusselatorModel = {
+  equilibrium: (c) => [c.a, c.b / c.a],
+  reaction: (u, v, c) => {
+    const q = u * u * v;
+    return [c.a - (c.b + 1) * u + q, c.b * u - q];
+  },
+  excite: false
+};
+const barkleyModel = {
+  equilibrium: () => [0, 0],
+  reaction: (u, v, c) => [
+    u * (1 - u) * (u - (v + c.b) / c.a) / c.epsilon,
+    u - v
+  ],
+  excite: true
+};
+function hash(i, seed) {
+  let n = Math.imul(i ^ seed, 73244475);
+  n = Math.imul(n ^ n >>> 16, 73244475);
+  return ((n ^ n >>> 16) >>> 0) / 4294967295;
+}
+const wrap = (value, size) => (value % size + size) % size;
+class ChemicalField {
+  constructor(options) {
+    this.options = options;
+    this.width = options.width;
+    this.height = options.height;
+    const n = this.width * this.height;
+    this.u = new Float32Array(n);
+    this.v = new Float32Array(n);
+    this.nextU = new Float32Array(n);
+    this.nextV = new Float32Array(n);
+    const c = options.config, [u, v] = options.model.equilibrium(c);
+    this.baseline = options.model.excite ? 1 : u;
+    for (let y = 0; y < this.height; y++)
+      for (let x = 0; x < this.width; x++) {
+        const i = y * this.width + x, noise = hash(i, c.seed) - 0.5;
+        if (options.model.excite) {
+          const split = 0.5 + (hash(0, c.seed) - 0.5) * 0.08;
+          this.u[i] = (x > this.width * 0.22 && x < this.width * 0.5 && y < this.height * split ? 1 : 0) * (1 - c.initialNoise * 0.5) + c.initialNoise * Math.max(0, noise);
+          this.v[i] = y >= this.height * split ? 0.5 : 0;
+        } else {
+          this.u[i] = u * (1 + c.initialNoise * noise);
+          this.v[i] = v * (1 + c.initialNoise * (hash(i, c.seed ^ 1367130551) - 0.5));
+        }
+      }
+    this.maxU = Math.max(1, u * 1.5);
+    this.maxV = Math.max(1, v * 1.5);
+    this.advance(c.warmupSteps);
+  }
+  options;
+  u;
+  v;
+  nextU;
+  nextV;
+  maxU = 1;
+  maxV = 1;
+  width;
+  height;
+  baseline;
+  advance(count) {
+    const { config: c, model } = this.options, w = this.width, h = this.height;
+    for (let step = 0; step < count; step++) {
+      const stiffness = model.excite ? 1 / c.epsilon : 1 + this.maxU * this.maxU + 2 * this.maxU * this.maxV + c.b;
+      const dt = Math.min(
+        0.05,
+        0.2 / Math.max(1e-3, c.diffusionU, c.diffusionV),
+        0.1 / stiffness
+      );
+      let maxU = 0, maxV = 0;
+      for (let y = 0; y < h; y++) {
+        const row = y * w, above = (y === 0 ? h - 1 : y - 1) * w, below = (y === h - 1 ? 0 : y + 1) * w;
+        for (let x = 0; x < w; x++) {
+          const i = row + x, left = row + (x === 0 ? w - 1 : x - 1), right = row + (x === w - 1 ? 0 : x + 1);
+          const u = this.u[i], v = this.v[i];
+          const lapU = this.u[left] + this.u[right] + this.u[above + x] + this.u[below + x] - 4 * u;
+          const lapV = this.v[left] + this.v[right] + this.v[above + x] + this.v[below + x] - 4 * v;
+          const [ru, rv] = model.reaction(u, v, c);
+          const nu = Math.max(0, u + dt * (c.diffusionU * lapU + ru));
+          const nv = Math.max(0, v + dt * (c.diffusionV * lapV + rv));
+          if (!Number.isFinite(nu) || !Number.isFinite(nv))
+            throw new RangeError(
+              "Reaction-diffusion parameters produced a nonfinite field."
+            );
+          this.nextU[i] = nu;
+          this.nextV[i] = nv;
+          maxU = Math.max(maxU, nu);
+          maxV = Math.max(maxV, nv);
+        }
+      }
+      [this.u, this.nextU] = [this.nextU, this.u];
+      [this.v, this.nextV] = [this.nextV, this.v];
+      this.maxU = maxU;
+      this.maxV = maxV;
+    }
+  }
+  paint(options) {
+    const { width: w, height: h } = this;
+    for (let y = 0; y < h; y++)
+      for (let x = 0; x < w; x++) {
+        const dx = wrap(x - options.u * w + w / 2, w) - w / 2, dy = wrap(y - options.v * h + h / 2, h) - h / 2;
+        const distance = Math.hypot(dx, dy) / options.radius;
+        if (distance >= 1) continue;
+        const dose = options.strength * (1 - distance * distance) ** 2;
+        const i = y * w + x;
+        this.u[i] = Math.max(
+          this.u[i],
+          this.options.model.excite ? dose : this.baseline * (1 + 2 * dose)
+        );
+      }
+    this.maxU = Math.max(this.maxU, this.baseline * 3);
+  }
+  stir(options) {
+    const w = this.width, h = this.height;
+    const sample = (field, x, y) => {
+      x = wrap(x, w);
+      y = wrap(y, h);
+      const ix = Math.floor(x), iy = Math.floor(y), fx = x - ix, fy = y - iy;
+      const x1 = (ix + 1) % w, y1 = (iy + 1) % h;
+      return (field[iy * w + ix] * (1 - fx) + field[iy * w + x1] * fx) * (1 - fy) + (field[y1 * w + ix] * (1 - fx) + field[y1 * w + x1] * fx) * fy;
+    };
+    for (let y = 0; y < h; y++)
+      for (let x = 0; x < w; x++) {
+        const dx = wrap(x - options.u * w + w / 2, w) - w / 2, dy = wrap(y - options.v * h + h / 2, h) - h / 2;
+        const t = Math.max(0, 1 - Math.hypot(dx, dy) / options.radius), weight = t * t * (3 - 2 * t);
+        const px = x - options.du * w * weight, py = y - options.dv * h * weight, i = y * w + x;
+        this.nextU[i] = sample(this.u, px, py);
+        this.nextV[i] = sample(this.v, px, py);
+      }
+    [this.u, this.nextU] = [this.nextU, this.u];
+    [this.v, this.nextV] = [this.nextV, this.v];
+  }
+  intensity(index) {
+    return this.options.model.excite ? this.u[index] : this.u[index] / (2 * this.baseline);
+  }
+  dispose() {
+    this.u = this.v = this.nextU = this.nextV = new Float32Array(0);
+  }
+}
+export {
+  ChemicalField,
+  barkleyModel,
+  brusselatorModel,
+  schnakenbergModel
+};
+//# sourceMappingURL=chemical-field.js.map

package/dist/features/textures/lib/chemical-field.js.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["../../../../src/features/textures/lib/chemical-field.ts"],
+  "sourcesContent": ["/** Private periodic finite-difference state; presentation adds a duplicate seam row/column. */\nexport interface ChemicalConfig {\n  readonly a: number;\n  readonly b: number;\n  readonly epsilon: number;\n  readonly diffusionU: number;\n  readonly diffusionV: number;\n  readonly seed: number;\n  readonly initialNoise: number;\n  readonly warmupSteps: number;\n  readonly speed: number;\n  readonly paused: boolean;\n}\ninterface Model {\n  readonly equilibrium: (c: ChemicalConfig) => readonly [number, number];\n  readonly reaction: (\n    u: number,\n    v: number,\n    c: ChemicalConfig,\n  ) => readonly [number, number];\n  readonly excite: boolean;\n}\nexport const schnakenbergModel: Model = {\n  equilibrium: (c) => [c.a + c.b, c.b / (c.a + c.b) ** 2],\n  reaction: (u, v, c) => {\n    const q = u * u * v;\n    return [c.a - u + q, c.b - q];\n  },\n  excite: false,\n};\nexport const brusselatorModel: Model = {\n  equilibrium: (c) => [c.a, c.b / c.a],\n  reaction: (u, v, c) => {\n    const q = u * u * v;\n    return [c.a - (c.b + 1) * u + q, c.b * u - q];\n  },\n  excite: false,\n};\nexport const barkleyModel: Model = {\n  equilibrium: () => [0, 0],\n  reaction: (u, v, c) => [\n    (u * (1 - u) * (u - (v + c.b) / c.a)) / c.epsilon,\n    u - v,\n  ],\n  excite: true,\n};\nfunction hash(i: number, seed: number) {\n  let n = Math.imul(i ^ seed, 0x45d9f3b);\n  n = Math.imul(n ^ (n >>> 16), 0x45d9f3b);\n  return ((n ^ (n >>> 16)) >>> 0) / 0xffffffff;\n}\nconst wrap = (value: number, size: number) => ((value % size) + size) % size;\nexport class ChemicalField {\n  u: Float32Array;\n  v: Float32Array;\n  private nextU: Float32Array;\n  private nextV: Float32Array;\n  private maxU = 1;\n  private maxV = 1;\n  readonly width: number;\n  readonly height: number;\n  readonly baseline: number;\n  constructor(\n    private readonly options: {\n      width: number;\n      height: number;\n      config: ChemicalConfig;\n      model: Model;\n    },\n  ) {\n    this.width = options.width;\n    this.height = options.height;\n    const n = this.width * this.height;\n    this.u = new Float32Array(n);\n    this.v = new Float32Array(n);\n    this.nextU = new Float32Array(n);\n    this.nextV = new Float32Array(n);\n    const c = options.config,\n      [u, v] = options.model.equilibrium(c);\n    this.baseline = options.model.excite ? 1 : u;\n    for (let y = 0; y < this.height; y++)\n      for (let x = 0; x < this.width; x++) {\n        const i = y * this.width + x,\n          noise = hash(i, c.seed) - 0.5;\n        if (options.model.excite) {\n          const split = 0.5 + (hash(0, c.seed) - 0.5) * 0.08;\n          this.u[i] =\n            (x > this.width * 0.22 &&\n            x < this.width * 0.5 &&\n            y < this.height * split\n              ? 1\n              : 0) *\n              (1 - c.initialNoise * 0.5) +\n            c.initialNoise * Math.max(0, noise);\n          this.v[i] = y >= this.height * split ? 0.5 : 0;\n        } else {\n          this.u[i] = u * (1 + c.initialNoise * noise);\n          this.v[i] =\n            v * (1 + c.initialNoise * (hash(i, c.seed ^ 0x517cc1b7) - 0.5));\n        }\n      }\n    this.maxU = Math.max(1, u * 1.5);\n    this.maxV = Math.max(1, v * 1.5);\n    this.advance(c.warmupSteps);\n  }\n  advance(count: number) {\n    const { config: c, model } = this.options,\n      w = this.width,\n      h = this.height;\n    for (let step = 0; step < count; step++) {\n      // Euler diffusion CFL and reaction stiffness limits. No user-facing time-step knob.\n      const stiffness = model.excite\n        ? 1 / c.epsilon\n        : 1 + this.maxU * this.maxU + 2 * this.maxU * this.maxV + c.b;\n      const dt = Math.min(\n        0.05,\n        0.2 / Math.max(0.001, c.diffusionU, c.diffusionV),\n        0.1 / stiffness,\n      );\n      let maxU = 0,\n        maxV = 0;\n      for (let y = 0; y < h; y++) {\n        const row = y * w,\n          above = (y === 0 ? h - 1 : y - 1) * w,\n          below = (y === h - 1 ? 0 : y + 1) * w;\n        for (let x = 0; x < w; x++) {\n          const i = row + x,\n            left = row + (x === 0 ? w - 1 : x - 1),\n            right = row + (x === w - 1 ? 0 : x + 1);\n          const u = this.u[i]!,\n            v = this.v[i]!;\n          const lapU =\n            this.u[left]! +\n            this.u[right]! +\n            this.u[above + x]! +\n            this.u[below + x]! -\n            4 * u;\n          const lapV =\n            this.v[left]! +\n            this.v[right]! +\n            this.v[above + x]! +\n            this.v[below + x]! -\n            4 * v;\n          const [ru, rv] = model.reaction(u, v, c);\n          const nu = Math.max(0, u + dt * (c.diffusionU * lapU + ru));\n          const nv = Math.max(0, v + dt * (c.diffusionV * lapV + rv));\n          if (!Number.isFinite(nu) || !Number.isFinite(nv))\n            throw new RangeError(\n              \"Reaction-diffusion parameters produced a nonfinite field.\",\n            );\n          this.nextU[i] = nu;\n          this.nextV[i] = nv;\n          maxU = Math.max(maxU, nu);\n          maxV = Math.max(maxV, nv);\n        }\n      }\n      [this.u, this.nextU] = [this.nextU, this.u];\n      [this.v, this.nextV] = [this.nextV, this.v];\n      this.maxU = maxU;\n      this.maxV = maxV;\n    }\n  }\n  paint(options: { u: number; v: number; radius: number; strength: number }) {\n    const { width: w, height: h } = this;\n    for (let y = 0; y < h; y++)\n      for (let x = 0; x < w; x++) {\n        const dx = wrap(x - options.u * w + w / 2, w) - w / 2,\n          dy = wrap(y - options.v * h + h / 2, h) - h / 2;\n        const distance = Math.hypot(dx, dy) / options.radius;\n        if (distance >= 1) continue;\n        const dose = options.strength * (1 - distance * distance) ** 2;\n        const i = y * w + x;\n        this.u[i] = Math.max(\n          this.u[i]!,\n          this.options.model.excite ? dose : this.baseline * (1 + 2 * dose),\n        );\n      }\n    this.maxU = Math.max(this.maxU, this.baseline * 3);\n  }\n  stir(options: {\n    u: number;\n    v: number;\n    du: number;\n    dv: number;\n    radius: number;\n  }) {\n    const w = this.width,\n      h = this.height;\n    const sample = (field: Float32Array, x: number, y: number) => {\n      x = wrap(x, w);\n      y = wrap(y, h);\n      const ix = Math.floor(x),\n        iy = Math.floor(y),\n        fx = x - ix,\n        fy = y - iy;\n      const x1 = (ix + 1) % w,\n        y1 = (iy + 1) % h;\n      return (\n        (field[iy * w + ix]! * (1 - fx) + field[iy * w + x1]! * fx) * (1 - fy) +\n        (field[y1 * w + ix]! * (1 - fx) + field[y1 * w + x1]! * fx) * fy\n      );\n    };\n    for (let y = 0; y < h; y++)\n      for (let x = 0; x < w; x++) {\n        const dx = wrap(x - options.u * w + w / 2, w) - w / 2,\n          dy = wrap(y - options.v * h + h / 2, h) - h / 2;\n        const t = Math.max(0, 1 - Math.hypot(dx, dy) / options.radius),\n          weight = t * t * (3 - 2 * t);\n        const px = x - options.du * w * weight,\n          py = y - options.dv * h * weight,\n          i = y * w + x;\n        this.nextU[i] = sample(this.u, px, py);\n        this.nextV[i] = sample(this.v, px, py);\n      }\n    [this.u, this.nextU] = [this.nextU, this.u];\n    [this.v, this.nextV] = [this.nextV, this.v];\n  }\n  intensity(index: number) {\n    return this.options.model.excite\n      ? this.u[index]!\n      : this.u[index]! / (2 * this.baseline);\n  }\n  dispose() {\n    this.u = this.v = this.nextU = this.nextV = new Float32Array(0);\n  }\n}\n"],
+  "mappings": "AAsBO,MAAM,oBAA2B;AAAA,EACtC,aAAa,CAAC,MAAM,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM,CAAC;AAAA,EACtD,UAAU,CAAC,GAAG,GAAG,MAAM;AACrB,UAAM,IAAI,IAAI,IAAI;AAClB,WAAO,CAAC,EAAE,IAAI,IAAI,GAAG,EAAE,IAAI,CAAC;AAAA,EAC9B;AAAA,EACA,QAAQ;AACV;AACO,MAAM,mBAA0B;AAAA,EACrC,aAAa,CAAC,MAAM,CAAC,EAAE,GAAG,EAAE,IAAI,EAAE,CAAC;AAAA,EACnC,UAAU,CAAC,GAAG,GAAG,MAAM;AACrB,UAAM,IAAI,IAAI,IAAI;AAClB,WAAO,CAAC,EAAE,KAAK,EAAE,IAAI,KAAK,IAAI,GAAG,EAAE,IAAI,IAAI,CAAC;AAAA,EAC9C;AAAA,EACA,QAAQ;AACV;AACO,MAAM,eAAsB;AAAA,EACjC,aAAa,MAAM,CAAC,GAAG,CAAC;AAAA,EACxB,UAAU,CAAC,GAAG,GAAG,MAAM;AAAA,IACpB,KAAK,IAAI,MAAM,KAAK,IAAI,EAAE,KAAK,EAAE,KAAM,EAAE;AAAA,IAC1C,IAAI;AAAA,EACN;AAAA,EACA,QAAQ;AACV;AACA,SAAS,KAAK,GAAW,MAAc;AACrC,MAAI,IAAI,KAAK,KAAK,IAAI,MAAM,QAAS;AACrC,MAAI,KAAK,KAAK,IAAK,MAAM,IAAK,QAAS;AACvC,WAAS,IAAK,MAAM,QAAS,KAAK;AACpC;AACA,MAAM,OAAO,CAAC,OAAe,UAAmB,QAAQ,OAAQ,QAAQ;AACjE,MAAM,cAAc;AAAA,EAUzB,YACmB,SAMjB;AANiB;AAOjB,SAAK,QAAQ,QAAQ;AACrB,SAAK,SAAS,QAAQ;AACtB,UAAM,IAAI,KAAK,QAAQ,KAAK;AAC5B,SAAK,IAAI,IAAI,aAAa,CAAC;AAC3B,SAAK,IAAI,IAAI,aAAa,CAAC;AAC3B,SAAK,QAAQ,IAAI,aAAa,CAAC;AAC/B,SAAK,QAAQ,IAAI,aAAa,CAAC;AAC/B,UAAM,IAAI,QAAQ,QAChB,CAAC,GAAG,CAAC,IAAI,QAAQ,MAAM,YAAY,CAAC;AACtC,SAAK,WAAW,QAAQ,MAAM,SAAS,IAAI;AAC3C,aAAS,IAAI,GAAG,IAAI,KAAK,QAAQ;AAC/B,eAAS,IAAI,GAAG,IAAI,KAAK,OAAO,KAAK;AACnC,cAAM,IAAI,IAAI,KAAK,QAAQ,GACzB,QAAQ,KAAK,GAAG,EAAE,IAAI,IAAI;AAC5B,YAAI,QAAQ,MAAM,QAAQ;AACxB,gBAAM,QAAQ,OAAO,KAAK,GAAG,EAAE,IAAI,IAAI,OAAO;AAC9C,eAAK,EAAE,CAAC,KACL,IAAI,KAAK,QAAQ,QAClB,IAAI,KAAK,QAAQ,OACjB,IAAI,KAAK,SAAS,QACd,IACA,MACD,IAAI,EAAE,eAAe,OACxB,EAAE,eAAe,KAAK,IAAI,GAAG,KAAK;AACpC,eAAK,EAAE,CAAC,IAAI,KAAK,KAAK,SAAS,QAAQ,MAAM;AAAA,QAC/C,OAAO;AACL,eAAK,EAAE,CAAC,IAAI,KAAK,IAAI,EAAE,eAAe;AACtC,eAAK,EAAE,CAAC,IACN,KAAK,IAAI,EAAE,gBAAgB,KAAK,GAAG,EAAE,OAAO,UAAU,IAAI;AAAA,QAC9D;AAAA,MACF;AACF,SAAK,OAAO,KAAK,IAAI,GAAG,IAAI,GAAG;AAC/B,SAAK,OAAO,KAAK,IAAI,GAAG,IAAI,GAAG;AAC/B,SAAK,QAAQ,EAAE,WAAW;AAAA,EAC5B;AAAA,EAzCmB;AAAA,EAVnB;AAAA,EACA;AAAA,EACQ;AAAA,EACA;AAAA,EACA,OAAO;AAAA,EACP,OAAO;AAAA,EACN;AAAA,EACA;AAAA,EACA;AAAA,EA4CT,QAAQ,OAAe;AACrB,UAAM,EAAE,QAAQ,GAAG,MAAM,IAAI,KAAK,SAChC,IAAI,KAAK,OACT,IAAI,KAAK;AACX,aAAS,OAAO,GAAG,OAAO,OAAO,QAAQ;AAEvC,YAAM,YAAY,MAAM,SACpB,IAAI,EAAE,UACN,IAAI,KAAK,OAAO,KAAK,OAAO,IAAI,KAAK,OAAO,KAAK,OAAO,EAAE;AAC9D,YAAM,KAAK,KAAK;AAAA,QACd;AAAA,QACA,MAAM,KAAK,IAAI,MAAO,EAAE,YAAY,EAAE,UAAU;AAAA,QAChD,MAAM;AAAA,MACR;AACA,UAAI,OAAO,GACT,OAAO;AACT,eAAS,IAAI,GAAG,IAAI,GAAG,KAAK;AAC1B,cAAM,MAAM,IAAI,GACd,SAAS,MAAM,IAAI,IAAI,IAAI,IAAI,KAAK,GACpC,SAAS,MAAM,IAAI,IAAI,IAAI,IAAI,KAAK;AACtC,iBAAS,IAAI,GAAG,IAAI,GAAG,KAAK;AAC1B,gBAAM,IAAI,MAAM,GACd,OAAO,OAAO,MAAM,IAAI,IAAI,IAAI,IAAI,IACpC,QAAQ,OAAO,MAAM,IAAI,IAAI,IAAI,IAAI;AACvC,gBAAM,IAAI,KAAK,EAAE,CAAC,GAChB,IAAI,KAAK,EAAE,CAAC;AACd,gBAAM,OACJ,KAAK,EAAE,IAAI,IACX,KAAK,EAAE,KAAK,IACZ,KAAK,EAAE,QAAQ,CAAC,IAChB,KAAK,EAAE,QAAQ,CAAC,IAChB,IAAI;AACN,gBAAM,OACJ,KAAK,EAAE,IAAI,IACX,KAAK,EAAE,KAAK,IACZ,KAAK,EAAE,QAAQ,CAAC,IAChB,KAAK,EAAE,QAAQ,CAAC,IAChB,IAAI;AACN,gBAAM,CAAC,IAAI,EAAE,IAAI,MAAM,SAAS,GAAG,GAAG,CAAC;AACvC,gBAAM,KAAK,KAAK,IAAI,GAAG,IAAI,MAAM,EAAE,aAAa,OAAO,GAAG;AAC1D,gBAAM,KAAK,KAAK,IAAI,GAAG,IAAI,MAAM,EAAE,aAAa,OAAO,GAAG;AAC1D,cAAI,CAAC,OAAO,SAAS,EAAE,KAAK,CAAC,OAAO,SAAS,EAAE;AAC7C,kBAAM,IAAI;AAAA,cACR;AAAA,YACF;AACF,eAAK,MAAM,CAAC,IAAI;AAChB,eAAK,MAAM,CAAC,IAAI;AAChB,iBAAO,KAAK,IAAI,MAAM,EAAE;AACxB,iBAAO,KAAK,IAAI,MAAM,EAAE;AAAA,QAC1B;AAAA,MACF;AACA,OAAC,KAAK,GAAG,KAAK,KAAK,IAAI,CAAC,KAAK,OAAO,KAAK,CAAC;AAC1C,OAAC,KAAK,GAAG,KAAK,KAAK,IAAI,CAAC,KAAK,OAAO,KAAK,CAAC;AAC1C,WAAK,OAAO;AACZ,WAAK,OAAO;AAAA,IACd;AAAA,EACF;AAAA,EACA,MAAM,SAAqE;AACzE,UAAM,EAAE,OAAO,GAAG,QAAQ,EAAE,IAAI;AAChC,aAAS,IAAI,GAAG,IAAI,GAAG;AACrB,eAAS,IAAI,GAAG,IAAI,GAAG,KAAK;AAC1B,cAAM,KAAK,KAAK,IAAI,QAAQ,IAAI,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAClD,KAAK,KAAK,IAAI,QAAQ,IAAI,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI;AAChD,cAAM,WAAW,KAAK,MAAM,IAAI,EAAE,IAAI,QAAQ;AAC9C,YAAI,YAAY,EAAG;AACnB,cAAM,OAAO,QAAQ,YAAY,IAAI,WAAW,aAAa;AAC7D,cAAM,IAAI,IAAI,IAAI;AAClB,aAAK,EAAE,CAAC,IAAI,KAAK;AAAA,UACf,KAAK,EAAE,CAAC;AAAA,UACR,KAAK,QAAQ,MAAM,SAAS,OAAO,KAAK,YAAY,IAAI,IAAI;AAAA,QAC9D;AAAA,MACF;AACF,SAAK,OAAO,KAAK,IAAI,KAAK,MAAM,KAAK,WAAW,CAAC;AAAA,EACnD;AAAA,EACA,KAAK,SAMF;AACD,UAAM,IAAI,KAAK,OACb,IAAI,KAAK;AACX,UAAM,SAAS,CAAC,OAAqB,GAAW,MAAc;AAC5D,UAAI,KAAK,GAAG,CAAC;AACb,UAAI,KAAK,GAAG,CAAC;AACb,YAAM,KAAK,KAAK,MAAM,CAAC,GACrB,KAAK,KAAK,MAAM,CAAC,GACjB,KAAK,IAAI,IACT,KAAK,IAAI;AACX,YAAM,MAAM,KAAK,KAAK,GACpB,MAAM,KAAK,KAAK;AAClB,cACG,MAAM,KAAK,IAAI,EAAE,KAAM,IAAI,MAAM,MAAM,KAAK,IAAI,EAAE,IAAK,OAAO,IAAI,OAClE,MAAM,KAAK,IAAI,EAAE,KAAM,IAAI,MAAM,MAAM,KAAK,IAAI,EAAE,IAAK,MAAM;AAAA,IAElE;AACA,aAAS,IAAI,GAAG,IAAI,GAAG;AACrB,eAAS,IAAI,GAAG,IAAI,GAAG,KAAK;AAC1B,cAAM,KAAK,KAAK,IAAI,QAAQ,IAAI,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,GAClD,KAAK,KAAK,IAAI,QAAQ,IAAI,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI;AAChD,cAAM,IAAI,KAAK,IAAI,GAAG,IAAI,KAAK,MAAM,IAAI,EAAE,IAAI,QAAQ,MAAM,GAC3D,SAAS,IAAI,KAAK,IAAI,IAAI;AAC5B,cAAM,KAAK,IAAI,QAAQ,KAAK,IAAI,QAC9B,KAAK,IAAI,QAAQ,KAAK,IAAI,QAC1B,IAAI,IAAI,IAAI;AACd,aAAK,MAAM,CAAC,IAAI,OAAO,KAAK,GAAG,IAAI,EAAE;AACrC,aAAK,MAAM,CAAC,IAAI,OAAO,KAAK,GAAG,IAAI,EAAE;AAAA,MACvC;AACF,KAAC,KAAK,GAAG,KAAK,KAAK,IAAI,CAAC,KAAK,OAAO,KAAK,CAAC;AAC1C,KAAC,KAAK,GAAG,KAAK,KAAK,IAAI,CAAC,KAAK,OAAO,KAAK,CAAC;AAAA,EAC5C;AAAA,EACA,UAAU,OAAe;AACvB,WAAO,KAAK,QAAQ,MAAM,SACtB,KAAK,EAAE,KAAK,IACZ,KAAK,EAAE,KAAK,KAAM,IAAI,KAAK;AAAA,EACjC;AAAA,EACA,UAAU;AACR,SAAK,IAAI,KAAK,IAAI,KAAK,QAAQ,KAAK,QAAQ,IAAI,aAAa,CAAC;AAAA,EAChE;AACF;",
+  "names": []
+}

package/dist/features/textures/lib/chemical-recipes.d.ts

@@ -0,0 +1,119 @@
+import { type TextureRecipe } from "@pibbl/core";
+import type { InkMessage } from "./definitions.js";
+/** Shared simulation timing, diffusion and initialization controls. See {@link SchnakenbergOptions}. */
+interface ChemicalOptions {
+    /** Unsigned 32-bit initialization seed; default 42. See {@link schnakenberg}. */
+    readonly seed?: number;
+    /** Activator diffusion per grid cell, finite 0–32. See {@link schnakenberg}. */
+    readonly diffusionU?: number;
+    /** Inhibitor diffusion per grid cell, finite 0–32. See {@link schnakenberg}. */
+    readonly diffusionV?: number;
+    /** Initial seeded disturbance amplitude, finite 0–1; default 0.1. See {@link schnakenberg}. */
+    readonly initialNoise?: number;
+    /** Initialization steps, integer 0–12000; default 4096. See {@link schnakenberg}. */
+    readonly warmupSteps?: number;
+    /** Multiplier of 480 numerical steps/second, finite 0–4; default 1. Zero idles. See {@link schnakenberg}. */
+    readonly speed?: number;
+    /** Initially suspend stepping; default false. See {@link InkMessage}. */
+    readonly paused?: boolean;
+}
+/** Schnakenberg production rates, diffusion and initialization. See {@link schnakenberg}. */
+export interface SchnakenbergOptions extends ChemicalOptions {
+    /** Activator source, finite 0.01–2; default 0.1. See {@link schnakenberg}. */
+    readonly a?: number;
+    /** Inhibitor source, finite 0.01–2; default 0.9. See {@link schnakenberg}. */
+    readonly b?: number;
+}
+/** Excitability, recovery, diffusion and initialization. See {@link barkley}. */
+export interface BarkleyOptions extends ChemicalOptions {
+    /** Excitability scale, finite 0.3–1.2; default 0.75. See {@link barkley}. */
+    readonly a?: number;
+    /** Excitation threshold offset, finite 0–0.2; default 0.02. See {@link barkley}. */
+    readonly b?: number;
+    /** Fast/slow time separation, finite 0.005–0.1; default 0.02. See {@link barkley}. */
+    readonly epsilon?: number;
+}
+/** Brusselator feed rates, diffusion and initialization. See {@link brusselator}. */
+export interface BrusselatorOptions extends ChemicalOptions {
+    /** Activator feed, finite 0.1–3; default 1. See {@link brusselator}. */
+    readonly a?: number;
+    /** Conversion rate, finite 0.1–5; default 1.8. See {@link brusselator}. */
+    readonly b?: number;
+}
+/** Immutable Schnakenberg starting configurations. See {@link SchnakenbergOptions}. */
+export declare const schnakenbergPresets: Readonly<{
+    /** Separated activator spots. See {@link schnakenberg}. */
+    spots: Readonly<{
+        a: 0.1;
+        b: 0.9;
+        diffusionU: 1;
+        diffusionV: 20;
+    }>;
+    /** Connected activator bands. See {@link schnakenberg}. */
+    labyrinths: Readonly<{
+        a: 0.01;
+        b: 1.7;
+        diffusionU: 1;
+        diffusionV: 20;
+        warmupSteps: 10000;
+    }>;
+}>;
+/** Immutable Barkley starting configurations. See {@link BarkleyOptions}. */
+export declare const barkleyPresets: Readonly<{
+    /** Broken excitation fronts form rotating spirals. See {@link barkley}. */
+    spirals: Readonly<{
+        a: 0.75;
+        b: 0.02;
+        epsilon: 0.02;
+        diffusionU: 1;
+        diffusionV: 0;
+    }>;
+    /** Broader propagating excitation fronts. See {@link barkley}. */
+    waves: Readonly<{
+        a: 0.9;
+        b: 0.01;
+        epsilon: 0.04;
+        diffusionU: 1;
+        diffusionV: 0;
+    }>;
+}>;
+/** Immutable Brusselator starting configurations. See {@link BrusselatorOptions}. */
+export declare const brusselatorPresets: Readonly<{
+    /** Differential diffusion supports stationary spots. See {@link brusselator}. */
+    spots: Readonly<{
+        a: 1;
+        b: 1.9;
+        diffusionU: 0.2;
+        diffusionV: 2;
+        warmupSteps: 10000;
+    }>;
+    /** Equal diffusion above the oscillation threshold. See {@link brusselator}. */
+    oscillations: Readonly<{
+        a: 1;
+        b: 3;
+        diffusionU: 0.3;
+        diffusionV: 0.3;
+    }>;
+}>;
+/**
+ * Creates a periodic Schnakenberg activator/inhibitor texture with paint, stir, pause and reset.
+ * @param options - Reaction rates, diffusion and initialization. See {@link SchnakenbergOptions}.
+ * @returns A mountable simulation recipe. See {@link TextureRecipe}, {@link InkMessage}.
+ * @see {@link schnakenbergPresets}
+ */
+export declare function schnakenberg(options?: SchnakenbergOptions): TextureRecipe<InkMessage>;
+/**
+ * Creates a periodic Barkley excitable texture initialized with a broken excitation front.
+ * @param options - Excitability, recovery, diffusion and initialization. See {@link BarkleyOptions}.
+ * @returns A mountable simulation recipe. See {@link TextureRecipe}, {@link InkMessage}.
+ * @see {@link barkleyPresets}
+ */
+export declare function barkley(options?: BarkleyOptions): TextureRecipe<InkMessage>;
+/**
+ * Creates a periodic Brusselator texture with paint, stir, pause and reset.
+ * @param options - Feed rates, diffusion and initialization. See {@link BrusselatorOptions}.
+ * @returns A mountable simulation recipe. See {@link TextureRecipe}, {@link InkMessage}.
+ * @see {@link brusselatorPresets}
+ */
+export declare function brusselator(options?: BrusselatorOptions): TextureRecipe<InkMessage>;
+export {};

package/dist/features/textures/lib/chemical-recipes.js

@@ -0,0 +1,275 @@
+import {
+  defineTexture
+} from "@pibbl/core";
+import { paletteLut, texturePalettes, validateSeed } from "./palette.js";
+import {
+  ChemicalField,
+  schnakenbergModel,
+  barkleyModel,
+  brusselatorModel
+} from "./chemical-field.js";
+const schnakenbergPresets = Object.freeze({
+  /** Separated activator spots. See {@link schnakenberg}. */
+  spots: Object.freeze({ a: 0.1, b: 0.9, diffusionU: 1, diffusionV: 20 }),
+  /** Connected activator bands. See {@link schnakenberg}. */
+  labyrinths: Object.freeze({
+    a: 0.01,
+    b: 1.7,
+    diffusionU: 1,
+    diffusionV: 20,
+    warmupSteps: 1e4
+  })
+});
+const barkleyPresets = Object.freeze({
+  /** Broken excitation fronts form rotating spirals. See {@link barkley}. */
+  spirals: Object.freeze({
+    a: 0.75,
+    b: 0.02,
+    epsilon: 0.02,
+    diffusionU: 1,
+    diffusionV: 0
+  }),
+  /** Broader propagating excitation fronts. See {@link barkley}. */
+  waves: Object.freeze({
+    a: 0.9,
+    b: 0.01,
+    epsilon: 0.04,
+    diffusionU: 1,
+    diffusionV: 0
+  })
+});
+const brusselatorPresets = Object.freeze({
+  /** Differential diffusion supports stationary spots. See {@link brusselator}. */
+  spots: Object.freeze({
+    a: 1,
+    b: 1.9,
+    diffusionU: 0.2,
+    diffusionV: 2,
+    warmupSteps: 1e4
+  }),
+  /** Equal diffusion above the oscillation threshold. See {@link brusselator}. */
+  oscillations: Object.freeze({ a: 1, b: 3, diffusionU: 0.3, diffusionV: 0.3 })
+});
+function bounded(name, value, min, max, integer = false) {
+  if (!Number.isFinite(value) || value < min || value > max || integer && !Number.isInteger(value))
+    throw new RangeError(
+      `${name} must be ${integer ? "an integer " : ""}within ${min}..${max}.`
+    );
+  return value;
+}
+function common(options) {
+  if (options.paused !== void 0 && typeof options.paused !== "boolean")
+    throw new TypeError("Chemical paused must be boolean.");
+  return {
+    seed: validateSeed(options.seed),
+    initialNoise: bounded(
+      "Chemical initialNoise",
+      options.initialNoise ?? 0.1,
+      0,
+      1
+    ),
+    warmupSteps: bounded(
+      "Chemical warmupSteps",
+      options.warmupSteps ?? 4096,
+      0,
+      12e3,
+      true
+    ),
+    speed: bounded("Chemical speed", options.speed ?? 1, 0, 4),
+    paused: options.paused ?? false
+  };
+}
+function diffusion(options, defaults) {
+  return {
+    diffusionU: bounded(
+      "Chemical diffusionU",
+      options.diffusionU ?? defaults.diffusionU,
+      0,
+      32
+    ),
+    diffusionV: bounded(
+      "Chemical diffusionV",
+      options.diffusionV ?? defaults.diffusionV,
+      0,
+      32
+    )
+  };
+}
+function schnakenbergConfig(options) {
+  return {
+    ...common(options),
+    ...diffusion(options, schnakenbergPresets.spots),
+    a: bounded("Schnakenberg a", options.a ?? 0.1, 0.01, 2),
+    b: bounded("Schnakenberg b", options.b ?? 0.9, 0.01, 2),
+    epsilon: 0.02
+  };
+}
+function barkleyConfig(options) {
+  return {
+    ...common(options),
+    ...diffusion(options, barkleyPresets.spirals),
+    a: bounded("Barkley a", options.a ?? 0.75, 0.3, 1.2),
+    b: bounded("Barkley b", options.b ?? 0.02, 0, 0.2),
+    epsilon: bounded("Barkley epsilon", options.epsilon ?? 0.02, 5e-3, 0.1)
+  };
+}
+function brusselatorConfig(options) {
+  return {
+    ...common(options),
+    ...diffusion(options, brusselatorPresets.spots),
+    a: bounded("Brusselator a", options.a ?? 1, 0.1, 3),
+    b: bounded("Brusselator b", options.b ?? 1.8, 0.1, 5),
+    epsilon: 0.02
+  };
+}
+function chemicalRecipe(model) {
+  const makeField = (context, config) => {
+    if (context.canvas.width < 3 || context.canvas.height < 3)
+      throw new RangeError(
+        "Chemical textures require at least 3\xD73 raster samples."
+      );
+    return new ChemicalField({
+      width: context.canvas.width - 1,
+      height: context.canvas.height - 1,
+      config,
+      model
+    });
+  };
+  const active = (state) => !state.paused && state.config.speed > 0;
+  const reset = (state, context) => {
+    const next = makeField(context, state.config);
+    state.field.dispose();
+    state.field = next;
+    state.elapsed = 0;
+    state.fresh = true;
+    context.invalidate();
+    if (active(state)) context.requestFrame();
+  };
+  return defineTexture({
+    defaultPlacement: { mode: "repeat" },
+    create(context, config) {
+      const state = {
+        config,
+        field: makeField(context, config),
+        paused: config.paused,
+        elapsed: 0,
+        fresh: true,
+        image: new ImageData(context.canvas.width, context.canvas.height),
+        lut: void 0,
+        paletteKey: ""
+      };
+      if (active(state)) context.requestFrame();
+      return state;
+    },
+    update(state, config, context) {
+      if (JSON.stringify(config) === JSON.stringify(state.config)) return;
+      const { speed: _speed, paused: _paused, ...previous } = state.config;
+      const { speed, paused, ...next } = config;
+      const changed = JSON.stringify(previous) !== JSON.stringify(next);
+      if (paused !== state.config.paused) state.paused = paused;
+      state.config = config;
+      state.elapsed = 0;
+      state.fresh = true;
+      if (changed) reset(state, context);
+      if (!state.paused && speed > 0) context.requestFrame();
+    },
+    receive(state, message, context) {
+      if (message.type === "pause") {
+        if (typeof message.paused !== "boolean")
+          throw new TypeError("Chemical pause requires a boolean.");
+        state.paused = message.paused;
+        state.elapsed = 0;
+        state.fresh = true;
+        if (active(state)) context.requestFrame();
+      } else if (message.type === "reset") reset(state, context);
+      else if (message.type === "paint" || message.type === "stir") {
+        const radius = bounded(
+          "Chemical brush radius",
+          message.radius ?? (message.type === "paint" ? 6 : 16),
+          0.01,
+          128
+        );
+        if (!Number.isFinite(message.u) || !Number.isFinite(message.v))
+          throw new RangeError("Chemical brush coordinates must be finite.");
+        if (message.type === "paint") {
+          const strength = bounded(
+            "Chemical brush strength",
+            message.strength ?? 0.8,
+            0,
+            1
+          );
+          if (message.u < 0 || message.u > 1 || message.v < 0 || message.v > 1)
+            return;
+          if (strength === 0) return;
+          state.field.paint({ ...message, radius, strength });
+        } else {
+          const du = bounded("Chemical stir du", message.du, -1, 1), dv = bounded("Chemical stir dv", message.dv, -1, 1);
+          if (message.u < 0 || message.u > 1 || message.v < 0 || message.v > 1)
+            return;
+          state.field.stir({ ...message, radius, du, dv });
+        }
+        context.invalidate();
+      } else throw new TypeError("Unknown chemical texture message.");
+    },
+    advance(state, frame, context) {
+      if (!active(state)) return;
+      if (!state.fresh) {
+        state.elapsed += Math.min(50, frame.delta) * state.config.speed * 0.48;
+        const count = Math.min(96, Math.floor(state.elapsed));
+        state.elapsed -= count;
+        if (count > 0) {
+          state.field.advance(count);
+          context.invalidate();
+        }
+      }
+      state.fresh = false;
+      context.requestFrame();
+    },
+    rasterize(state, context) {
+      const key = (context.palette ?? texturePalettes.lagoon).join();
+      if (key !== state.paletteKey) {
+        state.lut = paletteLut(context.palette ?? texturePalettes.lagoon);
+        state.paletteKey = key;
+      }
+      const { width, height } = context.canvas, data = state.image.data, lut = state.lut;
+      for (let y = 0; y < height; y++)
+        for (let x = 0; x < width; x++) {
+          const index = y % state.field.height * state.field.width + x % state.field.width;
+          const color = Math.round(
+            Math.max(0, Math.min(1, state.field.intensity(index))) * 255
+          ) * 4, i = (y * width + x) * 4;
+          data[i] = lut[color];
+          data[i + 1] = lut[color + 1];
+          data[i + 2] = lut[color + 2];
+          data[i + 3] = 255;
+        }
+      context.canvas.getContext("2d").putImageData(state.image, 0, 0);
+    },
+    dispose(state) {
+      state.field.dispose();
+      state.image = void 0;
+      state.lut = void 0;
+    }
+  });
+}
+const schnakenbergRecipe = /* @__PURE__ */ chemicalRecipe(schnakenbergModel);
+const barkleyRecipe = /* @__PURE__ */ chemicalRecipe(barkleyModel);
+const brusselatorRecipe = /* @__PURE__ */ chemicalRecipe(brusselatorModel);
+function schnakenberg(options = {}) {
+  return schnakenbergRecipe(schnakenbergConfig(options));
+}
+function barkley(options = {}) {
+  return barkleyRecipe(barkleyConfig(options));
+}
+function brusselator(options = {}) {
+  return brusselatorRecipe(brusselatorConfig(options));
+}
+export {
+  barkley,
+  barkleyPresets,
+  brusselator,
+  brusselatorPresets,
+  schnakenberg,
+  schnakenbergPresets
+};
+//# sourceMappingURL=chemical-recipes.js.map