diff-grok 1.2.0 → 1.3.1

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

@@ -1,150 +0,0 @@
-/**
- * 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 {CvodeMem, wrmsNorm, vLinearSum, vScale, CV_SUCCESS,
-  CV_RHSFUNC_FAIL, RHSFUNC_RECVR, CV_TOO_CLOSE,
-  HLB_FACTOR, HUB_FACTOR, H_BIAS, MAX_ITERS, TINY} 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: CvodeMem, tout: number): number {
-  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: number;
-
-  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: CvodeMem, tdist: number): number {
-  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: CvodeMem, hg: number): { retval: number; yddnrm: number } {
-  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};
-}

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

@@ -1 +0,0 @@
-{"version":3,"file":"cvode_io.js","sourceRoot":"","sources":["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,MAAA,MAAA,GAAG,CAAC,OAAO,0CAAE,GAAG,mCAAI,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/src/solver-tools/cvode/cvode_io.ts

@@ -1,258 +0,0 @@
-/**
- * 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 {
-  CvodeMem,
-  CV_SUCCESS,
-  CV_SS,
-  CV_SV,
-  ADAMS_Q_MAX,
-  BDF_Q_MAX,
-  CV_ADAMS,
-} from './common';
-
-// ============================================================================
-// Statistics interface
-// ============================================================================
-
-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;
-}
-
-// ============================================================================
-// Tolerance functions
-// ============================================================================
-
-/**
- * Set scalar relative and absolute tolerances.
- * Requires rtol >= 0 and atol >= 0.
- */
-export function cvodeSStolerances(mem: CvodeMem, rtol: number, atol: number): number {
-  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: CvodeMem, rtol: number, atol: Float64Array): number {
-  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: CvodeMem, mxstep: number): number {
-  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: CvodeMem, maxord: number): number {
-  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: CvodeMem, hmax: number): number {
-  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: CvodeMem, hmin: number): number {
-  mem.cv_hmin = Math.abs(hmin);
-  return CV_SUCCESS;
-}
-
-/**
- * Set the initial step size hint.
- */
-export function cvodeSetInitStep(mem: CvodeMem, hin: number): number {
-  mem.cv_hin = hin;
-  return CV_SUCCESS;
-}
-
-/**
- * Set a stop time. Integration will not proceed past tstop.
- */
-export function cvodeSetStopTime(mem: CvodeMem, tstop: number): number {
-  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: CvodeMem, data: any): number {
-  mem.cv_user_data = data;
-  return CV_SUCCESS;
-}
-
-/**
- * Enable or disable BDF stability limit detection.
- */
-export function cvodeSetStabLimDet(mem: CvodeMem, on: boolean): number {
-  mem.cv_sldeton = on;
-  return CV_SUCCESS;
-}
-
-// ============================================================================
-// Statistics getters
-// ============================================================================
-
-/** Return the total number of internal steps taken. */
-export function cvodeGetNumSteps(mem: CvodeMem): number {
-  return mem.cv_nst;
-}
-
-/** Return the total number of RHS evaluations. */
-export function cvodeGetNumRhsEvals(mem: CvodeMem): number {
-  return mem.cv_nfe;
-}
-
-/** Return the number of linear solver setup calls. */
-export function cvodeGetNumLinSolvSetups(mem: CvodeMem): number {
-  return mem.cv_nsetups;
-}
-
-/** Return the number of error test failures. */
-export function cvodeGetNumErrTestFails(mem: CvodeMem): number {
-  return mem.cv_netf;
-}
-
-/** Return the number of nonlinear solver iterations. */
-export function cvodeGetNumNonlinSolvIters(mem: CvodeMem): number {
-  return mem.cv_nni;
-}
-
-/** Return the number of nonlinear solver convergence failures. */
-export function cvodeGetNumNonlinSolvConvFails(mem: CvodeMem): number {
-  return mem.cv_ncfn;
-}
-
-/** Return the number of Jacobian evaluations (from the linear solver). */
-export function cvodeGetNumJacEvals(mem: CvodeMem): number {
-  return mem.cv_lmem?.nje ?? 0;
-}
-
-/** Return the total number of root-function evaluations. */
-export function cvodeGetNumGEvals(mem: CvodeMem): number {
-  return mem.cv_nge;
-}
-
-/** Return the method order used on the last successful step. */
-export function cvodeGetLastOrder(mem: CvodeMem): number {
-  return mem.cv_qu;
-}
-
-/** Return the step size used on the last successful step. */
-export function cvodeGetLastStep(mem: CvodeMem): number {
-  return mem.cv_hu;
-}
-
-/** Return the current internal time reached by the solver. */
-export function cvodeGetCurrentTime(mem: CvodeMem): number {
-  return mem.cv_tn;
-}
-
-/**
- * Return a snapshot of all integrator statistics as a plain object.
- */
-export function cvodeGetIntegratorStats(mem: CvodeMem): CvodeStats {
-  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),
-  };
-}

package/src/solver-tools/cvode/cvode_ls.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"cvode_ls.js","sourceRoot":"","sources":["cvode_ls.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,wDAAwD;AACxD,sCAAsC;AACtC,EAAE;AACF,uEAAuE;AACvE,yEAAyE;AAEzE,OAAO,EAAW,OAAO,EAAc,QAAQ,EAC7C,UAAU,EACM,aAAa,EAC7B,SAAS,EAAE,UAAU,EAAS,MAAM,UAAU,CAAC;AACjD,OAAO,EAAC,KAAK,EAAE,KAAK,EAAC,MAAM,gBAAgB,CAAC;AAE5C;;;GAGG;AACH,MAAM,UAAU,oBAAoB,CAAC,GAAa;IAChD,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IACnB,MAAM,IAAI,GAAG,IAAI,OAAO,EAAE,CAAC;IAE3B,6CAA6C;IAC7C,IAAI,CAAC,CAAC,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;IACtB,IAAI,CAAC,MAAM,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;IAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAChC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IACvC,CAAC;IACD,IAAI,CAAC,MAAM,GAAG,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC;IAEhC,eAAe;IACf,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;IAClB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;IAClB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;IACjB,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC;IACb,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC;IACf,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC;IACf,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC;IACd,IAAI,CAAC,IAAI,GAAG,SAAS,CAAC;IACtB,IAAI,CAAC,UAAU,GAAG,UAAU,CAAC;IAE7B,qBAAqB;IACrB,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC;IACnB,GAAG,CAAC,QAAQ,GAAG,cAAc,CAAC;IAC9B,GAAG,CAAC,SAAS,GAAG,SAAS,CAAC;IAC1B,GAAG,CAAC,SAAS,GAAG,SAAS,CAAC;IAC1B,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC;IAExB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,aAAa,CAAC,GAAa,EAAE,GAAsB;IACjE,MAAM,IAAI,GAAG,GAAG,CAAC,OAAO,CAAC;IACzB,IAAI,CAAC,IAAI;QAAE,OAAO,CAAC,CAAC,CAAC;IAErB,IAAI,GAAG,EAAE,CAAC;QACR,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,KAAK,GAAG,GAAG,CAAC;IACnB,CAAC;SAAM,CAAC;QACN,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;QAClB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;IACpB,CAAC;IACD,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;GAEG;AACH,SAAS,cAAc,CAAC,GAAa;IACnC,MAAM,IAAI,GAAG,GAAG,CAAC,OAAQ,CAAC;IAC1B,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC;IACb,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC;IACf,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC;IACf,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;IACjB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;;;;;;GAQG;AACH,SAAS,SAAS,CAAC,GAAa,EAAE,QAAgB,EAAE,KAAmB,EACrE,KAAmB,EAAE,OAA2B,EAChD,IAAkB,EAAE,IAAkB,EACtC,IAAkB;IAClB,MAAM,IAAI,GAAG,GAAG,CAAC,OAAQ,CAAC;IAC1B,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IACnB,IAAI,MAAc,CAAC;IAEnB,uCAAuC;IACvC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC;IACvB,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;IAEhB,IAAI,CAAC,IAAI,EAAE,CAAC;QACV,oBAAoB;QACpB,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,UAAU;YACjD,GAAG,GAAG,KAAK,CAAC;QAEd,kCAAkC;QAClC,IAAI,GAAG,CAAC,MAAM,IAAI,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI;YACtC,GAAG,GAAG,KAAK,CAAC;QAEd,iBAAiB;QACjB,IAAI,QAAQ,KAAK,aAAa;YAC5B,GAAG,GAAG,KAAK,CAAC;IAChB,CAAC;IAED,IAAI,GAAG,EAAE,CAAC;QACR,kDAAkD;QAClD,OAAO,CAAC,KAAK,GAAG,KAAK,CAAC;QAEtB,wCAAwC;QACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBACxB,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAEnD,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;QACtB,CAAC;IACH,CAAC;SAAM,CAAC;QACN,qBAAqB;QACrB,IAAI,CAAC,GAAG,EAAE,CAAC;QACX,IAAI,CAAC,KAAK,GAAG,GAAG,CAAC,MAAM,CAAC;QACxB,IAAI,CAAC,IAAI,GAAG,KAAK,CAAC;QAClB,OAAO,CAAC,KAAK,GAAG,IAAI,CAAC;QAErB,oBAAoB;QACpB,IAAI,IAAI,CAAC,KAAK;YACZ,MAAM,GAAG,SAAS,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,IAAI,CAAC,MAAM,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;;YAE1E,MAAM,GAAG,IAAI,CAAC,KAAM,CAAC,GAAG,CAAC,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,IAAI,CAAC,MAAM,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;QAG/E,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;YACf,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;YACjB,OAAO,CAAC,CAAC,CAAC;QACZ,CAAC;QACD,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;YACf,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;YACjB,OAAO,CAAC,CAAC;QACX,CAAC;QAED,uBAAuB;QACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBACxB,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAEnD,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;QACtB,CAAC;IACH,CAAC;IAED,kCAAkC;IAClC,MAAM,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;IAC3C,IAAI,IAAI,KAAK,CAAC,EAAE,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,OAAO,CAAC,CAAC,CAAC,cAAc;IAC1B,CAAC;IAED,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;GAEG;AACH,SAAS,SAAS,CAAC,GAAa,EAAE,CAAe,EAAE,MAAoB,EACrE,IAAkB,EAAE,IAAkB;IACtC,MAAM,IAAI,GAAG,GAAG,CAAC,OAAQ,CAAC;IAC1B,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IAEnB,KAAK,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IAEjC,sEAAsE;IACtE,IAAI,GAAG,CAAC,SAAS,KAAK,GAAG,EAAE,CAAC;QAC1B,MAAM,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,SAAS,CAAC,CAAC;QAC3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YACxB,CAAC,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;IACnB,CAAC;IAED,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;;;;GAMG;AACH,SAAS,SAAS,CAAC,GAAa,EAAE,CAAS,EAAE,CAAe,EAAE,EAAgB,EAC5E,CAAiB,EAAE,IAAkB,EAAE,IAAkB;IACzD,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IACnB,MAAM,IAAI,GAAG,GAAG,CAAC,OAAQ,CAAC;IAC1B,MAAM,MAAM,GAAG,sBAAsB,CAAC;IACtC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IAC/B,MAAM,SAAS,GAAG,QAAQ,CAAC,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC;IAC9C,MAAM,MAAM,GAAG,CAAC,SAAS,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;IAElG,+CAA+C;IAC/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,MAAM,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACrB,IAAI,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACrE,gEAAgE;QAEhE,eAAe;QACf,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACZ,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,GAAG,GAAG,CAAC;QACxB,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,iCAAiC;QAE1D,6BAA6B;QAC7B,MAAM,MAAM,GAAG,GAAG,CAAC,IAAK,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;QAC1D,IAAI,CAAC,KAAK,EAAE,CAAC;QACb,GAAG,CAAC,MAAM,EAAE,CAAC;QACb,IAAI,MAAM,KAAK,CAAC;YAAE,OAAO,MAAM,CAAC;QAEhC,wDAAwD;QACxD,MAAM,OAAO,GAAG,GAAG,GAAG,GAAG,CAAC;QAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YACxB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;IAC1C,CAAC;IAED,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;GAEG;AACH,SAAS,QAAQ,CAAC,GAAa;IAC7B,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC;AACrB,CAAC"}

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

@@ -1,249 +0,0 @@
-/**
- * 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 - Dense linear solver interface
-// Translated from SUNDIALS cvode_ls.c
-//
-// Instead of SUNDIALS' SUNLinearSolver / SUNMatrix abstraction, we use
-// direct Float64Array[] matrices and dgefa / dgesl from dense_linalg.ts.
-
-import {CvodeMem, CvLsMem, CvodeJacFn, wrmsNorm,
-  CV_SUCCESS, CV_RHSFUNC_FAIL, RHSFUNC_RECVR,
-  CV_NO_FAILURES, CV_FAIL_BAD_J, CV_FAIL_OTHER,
-  CVLS_MSBJ, CVLS_DGMAX, ONEPSM} from './common';
-import {dgefa, dgesl} from './dense_linalg';
-
-/**
- * cvodeSetLinearSolver - Initialize the dense linear solver.
- * Sets up CvLsMem, allocates matrix, sets function pointers on CvodeMem.
- */
-export function cvodeSetLinearSolver(mem: CvodeMem): number {
-  const N = mem.cv_N;
-  const lmem = new CvLsMem();
-
-  // Allocate dense matrix A (N x N, row-major)
-  lmem.A = new Array(N);
-  lmem.savedJ = new Array(N);
-  for (let i = 0; i < N; i++) {
-    lmem.A[i] = new Float64Array(N);
-    lmem.savedJ[i] = new Float64Array(N);
-  }
-  lmem.pivots = new Int32Array(N);
-
-  // Set defaults
-  lmem.jacDQ = true;
-  lmem.jacFn = null;
-  lmem.jbad = true;
-  lmem.nje = 0;
-  lmem.nfeDQ = 0;
-  lmem.nstlj = 0;
-  lmem.tnlj = 0;
-  lmem.msbj = CVLS_MSBJ;
-  lmem.dgmax_jbad = CVLS_DGMAX;
-
-  // Attach to CvodeMem
-  mem.cv_lmem = lmem;
-  mem.cv_linit = cvLsInitialize;
-  mem.cv_lsetup = cvLsSetup;
-  mem.cv_lsolve = cvLsSolve;
-  mem.cv_lfree = cvLsFree;
-
-  return CV_SUCCESS;
-}
-
-/**
- * cvodeSetJacFn - Set a user-supplied Jacobian function.
- */
-export function cvodeSetJacFn(mem: CvodeMem, jac: CvodeJacFn | null): number {
-  const lmem = mem.cv_lmem;
-  if (!lmem) return -1;
-
-  if (jac) {
-    lmem.jacDQ = false;
-    lmem.jacFn = jac;
-  } else {
-    lmem.jacDQ = true;
-    lmem.jacFn = null;
-  }
-  return CV_SUCCESS;
-}
-
-/**
- * cvLsInitialize - Initialize linear solver counters.
- */
-function cvLsInitialize(mem: CvodeMem): number {
-  const lmem = mem.cv_lmem!;
-  lmem.nje = 0;
-  lmem.nfeDQ = 0;
-  lmem.nstlj = 0;
-  lmem.jbad = true;
-  return CV_SUCCESS;
-}
-
-/**
- * cvLsSetup - Compute Jacobian and form/factor the linear system matrix A = I - gamma*J.
- *
- * Heuristic for when to recompute Jacobian:
- * - On first call (nstlj == 0)
- * - If convfail is CV_FAIL_BAD_J (NLS flagged Jacobian as bad)
- * - If enough steps since last eval (nst >= nstlj + msbj)
- * - If gamma changed significantly (|gamrat - 1| > dgmax_jbad)
- */
-function cvLsSetup(mem: CvodeMem, convfail: number, ypred: Float64Array,
-  fpred: Float64Array, jcurPtr: { value: boolean },
-  tmp1: Float64Array, tmp2: Float64Array,
-  tmp3: Float64Array): number {
-  const lmem = mem.cv_lmem!;
-  const N = mem.cv_N;
-  let retval: number;
-
-  // Determine if Jacobian needs updating
-  const jbad = lmem.jbad;
-  let jok = !jbad;
-
-  if (!jbad) {
-    // Check gamma ratio
-    if (Math.abs(mem.cv_gamrat - 1.0) > lmem.dgmax_jbad)
-      jok = false;
-
-    // Check steps since last Jacobian
-    if (mem.cv_nst >= lmem.nstlj + lmem.msbj)
-      jok = false;
-
-    // Check convfail
-    if (convfail === CV_FAIL_BAD_J)
-      jok = false;
-  }
-
-  if (jok) {
-    // Reuse saved Jacobian, just form A = I - gamma*J
-    jcurPtr.value = false;
-
-    // Copy savedJ to A and form I - gamma*J
-    for (let i = 0; i < N; i++) {
-      for (let j = 0; j < N; j++)
-        lmem.A[i][j] = -mem.cv_gamma * lmem.savedJ[i][j];
-
-      lmem.A[i][i] += 1.0;
-    }
-  } else {
-    // Recompute Jacobian
-    lmem.nje++;
-    lmem.nstlj = mem.cv_nst;
-    lmem.jbad = false;
-    jcurPtr.value = true;
-
-    // Evaluate Jacobian
-    if (lmem.jacDQ)
-      retval = cvLsDQJac(mem, mem.cv_tn, ypred, fpred, lmem.savedJ, tmp1, tmp2);
-    else
-      retval = lmem.jacFn!(mem.cv_tn, ypred, fpred, lmem.savedJ, mem.cv_user_data);
-
-
-    if (retval < 0) {
-      lmem.jbad = true;
-      return -1;
-    }
-    if (retval > 0) {
-      lmem.jbad = true;
-      return 1;
-    }
-
-    // Form A = I - gamma*J
-    for (let i = 0; i < N; i++) {
-      for (let j = 0; j < N; j++)
-        lmem.A[i][j] = -mem.cv_gamma * lmem.savedJ[i][j];
-
-      lmem.A[i][i] += 1.0;
-    }
-  }
-
-  // Factor A using LU decomposition
-  const info = dgefa(lmem.A, N, lmem.pivots);
-  if (info !== 0) {
-    lmem.jbad = true;
-    return 1; // recoverable
-  }
-
-  return CV_SUCCESS;
-}
-
-/**
- * cvLsSolve - Solve the linear system A*x = b using LU factors.
- */
-function cvLsSolve(mem: CvodeMem, b: Float64Array, weight: Float64Array,
-  ycur: Float64Array, fcur: Float64Array): number {
-  const lmem = mem.cv_lmem!;
-  const N = mem.cv_N;
-
-  dgesl(lmem.A, N, lmem.pivots, b);
-
-  // Apply gamma/gammap correction if gamma has changed since last setup
-  if (mem.cv_gamrat !== 1.0) {
-    const factor = 2.0 / (1.0 + mem.cv_gamrat);
-    for (let i = 0; i < N; i++)
-      b[i] *= factor;
-  }
-
-  return CV_SUCCESS;
-}
-
-/**
- * cvLsDQJac - Compute difference-quotient Jacobian approximation.
- *
- * For each column j:
- *   sigma_j = max(|y_j|, 1/ewt_j) * sqrt(uround)
- *   J[:,j] = (f(t, y + sigma_j*e_j) - f(t, y)) / sigma_j
- */
-function cvLsDQJac(mem: CvodeMem, t: number, y: Float64Array, fy: Float64Array,
-  J: Float64Array[], tmp1: Float64Array, tmp2: Float64Array): number {
-  const N = mem.cv_N;
-  const lmem = mem.cv_lmem!;
-  const uround = 2.2204460492503131e-16;
-  const srur = Math.sqrt(uround);
-  const fnorm_val = wrmsNorm(N, fy, mem.cv_ewt);
-  const minInc = (fnorm_val !== 0.0) ? (1000.0 * Math.abs(mem.cv_h) * uround * N * fnorm_val) : 1.0;
-
-  // Use tmp1 as perturbed y, tmp2 as perturbed f
-  for (let j = 0; j < N; j++) {
-    const yjsaved = y[j];
-    let inc = Math.max(srur * Math.abs(yjsaved), minInc / mem.cv_ewt[j]);
-    // Ensure inc has same sign (positive since we want to increase)
-
-    // Perturb y[j]
-    tmp1.set(y);
-    tmp1[j] = yjsaved + inc;
-    inc = tmp1[j] - yjsaved; // exact increment after rounding
-
-    // Evaluate f(t, y + inc*e_j)
-    const retval = mem.cv_f!(t, tmp1, tmp2, mem.cv_user_data);
-    lmem.nfeDQ++;
-    mem.cv_nfe++;
-    if (retval !== 0) return retval;
-
-    // Compute column j of Jacobian: (f_perturbed - f) / inc
-    const inc_inv = 1.0 / inc;
-    for (let i = 0; i < N; i++)
-      J[i][j] = (tmp2[i] - fy[i]) * inc_inv;
-  }
-
-  return CV_SUCCESS;
-}
-
-/**
- * cvLsFree - Free linear solver memory (no-op in TS, GC handles it).
- */
-function cvLsFree(mem: CvodeMem): void {
-  mem.cv_lmem = null;
-}

package/src/solver-tools/cvode/cvode_nls.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"cvode_nls.js","sourceRoot":"","sources":["cvode_nls.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,OAAO,EAAW,QAAQ,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,EACpD,UAAU,EAAE,MAAM,EAAY,eAAe,EAAE,cAAc,EAAE,cAAc,EAC7E,aAAa,EAAE,kBAAkB,EAAE,gBAAgB,EACnD,cAAc,EAAiB,aAAa,EAC5C,UAAU,EAAE,cAAc,EAAE,aAAa,EACzC,UAAU,EAAE,MAAM,EAAE,IAAI,EAAS,MAAM,UAAU,CAAC;AAEpD;;;;;;;;;GASG;AACH,MAAM,UAAU,KAAK,CAAC,GAAa,EAAE,KAAa;IAChD,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IACnB,IAAI,SAAkB,CAAC;IAEvB,6CAA6C;IAC7C,IAAI,GAAG,CAAC,SAAS,EAAE,CAAC;QAClB,GAAG,CAAC,QAAQ,GAAG,CAAC,KAAK,KAAK,UAAU,IAAI,KAAK,KAAK,aAAa,CAAC,CAAC,CAAC;YAChE,cAAc,CAAC,CAAC;YAChB,aAAa,CAAC;QAEhB,SAAS,GAAG,CAAC,KAAK,KAAK,cAAc,CAAC,IAAI,CAAC,KAAK,KAAK,aAAa,CAAC;YACvD,CAAC,GAAG,CAAC,MAAM,KAAK,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,uBAAuB,CAAC;YACnD,CAAC,GAAG,CAAC,MAAM,IAAI,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,OAAO,CAAC;YAC1C,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,eAAe,CAAC,CAAC;IACpE,CAAC;SAAM,CAAC;QACN,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC;QACnB,SAAS,GAAG,KAAK,CAAC;IACpB,CAAC;IAED,yCAAyC;IACzC,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IAE5B,wBAAwB;IACxB,GAAG,CAAC,WAAW,GAAG,KAAK,CAAC;IAExB,oCAAoC;IACpC,IAAI,IAAY,CAAC;IACjB,IAAI,GAAG,CAAC,MAAM,KAAK,MAAM;QACvB,IAAI,GAAG,iBAAiB,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC;;QAEzC,IAAI,GAAG,qBAAqB,CAAC,GAAG,CAAC,CAAC;IAGpC,IAAI,IAAI,KAAK,UAAU;QAAE,OAAO,IAAI,CAAC;IAErC,6CAA6C;IAC7C,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,OAAO,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAE7D,oCAAoC;IACpC,IAAI,CAAC,GAAG,CAAC,WAAW;QAClB,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,CAAC,EAAE,GAAG,CAAC,OAAO,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC;IAGtD,yBAAyB;IACzB,GAAG,CAAC,OAAO,GAAG,KAAK,CAAC;IAEpB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;;GAIG;AACH,SAAS,iBAAiB,CAAC,GAAa,EAAE,SAAkB;IAC1D,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IACnB,IAAI,MAAc,CAAC;IAEnB,qCAAqC;IACrC,IAAI,SAAS,IAAI,GAAG,CAAC,SAAS,EAAE,CAAC;QAC/B,2CAA2C;QAC3C,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,OAAO,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QAE7D,sCAAsC;QACtC,MAAM,GAAG,GAAG,CAAC,IAAK,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;QACxE,GAAG,CAAC,MAAM,EAAE,CAAC;QACb,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,eAAe,CAAC;QACvC,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,aAAa,CAAC;QAErC,cAAc;QACd,MAAM,OAAO,GAAG,EAAC,KAAK,EAAE,KAAK,EAAC,CAAC;QAC/B,MAAM,GAAG,GAAG,CAAC,SAAU,CAAC,GAAG,EAAE,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,QAAQ,EAAE,OAAO,EACxE,GAAG,CAAC,SAAS,EAAE,GAAG,CAAC,SAAS,EAAE,GAAG,CAAC,SAAS,CAAC,CAAC;QAC/C,GAAG,CAAC,UAAU,EAAE,CAAC;QAEjB,GAAG,CAAC,OAAO,GAAG,OAAO,CAAC,KAAK,CAAC;QAC5B,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC;QACpB,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,QAAQ,CAAC;QAC7B,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC;QACnB,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,MAAM,CAAC;QAE1B,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,cAAc,CAAC;QACtC,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,kBAAkB,CAAC;IAC5C,CAAC;IAED,wBAAwB;IACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,GAAG,CAAC,MAAM,EAAE,CAAC;QAEb,mBAAmB;QACnB,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,OAAO,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QAE7D,oBAAoB;QACpB,MAAM,GAAG,GAAG,CAAC,IAAK,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;QACxE,GAAG,CAAC,MAAM,EAAE,CAAC;QACb,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,eAAe,CAAC;QACvC,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,aAAa,CAAC;QAErC,uDAAuD;QACvD,iDAAiD;QACjD,MAAM,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC;QACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YACxB,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QAGxF,0EAA0E;QAC1E,6DAA6D;QAC7D,gCAAgC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAEzC,MAAM,GAAG,GAAG,CAAC,SAAU,CAAC,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC;QACpE,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,cAAc,CAAC;QACtC,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,kBAAkB,CAAC;QAE1C,mDAAmD;QACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAEnD,mBAAmB;QACnB,MAAM,GAAG,GAAG,QAAQ,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC;QACvC,MAAM,cAAc,GAAG,aAAa,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAEhE,IAAI,cAAc,KAAK,UAAU;YAAE,OAAO,UAAU,CAAC;QACrD,IAAI,cAAc,KAAK,kBAAkB,EAAE,CAAC;YAC1C,GAAG,CAAC,MAAM,EAAE,CAAC;YACb,OAAO,kBAAkB,CAAC;QAC5B,CAAC;QACD,uCAAuC;IACzC,CAAC;IAED,6CAA6C;IAC7C,GAAG,CAAC,MAAM,EAAE,CAAC;IACb,OAAO,kBAAkB,CAAC;AAC5B,CAAC;AAED;;;GAGG;AACH,SAAS,qBAAqB,CAAC,GAAa;IAC1C,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IACnB,IAAI,MAAc,CAAC;IAEnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,CAAC,EAAE,EAAE,CAAC;QACpC,GAAG,CAAC,MAAM,EAAE,CAAC;QAEb,mBAAmB;QACnB,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,OAAO,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QAE7D,sCAAsC;QACtC,MAAM,GAAG,GAAG,CAAC,IAAK,CAAC,GAAG,CAAC,KAAK,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;QACxE,GAAG,CAAC,MAAM,EAAE,CAAC;QACb,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,eAAe,CAAC;QACvC,IAAI,MAAM,GAAG,CAAC;YAAE,OAAO,aAAa,CAAC;QAErC,0CAA0C;QAC1C,mCAAmC;QACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YACxB,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAEjE,oBAAoB;QACpB,MAAM,OAAO,GAAG,GAAG,CAAC,SAAS,CAAC;QAC9B,MAAM,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,QAAQ,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC;QAE7C,iCAAiC;QACjC,MAAM,KAAK,GAAG,GAAG,CAAC,QAAQ,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YACxB,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;QAGzC,cAAc;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QAExD,mBAAmB;QACnB,MAAM,GAAG,GAAG,QAAQ,CAAC,CAAC,EAAE,KAAK,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC;QAC3C,MAAM,cAAc,GAAG,aAAa,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAEhE,IAAI,cAAc,KAAK,UAAU;YAAE,OAAO,UAAU,CAAC;QACrD,IAAI,cAAc,KAAK,kBAAkB,EAAE,CAAC;YAC1C,GAAG,CAAC,MAAM,EAAE,CAAC;YACb,OAAO,kBAAkB,CAAC;QAC5B,CAAC;QACD,uCAAuC;IACzC,CAAC;IAED,yBAAyB;IACzB,GAAG,CAAC,MAAM,EAAE,CAAC;IACb,OAAO,kBAAkB,CAAC;AAC5B,CAAC;AAED;;;;;;;;;GASG;AACH,SAAS,aAAa,CAAC,GAAa,EAAE,GAAW,EAAE,CAAS,EAAE,GAAW;IACvE,mCAAmC;IACnC,IAAI,CAAC,GAAG,CAAC;QACP,GAAG,CAAC,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,QAAQ,EAAE,GAAG,GAAG,GAAG,CAAC,OAAO,CAAC,CAAC;IAGpE,uBAAuB;IACvB,MAAM,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,QAAQ,CAAC,GAAG,GAAG,CAAC;IAErD,IAAI,IAAI,IAAI,GAAG,EAAE,CAAC;QAChB,YAAY;QACZ,IAAI,CAAC,KAAK,CAAC;YACT,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC;;YAEnB,GAAG,CAAC,QAAQ,GAAG,QAAQ,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,OAAO,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC;QAE7D,GAAG,CAAC,WAAW,GAAG,IAAI,CAAC;QACvB,OAAO,UAAU,CAAC;IACpB,CAAC;IAED,uBAAuB;IACvB,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,GAAG,IAAI,GAAG,GAAG,CAAC,OAAO;QACpC,OAAO,kBAAkB,CAAC;IAG5B,wBAAwB;IACxB,GAAG,CAAC,OAAO,GAAG,GAAG,CAAC;IAClB,OAAO,gBAAgB,CAAC;AAC1B,CAAC"}