@triforge/fluid-core 0.1.0

package/dist/MACGrid.d.ts

@@ -0,0 +1,58 @@
+/**
+ * MACGrid — Marker-And-Cell staggered velocity grid.
+ *
+ * Velocity components are stored at face centres (staggered):
+ *   u[i][j][k]  — X-velocity at face (i-½, j, k)
+ *   v[i][j][k]  — Y-velocity at face (i, j-½, k)
+ *   w[i][j][k]  — Z-velocity at face (i, j, k-½)
+ * Pressure p[i][j][k] is at cell centres.
+ *
+ * Grid cell (i,j,k) covers world-space region:
+ *   [ox + i*h, ox + (i+1)*h] × [oy + j*h, ...] × [oz + k*h, ...]
+ */
+export declare class MACGrid {
+    readonly nx: number;
+    readonly ny: number;
+    readonly nz: number;
+    readonly h: number;
+    readonly ox: number;
+    readonly oy: number;
+    readonly oz: number;
+    u: Float64Array;
+    v: Float64Array;
+    w: Float64Array;
+    uOld: Float64Array;
+    vOld: Float64Array;
+    wOld: Float64Array;
+    uW: Float64Array;
+    vW: Float64Array;
+    wW: Float64Array;
+    pressure: Float64Array;
+    cellType: Uint8Array;
+    constructor(nx: number, ny: number, nz: number, h: number, ox?: number, oy?: number, oz?: number);
+    /** Flat index for cell (i,j,k). */
+    idx(i: number, j: number, k: number): number;
+    /** Flat index for u-face (i,j,k) — grid is (nx+1)*ny*nz. */
+    uIdx(i: number, j: number, k: number): number;
+    /** Flat index for v-face (i,j,k) — grid is nx*(ny+1)*nz. */
+    vIdx(i: number, j: number, k: number): number;
+    /** Flat index for w-face (i,j,k) — grid is nx*ny*(nz+1). */
+    wIdx(i: number, j: number, k: number): number;
+    /** Trilinear interpolation of u-velocity at world position (wx,wy,wz). */
+    interpU(wx: number, wy: number, wz: number): number;
+    /** Trilinear interpolation of v-velocity at world position. */
+    interpV(wx: number, wy: number, wz: number): number;
+    /** Trilinear interpolation of w-velocity at world position. */
+    interpW(wx: number, wy: number, wz: number): number;
+    private _interp;
+    clearVelocities(): void;
+    saveVelocities(): void;
+    normalizeWeights(): void;
+    /** Mark border cells as solid, fluid cells from a predicate. */
+    markCells(isFluid: (i: number, j: number, k: number) => boolean): void;
+    /** Enforce solid (no-through) boundary: zero normal velocity at solid faces. */
+    enforceBoundary(): void;
+    /** Apply gravity to grid v velocities over timestep dt. */
+    applyGravity(gx: number, gy: number, gz: number, dt: number): void;
+}
+//# sourceMappingURL=MACGrid.d.ts.map

package/dist/MACGrid.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"MACGrid.d.ts","sourceRoot":"","sources":["../src/MACGrid.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AACH,qBAAa,OAAO;IAClB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAA;IACnB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAA;IACnB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAA;IACnB,QAAQ,CAAC,CAAC,EAAG,MAAM,CAAA;IACnB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAA;IACnB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAA;IACnB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAA;IAGnB,CAAC,EAAE,YAAY,CAAA;IACf,CAAC,EAAE,YAAY,CAAA;IACf,CAAC,EAAE,YAAY,CAAA;IAGf,IAAI,EAAE,YAAY,CAAA;IAClB,IAAI,EAAE,YAAY,CAAA;IAClB,IAAI,EAAE,YAAY,CAAA;IAGlB,EAAE,EAAE,YAAY,CAAA;IAChB,EAAE,EAAE,YAAY,CAAA;IAChB,EAAE,EAAE,YAAY,CAAA;IAGhB,QAAQ,EAAE,YAAY,CAAA;IAGtB,QAAQ,EAAE,UAAU,CAAA;gBAER,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,SAAI,EAAE,EAAE,SAAI,EAAE,EAAE,SAAI;IAoBjF,mCAAmC;IACnC,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAI5C,4DAA4D;IAC5D,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAI7C,4DAA4D;IAC5D,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAI7C,4DAA4D;IAC5D,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM;IAM7C,0EAA0E;IAC1E,OAAO,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM;IAKnD,+DAA+D;IAC/D,OAAO,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM;IAKnD,+DAA+D;IAC/D,OAAO,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM;IAKnD,OAAO,CAAC,OAAO;IAuBf,eAAe,IAAI,IAAI;IAKvB,cAAc,IAAI,IAAI;IAMtB,gBAAgB,IAAI,IAAI;IAMxB,gEAAgE;IAChE,SAAS,CAAC,OAAO,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,OAAO,GAAG,IAAI;IAUtE,gFAAgF;IAChF,eAAe,IAAI,IAAI;IAmBvB,2DAA2D;IAC3D,YAAY,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,IAAI;CASnE"}

package/dist/MACGrid.js

@@ -0,0 +1,154 @@
+/**
+ * MACGrid — Marker-And-Cell staggered velocity grid.
+ *
+ * Velocity components are stored at face centres (staggered):
+ *   u[i][j][k]  — X-velocity at face (i-½, j, k)
+ *   v[i][j][k]  — Y-velocity at face (i, j-½, k)
+ *   w[i][j][k]  — Z-velocity at face (i, j, k-½)
+ * Pressure p[i][j][k] is at cell centres.
+ *
+ * Grid cell (i,j,k) covers world-space region:
+ *   [ox + i*h, ox + (i+1)*h] × [oy + j*h, ...] × [oz + k*h, ...]
+ */
+export class MACGrid {
+    constructor(nx, ny, nz, h, ox = 0, oy = 0, oz = 0) {
+        this.nx = nx;
+        this.ny = ny;
+        this.nz = nz;
+        this.h = h;
+        this.ox = ox;
+        this.oy = oy;
+        this.oz = oz;
+        this.u = new Float64Array((nx + 1) * ny * nz);
+        this.v = new Float64Array(nx * (ny + 1) * nz);
+        this.w = new Float64Array(nx * ny * (nz + 1));
+        this.uOld = new Float64Array(this.u.length);
+        this.vOld = new Float64Array(this.v.length);
+        this.wOld = new Float64Array(this.w.length);
+        this.uW = new Float64Array(this.u.length);
+        this.vW = new Float64Array(this.v.length);
+        this.wW = new Float64Array(this.w.length);
+        this.pressure = new Float64Array(nx * ny * nz);
+        this.cellType = new Uint8Array(nx * ny * nz);
+    }
+    // ── Index helpers ─────────────────────────────────────────────────────────
+    /** Flat index for cell (i,j,k). */
+    idx(i, j, k) {
+        return i + this.nx * (j + this.ny * k);
+    }
+    /** Flat index for u-face (i,j,k) — grid is (nx+1)*ny*nz. */
+    uIdx(i, j, k) {
+        return i + (this.nx + 1) * (j + this.ny * k);
+    }
+    /** Flat index for v-face (i,j,k) — grid is nx*(ny+1)*nz. */
+    vIdx(i, j, k) {
+        return i + this.nx * (j + (this.ny + 1) * k);
+    }
+    /** Flat index for w-face (i,j,k) — grid is nx*ny*(nz+1). */
+    wIdx(i, j, k) {
+        return i + this.nx * (j + this.ny * k);
+    }
+    // ── Interpolation ─────────────────────────────────────────────────────────
+    /** Trilinear interpolation of u-velocity at world position (wx,wy,wz). */
+    interpU(wx, wy, wz) {
+        return this._interp(this.u, wx - this.ox, wy - this.oy - 0.5 * this.h, wz - this.oz - 0.5 * this.h, this.nx + 1, this.ny, this.nz);
+    }
+    /** Trilinear interpolation of v-velocity at world position. */
+    interpV(wx, wy, wz) {
+        return this._interp(this.v, wx - this.ox - 0.5 * this.h, wy - this.oy, wz - this.oz - 0.5 * this.h, this.nx, this.ny + 1, this.nz);
+    }
+    /** Trilinear interpolation of w-velocity at world position. */
+    interpW(wx, wy, wz) {
+        return this._interp(this.w, wx - this.ox - 0.5 * this.h, wy - this.oy - 0.5 * this.h, wz - this.oz, this.nx, this.ny, this.nz + 1);
+    }
+    _interp(arr, lx, ly, lz, Nx, Ny, Nz) {
+        const h = this.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)]));
+    }
+    // ── Utilities ─────────────────────────────────────────────────────────────
+    clearVelocities() {
+        this.u.fill(0);
+        this.v.fill(0);
+        this.w.fill(0);
+        this.uW.fill(0);
+        this.vW.fill(0);
+        this.wW.fill(0);
+    }
+    saveVelocities() {
+        this.uOld.set(this.u);
+        this.vOld.set(this.v);
+        this.wOld.set(this.w);
+    }
+    normalizeWeights() {
+        for (let i = 0; i < this.u.length; i++)
+            if (this.uW[i] > 0)
+                this.u[i] /= this.uW[i];
+        for (let i = 0; i < this.v.length; i++)
+            if (this.vW[i] > 0)
+                this.v[i] /= this.vW[i];
+        for (let i = 0; i < this.w.length; i++)
+            if (this.wW[i] > 0)
+                this.w[i] /= this.wW[i];
+    }
+    /** Mark border cells as solid, fluid cells from a predicate. */
+    markCells(isFluid) {
+        const { nx, ny, nz } = this;
+        for (let k = 0; k < nz; k++)
+            for (let j = 0; j < ny; j++)
+                for (let i = 0; i < nx; i++) {
+                    const solid = i === 0 || i === nx - 1 || j === 0 || j === ny - 1 || k === 0 || k === nz - 1;
+                    this.cellType[this.idx(i, j, k)] = solid ? 2 : (isFluid(i, j, k) ? 1 : 0);
+                }
+    }
+    /** Enforce solid (no-through) boundary: zero normal velocity at solid faces. */
+    enforceBoundary() {
+        const { nx, ny, nz } = this;
+        // Left/right walls (i=0 and i=nx)
+        for (let k = 0; k < nz; k++)
+            for (let j = 0; j < ny; j++) {
+                this.u[this.uIdx(0, j, k)] = 0;
+                this.u[this.uIdx(nx, j, k)] = 0;
+            }
+        // Bottom/top walls
+        for (let k = 0; k < nz; k++)
+            for (let i = 0; i < nx; i++) {
+                this.v[this.vIdx(i, 0, k)] = 0;
+                this.v[this.vIdx(i, ny, k)] = 0;
+            }
+        // Front/back walls
+        for (let j = 0; j < ny; j++)
+            for (let i = 0; i < nx; i++) {
+                this.w[this.wIdx(i, j, 0)] = 0;
+                this.w[this.wIdx(i, j, nz)] = 0;
+            }
+    }
+    /** Apply gravity to grid v velocities over timestep dt. */
+    applyGravity(gx, gy, gz, dt) {
+        const { nx, ny, nz } = this;
+        for (let k = 0; k < nz; k++)
+            for (let j = 0; j <= ny; j++)
+                for (let i = 0; i < nx; i++)
+                    this.v[this.vIdx(i, j, k)] += gy * dt;
+        for (let k = 0; k < nz; k++)
+            for (let j = 0; j < ny; j++)
+                for (let i = 0; i <= nx; i++)
+                    this.u[this.uIdx(i, j, k)] += gx * dt;
+        for (let k = 0; k <= nz; k++)
+            for (let j = 0; j < ny; j++)
+                for (let i = 0; i < nx; i++)
+                    this.w[this.wIdx(i, j, k)] += gz * dt;
+    }
+}
+//# sourceMappingURL=MACGrid.js.map

package/dist/MACGrid.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"MACGrid.js","sourceRoot":"","sources":["../src/MACGrid.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AACH,MAAM,OAAO,OAAO;IA8BlB,YAAY,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,CAAS,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC;QAC/E,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAA;QACxC,IAAI,CAAC,CAAC,GAAI,CAAC,CAAA;QACX,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAA;QAExC,IAAI,CAAC,CAAC,GAAU,IAAI,YAAY,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAA;QACpD,IAAI,CAAC,CAAC,GAAU,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAA;QACpD,IAAI,CAAC,CAAC,GAAU,IAAI,YAAY,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAA;QACpD,IAAI,CAAC,IAAI,GAAO,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAA;QAC/C,IAAI,CAAC,IAAI,GAAO,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAA;QAC/C,IAAI,CAAC,IAAI,GAAO,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAA;QAC/C,IAAI,CAAC,EAAE,GAAS,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAA;QAC/C,IAAI,CAAC,EAAE,GAAS,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAA;QAC/C,IAAI,CAAC,EAAE,GAAS,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAA;QAC/C,IAAI,CAAC,QAAQ,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAA;QAC9C,IAAI,CAAC,QAAQ,GAAG,IAAI,UAAU,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAA;IAC9C,CAAC;IAED,6EAA6E;IAE7E,mCAAmC;IACnC,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QACjC,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAA;IACxC,CAAC;IAED,4DAA4D;IAC5D,IAAI,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QAClC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAA;IAC9C,CAAC;IAED,4DAA4D;IAC5D,IAAI,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QAClC,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;IAC9C,CAAC;IAED,4DAA4D;IAC5D,IAAI,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS;QAClC,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAA;IACxC,CAAC;IAED,6EAA6E;IAE7E,0EAA0E;IAC1E,OAAO,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU;QACxC,OAAO,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,EAAE,EAAQ,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,EACpF,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAA;IACpD,CAAC;IAED,+DAA+D;IAC/D,OAAO,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU;QACxC,OAAO,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,EAAE,EAAQ,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,EACpF,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,CAAA;IACpD,CAAC;IAED,+DAA+D;IAC/D,OAAO,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU;QACxC,OAAO,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,EAAE,EAAE,GAAG,IAAI,CAAC,EAAE,EAC9E,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAA;IACpD,CAAC;IAEO,OAAO,CAAC,GAAiB,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EACrD,EAAU,EAAE,EAAU,EAAE,EAAU;QAChD,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAA;QAChB,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;QAE9C,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,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;IAE7E,eAAe;QACb,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;QAC9C,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;IACnD,CAAC;IAED,cAAc;QACZ,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;QACrB,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;QACrB,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAA;IACvB,CAAC;IAED,gBAAgB;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;gBAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAA;QACnF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;gBAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAA;QACnF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC;gBAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAA;IACrF,CAAC;IAED,gEAAgE;IAChE,SAAS,CAAC,OAAqD;QAC7D,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,IAAI,CAAA;QAC3B,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,KAAK,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,EAAE,GAAC,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,EAAE,GAAC,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,EAAE,GAAC,CAAC,CAAA;oBACrF,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;gBACvE,CAAC;IACH,CAAC;IAED,gFAAgF;IAChF,eAAe;QACb,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,IAAI,CAAA;QAC3B,kCAAkC;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBACzD,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAA;gBAC/B,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAA;YACjC,CAAC;QACD,mBAAmB;QACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBACzD,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAA;gBAC/B,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAA;YACjC,CAAC;QACD,mBAAmB;QACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBACzD,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAE,CAAC,GAAG,CAAC,CAAA;gBAC/B,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAA;YACjC,CAAC;IACH,CAAC;IAED,2DAA2D;IAC3D,YAAY,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU;QACzD,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,IAAI,CAAA;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE;gBAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;oBACpF,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,CAAA;QACvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;gBAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE;oBACpF,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,CAAA;QACvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE;YAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;gBAAE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;oBACpF,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,CAAA;IACzC,CAAC;CACF"}

package/dist/MarchingCubes.d.ts

@@ -0,0 +1,40 @@
+import { BufferGeometry } from 'three';
+export interface MarchingCubesOptions {
+    /** Grid resolution per side. Higher = more detail. Default 32. */
+    resolution?: number;
+    /** Iso-surface threshold. Default 0.5. */
+    isoLevel?: number;
+}
+/**
+ * MarchingCubes — extracts an iso-surface from a scalar field.
+ * Used to turn SPH particle density fields into a mesh.
+ * Blender: Fluid > Surface.
+ *
+ * Usage:
+ *   const mc = new MarchingCubes({ resolution: 32 })
+ *   mc.setField(myDensityFn)   // (x,y,z) => density
+ *   const geo = mc.extract()
+ */
+export declare class MarchingCubes {
+    parameters: {
+        resolution: number;
+        isoLevel: number;
+    };
+    private _field;
+    private _boundsMin;
+    private _boundsMax;
+    constructor(opts?: MarchingCubesOptions);
+    /**
+     * Set the scalar field function: (x,y,z) → density.
+     * For SPH, use a metaball kernel summed over nearby particles.
+     */
+    setField(fn: (x: number, y: number, z: number) => number): void;
+    /** Set world-space bounds of the sampling volume. */
+    setBounds(min: [number, number, number], max: [number, number, number]): void;
+    /**
+     * Sample the field and extract the iso-surface as a BufferGeometry.
+     * Returns null if no field is set.
+     */
+    extract(): BufferGeometry | null;
+}
+//# sourceMappingURL=MarchingCubes.d.ts.map

package/dist/MarchingCubes.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"MarchingCubes.d.ts","sourceRoot":"","sources":["../src/MarchingCubes.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,cAAc,EAAmB,MAAM,OAAO,CAAA;AAyTvD,MAAM,WAAW,oBAAoB;IACnC,kEAAkE;IAClE,UAAU,CAAC,EAAE,MAAM,CAAA;IACnB,0CAA0C;IAC1C,QAAQ,CAAC,EAAE,MAAM,CAAA;CAClB;AAED;;;;;;;;;GASG;AACH,qBAAa,aAAa;IACxB,UAAU,EAAE;QAAE,UAAU,EAAE,MAAM,CAAC;QAAC,QAAQ,EAAE,MAAM,CAAA;KAAE,CAAA;IAEpD,OAAO,CAAC,MAAM,CAA6D;IAC3E,OAAO,CAAC,UAAU,CAAqC;IACvD,OAAO,CAAC,UAAU,CAAqC;gBAE3C,IAAI,GAAE,oBAAyB;IAO3C;;;OAGG;IACH,QAAQ,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,MAAM,GAAG,IAAI;IAI/D,qDAAqD;IACrD,SAAS,CAAC,GAAG,EAAE,CAAC,MAAM,EAAC,MAAM,EAAC,MAAM,CAAC,EAAE,GAAG,EAAE,CAAC,MAAM,EAAC,MAAM,EAAC,MAAM,CAAC,GAAG,IAAI;IAKzE;;;OAGG;IACH,OAAO,IAAI,cAAc,GAAG,IAAI;CA8EjC"}