@johnhenry/math-sample 0.0.0-stage → 0.1.0

package/dist/complex.js

@@ -0,0 +1,236 @@
+import { assertCount, assertRange, Meter, progress } from "./budget.js";
+/** Sample `f` at the centre of each cell of a `width` x `height` grid over the rectangle. */
+export function sampleComplexRaster(f, reDomain, imDomain, options = {}) {
+    const width = options.width ?? 128;
+    const height = options.height ?? 128;
+    assertCount("width", width, 1);
+    assertCount("height", height, 1);
+    assertRange("re domain", reDomain);
+    assertRange("im domain", imDomain);
+    const meter = new Meter(options, options);
+    const re = new Float32Array(width * height).fill(Number.NaN);
+    const im = new Float32Array(width * height).fill(Number.NaN);
+    const abs = new Float32Array(width * height).fill(Number.NaN);
+    const arg = new Float32Array(width * height).fill(Number.NaN);
+    let rows = 0;
+    let reason;
+    for (let j = 0; j < height; j++) {
+        reason = meter.check(width);
+        if (reason)
+            break;
+        const y = imDomain.min + ((imDomain.max - imDomain.min) * (j + 0.5)) / height;
+        for (let i = 0; i < width; i++) {
+            const x = reDomain.min + ((reDomain.max - reDomain.min) * (i + 0.5)) / width;
+            const w = f(x, y);
+            const k = j * width + i;
+            if (Number.isFinite(w[0]) && Number.isFinite(w[1])) {
+                re[k] = w[0];
+                im[k] = w[1];
+                abs[k] = Math.hypot(w[0], w[1]);
+                arg[k] = Math.atan2(w[1], w[0]);
+            }
+        }
+        meter.spend(width);
+        rows++;
+    }
+    return {
+        width,
+        height: rows,
+        re: re.slice(0, width * rows),
+        im: im.slice(0, width * rows),
+        abs: abs.slice(0, width * rows),
+        arg: arg.slice(0, width * rows),
+        progress: progress(rows, height, reason),
+    };
+}
+/** The image under `f` of a rectangular grid of lines: where the grid lines go, drawn as the conformal-map picture. */
+export function conformalGrid(f, reDomain, imDomain, options = {}) {
+    const nRe = options.reLines ?? 11;
+    const nIm = options.imLines ?? 11;
+    const n = options.samplesPerLine ?? 200;
+    const bound = options.bound ?? 1e6;
+    assertCount("reLines", nRe, 1);
+    assertCount("imLines", nIm, 1);
+    assertCount("samplesPerLine", n);
+    assertRange("re domain", reDomain);
+    assertRange("im domain", imDomain);
+    const lines = [];
+    const trace = (kind, value, d) => {
+        const xs = [];
+        const ys = [];
+        let pending = false;
+        let px = Number.NaN;
+        let py = Number.NaN;
+        for (let k = 0; k < n; k++) {
+            const s = k === n - 1 ? d.max : d.min + ((d.max - d.min) * k) / (n - 1);
+            const w = kind === "re" ? f(value, s) : f(s, value);
+            const good = Number.isFinite(w[0]) && Number.isFinite(w[1]) && Math.hypot(w[0], w[1]) <= bound;
+            if (!good) {
+                pending = xs.length > 0;
+                px = Number.NaN;
+                continue;
+            }
+            const cut = options.maxStep !== undefined && !Number.isNaN(px) && Math.hypot(w[0] - px, w[1] - py) > options.maxStep;
+            if ((pending || cut) && xs.length > 0) {
+                xs.push(Number.NaN);
+                ys.push(Number.NaN);
+            }
+            xs.push(w[0]);
+            ys.push(w[1]);
+            px = w[0];
+            py = w[1];
+            pending = false;
+        }
+        return { kind, value, x: Float64Array.from(xs), y: Float64Array.from(ys) };
+    };
+    for (let i = 0; i < nRe; i++)
+        lines.push(trace("re", nRe === 1 ? reDomain.min : reDomain.min + ((reDomain.max - reDomain.min) * i) / (nRe - 1), imDomain));
+    for (let j = 0; j < nIm; j++)
+        lines.push(trace("im", nIm === 1 ? imDomain.min : imDomain.min + ((imDomain.max - imDomain.min) * j) / (nIm - 1), reDomain));
+    return lines;
+}
+const wrap = (d) => {
+    while (d > Math.PI)
+        d -= 2 * Math.PI;
+    while (d < -Math.PI)
+        d += 2 * Math.PI;
+    return d;
+};
+/**
+ * Zeros and poles of `f` inside a rectangle, by the argument principle: the winding number of `f` around each
+ * grid cell is (zeros - poles) inside it; non-zero cells are refined with Newton's method (on `f` for a
+ * zero, on `1/f` for a pole, scaled by the multiplicity). A zero and a pole inside one cell cancel and are missed:
+ * raise `grid`. The search grid is offset by an irrational fraction of a cell and extended past the
+ * rectangle, so round-number domains do not put roots on cell edges; results inside the rectangle (or on its edge) are
+ * kept and results within the tolerance are merged. `progress.rows` counts grid rows searched.
+ */
+export function findZerosAndPoles(f, reDomain, imDomain, options = {}) {
+    const grid = options.grid ?? 24;
+    assertCount("grid", grid, 1);
+    assertRange("re domain", reDomain);
+    assertRange("im domain", imDomain);
+    const size = Math.max(reDomain.max - reDomain.min, imDomain.max - imDomain.min);
+    const tol = (options.tolerance ?? 1e-10) * size;
+    const meter = new Meter(options, options);
+    const dx = (reDomain.max - reDomain.min) / grid;
+    const dy = (imDomain.max - imDomain.min) / grid;
+    const found = [];
+    const argAt = (x, y) => {
+        const w = f(x, y);
+        if (w[0] === 0 && w[1] === 0)
+            return Number.NaN; // an exact zero has no argument
+        return Math.atan2(w[1], w[0]);
+    };
+    // winding of f along segment p->q, subdividing until the phase moves < 1 rad per piece
+    const segment = (x0, y0, a0, x1, y1, a1, depth) => {
+        const d = wrap(a1 - a0);
+        // a minimum depth keeps a high-order root from aliasing a whole turn (3 theta moving 2 pi in one piece)
+        if ((depth >= 4 && Math.abs(d) < 1) || depth >= 14)
+            return d;
+        const xm = (x0 + x1) / 2;
+        const ym = (y0 + y1) / 2;
+        const am = argAt(xm, ym);
+        return segment(x0, y0, a0, xm, ym, am, depth + 1) + segment(xm, ym, am, x1, y1, a1, depth + 1);
+    };
+    const winding = (x0, y0, x1, y1) => {
+        const pts = [
+            [x0, y0],
+            [x1, y0],
+            [x1, y1],
+            [x0, y1],
+            [x0, y0],
+        ];
+        const args = pts.map(([x, y]) => argAt(x, y));
+        if (args.some((a) => !Number.isFinite(a)))
+            return Number.NaN;
+        let total = 0;
+        for (let k = 0; k < 4; k++) {
+            const p = pts[k];
+            const q = pts[k + 1];
+            total += segment(p[0], p[1], args[k], q[0], q[1], args[k + 1], 0);
+        }
+        return Math.round(total / (2 * Math.PI));
+    };
+    const newton = (x, y, kind, order, box) => {
+        const g = (a, b) => {
+            const w = f(a, b);
+            if (kind === "zero")
+                return [w[0], w[1]];
+            const m = w[0] * w[0] + w[1] * w[1];
+            return [w[0] / m, -w[1] / m]; // 1/f
+        };
+        let zx = x;
+        let zy = y;
+        for (let it = 0; it < 60; it++) {
+            const [gr, gi] = g(zx, zy);
+            if (!Number.isFinite(gr) || !Number.isFinite(gi))
+                return undefined;
+            const h = 1e-7 * Math.max(1, Math.abs(zx), Math.abs(zy));
+            const [ar, ai] = g(zx + h, zy);
+            const dr = (ar - gr) / h; // derivative of an analytic function = d/dx
+            const di = (ai - gi) / h;
+            const den = dr * dr + di * di;
+            if (!(den > 0) || !Number.isFinite(den))
+                return undefined;
+            // step = order * g / g'
+            const sr = (order * (gr * dr + gi * di)) / den;
+            const si = (order * (gi * dr - gr * di)) / den;
+            zx -= sr;
+            zy -= si;
+            if (Math.hypot(sr, si) < tol) {
+                const pad = 0.5 * Math.max(dx, dy);
+                if (zx < box[0] - pad || zx > box[1] + pad || zy < box[2] - pad || zy > box[3] + pad)
+                    return undefined;
+                return [zx, zy];
+            }
+            if (!Number.isFinite(zx) || !Number.isFinite(zy))
+                return undefined;
+        }
+        return undefined;
+    };
+    // The search grid is shifted by an irrational fraction of a cell and extended by one cell on each side, so
+    // roots and poles on a round-number domain (z^3 - z on [-2, 2]) sit inside cells, not on their edges or
+    // corners. Candidates outside the requested rectangle are discarded at the end.
+    const OFFSET = 0.3819660112501051;
+    const gx0 = reDomain.min - (1 - OFFSET) * dx;
+    const gy0 = imDomain.min - (1 - OFFSET) * dy;
+    const cells = grid + 1;
+    const slack = 1e-9 * size;
+    const insideRect = (x, y) => x >= reDomain.min - slack && x <= reDomain.max + slack && y >= imDomain.min - slack && y <= imDomain.max + slack;
+    let rows = 0;
+    let reason;
+    for (let j = 0; j < cells; j++) {
+        reason = meter.check(cells * 4);
+        if (reason)
+            break;
+        for (let i = 0; i < cells; i++) {
+            const x0 = gx0 + i * dx;
+            const y0 = gy0 + j * dy;
+            let w = winding(x0, y0, x0 + dx, y0 + dy);
+            let ex = 0;
+            if (Number.isNaN(w)) {
+                // f is zero or non-finite on this cell's boundary (an exact root/pole on an edge): retry slightly inflated
+                ex = 0.0123 * dx;
+                w = winding(x0 - ex, y0 - 0.0123 * dy, x0 + dx + ex, y0 + dy + 0.0123 * dy);
+            }
+            if (!Number.isFinite(w) || w === 0)
+                continue;
+            const kind = w > 0 ? "zero" : "pole";
+            const order = Math.abs(w);
+            const cx = x0 + dx / 2;
+            const cy = y0 + dy / 2;
+            const refined = newton(cx, cy, kind, order, [x0, x0 + dx, y0, y0 + dy]);
+            const [rx, ry] = refined ?? [cx, cy];
+            if (!insideRect(rx, ry))
+                continue;
+            if (found.some((p) => p.kind === kind && Math.hypot(p.re - rx, p.im - ry) < Math.max(tol * 1e3, 1e-7 * size)))
+                continue;
+            found.push({ re: rx, im: ry, kind, order, refined: refined !== undefined });
+        }
+        meter.spend(cells * 4);
+        rows++;
+    }
+    found.sort((a, b) => a.re - b.re || a.im - b.im);
+    return { items: found, progress: progress(rows, cells, reason) };
+}
+//# sourceMappingURL=complex.js.map

package/dist/complex.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"complex.js","sourceRoot":"","sources":["../src/complex.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,WAAW,EAAE,WAAW,EAAiC,KAAK,EAAiB,QAAQ,EAAE,MAAM,aAAa,CAAC;AA8BtH,6FAA6F;AAC7F,MAAM,UAAU,mBAAmB,CACjC,CAAY,EACZ,QAAgB,EAChB,QAAgB,EAChB,UAAgC,EAAE;IAElC,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,IAAI,GAAG,CAAC;IACnC,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,IAAI,GAAG,CAAC;IACrC,WAAW,CAAC,OAAO,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC;IAC/B,WAAW,CAAC,QAAQ,EAAE,MAAM,EAAE,CAAC,CAAC,CAAC;IACjC,WAAW,CAAC,WAAW,EAAE,QAAQ,CAAC,CAAC;IACnC,WAAW,CAAC,WAAW,EAAE,QAAQ,CAAC,CAAC;IACnC,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAC1C,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC7D,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC7D,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC9D,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC9D,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,MAAoC,CAAC;IACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QAC5B,IAAI,MAAM;YAAE,MAAM;QAClB,MAAM,CAAC,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC,CAAC,QAAQ,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QAC9E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC,CAAC,QAAQ,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC;YAC7E,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAClB,MAAM,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC;YACxB,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;gBACnD,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACb,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACb,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAClC,CAAC;QACH,CAAC;QACD,KAAK,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QACnB,IAAI,EAAE,CAAC;IACT,CAAC;IACD,OAAO;QACL,KAAK;QACL,MAAM,EAAE,IAAI;QACZ,EAAE,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,GAAG,IAAI,CAAC;QAC7B,EAAE,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,GAAG,IAAI,CAAC;QAC7B,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,GAAG,IAAI,CAAC;QAC/B,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,GAAG,IAAI,CAAC;QAC/B,QAAQ,EAAE,QAAQ,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC;KACzC,CAAC;AACJ,CAAC;AAwBD,uHAAuH;AACvH,MAAM,UAAU,aAAa,CAC3B,CAAY,EACZ,QAAgB,EAChB,QAAgB,EAChB,UAAgC,EAAE;IAElC,MAAM,GAAG,GAAG,OAAO,CAAC,OAAO,IAAI,EAAE,CAAC;IAClC,MAAM,GAAG,GAAG,OAAO,CAAC,OAAO,IAAI,EAAE,CAAC;IAClC,MAAM,CAAC,GAAG,OAAO,CAAC,cAAc,IAAI,GAAG,CAAC;IACxC,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,IAAI,GAAG,CAAC;IACnC,WAAW,CAAC,SAAS,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;IAC/B,WAAW,CAAC,SAAS,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;IAC/B,WAAW,CAAC,gBAAgB,EAAE,CAAC,CAAC,CAAC;IACjC,WAAW,CAAC,WAAW,EAAE,QAAQ,CAAC,CAAC;IACnC,WAAW,CAAC,WAAW,EAAE,QAAQ,CAAC,CAAC;IACnC,MAAM,KAAK,GAAoB,EAAE,CAAC;IAClC,MAAM,KAAK,GAAG,CAAC,IAAiB,EAAE,KAAa,EAAE,CAAS,EAAiB,EAAE;QAC3E,MAAM,EAAE,GAAa,EAAE,CAAC;QACxB,MAAM,EAAE,GAAa,EAAE,CAAC;QACxB,IAAI,OAAO,GAAG,KAAK,CAAC;QACpB,IAAI,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC;QACpB,IAAI,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC;QACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACxE,MAAM,CAAC,GAAG,IAAI,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;YACpD,MAAM,IAAI,GAAG,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC;YAC/F,IAAI,CAAC,IAAI,EAAE,CAAC;gBACV,OAAO,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;gBACxB,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC;gBAChB,SAAS;YACX,CAAC;YACD,MAAM,GAAG,GACP,OAAO,CAAC,OAAO,KAAK,SAAS,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,OAAO,CAAC,OAAO,CAAC;YAC3G,IAAI,CAAC,OAAO,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;gBACtC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;gBACpB,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;YACtB,CAAC;YACD,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACd,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACd,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACV,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACV,OAAO,GAAG,KAAK,CAAC;QAClB,CAAC;QACD,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,CAAC;IAC7E,CAAC,CAAC;IACF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;QAC1B,KAAK,CAAC,IAAI,CACR,KAAK,CAAC,IAAI,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,GAAG,CAAC,CAAC,QAAQ,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CACjH,CAAC;IACJ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;QAC1B,KAAK,CAAC,IAAI,CACR,KAAK,CAAC,IAAI,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,GAAG,CAAC,CAAC,QAAQ,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CACjH,CAAC;IACJ,OAAO,KAAK,CAAC;AACf,CAAC;AAwBD,MAAM,IAAI,GAAG,CAAC,CAAS,EAAE,EAAE;IACzB,OAAO,CAAC,GAAG,IAAI,CAAC,EAAE;QAAE,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IACrC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE;QAAE,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IACtC,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,UAAU,iBAAiB,CAC/B,CAAY,EACZ,QAAgB,EAChB,QAAgB,EAChB,UAA6B,EAAE;IAE/B,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI,IAAI,EAAE,CAAC;IAChC,WAAW,CAAC,MAAM,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC;IAC7B,WAAW,CAAC,WAAW,EAAE,QAAQ,CAAC,CAAC;IACnC,WAAW,CAAC,WAAW,EAAE,QAAQ,CAAC,CAAC;IACnC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,EAAE,QAAQ,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;IAChF,MAAM,GAAG,GAAG,CAAC,OAAO,CAAC,SAAS,IAAI,KAAK,CAAC,GAAG,IAAI,CAAC;IAChD,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAC1C,MAAM,EAAE,GAAG,CAAC,QAAQ,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;IAChD,MAAM,EAAE,GAAG,CAAC,QAAQ,CAAC,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;IAChD,MAAM,KAAK,GAAiB,EAAE,CAAC;IAE/B,MAAM,KAAK,GAAG,CAAC,CAAS,EAAE,CAAS,EAAU,EAAE;QAC7C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAClB,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC,gCAAgC;QACjF,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChC,CAAC,CAAC;IACF,uFAAuF;IACvF,MAAM,OAAO,GAAG,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,KAAa,EAAU,EAAE;QAChH,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACxB,wGAAwG;QACxG,IAAI,CAAC,KAAK,IAAI,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,KAAK,IAAI,EAAE;YAAE,OAAO,CAAC,CAAC;QAC7D,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;QACzB,OAAO,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC;IACjG,CAAC,CAAC;IACF,MAAM,OAAO,GAAG,CAAC,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAU,EAAE;QACzE,MAAM,GAAG,GAA4B;YACnC,CAAC,EAAE,EAAE,EAAE,CAAC;YACR,CAAC,EAAE,EAAE,EAAE,CAAC;YACR,CAAC,EAAE,EAAE,EAAE,CAAC;YACR,CAAC,EAAE,EAAE,EAAE,CAAC;YACR,CAAC,EAAE,EAAE,EAAE,CAAC;SACT,CAAC;QACF,MAAM,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC9C,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QAC7D,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,CAAqB,CAAC;YACrC,MAAM,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAqB,CAAC;YACzC,KAAK,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAW,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAW,EAAE,CAAC,CAAC,CAAC;QACxF,CAAC;QACD,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;IAC3C,CAAC,CAAC;IAEF,MAAM,MAAM,GAAG,CACb,CAAS,EACT,CAAS,EACT,IAAqB,EACrB,KAAa,EACb,GAAqC,EACrC,EAAE;QACF,MAAM,CAAC,GAAG,CAAC,CAAS,EAAE,CAAS,EAAoB,EAAE;YACnD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAClB,IAAI,IAAI,KAAK,MAAM;gBAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACzC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACpC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM;QACtC,CAAC,CAAC;QACF,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAC3B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;gBAAE,OAAO,SAAS,CAAC;YACnE,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YACzD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC;YAC/B,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,4CAA4C;YACtE,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;YACzB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;YAC9B,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC;gBAAE,OAAO,SAAS,CAAC;YAC1D,wBAAwB;YACxB,MAAM,EAAE,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC;YAC/C,MAAM,EAAE,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC;YAC/C,EAAE,IAAI,EAAE,CAAC;YACT,EAAE,IAAI,EAAE,CAAC;YACT,IAAI,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC;gBAC7B,MAAM,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;gBACnC,IAAI,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG;oBAAE,OAAO,SAAS,CAAC;gBACvG,OAAO,CAAC,EAAE,EAAE,EAAE,CAAU,CAAC;YAC3B,CAAC;YACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC;gBAAE,OAAO,SAAS,CAAC;QACrE,CAAC;QACD,OAAO,SAAS,CAAC;IACnB,CAAC,CAAC;IAEF,2GAA2G;IAC3G,wGAAwG;IACxG,gFAAgF;IAChF,MAAM,MAAM,GAAG,kBAAkB,CAAC;IAClC,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC;IAC7C,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC;IAC7C,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,CAAC;IACvB,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC;IAC1B,MAAM,UAAU,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAC1C,CAAC,IAAI,QAAQ,CAAC,GAAG,GAAG,KAAK,IAAI,CAAC,IAAI,QAAQ,CAAC,GAAG,GAAG,KAAK,IAAI,CAAC,IAAI,QAAQ,CAAC,GAAG,GAAG,KAAK,IAAI,CAAC,IAAI,QAAQ,CAAC,GAAG,GAAG,KAAK,CAAC;IACnH,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,MAAoC,CAAC;IACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAChC,IAAI,MAAM;YAAE,MAAM;QAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;YACxB,MAAM,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;YACxB,IAAI,CAAC,GAAG,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;YAC1C,IAAI,EAAE,GAAG,CAAC,CAAC;YACX,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;gBACpB,2GAA2G;gBAC3G,EAAE,GAAG,MAAM,GAAG,EAAE,CAAC;gBACjB,CAAC,GAAG,OAAO,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,MAAM,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,MAAM,GAAG,EAAE,CAAC,CAAC;YAC9E,CAAC;YACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC;gBAAE,SAAS;YAC7C,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC;YACrC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC1B,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YACvB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YACvB,MAAM,OAAO,GAAG,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC;YACxE,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,OAAO,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YACrC,IAAI,CAAC,UAAU,CAAC,EAAE,EAAE,EAAE,CAAC;gBAAE,SAAS;YAClC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,KAAK,IAAI,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,EAAE,IAAI,GAAG,IAAI,CAAC,CAAC;gBAC3G,SAAS;YACX,KAAK,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,OAAO,KAAK,SAAS,EAAE,CAAC,CAAC;QAC9E,CAAC;QACD,KAAK,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QACvB,IAAI,EAAE,CAAC;IACT,CAAC;IACD,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;IACjD,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,QAAQ,EAAE,QAAQ,CAAC,IAAI,EAAE,KAAK,EAAE,MAAM,CAAC,EAAE,CAAC;AACnE,CAAC"}

package/dist/expr.d.ts

@@ -0,0 +1,38 @@
+import { type Expr } from "@johnhenry/math";
+import type { Fn1 } from "./sample.ts";
+export interface CompileOptions {
+    /** The free variable that becomes the argument. Default `"x"`. */
+    variable?: string;
+    /** Fixed values for every other free variable (slider values, say). */
+    params?: Record<string, number>;
+    /**
+     * Read juxtaposition as multiplication in source text (`2x`, `3sin(x)`, `(x+1)(x-1)`), the way people
+     * type maths. Default `false`. Needs `@johnhenry/math` 0.1.0 or later; multi-letter variable names are
+     * unavailable when it is on (see `ParseOptions` in `@johnhenry/math`). Ignored when `expr` is already a parsed `Expr`.
+     */
+    implicitMultiplication?: boolean;
+}
+/**
+ * Compile an `@johnhenry/math` expression (source text or a parsed `Expr`)
+ * into a plain `(x) => number` the samplers take. The AST is walked once; the
+ * closure is cheap to call thousands of times. This is the only module that
+ * needs `@johnhenry/math` (an optional peer dependency), so a Worker can
+ * receive an expression string plus params, compile it here, and sample.
+ *
+ * Source text uses `@johnhenry/math`'s grammar: write `2*x`, or pass
+ * `implicitMultiplication: true` to accept `2x`.
+ */
+export declare function compileExpression(expr: Expr | string, options?: CompileOptions): Fn1;
+/**
+ * Taylor coefficients `a_k = f^(k)(x0) / k!` for `k = 0..order` of an expression. They are computed by truncated
+ * power-series ("jet") arithmetic, O(order^2) per node, so order 12 of `sin(x)*exp(x)/(1+x^2)` takes well under a
+ * millisecond. Supported: `+ - * /`, powers (integer or real exponent, or a variable exponent), and `exp ln sqrt cbrt sin
+ * cos tan sinh cosh tanh atan asin acos`. Any other node falls back to repeated symbolic differentiation with
+ * simplification after every derivative, which grows quickly: keep `order` at or below 8 there. A coefficient that is
+ * not finite at `x0` (a pole, `ln(0)`) is `NaN` or `Infinity`, and the partial sums built on it are gaps.
+ */
+export declare function taylorCoefficients(expr: Expr | string, x0: number, order: number, options?: {
+    variable?: string;
+    params?: Record<string, number>;
+}): Float64Array;
+//# sourceMappingURL=expr.d.ts.map

package/dist/expr.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"expr.d.ts","sourceRoot":"","sources":["../src/expr.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,IAAI,EAAY,MAAM,iBAAiB,CAAC;AACtD,OAAO,KAAK,EAAE,GAAG,EAAE,MAAM,aAAa,CAAC;AAEvC,MAAM,WAAW,cAAc;IAC7B,kEAAkE;IAClE,QAAQ,CAAC,EAAE,MAAM,CAAC;IAClB,uEAAuE;IACvE,MAAM,CAAC,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAChC;;;;OAIG;IACH,sBAAsB,CAAC,EAAE,OAAO,CAAC;CAClC;AAED;;;;;;;;;GASG;AACH,wBAAgB,iBAAiB,CAAC,IAAI,EAAE,IAAI,GAAG,MAAM,EAAE,OAAO,GAAE,cAAmB,GAAG,GAAG,CAYxF;AAqND;;;;;;;GAOG;AACH,wBAAgB,kBAAkB,CAChC,IAAI,EAAE,IAAI,GAAG,MAAM,EACnB,EAAE,EAAE,MAAM,EACV,KAAK,EAAE,MAAM,EACb,OAAO,GAAE;IAAE,QAAQ,CAAC,EAAE,MAAM,CAAC;IAAC,MAAM,CAAC,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAA;CAAO,GACnE,YAAY,CAkBd"}

package/dist/expr.js

@@ -0,0 +1,264 @@
+import { Symbolic } from "@johnhenry/math";
+/**
+ * Compile an `@johnhenry/math` expression (source text or a parsed `Expr`)
+ * into a plain `(x) => number` the samplers take. The AST is walked once; the
+ * closure is cheap to call thousands of times. This is the only module that
+ * needs `@johnhenry/math` (an optional peer dependency), so a Worker can
+ * receive an expression string plus params, compile it here, and sample.
+ *
+ * Source text uses `@johnhenry/math`'s grammar: write `2*x`, or pass
+ * `implicitMultiplication: true` to accept `2x`.
+ */
+export function compileExpression(expr, options = {}) {
+    const variable = options.variable ?? "x";
+    const tree = typeof expr === "string"
+        ? Symbolic.parse(expr, { implicitMultiplication: options.implicitMultiplication === true })
+        : expr;
+    const compiled = Symbolic.compile(tree);
+    const env = { ...options.params, [variable]: 0 };
+    return (x) => {
+        env[variable] = x;
+        return compiled(env);
+    };
+}
+class JetUnsupported extends Error {
+}
+function jmul(a, b) {
+    const n = a.length;
+    const c = new Float64Array(n);
+    for (let i = 0; i < n; i++) {
+        const ai = a[i];
+        if (ai === 0)
+            continue;
+        for (let j = 0; i + j < n; j++)
+            c[i + j] = c[i + j] + ai * b[j];
+    }
+    return c;
+}
+function jdiv(a, b) {
+    const n = a.length;
+    const c = new Float64Array(n);
+    for (let k = 0; k < n; k++) {
+        let acc = a[k];
+        for (let j = 1; j <= k; j++)
+            acc -= b[j] * c[k - j];
+        c[k] = acc / b[0];
+    }
+    return c;
+}
+const jconst = (v, n) => {
+    const j = new Float64Array(n);
+    j[0] = v;
+    return j;
+};
+function jexp(a) {
+    const n = a.length;
+    const e = new Float64Array(n);
+    e[0] = Math.exp(a[0]);
+    for (let k = 1; k < n; k++) {
+        let acc = 0;
+        for (let j = 1; j <= k; j++)
+            acc += j * a[j] * e[k - j];
+        e[k] = acc / k;
+    }
+    return e;
+}
+function jln(a) {
+    const n = a.length;
+    const l = new Float64Array(n);
+    l[0] = Math.log(a[0]);
+    for (let k = 1; k < n; k++) {
+        let acc = 0;
+        for (let j = 1; j < k; j++)
+            acc += j * l[j] * a[k - j];
+        l[k] = (a[k] - acc / k) / a[0];
+    }
+    return l;
+}
+/** sin/cos (hyperbolic = true for sinh/cosh) of a jet, together. */
+function jsincos(a, hyperbolic) {
+    const n = a.length;
+    const s = new Float64Array(n);
+    const c = new Float64Array(n);
+    s[0] = hyperbolic ? Math.sinh(a[0]) : Math.sin(a[0]);
+    c[0] = hyperbolic ? Math.cosh(a[0]) : Math.cos(a[0]);
+    for (let k = 1; k < n; k++) {
+        let ss = 0;
+        let cc = 0;
+        for (let j = 1; j <= k; j++) {
+            ss += j * a[j] * c[k - j];
+            cc += j * a[j] * s[k - j];
+        }
+        s[k] = ss / k;
+        c[k] = (hyperbolic ? cc : -cc) / k;
+    }
+    return [s, c];
+}
+function jpowReal(a, p) {
+    const n = a.length;
+    const f = new Float64Array(n);
+    const a0 = a[0];
+    f[0] = a0 ** p;
+    for (let k = 1; k < n; k++) {
+        let acc = 0;
+        for (let j = 1; j <= k; j++)
+            acc += (p * j - (k - j)) * a[j] * f[k - j];
+        f[k] = acc / (k * a0);
+    }
+    return f;
+}
+function jpowInt(a, m) {
+    let r = jconst(1, a.length);
+    let base = a;
+    let e = Math.abs(m);
+    while (e > 0) {
+        if (e & 1)
+            r = jmul(r, base);
+        base = jmul(base, base);
+        e = Math.floor(e / 2);
+    }
+    return m < 0 ? jdiv(jconst(1, a.length), r) : r;
+}
+/** Integrate a jet's derivative: given g = f', return f with f(0) = f0. */
+function jantiderivative(g, f0) {
+    const n = g.length;
+    const f = new Float64Array(n);
+    f[0] = f0;
+    for (let k = 1; k < n; k++)
+        f[k] = g[k - 1] / k;
+    return f;
+}
+function jderivative(a) {
+    const n = a.length;
+    const d = new Float64Array(n);
+    for (let k = 0; k + 1 < n; k++)
+        d[k] = (k + 1) * a[k + 1];
+    return d;
+}
+function jfunc(name, a) {
+    const n = a.length;
+    const one = jconst(1, n);
+    switch (name) {
+        case "exp":
+            return jexp(a);
+        case "ln":
+            return jln(a);
+        case "sqrt":
+            return jpowReal(a, 0.5);
+        case "cbrt":
+            return jpowReal(a, 1 / 3);
+        case "sin":
+            return jsincos(a, false)[0];
+        case "cos":
+            return jsincos(a, false)[1];
+        case "tan": {
+            const [s, c] = jsincos(a, false);
+            return jdiv(s, c);
+        }
+        case "sinh":
+            return jsincos(a, true)[0];
+        case "cosh":
+            return jsincos(a, true)[1];
+        case "tanh": {
+            const [s, c] = jsincos(a, true);
+            return jdiv(s, c);
+        }
+        case "atan": {
+            const d = jdiv(jderivative(a), jplus(one, jmul(a, a), 1));
+            return jantiderivative(d, Math.atan(a[0]));
+        }
+        case "asin":
+        case "acos": {
+            const inv = jpowReal(jplus(one, jmul(a, a), -1), -0.5);
+            const d = jmul(jderivative(a), inv);
+            return name === "asin"
+                ? jantiderivative(d, Math.asin(a[0]))
+                : jantiderivative(d.map((v) => -v), Math.acos(a[0]));
+        }
+        default:
+            throw new JetUnsupported(name);
+    }
+}
+function jplus(a, b, sign) {
+    const c = new Float64Array(a.length);
+    for (let k = 0; k < a.length; k++)
+        c[k] = a[k] + sign * b[k];
+    return c;
+}
+function evalJet(e, variable, x0, params, n) {
+    const go = (t) => {
+        switch (t.type) {
+            case "const":
+                return jconst(t.value, n);
+            case "var": {
+                if (t.name === variable) {
+                    const j = jconst(x0, n);
+                    if (n > 1)
+                        j[1] = 1;
+                    return j;
+                }
+                return jconst(params[t.name] ?? Number.NaN, n);
+            }
+            case "neg":
+                return go(t.arg).map((v) => -v);
+            case "add":
+                return jplus(go(t.left), go(t.right), 1);
+            case "sub":
+                return jplus(go(t.left), go(t.right), -1);
+            case "mul":
+                return jmul(go(t.left), go(t.right));
+            case "div":
+                return jdiv(go(t.left), go(t.right));
+            case "func":
+                return jfunc(t.name, go(t.arg));
+            case "pow": {
+                const b = go(t.base);
+                const ex = go(t.exp);
+                const constExp = ex.every((v, k) => k === 0 || v === 0);
+                if (constExp) {
+                    const p = ex[0];
+                    if (Number.isInteger(p) && Math.abs(p) <= 256)
+                        return jpowInt(b, p);
+                    return jpowReal(b, p);
+                }
+                return jexp(jmul(ex, jln(b)));
+            }
+            default:
+                throw new JetUnsupported(t.type);
+        }
+    };
+    return go(e);
+}
+/**
+ * Taylor coefficients `a_k = f^(k)(x0) / k!` for `k = 0..order` of an expression. They are computed by truncated
+ * power-series ("jet") arithmetic, O(order^2) per node, so order 12 of `sin(x)*exp(x)/(1+x^2)` takes well under a
+ * millisecond. Supported: `+ - * /`, powers (integer or real exponent, or a variable exponent), and `exp ln sqrt cbrt sin
+ * cos tan sinh cosh tanh atan asin acos`. Any other node falls back to repeated symbolic differentiation with
+ * simplification after every derivative, which grows quickly: keep `order` at or below 8 there. A coefficient that is
+ * not finite at `x0` (a pole, `ln(0)`) is `NaN` or `Infinity`, and the partial sums built on it are gaps.
+ */
+export function taylorCoefficients(expr, x0, order, options = {}) {
+    if (!Number.isInteger(order) || order < 0)
+        throw new RangeError("order must be a non-negative integer");
+    const variable = options.variable ?? "x";
+    const tree = typeof expr === "string" ? Symbolic.parse(expr) : expr;
+    try {
+        return evalJet(tree, variable, x0, options.params ?? {}, order + 1);
+    }
+    catch (err) {
+        if (!(err instanceof JetUnsupported))
+            throw err;
+    }
+    let term = tree;
+    const out = new Float64Array(order + 1);
+    let factorial = 1;
+    for (let k = 0; k <= order; k++) {
+        if (k > 0)
+            factorial *= k;
+        out[k] = Symbolic.evaluate(term, { ...options.params, [variable]: x0 }) / factorial;
+        if (k < order)
+            term = Symbolic.simplify(Symbolic.differentiate(term, variable));
+    }
+    return out;
+}
+//# sourceMappingURL=expr.js.map

package/dist/expr.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"expr.js","sourceRoot":"","sources":["../src/expr.ts"],"names":[],"mappings":"AAAA,OAAO,EAAa,QAAQ,EAAE,MAAM,iBAAiB,CAAC;AAgBtD;;;;;;;;;GASG;AACH,MAAM,UAAU,iBAAiB,CAAC,IAAmB,EAAE,UAA0B,EAAE;IACjF,MAAM,QAAQ,GAAG,OAAO,CAAC,QAAQ,IAAI,GAAG,CAAC;IACzC,MAAM,IAAI,GACR,OAAO,IAAI,KAAK,QAAQ;QACtB,CAAC,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,EAAE,EAAE,sBAAsB,EAAE,OAAO,CAAC,sBAAsB,KAAK,IAAI,EAAE,CAAC;QAC3F,CAAC,CAAC,IAAI,CAAC;IACX,MAAM,QAAQ,GAAG,QAAQ,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;IACxC,MAAM,GAAG,GAA2B,EAAE,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,CAAC;IACzE,OAAO,CAAC,CAAC,EAAE,EAAE;QACX,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;QAClB,OAAO,QAAQ,CAAC,GAAG,CAAC,CAAC;IACvB,CAAC,CAAC;AACJ,CAAC;AAID,MAAM,cAAe,SAAQ,KAAK;CAAG;AAErC,SAAS,IAAI,CAAC,CAAM,EAAE,CAAM;IAC1B,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAW,CAAC;QAC1B,IAAI,EAAE,KAAK,CAAC;YAAE,SAAS;QACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,GAAG,EAAE,GAAI,CAAC,CAAC,CAAC,CAAY,CAAC;IAC1F,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,IAAI,CAAC,CAAM,EAAE,CAAM;IAC1B,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAW,CAAC;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,IAAK,CAAC,CAAC,CAAC,CAAY,GAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,CAAC;QAC5E,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAI,CAAC,CAAC,CAAC,CAAY,CAAC;IAChC,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,MAAM,GAAG,CAAC,CAAS,EAAE,CAAS,EAAO,EAAE;IAC3C,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACT,OAAO,CAAC,CAAC;AACX,CAAC,CAAC;AAEF,SAAS,IAAI,CAAC,CAAM;IAClB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,IAAI,CAAC,GAAI,CAAC,CAAC,CAAC,CAAY,GAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,CAAC;QAChF,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;IACjB,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,GAAG,CAAC,CAAM;IACjB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,IAAI,CAAC,GAAI,CAAC,CAAC,CAAC,CAAY,GAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,CAAC;QAC/E,CAAC,CAAC,CAAC,CAAC,GAAG,CAAE,CAAC,CAAC,CAAC,CAAY,GAAG,GAAG,GAAG,CAAC,CAAC,GAAI,CAAC,CAAC,CAAC,CAAY,CAAC;IACzD,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,oEAAoE;AACpE,SAAS,OAAO,CAAC,CAAM,EAAE,UAAmB;IAC1C,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC;IACzE,CAAC,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC;IACzE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,EAAE,IAAI,CAAC,GAAI,CAAC,CAAC,CAAC,CAAY,GAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,CAAC;YAClD,EAAE,IAAI,CAAC,GAAI,CAAC,CAAC,CAAC,CAAY,GAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,CAAC;QACpD,CAAC;QACD,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACd,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;IACrC,CAAC;IACD,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AAChB,CAAC;AAED,SAAS,QAAQ,CAAC,CAAM,EAAE,CAAS;IACjC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAW,CAAC;IAC1B,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAI,CAAC,CAAC,CAAC,CAAY,GAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,CAAC;QAChG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;IACxB,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,OAAO,CAAC,CAAM,EAAE,CAAS;IAChC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC;IAC5B,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpB,OAAO,CAAC,GAAG,CAAC,EAAE,CAAC;QACb,IAAI,CAAC,GAAG,CAAC;YAAE,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;QAC7B,IAAI,GAAG,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QACxB,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACxB,CAAC;IACD,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAClD,CAAC;AAED,2EAA2E;AAC3E,SAAS,eAAe,CAAC,CAAM,EAAE,EAAU;IACzC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,CAAC,CAAC,GAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,GAAG,CAAC,CAAC;IAC5D,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,WAAW,CAAC,CAAM;IACzB,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAY,CAAC;IACtE,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,KAAK,CAAC,IAAY,EAAE,CAAM;IACjC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;IACnB,MAAM,GAAG,GAAG,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACzB,QAAQ,IAAI,EAAE,CAAC;QACb,KAAK,KAAK;YACR,OAAO,IAAI,CAAC,CAAC,CAAC,CAAC;QACjB,KAAK,IAAI;YACP,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC;QAChB,KAAK,MAAM;YACT,OAAO,QAAQ,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAC1B,KAAK,MAAM;YACT,OAAO,QAAQ,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5B,KAAK,KAAK;YACR,OAAO,OAAO,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9B,KAAK,KAAK;YACR,OAAO,OAAO,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC9B,KAAK,KAAK,CAAC,CAAC,CAAC;YACX,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;YACjC,OAAO,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACpB,CAAC;QACD,KAAK,MAAM;YACT,OAAO,OAAO,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,KAAK,MAAM;YACT,OAAO,OAAO,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7B,KAAK,MAAM,CAAC,CAAC,CAAC;YACZ,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;YAChC,OAAO,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACpB,CAAC;QACD,KAAK,MAAM,CAAC,CAAC,CAAC;YACZ,MAAM,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC1D,OAAO,eAAe,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC,CAAC;QACvD,CAAC;QACD,KAAK,MAAM,CAAC;QACZ,KAAK,MAAM,CAAC,CAAC,CAAC;YACZ,MAAM,GAAG,GAAG,QAAQ,CAAC,KAAK,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC;YACvD,MAAM,CAAC,GAAG,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;YACpC,OAAO,IAAI,KAAK,MAAM;gBACpB,CAAC,CAAC,eAAe,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC;gBAC/C,CAAC,CAAC,eAAe,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAQ,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC,CAAC,CAAC;QAC1E,CAAC;QACD;YACE,MAAM,IAAI,cAAc,CAAC,IAAI,CAAC,CAAC;IACnC,CAAC;AACH,CAAC;AAED,SAAS,KAAK,CAAC,CAAM,EAAE,CAAM,EAAE,IAAY;IACzC,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,CAAC,CAAC,GAAI,CAAC,CAAC,CAAC,CAAY,GAAG,IAAI,GAAI,CAAC,CAAC,CAAC,CAAY,CAAC;IACrF,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,OAAO,CAAC,CAAO,EAAE,QAAgB,EAAE,EAAU,EAAE,MAA8B,EAAE,CAAS;IAC/F,MAAM,EAAE,GAAG,CAAC,CAAO,EAAO,EAAE;QAC1B,QAAQ,CAAC,CAAC,IAAI,EAAE,CAAC;YACf,KAAK,OAAO;gBACV,OAAO,MAAM,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;YAC5B,KAAK,KAAK,CAAC,CAAC,CAAC;gBACX,IAAI,CAAC,CAAC,IAAI,KAAK,QAAQ,EAAE,CAAC;oBACxB,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;oBACxB,IAAI,CAAC,GAAG,CAAC;wBAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;oBACpB,OAAO,CAAC,CAAC;gBACX,CAAC;gBACD,OAAO,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YACjD,CAAC;YACD,KAAK,KAAK;gBACR,OAAO,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAQ,CAAC;YACzC,KAAK,KAAK;gBACR,OAAO,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;YAC3C,KAAK,KAAK;gBACR,OAAO,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC5C,KAAK,KAAK;gBACR,OAAO,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC;YACvC,KAAK,KAAK;gBACR,OAAO,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC;YACvC,KAAK,MAAM;gBACT,OAAO,KAAK,CAAC,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAClC,KAAK,KAAK,CAAC,CAAC,CAAC;gBACX,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;gBACrB,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACrB,MAAM,QAAQ,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;gBACxD,IAAI,QAAQ,EAAE,CAAC;oBACb,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAW,CAAC;oBAC1B,IAAI,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG;wBAAE,OAAO,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;oBACpE,OAAO,QAAQ,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACxB,CAAC;gBACD,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChC,CAAC;YACD;gBACE,MAAM,IAAI,cAAc,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QACrC,CAAC;IACH,CAAC,CAAC;IACF,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC;AACf,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,kBAAkB,CAChC,IAAmB,EACnB,EAAU,EACV,KAAa,EACb,UAAkE,EAAE;IAEpE,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,GAAG,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,sCAAsC,CAAC,CAAC;IACxG,MAAM,QAAQ,GAAG,OAAO,CAAC,QAAQ,IAAI,GAAG,CAAC;IACzC,MAAM,IAAI,GAAS,OAAO,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;IAC1E,IAAI,CAAC;QACH,OAAO,OAAO,CAAC,IAAI,EAAE,QAAQ,EAAE,EAAE,EAAE,OAAO,CAAC,MAAM,IAAI,EAAE,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC;IACtE,CAAC;IAAC,OAAO,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,CAAC,GAAG,YAAY,cAAc,CAAC;YAAE,MAAM,GAAG,CAAC;IAClD,CAAC;IACD,IAAI,IAAI,GAAS,IAAI,CAAC;IACtB,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IACxC,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAChC,IAAI,CAAC,GAAG,CAAC;YAAE,SAAS,IAAI,CAAC,CAAC;QAC1B,GAAG,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,QAAQ,CAAC,IAAI,EAAE,EAAE,GAAG,OAAO,CAAC,MAAM,EAAE,CAAC,QAAQ,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,SAAS,CAAC;QACpF,IAAI,CAAC,GAAG,KAAK;YAAE,IAAI,GAAG,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,aAAa,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC,CAAC;IAClF,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC"}

package/dist/field2d.d.ts

@@ -0,0 +1,58 @@
+import { type Budget, type Cancellable, type Progress } from "./budget.ts";
+import type { Domain, Samples } from "./sample.ts";
+export type Fn2 = (x: number, y: number) => number;
+export interface Grid2Options extends Budget, Cancellable {
+    /** Cells per side of the sampling grid (`grid + 1` nodes per side). Default 128. */
+    grid?: number | [number, number];
+}
+export interface ImplicitCurve {
+    /**
+     * Segments as rows: each segment is two rows (its endpoints) followed by one `NaN` separator row, so
+     * `Samples` gap rules hold and `runsOf` yields one run per segment (join them if you wish).
+     */
+    segments: Samples;
+    segmentCount: number;
+    progress: Progress;
+}
+/**
+ * The curve `f(x, y) = level` (default 0) by marching squares on a regular grid. A saddle cell is joined by its
+ * centre value. Cells with a non-finite corner are skipped (a hole, not a spike). With a budget,
+ * `progress.rows` counts grid rows finished.
+ */
+export declare function sampleImplicit(f: Fn2, xDomain: Domain, yDomain: Domain, options?: Grid2Options & {
+    level?: number;
+}): ImplicitCurve;
+export interface Contours {
+    levels: number[];
+    /** One entry per level, in `levels` order. */
+    curves: Array<{
+        level: number;
+        segments: Samples;
+        segmentCount: number;
+    }>;
+    progress: Progress;
+}
+/**
+ * Contour lines of `f` at the given `levels`, or `count` (default 9) levels evenly spread strictly between the
+ * grid's least and greatest values. One pass over the grid for all levels.
+ */
+export declare function sampleContours(f: Fn2, xDomain: Domain, yDomain: Domain, options?: Grid2Options & {
+    levels?: number[];
+    count?: number;
+}): Contours;
+export interface GradientField {
+    columns: number;
+    rows: number;
+    /** Arrow tails. */
+    x: Float64Array;
+    y: Float64Array;
+    /** Components; `NaN` where either component is non-finite. */
+    u: Float64Array;
+    v: Float64Array;
+    magnitude: Float64Array;
+    maxMagnitude: number;
+    progress: Progress;
+}
+/** Sample the vector field `(fx, fy)` (e.g. the two partial derivatives of a scalar field) on a grid including the edges. */
+export declare function sampleGradient(fx: Fn2, fy: Fn2, xDomain: Domain, yDomain: Domain, options?: Grid2Options): GradientField;
+//# sourceMappingURL=field2d.d.ts.map

package/dist/field2d.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"field2d.d.ts","sourceRoot":"","sources":["../src/field2d.ts"],"names":[],"mappings":"AAAA,OAAO,EAA4B,KAAK,MAAM,EAAE,KAAK,WAAW,EAAS,KAAK,QAAQ,EAAY,MAAM,aAAa,CAAC;AACtH,OAAO,KAAK,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,aAAa,CAAC;AAEnD,MAAM,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,MAAM,CAAC;AAEnD,MAAM,WAAW,YAAa,SAAQ,MAAM,EAAE,WAAW;IACvD,oFAAoF;IACpF,IAAI,CAAC,EAAE,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;CAClC;AAED,MAAM,WAAW,aAAa;IAC5B;;;OAGG;IACH,QAAQ,EAAE,OAAO,CAAC;IAClB,YAAY,EAAE,MAAM,CAAC;IACrB,QAAQ,EAAE,QAAQ,CAAC;CACpB;AA0GD;;;;GAIG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,GAAG,EACN,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,OAAO,GAAE,YAAY,GAAG;IAAE,KAAK,CAAC,EAAE,MAAM,CAAA;CAAO,GAC9C,aAAa,CAIf;AAED,MAAM,WAAW,QAAQ;IACvB,MAAM,EAAE,MAAM,EAAE,CAAC;IACjB,8CAA8C;IAC9C,MAAM,EAAE,KAAK,CAAC;QAAE,KAAK,EAAE,MAAM,CAAC;QAAC,QAAQ,EAAE,OAAO,CAAC;QAAC,YAAY,EAAE,MAAM,CAAA;KAAE,CAAC,CAAC;IAC1E,QAAQ,EAAE,QAAQ,CAAC;CACpB;AAED;;;GAGG;AACH,wBAAgB,cAAc,CAC5B,CAAC,EAAE,GAAG,EACN,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,OAAO,GAAE,YAAY,GAAG;IAAE,MAAM,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,KAAK,CAAC,EAAE,MAAM,CAAA;CAAO,GACjE,QAAQ,CA8BV;AAED,MAAM,WAAW,aAAa;IAC5B,OAAO,EAAE,MAAM,CAAC;IAChB,IAAI,EAAE,MAAM,CAAC;IACb,mBAAmB;IACnB,CAAC,EAAE,YAAY,CAAC;IAChB,CAAC,EAAE,YAAY,CAAC;IAChB,8DAA8D;IAC9D,CAAC,EAAE,YAAY,CAAC;IAChB,CAAC,EAAE,YAAY,CAAC;IAChB,SAAS,EAAE,YAAY,CAAC;IACxB,YAAY,EAAE,MAAM,CAAC;IACrB,QAAQ,EAAE,QAAQ,CAAC;CACpB;AAED,6HAA6H;AAC7H,wBAAgB,cAAc,CAC5B,EAAE,EAAE,GAAG,EACP,EAAE,EAAE,GAAG,EACP,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,OAAO,GAAE,YAAiB,GACzB,aAAa,CAiDf"}