@triforge/fluid-core 0.1.0

package/dist/SPHSimulator.d.ts

@@ -0,0 +1,139 @@
+import { BufferGeometry } from 'three';
+export interface SPHSimulatorOptions {
+    /** Number of fluid particles. Default 500. */
+    particleCount?: number;
+    /** Smoothing radius (h) — neighbourhood size. Blender: Fluid > Resolution. Default 0.3. */
+    smoothingRadius?: number;
+    /** Target rest density (ρ₀). Default 1000. */
+    restDensity?: number;
+    /** Pressure stiffness constant k. Higher = stiffer fluid. Default 200. */
+    pressureK?: number;
+    /** Viscosity coefficient μ. Blender: Fluid > Viscosity. Default 0.1. */
+    viscosity?: number;
+    /**
+     * Viscosity model.
+     * - `'standard'` — Müller 2003 viscosity Laplacian kernel (default)
+     * - `'xsph'`     — XSPH correction (Schechter & Bridson 2012): blends each
+     *                  particle's velocity toward its neighbours. More stable,
+     *                  less physically accurate. Good for thick fluids.
+     * - `'stiff'`    — Position-based viscosity: applies an extra velocity
+     *                  correction proportional to velocity divergence. Very stiff,
+     *                  best for near-incompressible high-viscosity fluids.
+     * Blender: Fluid > Viscosity Preset (thin/medium/thick water / honey / oil).
+     */
+    viscosityType?: 'standard' | 'xsph' | 'stiff';
+    /** Gravity [x,y,z]. Default [0,-9.8,0]. */
+    gravity?: [number, number, number];
+    /** Simulation substeps per frame. Default 4. */
+    substeps?: number;
+    /** Surface tension coefficient. Default 0.0728 (water). */
+    surfaceTension?: number;
+}
+/**
+ * SPHSimulator — Smoothed Particle Hydrodynamics fluid simulation.
+ * Blender: Fluid modifier (FLIP/SPH domain).
+ *
+ * Implements the Müller et al. 2003 SPH formulation:
+ *   - Poly6 kernel for density
+ *   - Spiky kernel gradient for pressure
+ *   - Viscosity kernel Laplacian for viscosity
+ *
+ * Workflow:
+ *   1. Construct with particle count + physics parameters
+ *   2. Call `setPositions(array)` to place particles in a volume
+ *   3. Call `step(dt)` each frame (dt ≈ 1/60)
+ *   4. Call `getGeometry()` to read back a THREE.BufferGeometry for rendering
+ *
+ * All physics parameters live in `parameters` for GSAP/st-keyframe:
+ *   gsap.to(sim.parameters, { viscosity: 0.5, duration: 1 })
+ */
+export declare class SPHSimulator {
+    /** All scalar simulation parameters — GSAP/st-keyframe compatible. */
+    parameters: {
+        smoothingRadius: number;
+        restDensity: number;
+        pressureK: number;
+        viscosity: number;
+        /** XSPH blend factor [0–1]. Only used when viscosityType='xsph'. Default 0.5. */
+        xsphFactor: number;
+        /** Stiff viscosity strength. Only used when viscosityType='stiff'. Default 0.5. */
+        stiffFactor: number;
+        gravityX: number;
+        gravityY: number;
+        gravityZ: number;
+        substeps: number;
+        surfaceTension: number;
+    };
+    /** Active viscosity model. Can be changed at runtime. */
+    viscosityType: 'standard' | 'xsph' | 'stiff';
+    readonly particleCount: number;
+    private _px;
+    private _py;
+    private _pz;
+    private _vx;
+    private _vy;
+    private _vz;
+    private _density;
+    private _pressure;
+    private _grid;
+    /** Simulation clock. */
+    time: number;
+    private _domainMin;
+    private _domainMax;
+    constructor(opts?: SPHSimulatorOptions);
+    private _initGrid;
+    /**
+     * Set particle positions from a flat array [x0,y0,z0, x1,y1,z1, ...].
+     */
+    setPositions(arr: ArrayLike<number>): void;
+    /**
+     * Set simulation domain (AABB boundary collider).
+     * Particles bounce off the domain walls.
+     */
+    setDomain(min: [number, number, number], max: [number, number, number]): void;
+    /**
+     * Advance the simulation by dt seconds.
+     * Recommended: pass clock.getDelta() clamped to 0.05 s.
+     */
+    step(dt: number): void;
+    /**
+     * Build and return a BufferGeometry with a 'position' attribute.
+     * Suitable for THREE.Points rendering.
+     */
+    getGeometry(): BufferGeometry;
+    /**
+     * Copy current particle positions into an existing BufferGeometry.
+     * More efficient than getGeometry() in the render loop.
+     */
+    updateGeometry(geo: BufferGeometry): void;
+    /** Per-particle density (read-only, updated after each step). */
+    getDensity(i: number): number;
+    /** Per-particle pressure (read-only). */
+    getPressure(i: number): number;
+    private _substep;
+    private _buildSpatialGrid;
+    private _cellKey;
+    private _neighbours;
+    private _poly6;
+    private _spikyGrad;
+    private _viscLaplacian;
+    private _computeDensityPressure;
+    private _integrate;
+    /**
+     * XSPH viscosity (Schechter & Bridson 2012).
+     * Blends each particle's velocity toward a weighted average of its neighbours.
+     * More stable than Müller viscosity, good for thick/honey-like fluids.
+     * Controlled by parameters.xsphFactor [0–1].
+     */
+    private _applyXSPH;
+    /**
+     * Stiff viscosity (position-based / divergence damping).
+     * Applies a velocity correction proportional to the local velocity divergence,
+     * strongly suppressing compressible motion.
+     * Best for near-incompressible high-viscosity simulations.
+     * Controlled by parameters.stiffFactor.
+     */
+    private _applyStiffViscosity;
+    private _enforceDomain;
+}
+//# sourceMappingURL=SPHSimulator.d.ts.map

package/dist/SPHSimulator.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SPHSimulator.d.ts","sourceRoot":"","sources":["../src/SPHSimulator.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAmB,MAAM,OAAO,CAAA;AAEvD,MAAM,WAAW,mBAAmB;IAClC,8CAA8C;IAC9C,aAAa,CAAC,EAAE,MAAM,CAAA;IACtB,2FAA2F;IAC3F,eAAe,CAAC,EAAE,MAAM,CAAA;IACxB,8CAA8C;IAC9C,WAAW,CAAC,EAAE,MAAM,CAAA;IACpB,0EAA0E;IAC1E,SAAS,CAAC,EAAE,MAAM,CAAA;IAClB,wEAAwE;IACxE,SAAS,CAAC,EAAE,MAAM,CAAA;IAClB;;;;;;;;;;OAUG;IACH,aAAa,CAAC,EAAE,UAAU,GAAG,MAAM,GAAG,OAAO,CAAA;IAC7C,2CAA2C;IAC3C,OAAO,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAA;IAClC,gDAAgD;IAChD,QAAQ,CAAC,EAAE,MAAM,CAAA;IACjB,2DAA2D;IAC3D,cAAc,CAAC,EAAE,MAAM,CAAA;CACxB;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,qBAAa,YAAY;IACvB,sEAAsE;IACtE,UAAU,EAAE;QACV,eAAe,EAAE,MAAM,CAAA;QACvB,WAAW,EAAM,MAAM,CAAA;QACvB,SAAS,EAAQ,MAAM,CAAA;QACvB,SAAS,EAAQ,MAAM,CAAA;QACvB,iFAAiF;QACjF,UAAU,EAAO,MAAM,CAAA;QACvB,mFAAmF;QACnF,WAAW,EAAM,MAAM,CAAA;QACvB,QAAQ,EAAS,MAAM,CAAA;QACvB,QAAQ,EAAS,MAAM,CAAA;QACvB,QAAQ,EAAS,MAAM,CAAA;QACvB,QAAQ,EAAS,MAAM,CAAA;QACvB,cAAc,EAAG,MAAM,CAAA;KACxB,CAAA;IAED,yDAAyD;IACzD,aAAa,EAAE,UAAU,GAAG,MAAM,GAAG,OAAO,CAAA;IAE5C,QAAQ,CAAC,aAAa,EAAE,MAAM,CAAA;IAG9B,OAAO,CAAC,GAAG,CAAc;IACzB,OAAO,CAAC,GAAG,CAAc;IACzB,OAAO,CAAC,GAAG,CAAc;IACzB,OAAO,CAAC,GAAG,CAAc;IACzB,OAAO,CAAC,GAAG,CAAc;IACzB,OAAO,CAAC,GAAG,CAAc;IACzB,OAAO,CAAC,QAAQ,CAAc;IAC9B,OAAO,CAAC,SAAS,CAAc;IAG/B,OAAO,CAAC,KAAK,CAAmC;IAEhD,wBAAwB;IACxB,IAAI,SAAI;IAGR,OAAO,CAAC,UAAU,CAAuC;IACzD,OAAO,CAAC,UAAU,CAAuC;gBAE7C,IAAI,GAAE,mBAAwB;IAqC1C,OAAO,CAAC,SAAS;IAkBjB;;OAEG;IACH,YAAY,CAAC,GAAG,EAAE,SAAS,CAAC,MAAM,CAAC,GAAG,IAAI;IAU1C;;;OAGG;IACH,SAAS,CACP,GAAG,EAAE,CAAC,MAAM,EAAC,MAAM,EAAC,MAAM,CAAC,EAC3B,GAAG,EAAE,CAAC,MAAM,EAAC,MAAM,EAAC,MAAM,CAAC,GAC1B,IAAI;IAOP;;;OAGG;IACH,IAAI,CAAC,EAAE,EAAE,MAAM,GAAG,IAAI;IAOtB;;;OAGG;IACH,WAAW,IAAI,cAAc;IAY7B;;;OAGG;IACH,cAAc,CAAC,GAAG,EAAE,cAAc,GAAG,IAAI;IASzC,iEAAiE;IACjE,UAAU,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAE7B,yCAAyC;IACzC,WAAW,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM;IAI9B,OAAO,CAAC,QAAQ;IAUhB,OAAO,CAAC,iBAAiB;IAezB,OAAO,CAAC,QAAQ;IAKhB,OAAO,CAAC,WAAW;IA0BnB,OAAO,CAAC,MAAM;IASd,OAAO,CAAC,UAAU;IAQlB,OAAO,CAAC,cAAc;IAMtB,OAAO,CAAC,uBAAuB;IAe/B,OAAO,CAAC,UAAU;IAqDlB;;;;;OAKG;IACH,OAAO,CAAC,UAAU;IA6BlB;;;;;;OAMG;IACH,OAAO,CAAC,oBAAoB;IAwC5B,OAAO,CAAC,cAAc;CAcvB"}

package/dist/SPHSimulator.js

@@ -0,0 +1,393 @@
+import { BufferGeometry, BufferAttribute } from 'three';
+/**
+ * SPHSimulator — Smoothed Particle Hydrodynamics fluid simulation.
+ * Blender: Fluid modifier (FLIP/SPH domain).
+ *
+ * Implements the Müller et al. 2003 SPH formulation:
+ *   - Poly6 kernel for density
+ *   - Spiky kernel gradient for pressure
+ *   - Viscosity kernel Laplacian for viscosity
+ *
+ * Workflow:
+ *   1. Construct with particle count + physics parameters
+ *   2. Call `setPositions(array)` to place particles in a volume
+ *   3. Call `step(dt)` each frame (dt ≈ 1/60)
+ *   4. Call `getGeometry()` to read back a THREE.BufferGeometry for rendering
+ *
+ * All physics parameters live in `parameters` for GSAP/st-keyframe:
+ *   gsap.to(sim.parameters, { viscosity: 0.5, duration: 1 })
+ */
+export class SPHSimulator {
+    constructor(opts = {}) {
+        // Spatial grid for O(n) neighbour search
+        this._grid = new Map();
+        /** Simulation clock. */
+        this.time = 0;
+        // Domain bounds (AABB collider)
+        this._domainMin = [-2, -2, -2];
+        this._domainMax = [2, 2, 2];
+        this.particleCount = Math.max(1, Math.round(opts.particleCount ?? 500));
+        const h = opts.smoothingRadius ?? 0.3;
+        const g = opts.gravity ?? [0, -9.8, 0];
+        this.viscosityType = opts.viscosityType ?? 'standard';
+        this.parameters = {
+            smoothingRadius: h,
+            restDensity: opts.restDensity ?? 1000,
+            pressureK: opts.pressureK ?? 200,
+            viscosity: opts.viscosity ?? 0.1,
+            xsphFactor: 0.5,
+            stiffFactor: 0.5,
+            gravityX: g[0],
+            gravityY: g[1],
+            gravityZ: g[2],
+            substeps: opts.substeps ?? 4,
+            surfaceTension: opts.surfaceTension ?? 0.0728,
+        };
+        this._px = new Float32Array(this.particleCount);
+        this._py = new Float32Array(this.particleCount);
+        this._pz = new Float32Array(this.particleCount);
+        this._vx = new Float32Array(this.particleCount);
+        this._vy = new Float32Array(this.particleCount);
+        this._vz = new Float32Array(this.particleCount);
+        this._density = new Float32Array(this.particleCount);
+        this._pressure = new Float32Array(this.particleCount);
+        // Default: arrange in a 3D grid block
+        this._initGrid();
+    }
+    // ── Initialisation ────────────────────────────────────────────────────────
+    _initGrid() {
+        const n = this.particleCount;
+        const cbrt = Math.ceil(Math.cbrt(n));
+        const spacing = this.parameters.smoothingRadius * 0.9;
+        let i = 0;
+        outer: for (let z = 0; z < cbrt; z++) {
+            for (let y = 0; y < cbrt; y++) {
+                for (let x = 0; x < cbrt; x++) {
+                    if (i >= n)
+                        break outer;
+                    this._px[i] = (x - cbrt / 2) * spacing;
+                    this._py[i] = (y - cbrt / 2) * spacing;
+                    this._pz[i] = (z - cbrt / 2) * spacing;
+                    i++;
+                }
+            }
+        }
+    }
+    /**
+     * Set particle positions from a flat array [x0,y0,z0, x1,y1,z1, ...].
+     */
+    setPositions(arr) {
+        const n = Math.min(this.particleCount, Math.floor(arr.length / 3));
+        for (let i = 0; i < n; i++) {
+            this._px[i] = arr[i * 3];
+            this._py[i] = arr[i * 3 + 1];
+            this._pz[i] = arr[i * 3 + 2];
+        }
+        this._vx.fill(0);
+        this._vy.fill(0);
+        this._vz.fill(0);
+    }
+    /**
+     * Set simulation domain (AABB boundary collider).
+     * Particles bounce off the domain walls.
+     */
+    setDomain(min, max) {
+        this._domainMin = [...min];
+        this._domainMax = [...max];
+    }
+    // ── Simulation ────────────────────────────────────────────────────────────
+    /**
+     * Advance the simulation by dt seconds.
+     * Recommended: pass clock.getDelta() clamped to 0.05 s.
+     */
+    step(dt) {
+        const sub = Math.max(1, Math.round(this.parameters.substeps));
+        const subDt = dt / sub;
+        for (let s = 0; s < sub; s++)
+            this._substep(subDt);
+        this.time += dt;
+    }
+    /**
+     * Build and return a BufferGeometry with a 'position' attribute.
+     * Suitable for THREE.Points rendering.
+     */
+    getGeometry() {
+        const positions = new Float32Array(this.particleCount * 3);
+        for (let i = 0; i < this.particleCount; i++) {
+            positions[i * 3] = this._px[i];
+            positions[i * 3 + 1] = this._py[i];
+            positions[i * 3 + 2] = this._pz[i];
+        }
+        const geo = new BufferGeometry();
+        geo.setAttribute('position', new BufferAttribute(positions, 3));
+        return geo;
+    }
+    /**
+     * Copy current particle positions into an existing BufferGeometry.
+     * More efficient than getGeometry() in the render loop.
+     */
+    updateGeometry(geo) {
+        const attr = geo.getAttribute('position');
+        const n = Math.min(attr.count, this.particleCount);
+        for (let i = 0; i < n; i++) {
+            attr.setXYZ(i, this._px[i], this._py[i], this._pz[i]);
+        }
+        attr.needsUpdate = true;
+    }
+    /** Per-particle density (read-only, updated after each step). */
+    getDensity(i) { return this._density[i]; }
+    /** Per-particle pressure (read-only). */
+    getPressure(i) { return this._pressure[i]; }
+    // ── Internals ─────────────────────────────────────────────────────────────
+    _substep(dt) {
+        this._buildSpatialGrid();
+        this._computeDensityPressure();
+        this._integrate(dt);
+        if (this.viscosityType === 'xsph')
+            this._applyXSPH();
+        if (this.viscosityType === 'stiff')
+            this._applyStiffViscosity(dt);
+        this._enforceDomain();
+    }
+    // Uniform spatial grid: maps cell key → particle indices
+    _buildSpatialGrid() {
+        const h = this.parameters.smoothingRadius;
+        this._grid.clear();
+        for (let i = 0; i < this.particleCount; i++) {
+            const key = this._cellKey(Math.floor(this._px[i] / h), Math.floor(this._py[i] / h), Math.floor(this._pz[i] / h));
+            let cell = this._grid.get(key);
+            if (!cell) {
+                cell = [];
+                this._grid.set(key, cell);
+            }
+            cell.push(i);
+        }
+    }
+    _cellKey(cx, cy, cz) {
+        // Simple hash: large primes, bit-truncated to 32-bit int range
+        return ((cx * 92837111) ^ (cy * 689287499) ^ (cz * 283923481)) | 0;
+    }
+    _neighbours(i, callback) {
+        const h = this.parameters.smoothingRadius;
+        const h2 = h * h;
+        const cx = Math.floor(this._px[i] / h);
+        const cy = Math.floor(this._py[i] / h);
+        const cz = Math.floor(this._pz[i] / h);
+        for (let dz = -1; dz <= 1; dz++) {
+            for (let dy = -1; dy <= 1; dy++) {
+                for (let dx = -1; dx <= 1; dx++) {
+                    const cell = this._grid.get(this._cellKey(cx + dx, cy + dy, cz + dz));
+                    if (!cell)
+                        continue;
+                    for (const j of cell) {
+                        if (j === i)
+                            continue;
+                        const rx = this._px[j] - this._px[i];
+                        const ry = this._py[j] - this._py[i];
+                        const rz = this._pz[j] - this._pz[i];
+                        const r2 = rx * rx + ry * ry + rz * rz;
+                        if (r2 < h2)
+                            callback(j, r2);
+                    }
+                }
+            }
+        }
+    }
+    // Müller 2003 Poly6 kernel
+    _poly6(r2) {
+        const h = this.parameters.smoothingRadius;
+        const h2 = h * h;
+        if (r2 > h2)
+            return 0;
+        const diff = h2 - r2;
+        return (315 / (64 * Math.PI * Math.pow(h, 9))) * diff * diff * diff;
+    }
+    // Spiky kernel gradient magnitude ∂W/∂r (scalar, apply direction separately)
+    _spikyGrad(r) {
+        const h = this.parameters.smoothingRadius;
+        if (r > h || r < 1e-6)
+            return 0;
+        const diff = h - r;
+        return -(45 / (Math.PI * Math.pow(h, 6))) * diff * diff;
+    }
+    // Viscosity kernel Laplacian
+    _viscLaplacian(r) {
+        const h = this.parameters.smoothingRadius;
+        if (r > h)
+            return 0;
+        return (45 / (Math.PI * Math.pow(h, 6))) * (h - r);
+    }
+    _computeDensityPressure() {
+        const { restDensity, pressureK } = this.parameters;
+        const particleMass = restDensity * Math.pow(this.parameters.smoothingRadius, 3) * 0.02;
+        for (let i = 0; i < this.particleCount; i++) {
+            let rho = this._poly6(0) * particleMass; // self contribution
+            this._neighbours(i, (_j, r2) => {
+                rho += particleMass * this._poly6(r2);
+            });
+            this._density[i] = Math.max(rho, restDensity * 0.001);
+            // Tait equation of state (water)
+            this._pressure[i] = Math.max(0, pressureK * (this._density[i] / restDensity - 1));
+        }
+    }
+    _integrate(dt) {
+        const { gravityX, gravityY, gravityZ, viscosity } = this.parameters;
+        const particleMass = this.parameters.restDensity *
+            Math.pow(this.parameters.smoothingRadius, 3) * 0.02;
+        for (let i = 0; i < this.particleCount; i++) {
+            let fx = 0, fy = 0, fz = 0;
+            const pi = this._pressure[i];
+            const rhoi = this._density[i];
+            this._neighbours(i, (j, r2) => {
+                const rx = this._px[j] - this._px[i];
+                const ry = this._py[j] - this._py[i];
+                const rz = this._pz[j] - this._pz[i];
+                const r = Math.sqrt(r2);
+                if (r < 1e-6)
+                    return;
+                const nx = rx / r, ny = ry / r, nz = rz / r;
+                const rhoj = this._density[j];
+                const pj = this._pressure[j];
+                const m = particleMass;
+                // Pressure force (symmetric formulation)
+                const spg = this._spikyGrad(r);
+                const pf = -m * (pi + pj) / (2 * rhoj) * spg;
+                fx += pf * nx;
+                fy += pf * ny;
+                fz += pf * nz;
+                // Standard viscosity (Müller 2003) — skipped for xsph/stiff modes
+                if (this.viscosityType === 'standard') {
+                    const visc = viscosity * m / rhoj * this._viscLaplacian(r);
+                    fx += visc * (this._vx[j] - this._vx[i]);
+                    fy += visc * (this._vy[j] - this._vy[i]);
+                    fz += visc * (this._vz[j] - this._vz[i]);
+                }
+            });
+            // Gravity
+            fx += gravityX * rhoi;
+            fy += gravityY * rhoi;
+            fz += gravityZ * rhoi;
+            // Euler integration
+            const invRho = 1 / rhoi;
+            this._vx[i] += (fx * invRho) * dt;
+            this._vy[i] += (fy * invRho) * dt;
+            this._vz[i] += (fz * invRho) * dt;
+            this._px[i] += this._vx[i] * dt;
+            this._py[i] += this._vy[i] * dt;
+            this._pz[i] += this._vz[i] * dt;
+        }
+    }
+    /**
+     * XSPH viscosity (Schechter & Bridson 2012).
+     * Blends each particle's velocity toward a weighted average of its neighbours.
+     * More stable than Müller viscosity, good for thick/honey-like fluids.
+     * Controlled by parameters.xsphFactor [0–1].
+     */
+    _applyXSPH() {
+        const h = this.parameters.smoothingRadius;
+        const h2 = h * h;
+        const factor = this.parameters.xsphFactor;
+        const particleMass = this.parameters.restDensity * Math.pow(h, 3) * 0.02;
+        // Temporary correction arrays
+        const cx = new Float32Array(this.particleCount);
+        const cy = new Float32Array(this.particleCount);
+        const cz = new Float32Array(this.particleCount);
+        for (let i = 0; i < this.particleCount; i++) {
+            this._neighbours(i, (j, r2) => {
+                const w = this._poly6(r2);
+                const avgRho = (this._density[i] + this._density[j]) * 0.5;
+                const wm = w * particleMass / avgRho;
+                cx[i] += wm * (this._vx[j] - this._vx[i]);
+                cy[i] += wm * (this._vy[j] - this._vy[i]);
+                cz[i] += wm * (this._vz[j] - this._vz[i]);
+            });
+        }
+        for (let i = 0; i < this.particleCount; i++) {
+            this._vx[i] += factor * cx[i];
+            this._vy[i] += factor * cy[i];
+            this._vz[i] += factor * cz[i];
+        }
+    }
+    /**
+     * Stiff viscosity (position-based / divergence damping).
+     * Applies a velocity correction proportional to the local velocity divergence,
+     * strongly suppressing compressible motion.
+     * Best for near-incompressible high-viscosity simulations.
+     * Controlled by parameters.stiffFactor.
+     */
+    _applyStiffViscosity(dt) {
+        const h = this.parameters.smoothingRadius;
+        const stiff = this.parameters.stiffFactor;
+        const particleMass = this.parameters.restDensity * Math.pow(h, 3) * 0.02;
+        const cx = new Float32Array(this.particleCount);
+        const cy = new Float32Array(this.particleCount);
+        const cz = new Float32Array(this.particleCount);
+        for (let i = 0; i < this.particleCount; i++) {
+            let divV = 0;
+            this._neighbours(i, (j, r2) => {
+                const r = Math.sqrt(r2);
+                if (r < 1e-6)
+                    return;
+                const rx = this._px[j] - this._px[i];
+                const ry = this._py[j] - this._py[i];
+                const rz = this._pz[j] - this._pz[i];
+                const dvx = this._vx[j] - this._vx[i];
+                const dvy = this._vy[j] - this._vy[i];
+                const dvz = this._vz[j] - this._vz[i];
+                // Divergence contribution: (v_ij · r_ij) / |r_ij|²
+                const dot = dvx * rx + dvy * ry + dvz * rz;
+                const spg = this._spikyGrad(r);
+                const coef = particleMass / this._density[j] * spg / r;
+                divV += dot / r2 * coef * (-1);
+                // Velocity correction: push toward incompressibility
+                cx[i] += coef * dot / r2 * rx;
+                cy[i] += coef * dot / r2 * ry;
+                cz[i] += coef * dot / r2 * rz;
+            });
+            void divV;
+        }
+        for (let i = 0; i < this.particleCount; i++) {
+            this._vx[i] -= stiff * dt * cx[i];
+            this._vy[i] -= stiff * dt * cy[i];
+            this._vz[i] -= stiff * dt * cz[i];
+        }
+    }
+    _enforceDomain() {
+        const restitution = 0.3;
+        const [minX, minY, minZ] = this._domainMin;
+        const [maxX, maxY, maxZ] = this._domainMax;
+        for (let i = 0; i < this.particleCount; i++) {
+            if (this._px[i] < minX) {
+                this._px[i] = minX;
+                if (this._vx[i] < 0)
+                    this._vx[i] *= -restitution;
+            }
+            if (this._px[i] > maxX) {
+                this._px[i] = maxX;
+                if (this._vx[i] > 0)
+                    this._vx[i] *= -restitution;
+            }
+            if (this._py[i] < minY) {
+                this._py[i] = minY;
+                if (this._vy[i] < 0)
+                    this._vy[i] *= -restitution;
+            }
+            if (this._py[i] > maxY) {
+                this._py[i] = maxY;
+                if (this._vy[i] > 0)
+                    this._vy[i] *= -restitution;
+            }
+            if (this._pz[i] < minZ) {
+                this._pz[i] = minZ;
+                if (this._vz[i] < 0)
+                    this._vz[i] *= -restitution;
+            }
+            if (this._pz[i] > maxZ) {
+                this._pz[i] = maxZ;
+                if (this._vz[i] > 0)
+                    this._vz[i] *= -restitution;
+            }
+        }
+    }
+}
+//# sourceMappingURL=SPHSimulator.js.map

package/dist/SPHSimulator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"SPHSimulator.js","sourceRoot":"","sources":["../src/SPHSimulator.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAE,eAAe,EAAE,MAAM,OAAO,CAAA;AAiCvD;;;;;;;;;;;;;;;;;GAiBG;AACH,MAAM,OAAO,YAAY;IA2CvB,YAAY,OAA4B,EAAE;QAV1C,yCAAyC;QACjC,UAAK,GAA0B,IAAI,GAAG,EAAE,CAAA;QAEhD,wBAAwB;QACxB,SAAI,GAAG,CAAC,CAAA;QAER,gCAAgC;QACxB,eAAU,GAA2B,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAA;QACjD,eAAU,GAA2B,CAAE,CAAC,EAAG,CAAC,EAAG,CAAC,CAAC,CAAA;QAGvD,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,aAAa,IAAI,GAAG,CAAC,CAAC,CAAA;QAEvE,MAAM,CAAC,GAAG,IAAI,CAAC,eAAe,IAAI,GAAG,CAAA;QACrC,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC,CAAA;QAEtC,IAAI,CAAC,aAAa,GAAG,IAAI,CAAC,aAAa,IAAI,UAAU,CAAA;QAErD,IAAI,CAAC,UAAU,GAAG;YAChB,eAAe,EAAE,CAAC;YAClB,WAAW,EAAM,IAAI,CAAC,WAAW,IAAO,IAAI;YAC5C,SAAS,EAAQ,IAAI,CAAC,SAAS,IAAS,GAAG;YAC3C,SAAS,EAAQ,IAAI,CAAC,SAAS,IAAS,GAAG;YAC3C,UAAU,EAAO,GAAG;YACpB,WAAW,EAAM,GAAG;YACpB,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,IAAU,CAAC;YACzC,cAAc,EAAG,IAAI,CAAC,cAAc,IAAI,MAAM;SAC/C,CAAA;QAED,IAAI,CAAC,GAAG,GAAS,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QACrD,IAAI,CAAC,GAAG,GAAS,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QACrD,IAAI,CAAC,GAAG,GAAS,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QACrD,IAAI,CAAC,GAAG,GAAS,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QACrD,IAAI,CAAC,GAAG,GAAS,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QACrD,IAAI,CAAC,GAAG,GAAS,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QACrD,IAAI,CAAC,QAAQ,GAAI,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QACrD,IAAI,CAAC,SAAS,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QAErD,sCAAsC;QACtC,IAAI,CAAC,SAAS,EAAE,CAAA;IAClB,CAAC;IAED,6EAA6E;IAErE,SAAS;QACf,MAAM,CAAC,GAAK,IAAI,CAAC,aAAa,CAAA;QAC9B,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAA;QACpC,MAAM,OAAO,GAAG,IAAI,CAAC,UAAU,CAAC,eAAe,GAAG,GAAG,CAAA;QACrD,IAAI,CAAC,GAAG,CAAC,CAAA;QACT,KAAK,EAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC9B,IAAI,CAAC,IAAI,CAAC;wBAAE,MAAM,KAAK,CAAA;oBACvB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAC,CAAC,CAAC,GAAG,OAAO,CAAA;oBACpC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAC,CAAC,CAAC,GAAG,OAAO,CAAA;oBACpC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAC,CAAC,CAAC,GAAG,OAAO,CAAA;oBACpC,CAAC,EAAE,CAAA;gBACL,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED;;OAEG;IACH,YAAY,CAAC,GAAsB;QACjC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,EAAE,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAA;QAClE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAC,CAAC,CAAC,CAAA;YACtB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,CAAA;YACxB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,CAAA;QAC1B,CAAC;QACD,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;IACtD,CAAC;IAED;;;OAGG;IACH,SAAS,CACP,GAA2B,EAC3B,GAA2B;QAE3B,IAAI,CAAC,UAAU,GAAG,CAAC,GAAG,GAAG,CAA2B,CAAA;QACpD,IAAI,CAAC,UAAU,GAAG,CAAC,GAAG,GAAG,CAA2B,CAAA;IACtD,CAAC;IAED,6EAA6E;IAE7E;;;OAGG;IACH,IAAI,CAAC,EAAU;QACb,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,UAAU,CAAC,QAAQ,CAAC,CAAC,CAAA;QAC7D,MAAM,KAAK,GAAG,EAAE,GAAG,GAAG,CAAA;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;YAAE,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAA;QAClD,IAAI,CAAC,IAAI,IAAI,EAAE,CAAA;IACjB,CAAC;IAED;;;OAGG;IACH,WAAW;QACT,MAAM,SAAS,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,GAAG,CAAC,CAAC,CAAA;QAC1D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,SAAS,CAAC,CAAC,GAAC,CAAC,CAAC,GAAK,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;YAC9B,SAAS,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;YAC9B,SAAS,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;QAChC,CAAC;QACD,MAAM,GAAG,GAAG,IAAI,cAAc,EAAE,CAAA;QAChC,GAAG,CAAC,YAAY,CAAC,UAAU,EAAE,IAAI,eAAe,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC,CAAA;QAC/D,OAAO,GAAG,CAAA;IACZ,CAAC;IAED;;;OAGG;IACH,cAAc,CAAC,GAAmB;QAChC,MAAM,IAAI,GAAG,GAAG,CAAC,YAAY,CAAC,UAAU,CAAoB,CAAA;QAC5D,MAAM,CAAC,GAAM,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,aAAa,CAAC,CAAA;QACrD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAA;QACvD,CAAC;QACD,IAAI,CAAC,WAAW,GAAG,IAAI,CAAA;IACzB,CAAC;IAED,iEAAiE;IACjE,UAAU,CAAC,CAAS,IAAY,OAAO,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAA,CAAC,CAAC;IAEzD,yCAAyC;IACzC,WAAW,CAAC,CAAS,IAAY,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAA,CAAC,CAAC;IAE3D,6EAA6E;IAErE,QAAQ,CAAC,EAAU;QACzB,IAAI,CAAC,iBAAiB,EAAE,CAAA;QACxB,IAAI,CAAC,uBAAuB,EAAE,CAAA;QAC9B,IAAI,CAAC,UAAU,CAAC,EAAE,CAAC,CAAA;QACnB,IAAI,IAAI,CAAC,aAAa,KAAK,MAAM;YAAG,IAAI,CAAC,UAAU,EAAE,CAAA;QACrD,IAAI,IAAI,CAAC,aAAa,KAAK,OAAO;YAAE,IAAI,CAAC,oBAAoB,CAAC,EAAE,CAAC,CAAA;QACjE,IAAI,CAAC,cAAc,EAAE,CAAA;IACvB,CAAC;IAED,yDAAyD;IACjD,iBAAiB;QACvB,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,eAAe,CAAA;QACzC,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,CAAA;QAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,MAAM,GAAG,GAAG,IAAI,CAAC,QAAQ,CACvB,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAC3B,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAC3B,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAC5B,CAAA;YACD,IAAI,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,CAAA;YAC9B,IAAI,CAAC,IAAI,EAAE,CAAC;gBAAC,IAAI,GAAG,EAAE,CAAC;gBAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAA;YAAC,CAAC;YACnD,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;QACd,CAAC;IACH,CAAC;IAEO,QAAQ,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU;QACjD,+DAA+D;QAC/D,OAAO,CAAC,CAAC,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,GAAG,SAAS,CAAC,GAAG,CAAC,EAAE,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAA;IACpE,CAAC;IAEO,WAAW,CAAC,CAAS,EAAE,QAAyC;QACtE,MAAM,CAAC,GAAI,IAAI,CAAC,UAAU,CAAC,eAAe,CAAA;QAC1C,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAA;QAChB,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QACtC,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QACtC,MAAM,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QAEtC,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAChC,KAAK,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;oBAChC,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAE,GAAC,EAAE,EAAE,EAAE,GAAC,EAAE,EAAE,EAAE,GAAC,EAAE,CAAC,CAAC,CAAA;oBAC/D,IAAI,CAAC,IAAI;wBAAE,SAAQ;oBACnB,KAAK,MAAM,CAAC,IAAI,IAAI,EAAE,CAAC;wBACrB,IAAI,CAAC,KAAK,CAAC;4BAAE,SAAQ;wBACrB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;wBACpC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;wBACpC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;wBACpC,MAAM,EAAE,GAAG,EAAE,GAAC,EAAE,GAAG,EAAE,GAAC,EAAE,GAAG,EAAE,GAAC,EAAE,CAAA;wBAChC,IAAI,EAAE,GAAG,EAAE;4BAAE,QAAQ,CAAC,CAAC,EAAE,EAAE,CAAC,CAAA;oBAC9B,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,2BAA2B;IACnB,MAAM,CAAC,EAAU;QACvB,MAAM,CAAC,GAAI,IAAI,CAAC,UAAU,CAAC,eAAe,CAAA;QAC1C,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAA;QAChB,IAAI,EAAE,GAAG,EAAE;YAAE,OAAO,CAAC,CAAA;QACrB,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAA;QACpB,OAAO,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,CAAA;IACrE,CAAC;IAED,6EAA6E;IACrE,UAAU,CAAC,CAAS;QAC1B,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,eAAe,CAAA;QACzC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,IAAI;YAAE,OAAO,CAAC,CAAA;QAC/B,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAA;QAClB,OAAO,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAA;IACzD,CAAC;IAED,6BAA6B;IACrB,cAAc,CAAC,CAAS;QAC9B,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC,eAAe,CAAA;QACzC,IAAI,CAAC,GAAG,CAAC;YAAE,OAAO,CAAC,CAAA;QACnB,OAAO,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;IACpD,CAAC;IAEO,uBAAuB;QAC7B,MAAM,EAAE,WAAW,EAAE,SAAS,EAAE,GAAG,IAAI,CAAC,UAAU,CAAA;QAClD,MAAM,YAAY,GAAG,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,UAAU,CAAC,eAAe,EAAE,CAAC,CAAC,GAAG,IAAI,CAAA;QAEtF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,IAAI,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,YAAY,CAAA,CAAE,oBAAoB;YAC7D,IAAI,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE;gBAC7B,GAAG,IAAI,YAAY,GAAG,IAAI,CAAC,MAAM,CAAC,EAAE,CAAC,CAAA;YACvC,CAAC,CAAC,CAAA;YACF,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAI,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,WAAW,GAAG,KAAK,CAAC,CAAA;YACtD,iCAAiC;YACjC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,SAAS,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,WAAW,GAAG,CAAC,CAAC,CAAC,CAAA;QACnF,CAAC;IACH,CAAC;IAEO,UAAU,CAAC,EAAU;QAC3B,MAAM,EAAE,QAAQ,EAAE,QAAQ,EAAE,QAAQ,EAAE,SAAS,EAAE,GAAG,IAAI,CAAC,UAAU,CAAA;QACnE,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,WAAW;YAC9C,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,UAAU,CAAC,eAAe,EAAE,CAAC,CAAC,GAAG,IAAI,CAAA;QAErD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAA;YAE1B,MAAM,EAAE,GAAK,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAA;YAC9B,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAA;YAE7B,IAAI,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;gBAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;gBACpC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;gBACpC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;gBACpC,MAAM,CAAC,GAAI,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAA;gBACxB,IAAI,CAAC,GAAG,IAAI;oBAAE,OAAM;gBAEpB,MAAM,EAAE,GAAG,EAAE,GAAC,CAAC,EAAE,EAAE,GAAG,EAAE,GAAC,CAAC,EAAE,EAAE,GAAG,EAAE,GAAC,CAAC,CAAA;gBACrC,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAA;gBAC7B,MAAM,EAAE,GAAK,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAA;gBAC9B,MAAM,CAAC,GAAM,YAAY,CAAA;gBAEzB,yCAAyC;gBACzC,MAAM,GAAG,GAAG,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAA;gBAC9B,MAAM,EAAE,GAAI,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,GAAG,CAAA;gBAC7C,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;gBAAE,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;gBAAE,EAAE,IAAI,EAAE,GAAG,EAAE,CAAA;gBAE7C,kEAAkE;gBAClE,IAAI,IAAI,CAAC,aAAa,KAAK,UAAU,EAAE,CAAC;oBACtC,MAAM,IAAI,GAAG,SAAS,GAAG,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAA;oBAC1D,EAAE,IAAI,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAA;oBACxC,EAAE,IAAI,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAA;oBACxC,EAAE,IAAI,IAAI,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAA;gBAC1C,CAAC;YACH,CAAC,CAAC,CAAA;YAEF,UAAU;YACV,EAAE,IAAI,QAAQ,GAAG,IAAI,CAAA;YACrB,EAAE,IAAI,QAAQ,GAAG,IAAI,CAAA;YACrB,EAAE,IAAI,QAAQ,GAAG,IAAI,CAAA;YAErB,oBAAoB;YACpB,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAA;YACvB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,EAAE,CAAA;YACjC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,EAAE,CAAA;YACjC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,EAAE,CAAA;YACjC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAA;YAC/B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAA;YAC/B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAA;QACjC,CAAC;IACH,CAAC;IAED;;;;;OAKG;IACK,UAAU;QAChB,MAAM,CAAC,GAAQ,IAAI,CAAC,UAAU,CAAC,eAAe,CAAA;QAC9C,MAAM,EAAE,GAAO,CAAC,GAAG,CAAC,CAAA;QACpB,MAAM,MAAM,GAAG,IAAI,CAAC,UAAU,CAAC,UAAU,CAAA;QACzC,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAA;QAExE,8BAA8B;QAC9B,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QAC/C,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QAC/C,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QAE/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,IAAI,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;gBAC5B,MAAM,CAAC,GAAQ,IAAI,CAAC,MAAM,CAAC,EAAE,CAAC,CAAA;gBAC9B,MAAM,MAAM,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAA;gBAC1D,MAAM,EAAE,GAAO,CAAC,GAAG,YAAY,GAAG,MAAM,CAAA;gBACxC,EAAE,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAA;gBACzC,EAAE,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAA;gBACzC,EAAE,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAA;YAC3C,CAAC,CAAC,CAAA;QACJ,CAAC;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC,CAAA;YAC7B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC,CAAA;YAC7B,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC,CAAA;QAC/B,CAAC;IACH,CAAC;IAED;;;;;;OAMG;IACK,oBAAoB,CAAC,EAAU;QACrC,MAAM,CAAC,GAAY,IAAI,CAAC,UAAU,CAAC,eAAe,CAAA;QAClD,MAAM,KAAK,GAAQ,IAAI,CAAC,UAAU,CAAC,WAAW,CAAA;QAC9C,MAAM,YAAY,GAAG,IAAI,CAAC,UAAU,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,IAAI,CAAA;QAExE,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QAC/C,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QAC/C,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,aAAa,CAAC,CAAA;QAE/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,IAAI,IAAI,GAAG,CAAC,CAAA;YACZ,IAAI,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE;gBAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAA;gBACvB,IAAI,CAAC,GAAG,IAAI;oBAAE,OAAM;gBACpB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;gBACpC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;gBACpC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;gBACpC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;gBACrC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;gBACrC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAA;gBACrC,mDAAmD;gBACnD,MAAM,GAAG,GAAI,GAAG,GAAC,EAAE,GAAG,GAAG,GAAC,EAAE,GAAG,GAAG,GAAC,EAAE,CAAA;gBACrC,MAAM,GAAG,GAAI,IAAI,CAAC,UAAU,CAAC,CAAC,CAAC,CAAA;gBAC/B,MAAM,IAAI,GAAG,YAAY,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAA;gBACtD,IAAI,IAAI,GAAG,GAAG,EAAE,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAA;gBAC9B,qDAAqD;gBACrD,EAAE,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAA;gBAC7B,EAAE,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAA;gBAC7B,EAAE,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAA;YAC/B,CAAC,CAAC,CAAA;YACF,KAAK,IAAI,CAAA;QACX,CAAC;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAA;YACjC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAA;YACjC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAA;QACnC,CAAC;IACH,CAAC;IAEO,cAAc;QACpB,MAAM,WAAW,GAAG,GAAG,CAAA;QACvB,MAAM,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC,UAAU,CAAA;QAC1C,MAAM,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,IAAI,CAAC,UAAU,CAAA;QAE1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5C,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,WAAW,CAAA;YAAC,CAAC;YAChG,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,WAAW,CAAA;YAAC,CAAC;YAChG,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,WAAW,CAAA;YAAC,CAAC;YAChG,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,WAAW,CAAA;YAAC,CAAC;YAChG,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,WAAW,CAAA;YAAC,CAAC;YAChG,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC;gBAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;gBAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;oBAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,WAAW,CAAA;YAAC,CAAC;QAClG,CAAC;IACH,CAAC;CACF"}

package/dist/index.d.ts

@@ -0,0 +1,10 @@
+export { FLIPSimulator } from './FLIPSimulator.js';
+export { MACGrid } from './MACGrid.js';
+export { SPHSimulator } from './SPHSimulator.js';
+export { FluidEmitter } from './FluidEmitter.js';
+export { MarchingCubes } from './MarchingCubes.js';
+export type { FLIPSimulatorOptions } from './FLIPSimulator.js';
+export type { SPHSimulatorOptions } from './SPHSimulator.js';
+export type { FluidEmitterOptions } from './FluidEmitter.js';
+export type { MarchingCubesOptions } from './MarchingCubes.js';
+//# sourceMappingURL=index.d.ts.map

package/dist/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,aAAa,EAAE,MAAO,oBAAoB,CAAA;AACnD,OAAO,EAAE,OAAO,EAAE,MAAa,cAAc,CAAA;AAC7C,OAAO,EAAE,YAAY,EAAE,MAAQ,mBAAmB,CAAA;AAClD,OAAO,EAAE,YAAY,EAAE,MAAQ,mBAAmB,CAAA;AAClD,OAAO,EAAE,aAAa,EAAE,MAAO,oBAAoB,CAAA;AAEnD,YAAY,EAAE,oBAAoB,EAAE,MAAM,oBAAoB,CAAA;AAC9D,YAAY,EAAE,mBAAmB,EAAE,MAAO,mBAAmB,CAAA;AAC7D,YAAY,EAAE,mBAAmB,EAAE,MAAO,mBAAmB,CAAA;AAC7D,YAAY,EAAE,oBAAoB,EAAE,MAAM,oBAAoB,CAAA"}