diff-grok 1.1.0 → 1.3.0

package/dist/src/solver-tools/cvode/cvode_class.js

@@ -0,0 +1,163 @@
+/**
+ * TypeScript port of SUNDIALS CVODE v7.5.0 (https://github.com/LLNL/sundials)
+ *
+ * Copyright (c) 2025, Lawrence Livermore National Security, University of
+ *   Maryland Baltimore County, and the SUNDIALS contributors.
+ * Copyright (c) 2013-2025, Lawrence Livermore National Security and
+ *   Southern Methodist University.
+ * Copyright (c) 2002-2013, Lawrence Livermore National Security.
+ * Copyright (c) 2026 Datagrok.
+ *
+ * Licensed under the BSD-3-Clause License. See THIRD_PARTY_LICENSES
+ * for the full license text and provenance chain of the original code.
+ */
+// CVODE TypeScript port - High-level Cvode class API
+// Follows the Lsoda class pattern from src/lsoda.ts
+import { cvodeGetDky, CV_ADAMS, CV_BDF, CV_NORMAL, CV_ONE_STEP, CV_ROOT_RETURN, } from './common';
+import { cvodeCreate, cvodeInit, cvode } from './cvode';
+import { cvodeSetLinearSolver, cvodeSetJacFn } from './cvode_ls';
+import { cvodeRootInit } from './cvode_root';
+import { cvodeGetIntegratorStats, cvodeSStolerances, cvodeSVtolerances, cvodeSetMaxNumSteps, cvodeSetMaxOrd, cvodeSetMaxStep, cvodeSetMinStep, cvodeSetInitStep, cvodeSetStopTime, cvodeSetUserData, } from './cvode_io';
+// ============================================================================
+// Cvode class
+// ============================================================================
+/**
+ * High-level CVODE solver class with a 0-based public interface.
+ *
+ * Usage:
+ *   const solver = new Cvode(f, neq, t0, y0, { lmm: 'bdf', rtol: 1e-6, atol: 1e-9 });
+ *   const result = solver.solve(tout);
+ */
+export class Cvode {
+    mem;
+    neq;
+    /** Output vector reused across calls (0-based, neq elements). */
+    yout;
+    /** Wrapped RHS stored to allow re-initialization. */
+    wrappedF;
+    constructor(f, neq, t0, y0, options) {
+        this.neq = neq;
+        this.yout = new Float64Array(neq);
+        const opts = options ?? {};
+        // --- 1. Create CvodeMem ---
+        const lmmValue = (opts.lmm === 'adams') ? CV_ADAMS : CV_BDF;
+        this.mem = cvodeCreate(lmmValue);
+        // --- 2. Wrap user's OdeFunction to CvodeRhsFn (void -> number) ---
+        this.wrappedF = (_t, _y, _ydot, _userData) => {
+            f(_t, _y, _ydot);
+            return 0;
+        };
+        // --- 3. Initialize ---
+        cvodeInit(this.mem, this.wrappedF, t0, y0);
+        // --- 4. Set tolerances ---
+        const rtol = opts.rtol ?? 1e-4;
+        const atol = opts.atol ?? 1e-6;
+        if (typeof atol === 'number')
+            cvodeSStolerances(this.mem, rtol, atol);
+        else
+            cvodeSVtolerances(this.mem, rtol, atol);
+        // --- 5. Set up dense linear solver (required for BDF; optional but useful for Adams) ---
+        cvodeSetLinearSolver(this.mem);
+        // --- 6. Set user-supplied Jacobian if provided ---
+        if (opts.jacFn) {
+            const userJac = opts.jacFn;
+            const wrappedJacFn = (t, y, fy, J, userData) => {
+                userJac(t, y, fy, J);
+                return 0;
+            };
+            cvodeSetJacFn(this.mem, wrappedJacFn);
+        }
+        // --- 7. Apply optional settings ---
+        if (opts.maxSteps !== undefined)
+            cvodeSetMaxNumSteps(this.mem, opts.maxSteps);
+        if (opts.maxOrder !== undefined)
+            cvodeSetMaxOrd(this.mem, opts.maxOrder);
+        if (opts.maxStep !== undefined)
+            cvodeSetMaxStep(this.mem, opts.maxStep);
+        if (opts.minStep !== undefined)
+            cvodeSetMinStep(this.mem, opts.minStep);
+        if (opts.initStep !== undefined)
+            cvodeSetInitStep(this.mem, opts.initStep);
+        if (opts.stopTime !== undefined)
+            cvodeSetStopTime(this.mem, opts.stopTime);
+        if (opts.userData !== undefined)
+            cvodeSetUserData(this.mem, opts.userData);
+        // --- 8. Set up rootfinding if requested ---
+        if (opts.rootFn && opts.nRootFns && opts.nRootFns > 0) {
+            const userRoot = opts.rootFn;
+            const wrappedRootFn = (t, y, gout, _userData) => {
+                userRoot(t, y, gout);
+                return 0;
+            };
+            cvodeRootInit(this.mem, opts.nRootFns, wrappedRootFn);
+        }
+    }
+    // ==========================================================================
+    // Public methods
+    // ==========================================================================
+    /**
+     * Integrate from the current time to tout using CV_NORMAL mode
+     * (integrator advances to tout and interpolates).
+     *
+     * @param tout - target output time
+     * @returns CvodeSolveResult with { t, y, flag, rootsFound? }
+     */
+    solve(tout) {
+        const { flag, t } = cvode(this.mem, tout, this.yout, CV_NORMAL);
+        return this._buildResult(t, flag);
+    }
+    /**
+     * Take a single internal step using CV_ONE_STEP mode.
+     * The solver advances by one internal step of whatever size it chooses.
+     *
+     * @returns CvodeSolveResult with { t, y, flag, rootsFound? }
+     */
+    step() {
+        const { flag, t } = cvode(this.mem, 1e30, this.yout, CV_ONE_STEP);
+        return this._buildResult(t, flag);
+    }
+    /**
+     * Dense output: compute the k-th derivative of the interpolating polynomial
+     * at time t using the current Nordsieck history array.
+     *
+     * @param t - interpolation time (must be within the last step interval)
+     * @param k - derivative order (0 = solution, 1 = first derivative, ...)
+     * @returns copy of the interpolated vector (length neq)
+     */
+    getDky(t, k) {
+        const dky = new Float64Array(this.neq);
+        const flag = cvodeGetDky(this.mem, t, k, dky);
+        if (flag !== 0)
+            throw new Error(`[Cvode] getDky failed: t=${t} out of range or k=${k} invalid`);
+        return dky;
+    }
+    /**
+     * Return current integrator statistics.
+     */
+    getStats() {
+        return cvodeGetIntegratorStats(this.mem);
+    }
+    /**
+     * Re-initialize the solver for a new initial-value problem with the same
+     * ODE function and options. Resets internal state and counters.
+     *
+     * @param t0 - new initial time
+     * @param y0 - new initial conditions (length neq)
+     */
+    reInit(t0, y0) {
+        cvodeInit(this.mem, this.wrappedF, t0, y0);
+    }
+    // ==========================================================================
+    // Private helpers
+    // ==========================================================================
+    _buildResult(t, flag) {
+        const y = new Float64Array(this.yout); // copy
+        const result = { t, y, flag };
+        if (flag === CV_ROOT_RETURN) {
+            // Return a copy of the root direction info array
+            result.rootsFound = new Int32Array(this.mem.cv_iroots);
+        }
+        return result;
+    }
+}
+//# sourceMappingURL=cvode_class.js.map

package/dist/src/solver-tools/cvode/cvode_class.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cvode_class.js","sourceRoot":"","sources":["../../../../src/solver-tools/cvode/cvode_class.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,qDAAqD;AACrD,oDAAoD;AAEpD,OAAO,EAKL,WAAW,EACX,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,WAAW,EACtB,cAAc,GACf,MAAM,UAAU,CAAC;AAElB,OAAO,EAAC,WAAW,EAAE,SAAS,EAAE,KAAK,EAAC,MAAM,SAAS,CAAC;AACtD,OAAO,EAAC,oBAAoB,EAAE,aAAa,EAAC,MAAM,YAAY,CAAC;AAC/D,OAAO,EAAC,aAAa,EAAC,MAAM,cAAc,CAAC;AAC3C,OAAO,EAEL,uBAAuB,EACvB,iBAAiB,EACjB,iBAAiB,EACjB,mBAAmB,EACnB,cAAc,EACd,eAAe,EACf,eAAe,EACf,gBAAgB,EAChB,gBAAgB,EAChB,gBAAgB,GACjB,MAAM,YAAY,CAAC;AA0DpB,+EAA+E;AAC/E,cAAc;AACd,+EAA+E;AAE/E;;;;;;GAMG;AACH,MAAM,OAAO,KAAK;IACR,GAAG,CAAW;IACd,GAAG,CAAS;IACpB,iEAAiE;IACzD,IAAI,CAAe;IAC3B,qDAAqD;IAC7C,QAAQ,CAAa;IAE7B,YACE,CAAc,EACd,GAAW,EACX,EAAU,EACV,EAAgB,EAChB,OAAsB;QAEtB,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;QAElC,MAAM,IAAI,GAAG,OAAO,IAAI,EAAE,CAAC;QAE3B,6BAA6B;QAC7B,MAAM,QAAQ,GAAG,CAAC,IAAI,CAAC,GAAG,KAAK,OAAO,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,MAAM,CAAC;QAC5D,IAAI,CAAC,GAAG,GAAG,WAAW,CAAC,QAAQ,CAAC,CAAC;QAEjC,oEAAoE;QACpE,IAAI,CAAC,QAAQ,GAAG,CAAC,EAAU,EAAE,EAAgB,EAAE,KAAmB,EAAE,SAAc,EAAU,EAAE;YAC5F,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC;YACjB,OAAO,CAAC,CAAC;QACX,CAAC,CAAC;QAEF,wBAAwB;QACxB,SAAS,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;QAE3C,4BAA4B;QAC5B,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC;QAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC;QAE/B,IAAI,OAAO,IAAI,KAAK,QAAQ;YAC1B,iBAAiB,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;;YAExC,iBAAiB,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAG1C,0FAA0F;QAC1F,oBAAoB,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAE/B,oDAAoD;QACpD,IAAI,IAAI,CAAC,KAAK,EAAE,CAAC;YACf,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC;YAC3B,MAAM,YAAY,GAAe,CAC/B,CAAS,EAAE,CAAe,EAAE,EAAgB,EAC5C,CAAiB,EAAE,QAAa,EACxB,EAAE;gBACV,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;gBACrB,OAAO,CAAC,CAAC;YACX,CAAC,CAAC;YACF,aAAa,CAAC,IAAI,CAAC,GAAG,EAAE,YAAY,CAAC,CAAC;QACxC,CAAC;QAED,qCAAqC;QACrC,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS;YAAE,mBAAmB,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC9E,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS;YAAE,cAAc,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QACzE,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAAE,eAAe,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;QACxE,IAAI,IAAI,CAAC,OAAO,KAAK,SAAS;YAAE,eAAe,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,OAAO,CAAC,CAAC;QACxE,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS;YAAE,gBAAgB,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC3E,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS;YAAE,gBAAgB,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC3E,IAAI,IAAI,CAAC,QAAQ,KAAK,SAAS;YAAE,gBAAgB,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QAE3E,6CAA6C;QAC7C,IAAI,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,QAAQ,IAAI,IAAI,CAAC,QAAQ,GAAG,CAAC,EAAE,CAAC;YACtD,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC;YAC7B,MAAM,aAAa,GAAgB,CACjC,CAAS,EAAE,CAAe,EAAE,IAAkB,EAAE,SAAc,EACtD,EAAE;gBACV,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;gBACrB,OAAO,CAAC,CAAC;YACX,CAAC,CAAC;YACF,aAAa,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,EAAE,aAAa,CAAC,CAAC;QACxD,CAAC;IACH,CAAC;IAED,6EAA6E;IAC7E,iBAAiB;IACjB,6EAA6E;IAE7E;;;;;;OAMG;IACH,KAAK,CAAC,IAAY;QAChB,MAAM,EAAC,IAAI,EAAE,CAAC,EAAC,GAAG,KAAK,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,CAAC,IAAI,EAAE,SAAS,CAAC,CAAC;QAC9D,OAAO,IAAI,CAAC,YAAY,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IACpC,CAAC;IAED;;;;;OAKG;IACH,IAAI;QACF,MAAM,EAAC,IAAI,EAAE,CAAC,EAAC,GAAG,KAAK,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,CAAC,IAAI,EAAE,WAAW,CAAC,CAAC;QAChE,OAAO,IAAI,CAAC,YAAY,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;IACpC,CAAC;IAED;;;;;;;OAOG;IACH,MAAM,CAAC,CAAS,EAAE,CAAS;QACzB,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACvC,MAAM,IAAI,GAAG,WAAW,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QAC9C,IAAI,IAAI,KAAK,CAAC;YACZ,MAAM,IAAI,KAAK,CAAC,4BAA4B,CAAC,sBAAsB,CAAC,UAAU,CAAC,CAAC;QAElF,OAAO,GAAG,CAAC;IACb,CAAC;IAED;;OAEG;IACH,QAAQ;QACN,OAAO,uBAAuB,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IAC3C,CAAC;IAED;;;;;;OAMG;IACH,MAAM,CAAC,EAAU,EAAE,EAAgB;QACjC,SAAS,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAC7C,CAAC;IAED,6EAA6E;IAC7E,kBAAkB;IAClB,6EAA6E;IAErE,YAAY,CAAC,CAAS,EAAE,IAAY;QAC1C,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,OAAO;QAC9C,MAAM,MAAM,GAAqB,EAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAC,CAAC;QAE9C,IAAI,IAAI,KAAK,cAAc,EAAE,CAAC;YAC5B,iDAAiD;YACjD,MAAM,CAAC,UAAU,GAAG,IAAI,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;QACzD,CAAC;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;CACF"}

package/dist/src/solver-tools/cvode/cvode_hin.d.ts

@@ -0,0 +1,24 @@
+/**
+ * TypeScript port of SUNDIALS CVODE v7.5.0 (https://github.com/LLNL/sundials)
+ *
+ * Copyright (c) 2025, Lawrence Livermore National Security, University of
+ *   Maryland Baltimore County, and the SUNDIALS contributors.
+ * Copyright (c) 2013-2025, Lawrence Livermore National Security and
+ *   Southern Methodist University.
+ * Copyright (c) 2002-2013, Lawrence Livermore National Security.
+ * Copyright (c) 2026 Datagrok.
+ *
+ * Licensed under the BSD-3-Clause License. See THIRD_PARTY_LICENSES
+ * for the full license text and provenance chain of the original code.
+ */
+import { CvodeMem } from './common';
+/**
+ * cvHin - Compute initial step size h0.
+ *
+ * The algorithm:
+ * 1. Compute a bound hub on |h| from ydot and tdist
+ * 2. Set hlb = HLB_FACTOR * tround  (lower bound)
+ * 3. Iterate to estimate ||y''|| and compute h candidates
+ * 4. h0 = H_BIAS * geometric mean of hlb and hub candidates
+ */
+export declare function cvHin(mem: CvodeMem, tout: number): number;

package/dist/src/solver-tools/cvode/cvode_hin.js

@@ -0,0 +1,119 @@
+/**
+ * TypeScript port of SUNDIALS CVODE v7.5.0 (https://github.com/LLNL/sundials)
+ *
+ * Copyright (c) 2025, Lawrence Livermore National Security, University of
+ *   Maryland Baltimore County, and the SUNDIALS contributors.
+ * Copyright (c) 2013-2025, Lawrence Livermore National Security and
+ *   Southern Methodist University.
+ * Copyright (c) 2002-2013, Lawrence Livermore National Security.
+ * Copyright (c) 2026 Datagrok.
+ *
+ * Licensed under the BSD-3-Clause License. See THIRD_PARTY_LICENSES
+ * for the full license text and provenance chain of the original code.
+ */
+// CVODE TypeScript port - Initial step size estimation
+// Ported from SUNDIALS CVODE cvode.c (cvHin, cvUpperBoundH0, cvYddNorm)
+import { wrmsNorm, vLinearSum, vScale, CV_SUCCESS, CV_RHSFUNC_FAIL, RHSFUNC_RECVR, CV_TOO_CLOSE, HLB_FACTOR, HUB_FACTOR, H_BIAS, MAX_ITERS } from './common';
+/**
+ * cvHin - Compute initial step size h0.
+ *
+ * The algorithm:
+ * 1. Compute a bound hub on |h| from ydot and tdist
+ * 2. Set hlb = HLB_FACTOR * tround  (lower bound)
+ * 3. Iterate to estimate ||y''|| and compute h candidates
+ * 4. h0 = H_BIAS * geometric mean of hlb and hub candidates
+ */
+export function cvHin(mem, tout) {
+    const N = mem.cv_N;
+    const sign = (tout - mem.cv_tn >= 0.0) ? 1.0 : -1.0;
+    const tdist = Math.abs(tout - mem.cv_tn);
+    const tround = Math.max(Math.abs(mem.cv_tn), Math.abs(tout)) * 2.2204460492503131e-16;
+    if (tdist < 2.0 * tround)
+        return CV_TOO_CLOSE;
+    // Set lower and upper bounds on h0, initial guess
+    const hlb = HLB_FACTOR * tround;
+    const hub = cvUpperBoundH0(mem, tdist);
+    // Initial guess - geometric mean
+    let hg = Math.sqrt(hlb * hub);
+    if (hub < 100.0 * hlb)
+        hg = 0.5 * hub;
+    // Iterate to estimate ydd norm
+    let hnew = hg;
+    let retval;
+    for (let count = 0; count < MAX_ITERS; count++) {
+        hg = hnew;
+        // Compute y''norm via finite difference
+        const result = cvYddNorm(mem, hg * sign);
+        retval = result.retval;
+        if (retval !== CV_SUCCESS)
+            return retval;
+        const yddnrm = result.yddnrm;
+        // If yddnrm is zero, use current hg
+        if (yddnrm * hub * hub > 2.0)
+            hnew = Math.sqrt(2.0 / yddnrm);
+        else
+            hnew = Math.sqrt(hg * hub);
+        // Refine: make hnew be a geometric mean biased toward smaller
+        hnew = H_BIAS * hnew + (1.0 - H_BIAS) * hg;
+        // Check convergence
+        if (hnew / hg > 2.0 || hnew / hg < 0.5)
+            continue;
+        // If converged, break
+        break;
+    }
+    // Bound hnew
+    hnew = Math.min(hnew, hub);
+    hnew = Math.max(hnew, hlb);
+    // Apply bounds on step size
+    if (mem.cv_hmax_inv > 0.0)
+        hnew = Math.min(hnew, 1.0 / mem.cv_hmax_inv);
+    if (mem.cv_hmin > 0.0)
+        hnew = Math.max(hnew, mem.cv_hmin);
+    // Bound by tdist
+    hnew = Math.min(hnew, tdist);
+    mem.cv_h = sign * hnew;
+    mem.cv_next_h = mem.cv_h;
+    mem.cv_hscale = mem.cv_h;
+    // Scale zn[1] = h * y'(t0)
+    vScale(mem.cv_h, mem.cv_zn[1], mem.cv_zn[1], N);
+    return CV_SUCCESS;
+}
+/**
+ * cvUpperBoundH0 - Compute upper bound on initial step size
+ * Based on ||y'|| weighted norm
+ */
+function cvUpperBoundH0(mem, tdist) {
+    const N = mem.cv_N;
+    // Get WRMS norm of y' (zn[1] = f(t0, y0) before scaling)
+    const hub_inv = wrmsNorm(N, mem.cv_zn[1], mem.cv_ewt);
+    // Bound based on tdist
+    let hub = HUB_FACTOR * tdist;
+    // Use smaller of two bounds
+    if (hub * hub_inv > 1.0)
+        hub = 1.0 / hub_inv;
+    return hub;
+}
+/**
+ * cvYddNorm - Estimate ||y''|| via finite difference
+ * y(t0 + hg) ~ y0 + hg*y'0
+ * y'(t0 + hg) = f(t0 + hg, y(t0 + hg))
+ * y'' ~ (y'(t0+hg) - y'(t0)) / hg
+ */
+function cvYddNorm(mem, hg) {
+    const N = mem.cv_N;
+    // y = zn[0] + hg * zn[1] (zn[1] still holds unscaled f at this point)
+    vLinearSum(hg, mem.cv_zn[1], 1.0, mem.cv_zn[0], mem.cv_y, N);
+    // Evaluate f(t0 + hg, y)
+    const retval = mem.cv_f(mem.cv_tn + hg, mem.cv_y, mem.cv_tempv, mem.cv_user_data);
+    mem.cv_nfe++;
+    if (retval < 0)
+        return { retval: CV_RHSFUNC_FAIL, yddnrm: 0 };
+    if (retval > 0)
+        return { retval: RHSFUNC_RECVR, yddnrm: 0 };
+    // tempv = (f(t0+hg, y) - zn[1]) / hg  (estimate of y'')
+    const hg_inv = 1.0 / hg;
+    vLinearSum(hg_inv, mem.cv_tempv, -hg_inv, mem.cv_zn[1], mem.cv_tempv, N);
+    const yddnrm = wrmsNorm(N, mem.cv_tempv, mem.cv_ewt);
+    return { retval: CV_SUCCESS, yddnrm };
+}
+//# sourceMappingURL=cvode_hin.js.map

package/dist/src/solver-tools/cvode/cvode_hin.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cvode_hin.js","sourceRoot":"","sources":["../../../../src/solver-tools/cvode/cvode_hin.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,uDAAuD;AACvD,wEAAwE;AAExE,OAAO,EAAW,QAAQ,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,EACxD,eAAe,EAAE,aAAa,EAAE,YAAY,EAC5C,UAAU,EAAE,UAAU,EAAE,MAAM,EAAE,SAAS,EAAO,MAAM,UAAU,CAAC;AAEnE;;;;;;;;GAQG;AACH,MAAM,UAAU,KAAK,CAAC,GAAa,EAAE,IAAY;IAC/C,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IACnB,MAAM,IAAI,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,KAAK,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IACpD,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC;IACzC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,sBAAsB,CAAC;IAEtF,IAAI,KAAK,GAAG,GAAG,GAAG,MAAM;QAAE,OAAO,YAAY,CAAC;IAE9C,kDAAkD;IAClD,MAAM,GAAG,GAAG,UAAU,GAAG,MAAM,CAAC;IAChC,MAAM,GAAG,GAAG,cAAc,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;IAEvC,iCAAiC;IACjC,IAAI,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;IAC9B,IAAI,GAAG,GAAG,KAAK,GAAG,GAAG;QACnB,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC;IAGjB,+BAA+B;IAC/B,IAAI,IAAI,GAAG,EAAE,CAAC;IACd,IAAI,MAAc,CAAC;IAEnB,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,SAAS,EAAE,KAAK,EAAE,EAAE,CAAC;QAC/C,EAAE,GAAG,IAAI,CAAC;QAEV,wCAAwC;QACxC,MAAM,MAAM,GAAG,SAAS,CAAC,GAAG,EAAE,EAAE,GAAG,IAAI,CAAC,CAAC;QACzC,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;QACvB,IAAI,MAAM,KAAK,UAAU;YAAE,OAAO,MAAM,CAAC;QAEzC,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;QAE7B,oCAAoC;QACpC,IAAI,MAAM,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG;YAC1B,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC;;YAE/B,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;QAG7B,8DAA8D;QAC9D,IAAI,GAAG,MAAM,GAAG,IAAI,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC;QAE3C,oBAAoB;QACpB,IAAI,IAAI,GAAG,EAAE,GAAG,GAAG,IAAI,IAAI,GAAG,EAAE,GAAG,GAAG;YAAE,SAAS;QACjD,sBAAsB;QACtB,MAAM;IACR,CAAC;IAED,aAAa;IACb,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;IAC3B,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;IAE3B,4BAA4B;IAC5B,IAAI,GAAG,CAAC,WAAW,GAAG,GAAG;QACvB,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,GAAG,GAAG,CAAC,WAAW,CAAC,CAAC;IAE/C,IAAI,GAAG,CAAC,OAAO,GAAG,GAAG;QACnB,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,OAAO,CAAC,CAAC;IAGrC,iBAAiB;IACjB,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;IAE7B,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC;IACvB,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,IAAI,CAAC;IACzB,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,IAAI,CAAC;IAEzB,2BAA2B;IAC3B,MAAM,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAEhD,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;GAGG;AACH,SAAS,cAAc,CAAC,GAAa,EAAE,KAAa;IAClD,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IAEnB,yDAAyD;IACzD,MAAM,OAAO,GAAG,QAAQ,CAAC,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC;IAEtD,uBAAuB;IACvB,IAAI,GAAG,GAAG,UAAU,GAAG,KAAK,CAAC;IAE7B,4BAA4B;IAC5B,IAAI,GAAG,GAAG,OAAO,GAAG,GAAG;QACrB,GAAG,GAAG,GAAG,GAAG,OAAO,CAAC;IAGtB,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;;GAKG;AACH,SAAS,SAAS,CAAC,GAAa,EAAE,EAAU;IAC1C,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IAEnB,sEAAsE;IACtE,UAAU,CAAC,EAAE,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAE7D,yBAAyB;IACzB,MAAM,MAAM,GAAG,GAAG,CAAC,IAAK,CAAC,GAAG,CAAC,KAAK,GAAG,EAAE,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IACnF,GAAG,CAAC,MAAM,EAAE,CAAC;IACb,IAAI,MAAM,GAAG,CAAC;QAAE,OAAO,EAAC,MAAM,EAAE,eAAe,EAAE,MAAM,EAAE,CAAC,EAAC,CAAC;IAC5D,IAAI,MAAM,GAAG,CAAC;QAAE,OAAO,EAAC,MAAM,EAAE,aAAa,EAAE,MAAM,EAAE,CAAC,EAAC,CAAC;IAE1D,wDAAwD;IACxD,MAAM,MAAM,GAAG,GAAG,GAAG,EAAE,CAAC;IACxB,UAAU,CAAC,MAAM,EAAE,GAAG,CAAC,QAAQ,EAAE,CAAC,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC;IAEzE,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,EAAE,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC;IAErD,OAAO,EAAC,MAAM,EAAE,UAAU,EAAE,MAAM,EAAC,CAAC;AACtC,CAAC"}

package/dist/src/solver-tools/cvode/cvode_io.d.ts

@@ -0,0 +1,98 @@
+/**
+ * TypeScript port of SUNDIALS CVODE v7.5.0 (https://github.com/LLNL/sundials)
+ *
+ * Copyright (c) 2025, Lawrence Livermore National Security, University of
+ *   Maryland Baltimore County, and the SUNDIALS contributors.
+ * Copyright (c) 2013-2025, Lawrence Livermore National Security and
+ *   Southern Methodist University.
+ * Copyright (c) 2002-2013, Lawrence Livermore National Security.
+ * Copyright (c) 2026 Datagrok.
+ *
+ * Licensed under the BSD-3-Clause License. See THIRD_PARTY_LICENSES
+ * for the full license text and provenance chain of the original code.
+ */
+import { CvodeMem } from './common';
+export interface CvodeStats {
+    nSteps: number;
+    nRhsEvals: number;
+    nLinSolvSetups: number;
+    nErrTestFails: number;
+    nNonlinSolvIters: number;
+    nNonlinSolvConvFails: number;
+    nJacEvals: number;
+    nGEvals: number;
+    lastOrder: number;
+    lastStep: number;
+    currentTime: number;
+}
+/**
+ * Set scalar relative and absolute tolerances.
+ * Requires rtol >= 0 and atol >= 0.
+ */
+export declare function cvodeSStolerances(mem: CvodeMem, rtol: number, atol: number): number;
+/**
+ * Set scalar relative tolerance and vector absolute tolerance.
+ * Requires rtol >= 0 and all atol[i] >= 0.
+ */
+export declare function cvodeSVtolerances(mem: CvodeMem, rtol: number, atol: Float64Array): number;
+/**
+ * Set the maximum number of internal steps per call to CVode.
+ */
+export declare function cvodeSetMaxNumSteps(mem: CvodeMem, mxstep: number): number;
+/**
+ * Set the maximum method order. Clamped to [1, qmax] where qmax is
+ * ADAMS_Q_MAX or BDF_Q_MAX depending on the linear multistep method.
+ */
+export declare function cvodeSetMaxOrd(mem: CvodeMem, maxord: number): number;
+/**
+ * Set the maximum step size. If hmax <= 0, no upper limit is imposed
+ * (cv_hmax_inv = 0). Otherwise cv_hmax_inv = 1 / hmax.
+ */
+export declare function cvodeSetMaxStep(mem: CvodeMem, hmax: number): number;
+/**
+ * Set the minimum step size (stored as absolute value).
+ */
+export declare function cvodeSetMinStep(mem: CvodeMem, hmin: number): number;
+/**
+ * Set the initial step size hint.
+ */
+export declare function cvodeSetInitStep(mem: CvodeMem, hin: number): number;
+/**
+ * Set a stop time. Integration will not proceed past tstop.
+ */
+export declare function cvodeSetStopTime(mem: CvodeMem, tstop: number): number;
+/**
+ * Attach user data that will be passed to the RHS, Jacobian, and root
+ * functions.
+ */
+export declare function cvodeSetUserData(mem: CvodeMem, data: any): number;
+/**
+ * Enable or disable BDF stability limit detection.
+ */
+export declare function cvodeSetStabLimDet(mem: CvodeMem, on: boolean): number;
+/** Return the total number of internal steps taken. */
+export declare function cvodeGetNumSteps(mem: CvodeMem): number;
+/** Return the total number of RHS evaluations. */
+export declare function cvodeGetNumRhsEvals(mem: CvodeMem): number;
+/** Return the number of linear solver setup calls. */
+export declare function cvodeGetNumLinSolvSetups(mem: CvodeMem): number;
+/** Return the number of error test failures. */
+export declare function cvodeGetNumErrTestFails(mem: CvodeMem): number;
+/** Return the number of nonlinear solver iterations. */
+export declare function cvodeGetNumNonlinSolvIters(mem: CvodeMem): number;
+/** Return the number of nonlinear solver convergence failures. */
+export declare function cvodeGetNumNonlinSolvConvFails(mem: CvodeMem): number;
+/** Return the number of Jacobian evaluations (from the linear solver). */
+export declare function cvodeGetNumJacEvals(mem: CvodeMem): number;
+/** Return the total number of root-function evaluations. */
+export declare function cvodeGetNumGEvals(mem: CvodeMem): number;
+/** Return the method order used on the last successful step. */
+export declare function cvodeGetLastOrder(mem: CvodeMem): number;
+/** Return the step size used on the last successful step. */
+export declare function cvodeGetLastStep(mem: CvodeMem): number;
+/** Return the current internal time reached by the solver. */
+export declare function cvodeGetCurrentTime(mem: CvodeMem): number;
+/**
+ * Return a snapshot of all integrator statistics as a plain object.
+ */
+export declare function cvodeGetIntegratorStats(mem: CvodeMem): CvodeStats;

package/dist/src/solver-tools/cvode/cvode_io.js

@@ -0,0 +1,202 @@
+/**
+ * TypeScript port of SUNDIALS CVODE v7.5.0 (https://github.com/LLNL/sundials)
+ *
+ * Copyright (c) 2025, Lawrence Livermore National Security, University of
+ *   Maryland Baltimore County, and the SUNDIALS contributors.
+ * Copyright (c) 2013-2025, Lawrence Livermore National Security and
+ *   Southern Methodist University.
+ * Copyright (c) 2002-2013, Lawrence Livermore National Security.
+ * Copyright (c) 2026 Datagrok.
+ *
+ * Licensed under the BSD-3-Clause License. See THIRD_PARTY_LICENSES
+ * for the full license text and provenance chain of the original code.
+ */
+// CVODE TypeScript port - Configuration and statistics I/O functions
+// Ported from SUNDIALS CVODE v7.5.0 (cvode_io.c)
+import { CV_SUCCESS, CV_SS, CV_SV, ADAMS_Q_MAX, BDF_Q_MAX, CV_ADAMS, } from './common';
+// ============================================================================
+// Tolerance functions
+// ============================================================================
+/**
+ * Set scalar relative and absolute tolerances.
+ * Requires rtol >= 0 and atol >= 0.
+ */
+export function cvodeSStolerances(mem, rtol, atol) {
+    if (rtol < 0) {
+        mem.cv_error = 'cvodeSStolerances: rtol must be >= 0';
+        return -1;
+    }
+    if (atol < 0) {
+        mem.cv_error = 'cvodeSStolerances: atol must be >= 0';
+        return -1;
+    }
+    mem.cv_reltol = rtol;
+    mem.cv_Sabstol = atol;
+    mem.cv_itol = CV_SS;
+    mem.cv_atolmin0 = (atol === 0);
+    return CV_SUCCESS;
+}
+/**
+ * Set scalar relative tolerance and vector absolute tolerance.
+ * Requires rtol >= 0 and all atol[i] >= 0.
+ */
+export function cvodeSVtolerances(mem, rtol, atol) {
+    if (rtol < 0) {
+        mem.cv_error = 'cvodeSVtolerances: rtol must be >= 0';
+        return -1;
+    }
+    for (let i = 0; i < atol.length; i++) {
+        if (atol[i] < 0) {
+            mem.cv_error = `cvodeSVtolerances: atol[${i}] must be >= 0`;
+            return -1;
+        }
+    }
+    mem.cv_reltol = rtol;
+    mem.cv_Vabstol = new Float64Array(atol);
+    mem.cv_itol = CV_SV;
+    let minAtol = Infinity;
+    for (let i = 0; i < atol.length; i++)
+        if (atol[i] < minAtol)
+            minAtol = atol[i];
+    mem.cv_atolmin0 = (minAtol === 0);
+    mem.cv_VabstolMallocDone = true;
+    return CV_SUCCESS;
+}
+// ============================================================================
+// Optional input functions
+// ============================================================================
+/**
+ * Set the maximum number of internal steps per call to CVode.
+ */
+export function cvodeSetMaxNumSteps(mem, mxstep) {
+    mem.cv_mxstep = mxstep;
+    return CV_SUCCESS;
+}
+/**
+ * Set the maximum method order. Clamped to [1, qmax] where qmax is
+ * ADAMS_Q_MAX or BDF_Q_MAX depending on the linear multistep method.
+ */
+export function cvodeSetMaxOrd(mem, maxord) {
+    const qmax = (mem.cv_lmm === CV_ADAMS) ? ADAMS_Q_MAX : BDF_Q_MAX;
+    if (maxord < 1)
+        maxord = 1;
+    if (maxord > qmax)
+        maxord = qmax;
+    mem.cv_qmax = maxord;
+    return CV_SUCCESS;
+}
+/**
+ * Set the maximum step size. If hmax <= 0, no upper limit is imposed
+ * (cv_hmax_inv = 0). Otherwise cv_hmax_inv = 1 / hmax.
+ */
+export function cvodeSetMaxStep(mem, hmax) {
+    if (hmax <= 0)
+        mem.cv_hmax_inv = 0;
+    else
+        mem.cv_hmax_inv = 1.0 / hmax;
+    return CV_SUCCESS;
+}
+/**
+ * Set the minimum step size (stored as absolute value).
+ */
+export function cvodeSetMinStep(mem, hmin) {
+    mem.cv_hmin = Math.abs(hmin);
+    return CV_SUCCESS;
+}
+/**
+ * Set the initial step size hint.
+ */
+export function cvodeSetInitStep(mem, hin) {
+    mem.cv_hin = hin;
+    return CV_SUCCESS;
+}
+/**
+ * Set a stop time. Integration will not proceed past tstop.
+ */
+export function cvodeSetStopTime(mem, tstop) {
+    mem.cv_tstop = tstop;
+    mem.cv_tstopset = true;
+    mem.cv_tstopinterp = true;
+    return CV_SUCCESS;
+}
+/**
+ * Attach user data that will be passed to the RHS, Jacobian, and root
+ * functions.
+ */
+export function cvodeSetUserData(mem, data) {
+    mem.cv_user_data = data;
+    return CV_SUCCESS;
+}
+/**
+ * Enable or disable BDF stability limit detection.
+ */
+export function cvodeSetStabLimDet(mem, on) {
+    mem.cv_sldeton = on;
+    return CV_SUCCESS;
+}
+// ============================================================================
+// Statistics getters
+// ============================================================================
+/** Return the total number of internal steps taken. */
+export function cvodeGetNumSteps(mem) {
+    return mem.cv_nst;
+}
+/** Return the total number of RHS evaluations. */
+export function cvodeGetNumRhsEvals(mem) {
+    return mem.cv_nfe;
+}
+/** Return the number of linear solver setup calls. */
+export function cvodeGetNumLinSolvSetups(mem) {
+    return mem.cv_nsetups;
+}
+/** Return the number of error test failures. */
+export function cvodeGetNumErrTestFails(mem) {
+    return mem.cv_netf;
+}
+/** Return the number of nonlinear solver iterations. */
+export function cvodeGetNumNonlinSolvIters(mem) {
+    return mem.cv_nni;
+}
+/** Return the number of nonlinear solver convergence failures. */
+export function cvodeGetNumNonlinSolvConvFails(mem) {
+    return mem.cv_ncfn;
+}
+/** Return the number of Jacobian evaluations (from the linear solver). */
+export function cvodeGetNumJacEvals(mem) {
+    return mem.cv_lmem?.nje ?? 0;
+}
+/** Return the total number of root-function evaluations. */
+export function cvodeGetNumGEvals(mem) {
+    return mem.cv_nge;
+}
+/** Return the method order used on the last successful step. */
+export function cvodeGetLastOrder(mem) {
+    return mem.cv_qu;
+}
+/** Return the step size used on the last successful step. */
+export function cvodeGetLastStep(mem) {
+    return mem.cv_hu;
+}
+/** Return the current internal time reached by the solver. */
+export function cvodeGetCurrentTime(mem) {
+    return mem.cv_tn;
+}
+/**
+ * Return a snapshot of all integrator statistics as a plain object.
+ */
+export function cvodeGetIntegratorStats(mem) {
+    return {
+        nSteps: cvodeGetNumSteps(mem),
+        nRhsEvals: cvodeGetNumRhsEvals(mem),
+        nLinSolvSetups: cvodeGetNumLinSolvSetups(mem),
+        nErrTestFails: cvodeGetNumErrTestFails(mem),
+        nNonlinSolvIters: cvodeGetNumNonlinSolvIters(mem),
+        nNonlinSolvConvFails: cvodeGetNumNonlinSolvConvFails(mem),
+        nJacEvals: cvodeGetNumJacEvals(mem),
+        nGEvals: cvodeGetNumGEvals(mem),
+        lastOrder: cvodeGetLastOrder(mem),
+        lastStep: cvodeGetLastStep(mem),
+        currentTime: cvodeGetCurrentTime(mem),
+    };
+}
+//# sourceMappingURL=cvode_io.js.map

package/dist/src/solver-tools/cvode/cvode_io.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"cvode_io.js","sourceRoot":"","sources":["../../../../src/solver-tools/cvode/cvode_io.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,qEAAqE;AACrE,iDAAiD;AAEjD,OAAO,EAEL,UAAU,EACV,KAAK,EACL,KAAK,EACL,WAAW,EACX,SAAS,EACT,QAAQ,GACT,MAAM,UAAU,CAAC;AAoBlB,+EAA+E;AAC/E,sBAAsB;AACtB,+EAA+E;AAE/E;;;GAGG;AACH,MAAM,UAAU,iBAAiB,CAAC,GAAa,EAAE,IAAY,EAAE,IAAY;IACzE,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;QACb,GAAG,CAAC,QAAQ,GAAG,sCAAsC,CAAC;QACtD,OAAO,CAAC,CAAC,CAAC;IACZ,CAAC;IACD,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;QACb,GAAG,CAAC,QAAQ,GAAG,sCAAsC,CAAC;QACtD,OAAO,CAAC,CAAC,CAAC;IACZ,CAAC;IAED,GAAG,CAAC,SAAS,GAAG,IAAI,CAAC;IACrB,GAAG,CAAC,UAAU,GAAG,IAAI,CAAC;IACtB,GAAG,CAAC,OAAO,GAAG,KAAK,CAAC;IACpB,GAAG,CAAC,WAAW,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC;IAE/B,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,iBAAiB,CAAC,GAAa,EAAE,IAAY,EAAE,IAAkB;IAC/E,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;QACb,GAAG,CAAC,QAAQ,GAAG,sCAAsC,CAAC;QACtD,OAAO,CAAC,CAAC,CAAC;IACZ,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC;YAChB,GAAG,CAAC,QAAQ,GAAG,2BAA2B,CAAC,gBAAgB,CAAC;YAC5D,OAAO,CAAC,CAAC,CAAC;QACZ,CAAC;IACH,CAAC;IAED,GAAG,CAAC,SAAS,GAAG,IAAI,CAAC;IACrB,GAAG,CAAC,UAAU,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC;IACxC,GAAG,CAAC,OAAO,GAAG,KAAK,CAAC;IAEpB,IAAI,OAAO,GAAG,QAAQ,CAAC;IACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE;QAClC,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO;YAAE,OAAO,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IAE3C,GAAG,CAAC,WAAW,GAAG,CAAC,OAAO,KAAK,CAAC,CAAC,CAAC;IAClC,GAAG,CAAC,oBAAoB,GAAG,IAAI,CAAC;IAEhC,OAAO,UAAU,CAAC;AACpB,CAAC;AAED,+EAA+E;AAC/E,2BAA2B;AAC3B,+EAA+E;AAE/E;;GAEG;AACH,MAAM,UAAU,mBAAmB,CAAC,GAAa,EAAE,MAAc;IAC/D,GAAG,CAAC,SAAS,GAAG,MAAM,CAAC;IACvB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,cAAc,CAAC,GAAa,EAAE,MAAc;IAC1D,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,MAAM,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,SAAS,CAAC;IACjE,IAAI,MAAM,GAAG,CAAC;QAAE,MAAM,GAAG,CAAC,CAAC;IAC3B,IAAI,MAAM,GAAG,IAAI;QAAE,MAAM,GAAG,IAAI,CAAC;IACjC,GAAG,CAAC,OAAO,GAAG,MAAM,CAAC;IACrB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,eAAe,CAAC,GAAa,EAAE,IAAY;IACzD,IAAI,IAAI,IAAI,CAAC;QACX,GAAG,CAAC,WAAW,GAAG,CAAC,CAAC;;QAEpB,GAAG,CAAC,WAAW,GAAG,GAAG,GAAG,IAAI,CAAC;IAE/B,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,eAAe,CAAC,GAAa,EAAE,IAAY;IACzD,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAC7B,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,gBAAgB,CAAC,GAAa,EAAE,GAAW;IACzD,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC;IACjB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,gBAAgB,CAAC,GAAa,EAAE,KAAa;IAC3D,GAAG,CAAC,QAAQ,GAAG,KAAK,CAAC;IACrB,GAAG,CAAC,WAAW,GAAG,IAAI,CAAC;IACvB,GAAG,CAAC,cAAc,GAAG,IAAI,CAAC;IAC1B,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,gBAAgB,CAAC,GAAa,EAAE,IAAS;IACvD,GAAG,CAAC,YAAY,GAAG,IAAI,CAAC;IACxB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,kBAAkB,CAAC,GAAa,EAAE,EAAW;IAC3D,GAAG,CAAC,UAAU,GAAG,EAAE,CAAC;IACpB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED,+EAA+E;AAC/E,qBAAqB;AACrB,+EAA+E;AAE/E,uDAAuD;AACvD,MAAM,UAAU,gBAAgB,CAAC,GAAa;IAC5C,OAAO,GAAG,CAAC,MAAM,CAAC;AACpB,CAAC;AAED,kDAAkD;AAClD,MAAM,UAAU,mBAAmB,CAAC,GAAa;IAC/C,OAAO,GAAG,CAAC,MAAM,CAAC;AACpB,CAAC;AAED,sDAAsD;AACtD,MAAM,UAAU,wBAAwB,CAAC,GAAa;IACpD,OAAO,GAAG,CAAC,UAAU,CAAC;AACxB,CAAC;AAED,gDAAgD;AAChD,MAAM,UAAU,uBAAuB,CAAC,GAAa;IACnD,OAAO,GAAG,CAAC,OAAO,CAAC;AACrB,CAAC;AAED,wDAAwD;AACxD,MAAM,UAAU,0BAA0B,CAAC,GAAa;IACtD,OAAO,GAAG,CAAC,MAAM,CAAC;AACpB,CAAC;AAED,kEAAkE;AAClE,MAAM,UAAU,8BAA8B,CAAC,GAAa;IAC1D,OAAO,GAAG,CAAC,OAAO,CAAC;AACrB,CAAC;AAED,0EAA0E;AAC1E,MAAM,UAAU,mBAAmB,CAAC,GAAa;IAC/C,OAAO,GAAG,CAAC,OAAO,EAAE,GAAG,IAAI,CAAC,CAAC;AAC/B,CAAC;AAED,4DAA4D;AAC5D,MAAM,UAAU,iBAAiB,CAAC,GAAa;IAC7C,OAAO,GAAG,CAAC,MAAM,CAAC;AACpB,CAAC;AAED,gEAAgE;AAChE,MAAM,UAAU,iBAAiB,CAAC,GAAa;IAC7C,OAAO,GAAG,CAAC,KAAK,CAAC;AACnB,CAAC;AAED,6DAA6D;AAC7D,MAAM,UAAU,gBAAgB,CAAC,GAAa;IAC5C,OAAO,GAAG,CAAC,KAAK,CAAC;AACnB,CAAC;AAED,8DAA8D;AAC9D,MAAM,UAAU,mBAAmB,CAAC,GAAa;IAC/C,OAAO,GAAG,CAAC,KAAK,CAAC;AACnB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,uBAAuB,CAAC,GAAa;IACnD,OAAO;QACL,MAAM,EAAE,gBAAgB,CAAC,GAAG,CAAC;QAC7B,SAAS,EAAE,mBAAmB,CAAC,GAAG,CAAC;QACnC,cAAc,EAAE,wBAAwB,CAAC,GAAG,CAAC;QAC7C,aAAa,EAAE,uBAAuB,CAAC,GAAG,CAAC;QAC3C,gBAAgB,EAAE,0BAA0B,CAAC,GAAG,CAAC;QACjD,oBAAoB,EAAE,8BAA8B,CAAC,GAAG,CAAC;QACzD,SAAS,EAAE,mBAAmB,CAAC,GAAG,CAAC;QACnC,OAAO,EAAE,iBAAiB,CAAC,GAAG,CAAC;QAC/B,SAAS,EAAE,iBAAiB,CAAC,GAAG,CAAC;QACjC,QAAQ,EAAE,gBAAgB,CAAC,GAAG,CAAC;QAC/B,WAAW,EAAE,mBAAmB,CAAC,GAAG,CAAC;KACtC,CAAC;AACJ,CAAC"}

package/dist/src/solver-tools/cvode/cvode_ls.d.ts

@@ -0,0 +1,23 @@
+/**
+ * TypeScript port of SUNDIALS CVODE v7.5.0 (https://github.com/LLNL/sundials)
+ *
+ * Copyright (c) 2025, Lawrence Livermore National Security, University of
+ *   Maryland Baltimore County, and the SUNDIALS contributors.
+ * Copyright (c) 2013-2025, Lawrence Livermore National Security and
+ *   Southern Methodist University.
+ * Copyright (c) 2002-2013, Lawrence Livermore National Security.
+ * Copyright (c) 2026 Datagrok.
+ *
+ * Licensed under the BSD-3-Clause License. See THIRD_PARTY_LICENSES
+ * for the full license text and provenance chain of the original code.
+ */
+import { CvodeMem, CvodeJacFn } from './common';
+/**
+ * cvodeSetLinearSolver - Initialize the dense linear solver.
+ * Sets up CvLsMem, allocates matrix, sets function pointers on CvodeMem.
+ */
+export declare function cvodeSetLinearSolver(mem: CvodeMem): number;
+/**
+ * cvodeSetJacFn - Set a user-supplied Jacobian function.
+ */
+export declare function cvodeSetJacFn(mem: CvodeMem, jac: CvodeJacFn | null): number;