@johnhenry/math-sample 0.0.1 → 0.1.0

package/dist/expr.js

@@ -21,4 +21,244 @@ export function compileExpression(expr, options = {}) {
         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

@@ -1 +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"}
+{"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"}

package/dist/field2d.js

@@ -0,0 +1,199 @@
+import { assertCount, assertRange, Meter, progress } from "./budget.js";
+function dims(g) {
+    const [nx, ny] = typeof g === "number" ? [g, g] : (g ?? [128, 128]);
+    assertCount("grid", nx, 1);
+    assertCount("grid", ny, 1);
+    return [nx, ny];
+}
+function marching(f, xDomain, yDomain, levels, options) {
+    const [nx, ny] = dims(options.grid);
+    assertRange("x domain", xDomain);
+    assertRange("y domain", yDomain);
+    const meter = new Meter(options, options);
+    const hx = (xDomain.max - xDomain.min) / nx;
+    const hy = (yDomain.max - yDomain.min) / ny;
+    const segs = levels.map(() => []);
+    let prev = new Float64Array(nx + 1);
+    let cur = new Float64Array(nx + 1);
+    const row = (j, out) => {
+        const y = yDomain.min + j * hy;
+        for (let i = 0; i <= nx; i++)
+            out[i] = f(xDomain.min + i * hx, y);
+    };
+    let reason;
+    row(0, prev);
+    meter.spend(nx + 1);
+    let done = 0;
+    for (let j = 0; j < ny; j++) {
+        reason = meter.check(nx + 1);
+        if (reason)
+            break;
+        row(j + 1, cur);
+        meter.spend(nx + 1);
+        const y0 = yDomain.min + j * hy;
+        for (let i = 0; i < nx; i++) {
+            const x0 = xDomain.min + i * hx;
+            const v = [prev[i], prev[i + 1], cur[i + 1], cur[i]]; // corners ccw from bottom-left
+            if (!v.every(Number.isFinite))
+                continue; // holes: no segments through undefined cells
+            const cx = [x0, x0 + hx, x0 + hx, x0];
+            const cy = [y0, y0, y0 + hy, y0 + hy];
+            levels.forEach((level, li) => {
+                const s = v.map((q) => q - level);
+                // edge e joins corner e and corner e+1
+                const cross = (e) => {
+                    const a = s[e];
+                    const b = s[(e + 1) % 4];
+                    if (a < 0 === b < 0)
+                        return undefined; // a === 0 counts as >= 0 on the positive side
+                    const t = a / (a - b);
+                    const ea = e;
+                    const eb = (e + 1) % 4;
+                    return [
+                        cx[ea] + t * (cx[eb] - cx[ea]),
+                        cy[ea] + t * (cy[eb] - cy[ea]),
+                    ];
+                };
+                const pts = [0, 1, 2, 3].map(cross);
+                const idx = [0, 1, 2, 3].filter((e) => pts[e]);
+                const out = segs[li];
+                const emit = (p, q) => out.push(p[0], p[1], q[0], q[1]);
+                if (idx.length === 2)
+                    emit(pts[idx[0]], pts[idx[1]]);
+                else if (idx.length === 4) {
+                    // saddle: the cell centre decides which corners are connected
+                    const centre = (s[0] + s[1] + s[2] + s[3]) / 4;
+                    const p = pts;
+                    const bottomLeftPositive = s[0] >= 0;
+                    const joinsCorner0 = centre >= 0 === bottomLeftPositive; // centre agrees with corner 0: corners 0 and 2 are linked through the middle, so the curve cuts off corners 1 and 3
+                    if (joinsCorner0) {
+                        emit(p[0], p[1]);
+                        emit(p[2], p[3]);
+                    }
+                    else {
+                        emit(p[3], p[0]);
+                        emit(p[1], p[2]);
+                    }
+                }
+            });
+        }
+        [prev, cur] = [cur, prev];
+        done++;
+    }
+    return { segs, progress: progress(done, ny, reason) };
+}
+function toSamples(flat) {
+    const n = flat.length / 4;
+    const x = new Float64Array(n === 0 ? 0 : n * 3 - 1);
+    const y = new Float64Array(x.length);
+    for (let s = 0; s < n; s++) {
+        const o = s * 3;
+        x[o] = flat[4 * s];
+        y[o] = flat[4 * s + 1];
+        x[o + 1] = flat[4 * s + 2];
+        y[o + 1] = flat[4 * s + 3];
+        if (s < n - 1) {
+            x[o + 2] = Number.NaN;
+            y[o + 2] = Number.NaN;
+        }
+    }
+    return { x, y };
+}
+/**
+ * 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 function sampleImplicit(f, xDomain, yDomain, options = {}) {
+    const { segs, progress: p } = marching(f, xDomain, yDomain, [options.level ?? 0], options);
+    const flat = segs[0];
+    return { segments: toSamples(flat), segmentCount: flat.length / 4, progress: p };
+}
+/**
+ * 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 function sampleContours(f, xDomain, yDomain, options = {}) {
+    let levels = options.levels;
+    if (!levels) {
+        const [nx, ny] = dims(options.grid);
+        let lo = Infinity;
+        let hi = -Infinity;
+        for (let j = 0; j <= ny; j++)
+            for (let i = 0; i <= nx; i++) {
+                const v = f(xDomain.min + ((xDomain.max - xDomain.min) * i) / nx, yDomain.min + ((yDomain.max - yDomain.min) * j) / ny);
+                if (Number.isFinite(v)) {
+                    if (v < lo)
+                        lo = v;
+                    if (v > hi)
+                        hi = v;
+                }
+            }
+        const count = options.count ?? 9;
+        levels = lo < hi ? Array.from({ length: count }, (_, k) => lo + ((hi - lo) * (k + 1)) / (count + 1)) : [];
+    }
+    const { segs, progress: p } = marching(f, xDomain, yDomain, levels, options);
+    return {
+        levels,
+        curves: levels.map((level, k) => ({
+            level,
+            segments: toSamples(segs[k]),
+            segmentCount: segs[k].length / 4,
+        })),
+        progress: p,
+    };
+}
+/** Sample the vector field `(fx, fy)` (e.g. the two partial derivatives of a scalar field) on a grid including the edges. */
+export function sampleGradient(fx, fy, xDomain, yDomain, options = {}) {
+    const [cols, rws] = typeof options.grid === "number" ? [options.grid, options.grid] : (options.grid ?? [16, 16]);
+    assertCount("grid", cols);
+    assertCount("grid", rws);
+    assertRange("x domain", xDomain);
+    assertRange("y domain", yDomain);
+    const meter = new Meter(options, options);
+    const n = cols * rws;
+    const X = new Float64Array(n);
+    const Y = new Float64Array(n);
+    const U = new Float64Array(n).fill(Number.NaN);
+    const V = new Float64Array(n).fill(Number.NaN);
+    const M = new Float64Array(n).fill(Number.NaN);
+    let done = 0;
+    let reason;
+    let max = 0;
+    for (let j = 0; j < rws; j++) {
+        reason = meter.check(cols);
+        if (reason)
+            break;
+        const y = j === rws - 1 ? yDomain.max : yDomain.min + ((yDomain.max - yDomain.min) * j) / (rws - 1);
+        for (let i = 0; i < cols; i++) {
+            const x = i === cols - 1 ? xDomain.max : xDomain.min + ((xDomain.max - xDomain.min) * i) / (cols - 1);
+            const k = j * cols + i;
+            X[k] = x;
+            Y[k] = y;
+            const a = fx(x, y);
+            const b = fy(x, y);
+            if (Number.isFinite(a) && Number.isFinite(b)) {
+                U[k] = a;
+                V[k] = b;
+                M[k] = Math.hypot(a, b);
+                if (M[k] > max)
+                    max = M[k];
+            }
+        }
+        meter.spend(cols);
+        done++;
+    }
+    const m = cols * done;
+    return {
+        columns: cols,
+        rows: done,
+        x: X.slice(0, m),
+        y: Y.slice(0, m),
+        u: U.slice(0, m),
+        v: V.slice(0, m),
+        magnitude: M.slice(0, m),
+        maxMagnitude: max,
+        progress: progress(done, rws, reason),
+    };
+}
+//# sourceMappingURL=field2d.js.map

package/dist/field2d.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"field2d.js","sourceRoot":"","sources":["../src/field2d.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,WAAW,EAAE,WAAW,EAAiC,KAAK,EAAiB,QAAQ,EAAE,MAAM,aAAa,CAAC;AAoBtH,SAAS,IAAI,CAAC,CAAwC;IACpD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,OAAO,CAAC,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;IACpE,WAAW,CAAC,MAAM,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IAC3B,WAAW,CAAC,MAAM,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IAC3B,OAAO,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;AAClB,CAAC;AAED,SAAS,QAAQ,CACf,CAAM,EACN,OAAe,EACf,OAAe,EACf,MAAgB,EAChB,OAAqB;IAErB,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;IACpC,WAAW,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IACjC,WAAW,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAC1C,MAAM,EAAE,GAAG,CAAC,OAAO,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;IAC5C,MAAM,EAAE,GAAG,CAAC,OAAO,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC;IAC5C,MAAM,IAAI,GAAe,MAAM,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;IAC9C,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpC,IAAI,GAAG,GAAG,IAAI,YAAY,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACnC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,GAAiB,EAAE,EAAE;QAC3C,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,CAAC;IACpE,CAAC,CAAC;IACF,IAAI,MAAoC,CAAC;IACzC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IACb,KAAK,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IACpB,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC7B,IAAI,MAAM;YAAE,MAAM;QAClB,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;QAChB,KAAK,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QACpB,MAAM,EAAE,GAAG,OAAO,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,EAAE,GAAG,OAAO,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;YAChC,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAa,CAAC,CAAC,+BAA+B;YACjG,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;gBAAE,SAAS,CAAC,6CAA6C;YACtF,MAAM,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;YACtC,MAAM,CAAC,OAAO,CAAC,CAAC,KAAK,EAAE,EAAE,EAAE,EAAE;gBAC3B,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;gBAClC,uCAAuC;gBACvC,MAAM,KAAK,GAAG,CAAC,CAAS,EAAgC,EAAE;oBACxD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAW,CAAC;oBACzB,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAW,CAAC;oBACnC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC;wBAAE,OAAO,SAAS,CAAC,CAAC,8CAA8C;oBACrF,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;oBACtB,MAAM,EAAE,GAAG,CAAC,CAAC;oBACb,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;oBACvB,OAAO;wBACJ,EAAE,CAAC,EAAE,CAAY,GAAG,CAAC,GAAG,CAAE,EAAE,CAAC,EAAE,CAAY,GAAI,EAAE,CAAC,EAAE,CAAY,CAAC;wBACjE,EAAE,CAAC,EAAE,CAAY,GAAG,CAAC,GAAG,CAAE,EAAE,CAAC,EAAE,CAAY,GAAI,EAAE,CAAC,EAAE,CAAY,CAAC;qBACnE,CAAC;gBACJ,CAAC,CAAC;gBACF,MAAM,GAAG,GAAwC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;gBACzE,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/C,MAAM,GAAG,GAAG,IAAI,CAAC,EAAE,CAAa,CAAC;gBACjC,MAAM,IAAI,GAAG,CAAC,CAAmB,EAAE,CAAmB,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5F,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC;oBAClB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAW,CAAqB,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAW,CAAqB,CAAC,CAAC;qBACxF,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;oBAC1B,8DAA8D;oBAC9D,MAAM,MAAM,GAAG,CAAE,CAAC,CAAC,CAAC,CAAY,GAAI,CAAC,CAAC,CAAC,CAAY,GAAI,CAAC,CAAC,CAAC,CAAY,GAAI,CAAC,CAAC,CAAC,CAAY,CAAC,GAAG,CAAC,CAAC;oBAC/F,MAAM,CAAC,GAAG,GAA8B,CAAC;oBACzC,MAAM,kBAAkB,GAAI,CAAC,CAAC,CAAC,CAAY,IAAI,CAAC,CAAC;oBACjD,MAAM,YAAY,GAAG,MAAM,IAAI,CAAC,KAAK,kBAAkB,CAAC,CAAC,oHAAoH;oBAC7K,IAAI,YAAY,EAAE,CAAC;wBACjB,IAAI,CAAC,CAAC,CAAC,CAAC,CAAqB,EAAE,CAAC,CAAC,CAAC,CAAqB,CAAC,CAAC;wBACzD,IAAI,CAAC,CAAC,CAAC,CAAC,CAAqB,EAAE,CAAC,CAAC,CAAC,CAAqB,CAAC,CAAC;oBAC3D,CAAC;yBAAM,CAAC;wBACN,IAAI,CAAC,CAAC,CAAC,CAAC,CAAqB,EAAE,CAAC,CAAC,CAAC,CAAqB,CAAC,CAAC;wBACzD,IAAI,CAAC,CAAC,CAAC,CAAC,CAAqB,EAAE,CAAC,CAAC,CAAC,CAAqB,CAAC,CAAC;oBAC3D,CAAC;gBACH,CAAC;YACH,CAAC,CAAC,CAAC;QACL,CAAC;QACD,CAAC,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;QAC1B,IAAI,EAAE,CAAC;IACT,CAAC;IACD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,QAAQ,CAAC,IAAI,EAAE,EAAE,EAAE,MAAM,CAAC,EAAE,CAAC;AACxD,CAAC;AAED,SAAS,SAAS,CAAC,IAAc;IAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;IAC1B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACpD,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAChB,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAW,CAAC;QAC7B,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAW,CAAC;QACjC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAW,CAAC;QACrC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAW,CAAC;QACrC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC;YACd,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;YACtB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;QACxB,CAAC;IACH,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAClB,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,cAAc,CAC5B,CAAM,EACN,OAAe,EACf,OAAe,EACf,UAA6C,EAAE;IAE/C,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,CAAC,EAAE,GAAG,QAAQ,CAAC,CAAC,EAAE,OAAO,EAAE,OAAO,EAAE,CAAC,OAAO,CAAC,KAAK,IAAI,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;IAC3F,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,CAAa,CAAC;IACjC,OAAO,EAAE,QAAQ,EAAE,SAAS,CAAC,IAAI,CAAC,EAAE,YAAY,EAAE,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,CAAC;AACnF,CAAC;AASD;;;GAGG;AACH,MAAM,UAAU,cAAc,CAC5B,CAAM,EACN,OAAe,EACf,OAAe,EACf,UAAgE,EAAE;IAElE,IAAI,MAAM,GAAG,OAAO,CAAC,MAAM,CAAC;IAC5B,IAAI,CAAC,MAAM,EAAE,CAAC;QACZ,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;QACpC,IAAI,EAAE,GAAG,QAAQ,CAAC;QAClB,IAAI,EAAE,GAAG,CAAC,QAAQ,CAAC;QACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE;YAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,MAAM,CAAC,GAAG,CAAC,CACT,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EACpD,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CACrD,CAAC;gBACF,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;oBACvB,IAAI,CAAC,GAAG,EAAE;wBAAE,EAAE,GAAG,CAAC,CAAC;oBACnB,IAAI,CAAC,GAAG,EAAE;wBAAE,EAAE,GAAG,CAAC,CAAC;gBACrB,CAAC;YACH,CAAC;QACH,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,IAAI,CAAC,CAAC;QACjC,MAAM,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;IAC5G,CAAC;IACD,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,CAAC,EAAE,GAAG,QAAQ,CAAC,CAAC,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;IAC7E,OAAO;QACL,MAAM;QACN,MAAM,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;YAChC,KAAK;YACL,QAAQ,EAAE,SAAS,CAAC,IAAI,CAAC,CAAC,CAAa,CAAC;YACxC,YAAY,EAAG,IAAI,CAAC,CAAC,CAAc,CAAC,MAAM,GAAG,CAAC;SAC/C,CAAC,CAAC;QACH,QAAQ,EAAE,CAAC;KACZ,CAAC;AACJ,CAAC;AAgBD,6HAA6H;AAC7H,MAAM,UAAU,cAAc,CAC5B,EAAO,EACP,EAAO,EACP,OAAe,EACf,OAAe,EACf,UAAwB,EAAE;IAE1B,MAAM,CAAC,IAAI,EAAE,GAAG,CAAC,GAAG,OAAO,OAAO,CAAC,IAAI,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,IAAI,EAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,IAAI,IAAI,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IACjH,WAAW,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IAC1B,WAAW,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;IACzB,WAAW,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IACjC,WAAW,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAC1C,MAAM,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC;IACrB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAC9B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC/C,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC/C,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC/C,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,MAAoC,CAAC;IACzC,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC3B,IAAI,MAAM;YAAE,MAAM;QAClB,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;QACpG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAC9B,MAAM,CAAC,GAAG,CAAC,KAAK,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YACtG,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;YACvB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACT,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACT,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACnB,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACnB,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;gBAC7C,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACT,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACT,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACxB,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG;oBAAE,GAAG,GAAG,CAAC,CAAC,CAAC,CAAW,CAAC;YACvC,CAAC;QACH,CAAC;QACD,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAClB,IAAI,EAAE,CAAC;IACT,CAAC;IACD,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC;IACtB,OAAO;QACL,OAAO,EAAE,IAAI;QACb,IAAI,EAAE,IAAI;QACV,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;QAChB,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;QAChB,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;QAChB,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;QAChB,SAAS,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;QACxB,YAAY,EAAE,GAAG;QACjB,QAAQ,EAAE,QAAQ,CAAC,IAAI,EAAE,GAAG,EAAE,MAAM,CAAC;KACtC,CAAC;AACJ,CAAC"}

package/dist/index.d.ts

@@ -1,4 +1,10 @@
+export type { Budget, Cancellable, Progress } from "./budget.ts";
+export { type ComplexFn, type ComplexRaster, type ComplexRasterOptions, type ConformalGridOptions, type ConformalLine, conformalGrid, findZerosAndPoles, sampleComplexRaster, type ZeroOrPole, type ZerosPoles, type ZerosPolesOptions, } from "./complex.ts";
+export { type Contours, type GradientField, type Grid2Options, type ImplicitCurve, sampleContours, sampleGradient, sampleImplicit, } from "./field2d.ts";
+export { type Curve3, type CurveFn3, type CurveOptions, curve3Runs, type Fn2, type Mesh, type MeshOptions, type ParametricSurfaceFn, sampleParametricSurface, sampleSpaceCurve, sampleSurface, sampleVectorField3, type VectorField3, type VectorField3Options, type VectorFn3, } from "./mesh.ts";
+export { type AdaptiveOdeOptions, classifyFixedPoint, type Field2, type FixedPoint, type FixedPointType, integrateAdaptive, integrateRk4, jacobian2, type OdeOptions, type PhasePortrait, type PhasePortraitOptions, type Rk4Options, type Slope, type SlopeField, type SlopeFieldOptions, type SolutionCurveOptions, type StopReason, type System, samplePhasePortrait, sampleSlopeField, sampleSolutionCurve, type Trajectory, } from "./ode.ts";
 export { type AdaptiveOptions, type CancelSignal, type Domain, type Fn1, findCrossings, findDiscontinuities, findIntersections, findRoots, type Point, type RootOptions, runsOf, SampleCancelledError, type SampleOptions, type Samples, sampleAdaptive, sampleUniform, toPoints, type YRange, } from "./sample.ts";
+export { evaluateFourier, evaluateTaylor, type FourierCoefficients, fourierCoefficients, type NumericSeries, type PartialSumCurves, partialSums, sampleFourierPartialSums, sampleSeriesPartialSums, sampleTaylorPartialSums, } from "./series.ts";
 export { type AxesOptions, axesToSvg, type PathOptions, samplesToSvgPath, svgDocument, type XY } from "./svg.ts";
 export { formatTick, niceTicks, tickStep } from "./ticks.ts";
 export { backingStoreSize, fitViewport, panViewport, pinchZoomFactor, toDataX, toDataY, toScreenX, toScreenY, type Viewport, viewportFromAnchor, wheelZoomFactor, zoomViewport, } from "./viewport.ts";

package/dist/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EACL,KAAK,eAAe,EACpB,KAAK,YAAY,EACjB,KAAK,MAAM,EACX,KAAK,GAAG,EACR,aAAa,EACb,mBAAmB,EACnB,iBAAiB,EACjB,SAAS,EACT,KAAK,KAAK,EACV,KAAK,WAAW,EAChB,MAAM,EACN,oBAAoB,EACpB,KAAK,aAAa,EAClB,KAAK,OAAO,EACZ,cAAc,EACd,aAAa,EACb,QAAQ,EACR,KAAK,MAAM,GACZ,MAAM,aAAa,CAAC;AACrB,OAAO,EAAE,KAAK,WAAW,EAAE,SAAS,EAAE,KAAK,WAAW,EAAE,gBAAgB,EAAE,WAAW,EAAE,KAAK,EAAE,EAAE,MAAM,UAAU,CAAC;AACjH,OAAO,EAAE,UAAU,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAC7D,OAAO,EACL,gBAAgB,EAChB,WAAW,EACX,WAAW,EACX,eAAe,EACf,OAAO,EACP,OAAO,EACP,SAAS,EACT,SAAS,EACT,KAAK,QAAQ,EACb,kBAAkB,EAClB,eAAe,EACf,YAAY,GACb,MAAM,eAAe,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,YAAY,EAAE,MAAM,EAAE,WAAW,EAAE,QAAQ,EAAE,MAAM,aAAa,CAAC;AACjE,OAAO,EACL,KAAK,SAAS,EACd,KAAK,aAAa,EAClB,KAAK,oBAAoB,EACzB,KAAK,oBAAoB,EACzB,KAAK,aAAa,EAClB,aAAa,EACb,iBAAiB,EACjB,mBAAmB,EACnB,KAAK,UAAU,EACf,KAAK,UAAU,EACf,KAAK,iBAAiB,GACvB,MAAM,cAAc,CAAC;AACtB,OAAO,EACL,KAAK,QAAQ,EACb,KAAK,aAAa,EAClB,KAAK,YAAY,EACjB,KAAK,aAAa,EAClB,cAAc,EACd,cAAc,EACd,cAAc,GACf,MAAM,cAAc,CAAC;AACtB,OAAO,EACL,KAAK,MAAM,EACX,KAAK,QAAQ,EACb,KAAK,YAAY,EACjB,UAAU,EACV,KAAK,GAAG,EACR,KAAK,IAAI,EACT,KAAK,WAAW,EAChB,KAAK,mBAAmB,EACxB,uBAAuB,EACvB,gBAAgB,EAChB,aAAa,EACb,kBAAkB,EAClB,KAAK,YAAY,EACjB,KAAK,mBAAmB,EACxB,KAAK,SAAS,GACf,MAAM,WAAW,CAAC;AACnB,OAAO,EACL,KAAK,kBAAkB,EACvB,kBAAkB,EAClB,KAAK,MAAM,EACX,KAAK,UAAU,EACf,KAAK,cAAc,EACnB,iBAAiB,EACjB,YAAY,EACZ,SAAS,EACT,KAAK,UAAU,EACf,KAAK,aAAa,EAClB,KAAK,oBAAoB,EACzB,KAAK,UAAU,EACf,KAAK,KAAK,EACV,KAAK,UAAU,EACf,KAAK,iBAAiB,EACtB,KAAK,oBAAoB,EACzB,KAAK,UAAU,EACf,KAAK,MAAM,EACX,mBAAmB,EACnB,gBAAgB,EAChB,mBAAmB,EACnB,KAAK,UAAU,GAChB,MAAM,UAAU,CAAC;AAClB,OAAO,EACL,KAAK,eAAe,EACpB,KAAK,YAAY,EACjB,KAAK,MAAM,EACX,KAAK,GAAG,EACR,aAAa,EACb,mBAAmB,EACnB,iBAAiB,EACjB,SAAS,EACT,KAAK,KAAK,EACV,KAAK,WAAW,EAChB,MAAM,EACN,oBAAoB,EACpB,KAAK,aAAa,EAClB,KAAK,OAAO,EACZ,cAAc,EACd,aAAa,EACb,QAAQ,EACR,KAAK,MAAM,GACZ,MAAM,aAAa,CAAC;AACrB,OAAO,EACL,eAAe,EACf,cAAc,EACd,KAAK,mBAAmB,EACxB,mBAAmB,EACnB,KAAK,aAAa,EAClB,KAAK,gBAAgB,EACrB,WAAW,EACX,wBAAwB,EACxB,uBAAuB,EACvB,uBAAuB,GACxB,MAAM,aAAa,CAAC;AACrB,OAAO,EAAE,KAAK,WAAW,EAAE,SAAS,EAAE,KAAK,WAAW,EAAE,gBAAgB,EAAE,WAAW,EAAE,KAAK,EAAE,EAAE,MAAM,UAAU,CAAC;AACjH,OAAO,EAAE,UAAU,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAC7D,OAAO,EACL,gBAAgB,EAChB,WAAW,EACX,WAAW,EACX,eAAe,EACf,OAAO,EACP,OAAO,EACP,SAAS,EACT,SAAS,EACT,KAAK,QAAQ,EACb,kBAAkB,EAClB,eAAe,EACf,YAAY,GACb,MAAM,eAAe,CAAC"}

package/dist/index.js

@@ -1,4 +1,9 @@
+export { conformalGrid, findZerosAndPoles, sampleComplexRaster, } from "./complex.js";
+export { sampleContours, sampleGradient, sampleImplicit, } from "./field2d.js";
+export { curve3Runs, sampleParametricSurface, sampleSpaceCurve, sampleSurface, sampleVectorField3, } from "./mesh.js";
+export { classifyFixedPoint, integrateAdaptive, integrateRk4, jacobian2, samplePhasePortrait, sampleSlopeField, sampleSolutionCurve, } from "./ode.js";
 export { findCrossings, findDiscontinuities, findIntersections, findRoots, runsOf, SampleCancelledError, sampleAdaptive, sampleUniform, toPoints, } from "./sample.js";
+export { evaluateFourier, evaluateTaylor, fourierCoefficients, partialSums, sampleFourierPartialSums, sampleSeriesPartialSums, sampleTaylorPartialSums, } from "./series.js";
 export { axesToSvg, samplesToSvgPath, svgDocument } from "./svg.js";
 export { formatTick, niceTicks, tickStep } from "./ticks.js";
 export { backingStoreSize, fitViewport, panViewport, pinchZoomFactor, toDataX, toDataY, toScreenX, toScreenY, viewportFromAnchor, wheelZoomFactor, zoomViewport, } from "./viewport.js";

package/dist/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAKL,aAAa,EACb,mBAAmB,EACnB,iBAAiB,EACjB,SAAS,EAGT,MAAM,EACN,oBAAoB,EAGpB,cAAc,EACd,aAAa,EACb,QAAQ,GAET,MAAM,aAAa,CAAC;AACrB,OAAO,EAAoB,SAAS,EAAoB,gBAAgB,EAAE,WAAW,EAAW,MAAM,UAAU,CAAC;AACjH,OAAO,EAAE,UAAU,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAC7D,OAAO,EACL,gBAAgB,EAChB,WAAW,EACX,WAAW,EACX,eAAe,EACf,OAAO,EACP,OAAO,EACP,SAAS,EACT,SAAS,EAET,kBAAkB,EAClB,eAAe,EACf,YAAY,GACb,MAAM,eAAe,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AACA,OAAO,EAML,aAAa,EACb,iBAAiB,EACjB,mBAAmB,GAIpB,MAAM,cAAc,CAAC;AACtB,OAAO,EAKL,cAAc,EACd,cAAc,EACd,cAAc,GACf,MAAM,cAAc,CAAC;AACtB,OAAO,EAIL,UAAU,EAKV,uBAAuB,EACvB,gBAAgB,EAChB,aAAa,EACb,kBAAkB,GAInB,MAAM,WAAW,CAAC;AACnB,OAAO,EAEL,kBAAkB,EAIlB,iBAAiB,EACjB,YAAY,EACZ,SAAS,EAWT,mBAAmB,EACnB,gBAAgB,EAChB,mBAAmB,GAEpB,MAAM,UAAU,CAAC;AAClB,OAAO,EAKL,aAAa,EACb,mBAAmB,EACnB,iBAAiB,EACjB,SAAS,EAGT,MAAM,EACN,oBAAoB,EAGpB,cAAc,EACd,aAAa,EACb,QAAQ,GAET,MAAM,aAAa,CAAC;AACrB,OAAO,EACL,eAAe,EACf,cAAc,EAEd,mBAAmB,EAGnB,WAAW,EACX,wBAAwB,EACxB,uBAAuB,EACvB,uBAAuB,GACxB,MAAM,aAAa,CAAC;AACrB,OAAO,EAAoB,SAAS,EAAoB,gBAAgB,EAAE,WAAW,EAAW,MAAM,UAAU,CAAC;AACjH,OAAO,EAAE,UAAU,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAC7D,OAAO,EACL,gBAAgB,EAChB,WAAW,EACX,WAAW,EACX,eAAe,EACf,OAAO,EACP,OAAO,EACP,SAAS,EACT,SAAS,EAET,kBAAkB,EAClB,eAAe,EACf,YAAY,GACb,MAAM,eAAe,CAAC"}