@triforge/fluid-core 0.1.0

package/dist/FLIPSimulator.d.ts

@@ -0,0 +1,98 @@
+import { MACGrid } from './MACGrid.js';
+export interface FLIPSimulatorOptions {
+    /**
+     * Grid resolution along each axis.
+     * Higher = more detail, higher cost. Default 24.
+     * Blender: Fluid > Resolution Divisions.
+     */
+    resolution?: number;
+    /**
+     * World-space size of the fluid domain cube.
+     * Default 2 (a 2×2×2 metre box).
+     */
+    domainSize?: number;
+    /**
+     * World-space origin (bottom-left-front corner of the domain).
+     * Default [−1, −1, −1].
+     */
+    origin?: [number, number, number];
+    /**
+     * Number of fluid particles. Default 2000.
+     * Blender: Fluid > Particle Count.
+     */
+    particleCount?: number;
+    /**
+     * FLIP blend factor [0–1].
+     * 1 = pure FLIP (energetic, noisy), 0 = pure PIC (damped, stable).
+     * 0.95 matches Blender's default. Blender: Fluid > FLIP Ratio.
+     */
+    flipRatio?: number;
+    /** Gravity vector [x,y,z]. Default [0,−9.8,0]. */
+    gravity?: [number, number, number];
+    /** Pressure solver maximum iterations. Default 80. */
+    solverMaxIter?: number;
+    /** Pressure solver convergence tolerance. Default 1e-4. */
+    solverTolerance?: number;
+    /** Simulation substeps per frame. Default 3. */
+    substeps?: number;
+}
+/**
+ * FLIPSimulator — full FLIP (Fluid-Implicit-Particle) fluid simulation.
+ *
+ * Implements the Zhu & Bridson 2005 FLIP algorithm:
+ *   1. Transfer particle velocities → MAC grid (P2G, trilinear splat)
+ *   2. Apply body forces (gravity) on grid
+ *   3. Enforce solid boundary conditions
+ *   4. Solve pressure Poisson equation (conjugate gradient)
+ *   5. Subtract pressure gradient from grid velocity (projection)
+ *   6. Transfer grid velocity change back to particles (G2P, FLIP blend)
+ *   7. Advect particles through the corrected velocity field
+ *
+ * All parameters live in `parameters` for GSAP/st-keyframe animation.
+ *
+ * ```typescript
+ * const sim = new FLIPSimulator({ resolution: 20, particleCount: 1500 })
+ * sim.fillBox(-0.5, -0.5, -0.5, 0.5, 0.5, 0.5)
+ *
+ * // Each frame:
+ * sim.step(dt)
+ * const pts = sim.getPositions()  // Float32Array  [x,y,z, x,y,z, ...]
+ * ```
+ */
+export declare class FLIPSimulator {
+    parameters: Record<string, number>;
+    readonly particleCount: number;
+    private readonly _grid;
+    private readonly _particles;
+    private readonly _b;
+    private readonly _x;
+    private readonly _nx;
+    private readonly _ny;
+    private readonly _nz;
+    constructor(opts?: FLIPSimulatorOptions);
+    /**
+     * Fill a world-space axis-aligned box with fluid particles.
+     * Particles are distributed using a jittered grid to avoid clumping.
+     * Multiple calls accumulate particles (capped at particleCount).
+     */
+    fillBox(x0: number, y0: number, z0: number, x1: number, y1: number, z1: number): void;
+    /** Advance the simulation by `dt` seconds (call once per frame, e.g. dt=1/60). */
+    step(dt: number): void;
+    private _substep;
+    private _p2g;
+    private _splatToFace;
+    private _pressureSolve;
+    private _g2p;
+    private _interpOld;
+    private _advect;
+    private _clampParticles;
+    /** Returns particle positions as a flat Float32Array [x,y,z, x,y,z, ...]. */
+    getPositions(): Float32Array;
+    /** Returns particle velocities as a flat Float32Array [vx,vy,vz, ...]. */
+    getVelocities(): Float32Array;
+    /** Number of live particles. */
+    get liveCount(): number;
+    /** The underlying MAC grid (for advanced use / surface reconstruction). */
+    get grid(): MACGrid;
+}
+//# sourceMappingURL=FLIPSimulator.d.ts.map

package/dist/FLIPSimulator.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"FLIPSimulator.d.ts","sourceRoot":"","sources":["../src/FLIPSimulator.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAA;AAEtC,MAAM,WAAW,oBAAoB;IACnC;;;;OAIG;IACH,UAAU,CAAC,EAAE,MAAM,CAAA;IAEnB;;;OAGG;IACH,UAAU,CAAC,EAAE,MAAM,CAAA;IAEnB;;;OAGG;IACH,MAAM,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAA;IAEjC;;;OAGG;IACH,aAAa,CAAC,EAAE,MAAM,CAAA;IAEtB;;;;OAIG;IACH,SAAS,CAAC,EAAE,MAAM,CAAA;IAElB,kDAAkD;IAClD,OAAO,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAA;IAElC,sDAAsD;IACtD,aAAa,CAAC,EAAE,MAAM,CAAA;IAEtB,2DAA2D;IAC3D,eAAe,CAAC,EAAE,MAAM,CAAA;IAExB,gDAAgD;IAChD,QAAQ,CAAC,EAAE,MAAM,CAAA;CAClB;AAQD;;;;;;;;;;;;;;;;;;;;;;GAsBG;AACH,qBAAa,aAAa;IACxB,UAAU,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAA;IAElC,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAA;IAC9B,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAS;IAC/B,OAAO,CAAC,QAAQ,CAAC,UAAU,CAAgB;IAG3C,OAAO,CAAC,QAAQ,CAAC,EAAE,CAAe;IAClC,OAAO,CAAC,QAAQ,CAAC,EAAE,CAAe;IAClC,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAQ;IAC5B,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAQ;IAC5B,OAAO,CAAC,QAAQ,CAAC,GAAG,CAAQ;gBAEhB,IAAI,GAAE,oBAAyB;IA8B3C;;;;OAIG;IACH,OAAO,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAClC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,IAAI;IAqBjD,kFAAkF;IAClF,IAAI,CAAC,EAAE,EAAE,MAAM,GAAG,IAAI;IAMtB,OAAO,CAAC,QAAQ;IAwChB,OAAO,CAAC,IAAI;IAsBZ,OAAO,CAAC,YAAY;IA0BpB,OAAO,CAAC,cAAc;IAgHtB,OAAO,CAAC,IAAI;IAwBZ,OAAO,CAAC,UAAU;IAoBlB,OAAO,CAAC,OAAO;IAQf,OAAO,CAAC,eAAe;IAmBvB,6EAA6E;IAC7E,YAAY,IAAI,YAAY;IAU5B,0EAA0E;IAC1E,aAAa,IAAI,YAAY;IAU7B,gCAAgC;IAChC,IAAI,SAAS,IAAI,MAAM,CAAkC;IAEzD,2EAA2E;IAC3E,IAAI,IAAI,IAAI,OAAO,CAAsB;CAC1C"}

package/dist/FLIPSimulator.js

@@ -0,0 +1,404 @@
+import pcg from 'conjugate-gradient';
+import { MACGrid } from './MACGrid.js';
+/**
+ * FLIPSimulator — full FLIP (Fluid-Implicit-Particle) fluid simulation.
+ *
+ * Implements the Zhu & Bridson 2005 FLIP algorithm:
+ *   1. Transfer particle velocities → MAC grid (P2G, trilinear splat)
+ *   2. Apply body forces (gravity) on grid
+ *   3. Enforce solid boundary conditions
+ *   4. Solve pressure Poisson equation (conjugate gradient)
+ *   5. Subtract pressure gradient from grid velocity (projection)
+ *   6. Transfer grid velocity change back to particles (G2P, FLIP blend)
+ *   7. Advect particles through the corrected velocity field
+ *
+ * All parameters live in `parameters` for GSAP/st-keyframe animation.
+ *
+ * ```typescript
+ * const sim = new FLIPSimulator({ resolution: 20, particleCount: 1500 })
+ * sim.fillBox(-0.5, -0.5, -0.5, 0.5, 0.5, 0.5)
+ *
+ * // Each frame:
+ * sim.step(dt)
+ * const pts = sim.getPositions()  // Float32Array  [x,y,z, x,y,z, ...]
+ * ```
+ */
+export class FLIPSimulator {
+    constructor(opts = {}) {
+        const res = opts.resolution ?? 24;
+        const domainSize = opts.domainSize ?? 2;
+        const origin = opts.origin ?? [-1, -1, -1];
+        const h = domainSize / res;
+        this._nx = res;
+        this._ny = res;
+        this._nz = res;
+        this._grid = new MACGrid(res, res, res, h, origin[0], origin[1], origin[2]);
+        this.particleCount = opts.particleCount ?? 2000;
+        this._particles = [];
+        this._b = new Float64Array(res * res * res);
+        this._x = new Float64Array(res * res * res);
+        const g = opts.gravity ?? [0, -9.8, 0];
+        this.parameters = {
+            flipRatio: opts.flipRatio ?? 0.95,
+            gravityX: g[0],
+            gravityY: g[1],
+            gravityZ: g[2],
+            substeps: opts.substeps ?? 3,
+            solverMaxIter: opts.solverMaxIter ?? 80,
+            solverTolerance: opts.solverTolerance ?? 1e-4,
+        };
+    }
+    // ── Particle placement ────────────────────────────────────────────────────
+    /**
+     * Fill a world-space axis-aligned box with fluid particles.
+     * Particles are distributed using a jittered grid to avoid clumping.
+     * Multiple calls accumulate particles (capped at particleCount).
+     */
+    fillBox(x0, y0, z0, x1, y1, z1) {
+        this._particles.length = 0;
+        const h = this._grid.h;
+        const pph = 2; // particles per cell per axis (2³ = 8 per cell, standard FLIP)
+        const sp = h / pph;
+        for (let z = z0; z < z1; z += sp)
+            for (let y = y0; y < y1; y += sp)
+                for (let x = x0; x < x1; x += sp) {
+                    if (this._particles.length >= this.particleCount)
+                        break;
+                    this._particles.push({
+                        x: x + sp * Math.random(),
+                        y: y + sp * Math.random(),
+                        z: z + sp * Math.random(),
+                        vx: 0, vy: 0, vz: 0,
+                    });
+                }
+    }
+    // ── Main step ─────────────────────────────────────────────────────────────
+    /** Advance the simulation by `dt` seconds (call once per frame, e.g. dt=1/60). */
+    step(dt) {
+        const substeps = Math.max(1, Math.round(this.parameters.substeps));
+        const subDt = dt / substeps;
+        for (let s = 0; s < substeps; s++)
+            this._substep(subDt);
+    }
+    _substep(dt) {
+        const grid = this._grid;
+        const { nx, ny, nz, h } = grid;
+        // 1. Mark cells as fluid / air based on particle positions
+        const fluidSet = new Set();
+        for (const p of this._particles) {
+            const i = Math.floor((p.x - grid.ox) / h);
+            const j = Math.floor((p.y - grid.oy) / h);
+            const k = Math.floor((p.z - grid.oz) / h);
+            if (i >= 1 && i < nx - 1 && j >= 1 && j < ny - 1 && k >= 1 && k < nz - 1)
+                fluidSet.add(grid.idx(i, j, k));
+        }
+        grid.markCells((i, j, k) => fluidSet.has(grid.idx(i, j, k)));
+        // 2. P2G — splat particle velocities to grid faces
+        grid.clearVelocities();
+        this._p2g();
+        grid.normalizeWeights();
+        grid.saveVelocities();
+        // 3. Body forces + boundary
+        grid.applyGravity(this.parameters.gravityX, this.parameters.gravityY, this.parameters.gravityZ, dt);
+        grid.enforceBoundary();
+        // 4. Pressure projection
+        this._pressureSolve(dt);
+        // 5. G2P — transfer corrected grid velocities back to particles (FLIP blend)
+        this._g2p();
+        // 6. Advect particles
+        this._advect(dt);
+        // 7. Clamp particles inside domain
+        this._clampParticles();
+    }
+    // ── P2G (particle → grid) ─────────────────────────────────────────────────
+    _p2g() {
+        const grid = this._grid;
+        const { nx, ny, nz, h, ox, oy, oz } = grid;
+        for (const p of this._particles) {
+            // U faces (i-½, j, k) — offset by (0, h/2, h/2)
+            this._splatToFace(grid.u, grid.uW, p.vx, p.x - ox, p.y - oy - 0.5 * h, p.z - oz - 0.5 * h, nx + 1, ny, nz, h);
+            // V faces (i, j-½, k) — offset by (h/2, 0, h/2)
+            this._splatToFace(grid.v, grid.vW, p.vy, p.x - ox - 0.5 * h, p.y - oy, p.z - oz - 0.5 * h, nx, ny + 1, nz, h);
+            // W faces (i, j, k-½) — offset by (h/2, h/2, 0)
+            this._splatToFace(grid.w, grid.wW, p.vz, p.x - ox - 0.5 * h, p.y - oy - 0.5 * h, p.z - oz, nx, ny, nz + 1, h);
+        }
+    }
+    _splatToFace(arr, wArr, val, lx, ly, lz, Nx, Ny, Nz, h) {
+        let fi = lx / h, fj = ly / h, fk = lz / h;
+        fi = Math.max(0, Math.min(Nx - 1.0001, fi));
+        fj = Math.max(0, Math.min(Ny - 1.0001, fj));
+        fk = Math.max(0, Math.min(Nz - 1.0001, fk));
+        const i0 = Math.floor(fi), i1 = i0 + 1;
+        const j0 = Math.floor(fj), j1 = j0 + 1;
+        const k0 = Math.floor(fk), k1 = k0 + 1;
+        const tx = fi - i0, ty = fj - j0, tz = fk - k0;
+        const idx = (i, j, k) => i + Nx * (j + Ny * k);
+        const w000 = (1 - tx) * (1 - ty) * (1 - tz);
+        arr[idx(i0, j0, k0)] += w000 * val;
+        wArr[idx(i0, j0, k0)] += w000;
+        const w100 = tx * (1 - ty) * (1 - tz);
+        arr[idx(i1, j0, k0)] += w100 * val;
+        wArr[idx(i1, j0, k0)] += w100;
+        const w010 = (1 - tx) * ty * (1 - tz);
+        arr[idx(i0, j1, k0)] += w010 * val;
+        wArr[idx(i0, j1, k0)] += w010;
+        const w110 = tx * ty * (1 - tz);
+        arr[idx(i1, j1, k0)] += w110 * val;
+        wArr[idx(i1, j1, k0)] += w110;
+        const w001 = (1 - tx) * (1 - ty) * tz;
+        arr[idx(i0, j0, k1)] += w001 * val;
+        wArr[idx(i0, j0, k1)] += w001;
+        const w101 = tx * (1 - ty) * tz;
+        arr[idx(i1, j0, k1)] += w101 * val;
+        wArr[idx(i1, j0, k1)] += w101;
+        const w011 = (1 - tx) * ty * tz;
+        arr[idx(i0, j1, k1)] += w011 * val;
+        wArr[idx(i0, j1, k1)] += w011;
+        const w111 = tx * ty * tz;
+        arr[idx(i1, j1, k1)] += w111 * val;
+        wArr[idx(i1, j1, k1)] += w111;
+    }
+    // ── Pressure projection ───────────────────────────────────────────────────
+    _pressureSolve(dt) {
+        const grid = this._grid;
+        const { nx, ny, nz, h } = grid;
+        const N = nx * ny * nz;
+        const invH2 = 1.0 / (h * h);
+        // Build RHS: b[i] = div(u)[i] / dt  (standard incompressible projection)
+        this._b.fill(0);
+        for (let k = 1; k < nz - 1; k++)
+            for (let j = 1; j < ny - 1; j++)
+                for (let i = 1; i < nx - 1; i++) {
+                    const ci = grid.idx(i, j, k);
+                    if (grid.cellType[ci] !== 1)
+                        continue;
+                    const div = (grid.u[grid.uIdx(i + 1, j, k)] - grid.u[grid.uIdx(i, j, k)]) / h
+                        + (grid.v[grid.vIdx(i, j + 1, k)] - grid.v[grid.vIdx(i, j, k)]) / h
+                        + (grid.w[grid.wIdx(i, j, k + 1)] - grid.w[grid.wIdx(i, j, k)]) / h;
+                    this._b[ci] = div / dt;
+                }
+        this._x.fill(0);
+        const cellType = grid.cellType;
+        const idxFn = (i, j, k) => grid.idx(i, j, k);
+        // A = discrete positive-definite Laplacian  (1/h² per neighbour pair)
+        const A = {
+            rowCount: N,
+            get: (row, _col) => {
+                if (cellType[row] !== 1)
+                    return 1;
+                const kk = Math.floor(row / (nx * ny));
+                const rem = row - kk * nx * ny;
+                const jj = Math.floor(rem / nx);
+                const ii = rem - jj * nx;
+                let d = 0;
+                if (ii > 0 && cellType[idxFn(ii - 1, jj, kk)] !== 2)
+                    d += invH2;
+                if (ii < nx - 1 && cellType[idxFn(ii + 1, jj, kk)] !== 2)
+                    d += invH2;
+                if (jj > 0 && cellType[idxFn(ii, jj - 1, kk)] !== 2)
+                    d += invH2;
+                if (jj < ny - 1 && cellType[idxFn(ii, jj + 1, kk)] !== 2)
+                    d += invH2;
+                if (kk > 0 && cellType[idxFn(ii, jj, kk - 1)] !== 2)
+                    d += invH2;
+                if (kk < nz - 1 && cellType[idxFn(ii, jj, kk + 1)] !== 2)
+                    d += invH2;
+                return d > 0 ? d : 1;
+            },
+            apply: (xVec, out) => {
+                out.fill(0);
+                for (let k2 = 0; k2 < nz; k2++)
+                    for (let j2 = 0; j2 < ny; j2++)
+                        for (let i2 = 0; i2 < nx; i2++) {
+                            const ci = idxFn(i2, j2, k2);
+                            if (cellType[ci] !== 1) {
+                                out[ci] = xVec[ci];
+                                continue;
+                            }
+                            let sum = 0, cnt = 0;
+                            if (i2 > 0 && cellType[idxFn(i2 - 1, j2, k2)] !== 2) {
+                                sum += xVec[idxFn(i2 - 1, j2, k2)];
+                                cnt++;
+                            }
+                            if (i2 < nx - 1 && cellType[idxFn(i2 + 1, j2, k2)] !== 2) {
+                                sum += xVec[idxFn(i2 + 1, j2, k2)];
+                                cnt++;
+                            }
+                            if (j2 > 0 && cellType[idxFn(i2, j2 - 1, k2)] !== 2) {
+                                sum += xVec[idxFn(i2, j2 - 1, k2)];
+                                cnt++;
+                            }
+                            if (j2 < ny - 1 && cellType[idxFn(i2, j2 + 1, k2)] !== 2) {
+                                sum += xVec[idxFn(i2, j2 + 1, k2)];
+                                cnt++;
+                            }
+                            if (k2 > 0 && cellType[idxFn(i2, j2, k2 - 1)] !== 2) {
+                                sum += xVec[idxFn(i2, j2, k2 - 1)];
+                                cnt++;
+                            }
+                            if (k2 < nz - 1 && cellType[idxFn(i2, j2, k2 + 1)] !== 2) {
+                                sum += xVec[idxFn(i2, j2, k2 + 1)];
+                                cnt++;
+                            }
+                            out[ci] = invH2 * (cnt * xVec[ci] - sum);
+                        }
+            },
+        };
+        try {
+            pcg(A, this._b, this._x, this.parameters.solverTolerance, Math.round(this.parameters.solverMaxIter));
+        }
+        catch (_) { /* degenerate — keep zero pressure */ }
+        // Guard against NaN (degenerate geometry — no fluid cells, etc.)
+        for (let i = 0; i < this._x.length; i++) {
+            if (!isFinite(this._x[i])) {
+                this._x.fill(0);
+                break;
+            }
+        }
+        // Subtract pressure gradient — loop over faces (each face updated exactly once)
+        // U-faces: (i, j, k) for i in [1, nx-1]
+        for (let k = 0; k < nz; k++)
+            for (let j = 0; j < ny; j++)
+                for (let i = 1; i < nx; i++) {
+                    const cL = grid.idx(i - 1, j, k);
+                    const cR = grid.idx(i, j, k);
+                    if ((cellType[cL] === 1 || cellType[cR] === 1) &&
+                        cellType[cL] !== 2 && cellType[cR] !== 2) {
+                        grid.u[grid.uIdx(i, j, k)] -= dt * (this._x[cR] - this._x[cL]) / h;
+                    }
+                }
+        // V-faces: (i, j, k) for j in [1, ny-1]
+        for (let k = 0; k < nz; k++)
+            for (let j = 1; j < ny; j++)
+                for (let i = 0; i < nx; i++) {
+                    const cB = grid.idx(i, j - 1, k);
+                    const cT = grid.idx(i, j, k);
+                    if ((cellType[cB] === 1 || cellType[cT] === 1) &&
+                        cellType[cB] !== 2 && cellType[cT] !== 2) {
+                        grid.v[grid.vIdx(i, j, k)] -= dt * (this._x[cT] - this._x[cB]) / h;
+                    }
+                }
+        // W-faces: (i, j, k) for k in [1, nz-1]
+        for (let k = 1; k < nz; k++)
+            for (let j = 0; j < ny; j++)
+                for (let i = 0; i < nx; i++) {
+                    const cN = grid.idx(i, j, k - 1);
+                    const cF = grid.idx(i, j, k);
+                    if ((cellType[cN] === 1 || cellType[cF] === 1) &&
+                        cellType[cN] !== 2 && cellType[cF] !== 2) {
+                        grid.w[grid.wIdx(i, j, k)] -= dt * (this._x[cF] - this._x[cN]) / h;
+                    }
+                }
+        grid.enforceBoundary();
+    }
+    // ── G2P (grid → particle) ─────────────────────────────────────────────────
+    _g2p() {
+        const grid = this._grid;
+        const ratio = this.parameters.flipRatio;
+        for (const p of this._particles) {
+            // New grid velocity at particle position
+            const nu = grid.interpU(p.x, p.y, p.z);
+            const nv = grid.interpV(p.x, p.y, p.z);
+            const nw = grid.interpW(p.x, p.y, p.z);
+            // Old grid velocity at particle position (before pressure solve)
+            const ou = this._interpOld(grid.uOld, p.x - grid.ox, p.y - grid.oy - 0.5 * grid.h, p.z - grid.oz - 0.5 * grid.h, grid.nx + 1, grid.ny, grid.nz);
+            const ov = this._interpOld(grid.vOld, p.x - grid.ox - 0.5 * grid.h, p.y - grid.oy, p.z - grid.oz - 0.5 * grid.h, grid.nx, grid.ny + 1, grid.nz);
+            const ow = this._interpOld(grid.wOld, p.x - grid.ox - 0.5 * grid.h, p.y - grid.oy - 0.5 * grid.h, p.z - grid.oz, grid.nx, grid.ny, grid.nz + 1);
+            // FLIP: particle velocity += (new - old) grid velocity
+            // PIC:  particle velocity  = new grid velocity
+            // Blend
+            p.vx = ratio * (p.vx + nu - ou) + (1 - ratio) * nu;
+            p.vy = ratio * (p.vy + nv - ov) + (1 - ratio) * nv;
+            p.vz = ratio * (p.vz + nw - ow) + (1 - ratio) * nw;
+        }
+    }
+    _interpOld(arr, lx, ly, lz, Nx, Ny, Nz) {
+        const h = this._grid.h;
+        let fi = lx / h, fj = ly / h, fk = lz / h;
+        fi = Math.max(0, Math.min(Nx - 1.0001, fi));
+        fj = Math.max(0, Math.min(Ny - 1.0001, fj));
+        fk = Math.max(0, Math.min(Nz - 1.0001, fk));
+        const i0 = Math.floor(fi), i1 = i0 + 1;
+        const j0 = Math.floor(fj), j1 = j0 + 1;
+        const k0 = Math.floor(fk), k1 = k0 + 1;
+        const tx = fi - i0, ty = fj - j0, tz = fk - k0;
+        const idx = (i, j, k) => i + Nx * (j + Ny * k);
+        return (1 - tz) * ((1 - ty) * ((1 - tx) * arr[idx(i0, j0, k0)] + tx * arr[idx(i1, j0, k0)]) +
+            ty * ((1 - tx) * arr[idx(i0, j1, k0)] + tx * arr[idx(i1, j1, k0)])) +
+            tz * ((1 - ty) * ((1 - tx) * arr[idx(i0, j0, k1)] + tx * arr[idx(i1, j0, k1)]) +
+                ty * ((1 - tx) * arr[idx(i0, j1, k1)] + tx * arr[idx(i1, j1, k1)]));
+    }
+    // ── Advection ─────────────────────────────────────────────────────────────
+    _advect(dt) {
+        for (const p of this._particles) {
+            p.x += p.vx * dt;
+            p.y += p.vy * dt;
+            p.z += p.vz * dt;
+        }
+    }
+    _clampParticles() {
+        const { ox, oy, oz, nx, ny, nz, h } = this._grid;
+        const margin = h * 0.5;
+        const xMin = ox + margin, xMax = ox + nx * h - margin;
+        const yMin = oy + margin, yMax = oy + ny * h - margin;
+        const zMin = oz + margin, zMax = oz + nz * h - margin;
+        for (const p of this._particles) {
+            if (p.x < xMin) {
+                p.x = xMin;
+                if (p.vx < 0)
+                    p.vx = 0;
+            }
+            if (p.x > xMax) {
+                p.x = xMax;
+                if (p.vx > 0)
+                    p.vx = 0;
+            }
+            if (p.y < yMin) {
+                p.y = yMin;
+                if (p.vy < 0)
+                    p.vy = 0;
+            }
+            if (p.y > yMax) {
+                p.y = yMax;
+                if (p.vy > 0)
+                    p.vy = 0;
+            }
+            if (p.z < zMin) {
+                p.z = zMin;
+                if (p.vz < 0)
+                    p.vz = 0;
+            }
+            if (p.z > zMax) {
+                p.z = zMax;
+                if (p.vz > 0)
+                    p.vz = 0;
+            }
+        }
+    }
+    // ── Output ────────────────────────────────────────────────────────────────
+    /** Returns particle positions as a flat Float32Array [x,y,z, x,y,z, ...]. */
+    getPositions() {
+        const out = new Float32Array(this._particles.length * 3);
+        for (let i = 0; i < this._particles.length; i++) {
+            out[i * 3] = this._particles[i].x;
+            out[i * 3 + 1] = this._particles[i].y;
+            out[i * 3 + 2] = this._particles[i].z;
+        }
+        return out;
+    }
+    /** Returns particle velocities as a flat Float32Array [vx,vy,vz, ...]. */
+    getVelocities() {
+        const out = new Float32Array(this._particles.length * 3);
+        for (let i = 0; i < this._particles.length; i++) {
+            out[i * 3] = this._particles[i].vx;
+            out[i * 3 + 1] = this._particles[i].vy;
+            out[i * 3 + 2] = this._particles[i].vz;
+        }
+        return out;
+    }
+    /** Number of live particles. */
+    get liveCount() { return this._particles.length; }
+    /** The underlying MAC grid (for advanced use / surface reconstruction). */
+    get grid() { return this._grid; }
+}
+//# sourceMappingURL=FLIPSimulator.js.map

package/dist/FLIPSimulator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"FLIPSimulator.js","sourceRoot":"","sources":["../src/FLIPSimulator.ts"],"names":[],"mappings":"AAAA,OAAO,GAAG,MAAM,oBAAoB,CAAA;AACpC,OAAO,EAAE,OAAO,EAAE,MAAM,cAAc,CAAA;AAsDtC;;;;;;;;;;;;;;;;;;;;;;GAsBG;AACH,MAAM,OAAO,aAAa;IAcxB,YAAY,OAA6B,EAAE;QACzC,MAAM,GAAG,GAAU,IAAI,CAAC,UAAU,IAAO,EAAE,CAAA;QAC3C,MAAM,UAAU,GAAG,IAAI,CAAC,UAAU,IAAO,CAAC,CAAA;QAC1C,MAAM,MAAM,GAAO,IAAI,CAAC,MAAM,IAAW,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAA;QACrD,MAAM,CAAC,GAAY,UAAU,GAAG,GAAG,CAAA;QAEnC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAA;QAC9C,IAAI,CAAC,KAAK,GAAG,IAAI,OAAO,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAA;QAE3E,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,aAAa,IAAI,IAAI,CAAA;QAC/C,IAAI,CAAC,UAAU,GAAM,EAAE,CAAA;QAEvB,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAA;QAC3C,IAAI,CAAC,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAA;QAE3C,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAA;QAEtC,IAAI,CAAC,UAAU,GAAG;YAChB,SAAS,EAAQ,IAAI,CAAC,SAAS,IAAU,IAAI;YAC7C,QAAQ,EAAS,CAAC,CAAC,CAAC,CAAC;YACrB,QAAQ,EAAS,CAAC,CAAC,CAAC,CAAC;YACrB,QAAQ,EAAS,CAAC,CAAC,CAAC,CAAC;YACrB,QAAQ,EAAS,IAAI,CAAC,QAAQ,IAAW,CAAC;YAC1C,aAAa,EAAI,IAAI,CAAC,aAAa,IAAM,EAAE;YAC3C,eAAe,EAAE,IAAI,CAAC,eAAe,IAAI,IAAI;SAC9C,CAAA;IACH,CAAC;IAED,6EAA6E;IAE7E;;;;OAIG;IACH,OAAO,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAClC,EAAU,EAAE,EAAU,EAAE,EAAU;QACxC,IAAI,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAA;QAC1B,MAAM,CAAC,GAAK,IAAI,CAAC,KAAK,CAAC,CAAC,CAAA;QACxB,MAAM,GAAG,GAAG,CAAC,CAAA,CAAG,+DAA+D;QAC/E,MAAM,EAAE,GAAI,CAAC,GAAG,GAAG,CAAA;QAEnB,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE;YAChC,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE;gBAChC,KAAK,IAAI,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC;oBACjC,IAAI,IAAI,CAAC,UAAU,CAAC,MAAM,IAAI,IAAI,CAAC,aAAa;wBAAE,MAAK;oBACvD,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC;wBACnB,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,MAAM,EAAE;wBACzB,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,MAAM,EAAE;wBACzB,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,IAAI,CAAC,MAAM,EAAE;wBACzB,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC;qBACpB,CAAC,CAAA;gBACJ,CAAC;IACH,CAAC;IAED,6EAA6E;IAE7E,kFAAkF;IAClF,IAAI,CAAC,EAAU;QACb,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,UAAU,CAAC,QAAQ,CAAC,CAAC,CAAA;QAClE,MAAM,KAAK,GAAM,EAAE,GAAG,QAAQ,CAAA;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAA;IACzD,CAAC;IAEO,QAAQ,CAAC,EAAU;QACzB,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAA;QACvB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,IAAI,CAAA;QAE9B,2DAA2D;QAC3D,MAAM,QAAQ,GAAG,IAAI,GAAG,EAAU,CAAA;QAClC,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YAChC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAA;YACzC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAA;YACzC,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAA;YACzC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,EAAE,GAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,EAAE,GAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,EAAE,GAAC,CAAC;gBAChE,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAA;QACnC,CAAC;QACD,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAA;QAE5D,mDAAmD;QACnD,IAAI,CAAC,eAAe,EAAE,CAAA;QACtB,IAAI,CAAC,IAAI,EAAE,CAAA;QACX,IAAI,CAAC,gBAAgB,EAAE,CAAA;QACvB,IAAI,CAAC,cAAc,EAAE,CAAA;QAErB,4BAA4B;QAC5B,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,UAAU,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,CAAC,QAAQ,EAAE,IAAI,CAAC,UAAU,CAAC,QAAQ,EAAE,EAAE,CAAC,CAAA;QACnG,IAAI,CAAC,eAAe,EAAE,CAAA;QAEtB,yBAAyB;QACzB,IAAI,CAAC,cAAc,CAAC,EAAE,CAAC,CAAA;QAEvB,6EAA6E;QAC7E,IAAI,CAAC,IAAI,EAAE,CAAA;QAEX,sBAAsB;QACtB,IAAI,CAAC,OAAO,CAAC,EAAE,CAAC,CAAA;QAEhB,mCAAmC;QACnC,IAAI,CAAC,eAAe,EAAE,CAAA;IACxB,CAAC;IAED,6EAA6E;IAErE,IAAI;QACV,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAA;QACvB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,IAAI,CAAA;QAE1C,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YAChC,gDAAgD;YAChD,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EACrC,CAAC,CAAC,CAAC,GAAG,EAAE,EAAU,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAC,CAAC,EACpD,EAAE,GAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAA;YAElB,gDAAgD;YAChD,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EACrC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,EAAU,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAC,CAAC,EACpD,EAAE,EAAE,EAAE,GAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAA;YAElB,gDAAgD;YAChD,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EACrC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,EAC5C,EAAE,EAAE,EAAE,EAAE,EAAE,GAAC,CAAC,EAAE,CAAC,CAAC,CAAA;QACpB,CAAC;IACH,CAAC;IAEO,YAAY,CAAC,GAAiB,EAAE,IAAkB,EAAE,GAAW,EAClD,EAAU,EAAE,EAAU,EAAE,EAAU,EAClC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,CAAS;QAChE,IAAI,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAA;QACzC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,CAAC,CAAA;QAC3C,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,CAAC,CAAA;QAC3C,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,CAAC,CAAA;QAE3C,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAA;QACtC,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAA;QACtC,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAA;QACtC,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,CAAA;QAC9C,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAA;QAEtE,MAAM,IAAI,GAAG,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,GAAC,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAA;QAC9F,MAAM,IAAI,GAAM,EAAE,GAAE,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,GAAC,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAA;QAC9F,MAAM,IAAI,GAAG,CAAC,CAAC,GAAC,EAAE,CAAC,GAAI,EAAE,GAAE,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,GAAC,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAA;QAC9F,MAAM,IAAI,GAAM,EAAE,GAAK,EAAE,GAAE,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC;QAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,GAAC,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAA;QAC9F,MAAM,IAAI,GAAG,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAI,EAAE,CAAE;QAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,GAAC,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAA;QAC9F,MAAM,IAAI,GAAM,EAAE,GAAE,CAAC,CAAC,GAAC,EAAE,CAAC,GAAI,EAAE,CAAE;QAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,GAAC,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAA;QAC9F,MAAM,IAAI,GAAG,CAAC,CAAC,GAAC,EAAE,CAAC,GAAI,EAAE,GAAK,EAAE,CAAE;QAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,GAAC,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAA;QAC9F,MAAM,IAAI,GAAM,EAAE,GAAK,EAAE,GAAK,EAAE,CAAE;QAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,GAAC,GAAG,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,IAAI,IAAI,CAAA;IAChG,CAAC;IAED,6EAA6E;IAErE,cAAc,CAAC,EAAU;QAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAA;QACvB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,IAAI,CAAA;QAC9B,MAAM,CAAC,GAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAA;QACzB,MAAM,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QAE3B,yEAAyE;QACzE,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;QACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,GAAC,CAAC,EAAE,CAAC,EAAE;YAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,GAAC,CAAC,EAAE,CAAC,EAAE;gBAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,GAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAA;oBAC5B,IAAI,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC;wBAAE,SAAQ;oBACrC,MAAM,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;0BAC3D,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAC,CAAC,GAAC,CAAC,EAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;0BAC3D,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,GAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAA;oBACvE,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,EAAE,CAAA;gBACxB,CAAC;QAED,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;QAEf,MAAM,QAAQ,GAAG,IAAI,CAAC,QAAQ,CAAA;QAC9B,MAAM,KAAK,GAAM,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAA;QAEvE,sEAAsE;QACtE,MAAM,CAAC,GAAG;YACR,QAAQ,EAAE,CAAC;YACX,GAAG,EAAE,CAAC,GAAW,EAAE,IAAY,EAAU,EAAE;gBACzC,IAAI,QAAQ,CAAC,GAAG,CAAC,KAAK,CAAC;oBAAE,OAAO,CAAC,CAAA;gBACjC,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAA;gBACtC,MAAM,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAA;gBAC9B,MAAM,EAAE,GAAI,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,EAAE,CAAC,CAAA;gBAChC,MAAM,EAAE,GAAI,GAAG,GAAG,EAAE,GAAG,EAAE,CAAA;gBACzB,IAAI,CAAC,GAAG,CAAC,CAAA;gBACT,IAAI,EAAE,GAAG,CAAC,IAAO,QAAQ,CAAC,KAAK,CAAC,EAAE,GAAC,CAAC,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,KAAK,CAAC;oBAAE,CAAC,IAAI,KAAK,CAAA;gBAC9D,IAAI,EAAE,GAAG,EAAE,GAAC,CAAC,IAAI,QAAQ,CAAC,KAAK,CAAC,EAAE,GAAC,CAAC,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,KAAK,CAAC;oBAAE,CAAC,IAAI,KAAK,CAAA;gBAC9D,IAAI,EAAE,GAAG,CAAC,IAAO,QAAQ,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,GAAC,CAAC,EAAC,EAAE,CAAC,CAAC,KAAK,CAAC;oBAAE,CAAC,IAAI,KAAK,CAAA;gBAC9D,IAAI,EAAE,GAAG,EAAE,GAAC,CAAC,IAAI,QAAQ,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,GAAC,CAAC,EAAC,EAAE,CAAC,CAAC,KAAK,CAAC;oBAAE,CAAC,IAAI,KAAK,CAAA;gBAC9D,IAAI,EAAE,GAAG,CAAC,IAAO,QAAQ,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,GAAC,CAAC,CAAC,CAAC,KAAK,CAAC;oBAAE,CAAC,IAAI,KAAK,CAAA;gBAC9D,IAAI,EAAE,GAAG,EAAE,GAAC,CAAC,IAAI,QAAQ,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,GAAC,CAAC,CAAC,CAAC,KAAK,CAAC;oBAAE,CAAC,IAAI,KAAK,CAAA;gBAC9D,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YACtB,CAAC;YACD,KAAK,EAAE,CAAC,IAAkB,EAAE,GAAiB,EAAQ,EAAE;gBACrD,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;gBACX,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE;oBAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE;wBAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;4BAC/B,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAA;4BAC5B,IAAI,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC;gCAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;gCAAC,SAAQ;4BAAC,CAAC;4BACxD,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAA;4BACpB,IAAI,EAAE,GAAG,CAAC,IAAO,QAAQ,CAAC,KAAK,CAAC,EAAE,GAAC,CAAC,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;gCAAC,GAAG,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,GAAC,CAAC,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;gCAAC,GAAG,EAAE,CAAA;4BAAC,CAAC;4BAC7F,IAAI,EAAE,GAAG,EAAE,GAAC,CAAC,IAAI,QAAQ,CAAC,KAAK,CAAC,EAAE,GAAC,CAAC,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;gCAAC,GAAG,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,GAAC,CAAC,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;gCAAC,GAAG,EAAE,CAAA;4BAAC,CAAC;4BAC7F,IAAI,EAAE,GAAG,CAAC,IAAO,QAAQ,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,GAAC,CAAC,EAAC,EAAE,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;gCAAC,GAAG,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,GAAC,CAAC,EAAC,EAAE,CAAC,CAAC,CAAC;gCAAC,GAAG,EAAE,CAAA;4BAAC,CAAC;4BAC7F,IAAI,EAAE,GAAG,EAAE,GAAC,CAAC,IAAI,QAAQ,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,GAAC,CAAC,EAAC,EAAE,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;gCAAC,GAAG,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,GAAC,CAAC,EAAC,EAAE,CAAC,CAAC,CAAC;gCAAC,GAAG,EAAE,CAAA;4BAAC,CAAC;4BAC7F,IAAI,EAAE,GAAG,CAAC,IAAO,QAAQ,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,GAAC,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;gCAAC,GAAG,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,GAAC,CAAC,CAAC,CAAC,CAAC;gCAAC,GAAG,EAAE,CAAA;4BAAC,CAAC;4BAC7F,IAAI,EAAE,GAAG,EAAE,GAAC,CAAC,IAAI,QAAQ,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,GAAC,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;gCAAC,GAAG,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,GAAC,CAAC,CAAC,CAAC,CAAC;gCAAC,GAAG,EAAE,CAAA;4BAAC,CAAC;4BAC7F,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAA;wBAC1C,CAAC;YACH,CAAC;SACF,CAAA;QAED,IAAI,CAAC;YACH,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,UAAU,CAAC,eAAe,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,UAAU,CAAC,aAAa,CAAC,CAAC,CAAA;QACtG,CAAC;QAAC,OAAO,CAAC,EAAE,CAAC,CAAC,qCAAqC,CAAC,CAAC;QAErD,iEAAiE;QACjE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACxC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;gBAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gBAAC,MAAK;YAAC,CAAC;QACvD,CAAC;QAED,gFAAgF;QAChF,wCAAwC;QACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAA;oBAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAI,CAAC,EAAE,CAAC,CAAC,CAAA;oBAC9B,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;wBACzC,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC;wBAC9C,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAA;oBAClE,CAAC;gBACH,CAAC;QAED,wCAAwC;QACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAC,CAAC,EAAE,CAAC,CAAC,CAAA;oBAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAI,CAAC,CAAC,CAAA;oBAC9B,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;wBACzC,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC;wBAC9C,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAA;oBAClE,CAAC;gBACH,CAAC;QAED,wCAAwC;QACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;gBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,GAAC,CAAC,CAAC,CAAA;oBAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAG,CAAA;oBAC9B,IAAI,CAAC,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;wBACzC,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,QAAQ,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC;wBAC9C,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAA;oBAClE,CAAC;gBACH,CAAC;QAED,IAAI,CAAC,eAAe,EAAE,CAAA;IACxB,CAAC;IAED,6EAA6E;IAErE,IAAI;QACV,MAAM,IAAI,GAAI,IAAI,CAAC,KAAK,CAAA;QACxB,MAAM,KAAK,GAAG,IAAI,CAAC,UAAU,CAAC,SAAS,CAAA;QAEvC,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YAChC,yCAAyC;YACzC,MAAM,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAA;YACtC,MAAM,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAA;YACtC,MAAM,EAAE,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAA;YAEtC,iEAAiE;YACjE,MAAM,EAAE,GAAG,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,EAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,GAAC,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAA;YACjJ,MAAM,EAAE,GAAG,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,EAAU,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,GAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAA;YACjJ,MAAM,EAAE,GAAG,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,EAAU,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,GAAC,CAAC,CAAC,CAAA;YAEjJ,uDAAuD;YACvD,+CAA+C;YAC/C,QAAQ;YACR,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,EAAE,CAAA;YAClD,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,EAAE,CAAA;YAClD,CAAC,CAAC,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,EAAE,CAAA;QACpD,CAAC;IACH,CAAC;IAEO,UAAU,CAAC,GAAiB,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EACrD,EAAU,EAAE,EAAU,EAAE,EAAU;QACnD,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAA;QACtB,IAAI,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAA;QACzC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,CAAC,CAAA;QAC3C,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,CAAC,CAAA;QAC3C,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,MAAM,EAAE,EAAE,CAAC,CAAC,CAAA;QAC3C,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAA;QACtC,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAA;QACtC,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC,CAAA;QACtC,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,CAAA;QAC9C,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAA;QACtE,OAAO,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,GAAG,EAAE,GAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;YACvD,EAAE,GAAE,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,GAAG,EAAE,GAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC/D,EAAE,GAAC,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,GAAG,EAAE,GAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC;gBACvD,EAAE,GAAE,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,GAAG,EAAE,GAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAC,EAAE,EAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAA;IAC5E,CAAC;IAED,6EAA6E;IAErE,OAAO,CAAC,EAAU;QACxB,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YAChC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,EAAE,CAAA;YAChB,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,EAAE,CAAA;YAChB,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,EAAE,CAAA;QAClB,CAAC;IACH,CAAC;IAEO,eAAe;QACrB,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAA;QAChD,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,CAAA;QACtB,MAAM,IAAI,GAAG,EAAE,GAAG,MAAM,EAAE,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,MAAM,CAAA;QACrD,MAAM,IAAI,GAAG,EAAE,GAAG,MAAM,EAAE,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,MAAM,CAAA;QACrD,MAAM,IAAI,GAAG,EAAE,GAAG,MAAM,EAAE,IAAI,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,GAAG,MAAM,CAAA;QAErD,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,UAAU,EAAE,CAAC;YAChC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC;oBAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAA;YAAC,CAAC;YACtD,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC;oBAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAA;YAAC,CAAC;YACtD,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC;oBAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAA;YAAC,CAAC;YACtD,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC;oBAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAA;YAAC,CAAC;YACtD,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC;oBAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAA;YAAC,CAAC;YACtD,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC;oBAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAA;YAAC,CAAC;QACxD,CAAC;IACH,CAAC;IAED,6EAA6E;IAE7E,6EAA6E;IAC7E,YAAY;QACV,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;QACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAChD,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,GAAK,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YACjC,GAAG,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YACjC,GAAG,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QACnC,CAAC;QACD,OAAO,GAAG,CAAA;IACZ,CAAC;IAED,0EAA0E;IAC1E,aAAa;QACX,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;QACxD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAChD,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,GAAK,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAA;YAClC,GAAG,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAA;YAClC,GAAG,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAA;QACpC,CAAC;QACD,OAAO,GAAG,CAAA;IACZ,CAAC;IAED,gCAAgC;IAChC,IAAI,SAAS,KAAa,OAAO,IAAI,CAAC,UAAU,CAAC,MAAM,CAAA,CAAC,CAAC;IAEzD,2EAA2E;IAC3E,IAAI,IAAI,KAAc,OAAO,IAAI,CAAC,KAAK,CAAA,CAAC,CAAC;CAC1C"}

package/dist/FluidEmitter.d.ts

@@ -0,0 +1,35 @@
+/**
+ * FluidEmitter — defines a region that spawns fluid particles.
+ * Blender: Fluid > Emitter object.
+ *
+ * Used to spawn batches of particles from SPHSimulator.spawnFrom(emitter).
+ */
+export interface FluidEmitterOptions {
+    /** World-space position [x,y,z]. Default [0,0,0]. */
+    position?: [number, number, number];
+    /** Spawn radius. Default 0.5. */
+    radius?: number;
+    /** Initial velocity of spawned particles [x,y,z]. Default [0,0,0]. */
+    velocity?: [number, number, number];
+    /** Jitter applied to initial velocity. Default 0.1. */
+    jitter?: number;
+}
+export declare class FluidEmitter {
+    parameters: {
+        posX: number;
+        posY: number;
+        posZ: number;
+        radius: number;
+        velX: number;
+        velY: number;
+        velZ: number;
+        jitter: number;
+    };
+    enabled: boolean;
+    constructor(opts?: FluidEmitterOptions);
+    /** Sample a random spawn position within the emitter sphere. */
+    samplePosition(): [number, number, number];
+    /** Sample initial velocity with jitter. */
+    sampleVelocity(): [number, number, number];
+}
+//# sourceMappingURL=FluidEmitter.d.ts.map

package/dist/FluidEmitter.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"FluidEmitter.d.ts","sourceRoot":"","sources":["../src/FluidEmitter.ts"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,MAAM,WAAW,mBAAmB;IAClC,qDAAqD;IACrD,QAAQ,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAA;IACnC,iCAAiC;IACjC,MAAM,CAAC,EAAE,MAAM,CAAA;IACf,sEAAsE;IACtE,QAAQ,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAA;IACnC,uDAAuD;IACvD,MAAM,CAAC,EAAE,MAAM,CAAA;CAChB;AAED,qBAAa,YAAY;IACvB,UAAU,EAAE;QACV,IAAI,EAAM,MAAM,CAAA;QAChB,IAAI,EAAM,MAAM,CAAA;QAChB,IAAI,EAAM,MAAM,CAAA;QAChB,MAAM,EAAI,MAAM,CAAA;QAChB,IAAI,EAAM,MAAM,CAAA;QAChB,IAAI,EAAM,MAAM,CAAA;QAChB,IAAI,EAAM,MAAM,CAAA;QAChB,MAAM,EAAI,MAAM,CAAA;KACjB,CAAA;IAED,OAAO,UAAO;gBAEF,IAAI,GAAE,mBAAwB;IAW1C,gEAAgE;IAChE,cAAc,IAAI,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC;IAY1C,2CAA2C;IAC3C,cAAc,IAAI,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC;CAQ3C"}

package/dist/FluidEmitter.js

@@ -0,0 +1,35 @@
+export class FluidEmitter {
+    constructor(opts = {}) {
+        this.enabled = true;
+        const p = opts.position ?? [0, 0, 0];
+        const v = opts.velocity ?? [0, 0, 0];
+        this.parameters = {
+            posX: p[0], posY: p[1], posZ: p[2],
+            radius: opts.radius ?? 0.5,
+            velX: v[0], velY: v[1], velZ: v[2],
+            jitter: opts.jitter ?? 0.1,
+        };
+    }
+    /** Sample a random spawn position within the emitter sphere. */
+    samplePosition() {
+        const r = this.parameters.radius * Math.cbrt(Math.random());
+        const theta = Math.random() * Math.PI * 2;
+        const phi = Math.acos(2 * Math.random() - 1);
+        const sinPhi = Math.sin(phi);
+        return [
+            this.parameters.posX + r * sinPhi * Math.cos(theta),
+            this.parameters.posY + r * sinPhi * Math.sin(theta),
+            this.parameters.posZ + r * Math.cos(phi),
+        ];
+    }
+    /** Sample initial velocity with jitter. */
+    sampleVelocity() {
+        const j = this.parameters.jitter;
+        return [
+            this.parameters.velX + (Math.random() - 0.5) * j * 2,
+            this.parameters.velY + (Math.random() - 0.5) * j * 2,
+            this.parameters.velZ + (Math.random() - 0.5) * j * 2,
+        ];
+    }
+}
+//# sourceMappingURL=FluidEmitter.js.map

package/dist/FluidEmitter.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"FluidEmitter.js","sourceRoot":"","sources":["../src/FluidEmitter.ts"],"names":[],"mappings":"AAiBA,MAAM,OAAO,YAAY;IAcvB,YAAY,OAA4B,EAAE;QAF1C,YAAO,GAAG,IAAI,CAAA;QAGZ,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAA;QACpC,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAA;QACpC,IAAI,CAAC,UAAU,GAAG;YAChB,IAAI,EAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC;YACpC,MAAM,EAAE,IAAI,CAAC,MAAM,IAAI,GAAG;YAC1B,IAAI,EAAI,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC;YACpC,MAAM,EAAE,IAAI,CAAC,MAAM,IAAI,GAAG;SAC3B,CAAA;IACH,CAAC;IAED,gEAAgE;IAChE,cAAc;QACZ,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,CAAA;QAC3D,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAA;QACzC,MAAM,GAAG,GAAK,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAA;QAC9C,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;QAC5B,OAAO;YACL,IAAI,CAAC,UAAU,CAAC,IAAI,GAAG,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;YACnD,IAAI,CAAC,UAAU,CAAC,IAAI,GAAG,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;YACnD,IAAI,CAAC,UAAU,CAAC,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC;SACzC,CAAA;IACH,CAAC;IAED,2CAA2C;IAC3C,cAAc;QACZ,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,MAAM,CAAA;QAChC,OAAO;YACL,IAAI,CAAC,UAAU,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC;YACpD,IAAI,CAAC,UAAU,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC;YACpD,IAAI,CAAC,UAAU,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,MAAM,EAAE,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC;SACrD,CAAA;IACH,CAAC;CACF"}