diff-grok 1.1.0 → 1.2.0

package/src/solver-tools/cvode/cvode_adams.js

@@ -0,0 +1,124 @@
+/**
+ * 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 { L_MAX } from './common';
+/**
+ * cvAltSum - Alternating sum: sum_{i=0}^{iend} (-1)^i * a[i] / (i + k)
+ */
+function cvAltSum(iend, a, k) {
+    if (iend < 0)
+        return 0.0;
+    let sum = 0.0;
+    let sign = 1;
+    for (let i = 0; i <= iend; i++) {
+        sum += sign * (a[i] / (i + k));
+        sign = -sign;
+    }
+    return sum;
+}
+/**
+ * cvAdamsStart - Generate product polynomial coefficients for Adams method
+ * Returns hsum.
+ */
+function cvAdamsStart(mem, m) {
+    let hsum = mem.cv_h;
+    m[0] = 1.0;
+    for (let i = 1; i <= mem.cv_q; i++)
+        m[i] = 0.0;
+    for (let j = 1; j < mem.cv_q; j++) {
+        if (j === mem.cv_q - 1 && mem.cv_qwait === 1) {
+            const sum = cvAltSum(mem.cv_q - 2, m, 2);
+            mem.cv_tq[1] = mem.cv_q * sum / m[mem.cv_q - 2];
+        }
+        const xi_inv = mem.cv_h / hsum;
+        for (let i = j; i >= 1; i--)
+            m[i] += m[i - 1] * xi_inv;
+        hsum += mem.cv_tau[j];
+    }
+    return hsum;
+}
+/**
+ * cvAdamsFinish - Complete Adams l and tq computation
+ */
+function cvAdamsFinish(mem, m, M, hsum) {
+    const M0_inv = 1.0 / M[0];
+    mem.cv_l[0] = 1.0;
+    for (let i = 1; i <= mem.cv_q; i++)
+        mem.cv_l[i] = M0_inv * (m[i - 1] / i);
+    const xi = hsum / mem.cv_h;
+    const xi_inv = 1.0 / xi;
+    mem.cv_tq[2] = M[1] * M0_inv / xi;
+    mem.cv_tq[5] = xi / mem.cv_l[mem.cv_q];
+    if (mem.cv_qwait === 1) {
+        for (let i = mem.cv_q; i >= 1; i--)
+            m[i] += m[i - 1] * xi_inv;
+        M[2] = cvAltSum(mem.cv_q, m, 2);
+        mem.cv_tq[3] = M[2] * M0_inv / mem.cv_L;
+    }
+    mem.cv_tq[4] = mem.cv_nlscoef / mem.cv_tq[2];
+}
+/**
+ * cvSetAdams - Compute Adams method coefficients l[] and tq[]
+ */
+export function cvSetAdams(mem) {
+    const m = new Float64Array(L_MAX);
+    const M = new Float64Array(3);
+    if (mem.cv_q === 1) {
+        mem.cv_l[0] = 1.0;
+        mem.cv_l[1] = 1.0;
+        mem.cv_tq[1] = 1.0;
+        mem.cv_tq[5] = 1.0;
+        mem.cv_tq[2] = 0.5;
+        mem.cv_tq[3] = 1.0 / 12.0;
+        mem.cv_tq[4] = mem.cv_nlscoef / mem.cv_tq[2];
+        return;
+    }
+    const hsum = cvAdamsStart(mem, m);
+    M[0] = cvAltSum(mem.cv_q - 1, m, 1);
+    M[1] = cvAltSum(mem.cv_q - 1, m, 2);
+    cvAdamsFinish(mem, m, M, hsum);
+}
+/**
+ * cvAdjustAdams - Adjust history array for Adams order change
+ */
+export function cvAdjustAdams(mem, deltaq) {
+    const N = mem.cv_N;
+    if (deltaq === 1) {
+        // On order increase, set new column of zn to zero
+        mem.cv_zn[mem.cv_L].fill(0.0, 0, N);
+        return;
+    }
+    // Order decrease: adjust zn[j] by multiples of zn[q]
+    for (let i = 0; i <= mem.cv_qmax; i++)
+        mem.cv_l[i] = 0.0;
+    mem.cv_l[1] = 1.0;
+    let hsum = 0.0;
+    for (let j = 1; j <= mem.cv_q - 2; j++) {
+        hsum += mem.cv_tau[j];
+        const xi = hsum / mem.cv_hscale;
+        for (let i = j + 1; i >= 1; i--)
+            mem.cv_l[i] = mem.cv_l[i] * xi + mem.cv_l[i - 1];
+    }
+    for (let j = 1; j <= mem.cv_q - 2; j++)
+        mem.cv_l[j + 1] = mem.cv_q * (mem.cv_l[j] / (j + 1));
+    if (mem.cv_q > 2) {
+        const znq = mem.cv_zn[mem.cv_q];
+        for (let j = 2; j < mem.cv_q; j++) {
+            const c = -mem.cv_l[j];
+            const znj = mem.cv_zn[j];
+            for (let i = 0; i < N; i++)
+                znj[i] += c * znq[i];
+        }
+    }
+}
+//# sourceMappingURL=cvode_adams.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"cvode_adams.js","sourceRoot":"","sources":["cvode_adams.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,OAAO,EAAW,KAAK,EAAC,MAAM,UAAU,CAAC;AAEzC;;GAEG;AACH,SAAS,QAAQ,CAAC,IAAY,EAAE,CAAe,EAAE,CAAS;IACxD,IAAI,IAAI,GAAG,CAAC;QAAE,OAAO,GAAG,CAAC;IACzB,IAAI,GAAG,GAAG,GAAG,CAAC;IACd,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,GAAG,IAAI,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC/B,IAAI,GAAG,CAAC,IAAI,CAAC;IACf,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;GAGG;AACH,SAAS,YAAY,CAAC,GAAa,EAAE,CAAe;IAClD,IAAI,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;IACpB,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAE/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,QAAQ,KAAK,CAAC,EAAE,CAAC;YAC7C,MAAM,GAAG,GAAG,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACzC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QAClD,CAAC;QACD,MAAM,MAAM,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YACzB,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QAE5B,IAAI,IAAI,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IACxB,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED;;GAEG;AACH,SAAS,aAAa,CAAC,GAAa,EAAE,CAAe,EAAE,CAAe,EAAE,IAAY;IAClF,MAAM,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAE1B,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE;QAChC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAGxC,MAAM,EAAE,GAAG,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;IAC3B,MAAM,MAAM,GAAG,GAAG,GAAG,EAAE,CAAC;IAExB,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,EAAE,CAAC;IAClC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAEvC,IAAI,GAAG,CAAC,QAAQ,KAAK,CAAC,EAAE,CAAC;QACvB,KAAK,IAAI,CAAC,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAChC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QAE5B,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAChC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,GAAG,CAAC,IAAI,CAAC;IAC1C,CAAC;IAED,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,UAAU,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;AAC/C,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,UAAU,CAAC,GAAa;IACtC,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC;IAClC,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAE9B,IAAI,GAAG,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;QACnB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAClB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAClB,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACnB,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACnB,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACnB,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC;QAC1B,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,UAAU,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QAC7C,OAAO;IACT,CAAC;IAED,MAAM,IAAI,GAAG,YAAY,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpC,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACpC,aAAa,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;AACjC,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,aAAa,CAAC,GAAa,EAAE,MAAc;IACzD,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IAEnB,IAAI,MAAM,KAAK,CAAC,EAAE,CAAC;QACjB,kDAAkD;QAClD,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACpC,OAAO;IACT,CAAC;IAED,qDAAqD;IACrD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,OAAO,EAAE,CAAC,EAAE;QAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACzD,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAClB,IAAI,IAAI,GAAG,GAAG,CAAC;IAEf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,IAAI,IAAI,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACtB,MAAM,EAAE,GAAG,IAAI,GAAG,GAAG,CAAC,SAAS,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAC7B,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACrD,CAAC;IAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,EAAE;QACpC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAGvD,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;QACjB,MAAM,GAAG,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACvB,MAAM,GAAG,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACnD,CAAC;IACH,CAAC;AACH,CAAC"}

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

@@ -0,0 +1,142 @@
+/**
+ * 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, L_MAX} from './common';
+
+/**
+ * cvAltSum - Alternating sum: sum_{i=0}^{iend} (-1)^i * a[i] / (i + k)
+ */
+function cvAltSum(iend: number, a: Float64Array, k: number): number {
+  if (iend < 0) return 0.0;
+  let sum = 0.0;
+  let sign = 1;
+  for (let i = 0; i <= iend; i++) {
+    sum += sign * (a[i] / (i + k));
+    sign = -sign;
+  }
+  return sum;
+}
+
+/**
+ * cvAdamsStart - Generate product polynomial coefficients for Adams method
+ * Returns hsum.
+ */
+function cvAdamsStart(mem: CvodeMem, m: Float64Array): number {
+  let hsum = mem.cv_h;
+  m[0] = 1.0;
+  for (let i = 1; i <= mem.cv_q; i++) m[i] = 0.0;
+
+  for (let j = 1; j < mem.cv_q; j++) {
+    if (j === mem.cv_q - 1 && mem.cv_qwait === 1) {
+      const sum = cvAltSum(mem.cv_q - 2, m, 2);
+      mem.cv_tq[1] = mem.cv_q * sum / m[mem.cv_q - 2];
+    }
+    const xi_inv = mem.cv_h / hsum;
+    for (let i = j; i >= 1; i--)
+      m[i] += m[i - 1] * xi_inv;
+
+    hsum += mem.cv_tau[j];
+  }
+  return hsum;
+}
+
+/**
+ * cvAdamsFinish - Complete Adams l and tq computation
+ */
+function cvAdamsFinish(mem: CvodeMem, m: Float64Array, M: Float64Array, hsum: number): void {
+  const M0_inv = 1.0 / M[0];
+
+  mem.cv_l[0] = 1.0;
+  for (let i = 1; i <= mem.cv_q; i++)
+    mem.cv_l[i] = M0_inv * (m[i - 1] / i);
+
+
+  const xi = hsum / mem.cv_h;
+  const xi_inv = 1.0 / xi;
+
+  mem.cv_tq[2] = M[1] * M0_inv / xi;
+  mem.cv_tq[5] = xi / mem.cv_l[mem.cv_q];
+
+  if (mem.cv_qwait === 1) {
+    for (let i = mem.cv_q; i >= 1; i--)
+      m[i] += m[i - 1] * xi_inv;
+
+    M[2] = cvAltSum(mem.cv_q, m, 2);
+    mem.cv_tq[3] = M[2] * M0_inv / mem.cv_L;
+  }
+
+  mem.cv_tq[4] = mem.cv_nlscoef / mem.cv_tq[2];
+}
+
+/**
+ * cvSetAdams - Compute Adams method coefficients l[] and tq[]
+ */
+export function cvSetAdams(mem: CvodeMem): void {
+  const m = new Float64Array(L_MAX);
+  const M = new Float64Array(3);
+
+  if (mem.cv_q === 1) {
+    mem.cv_l[0] = 1.0;
+    mem.cv_l[1] = 1.0;
+    mem.cv_tq[1] = 1.0;
+    mem.cv_tq[5] = 1.0;
+    mem.cv_tq[2] = 0.5;
+    mem.cv_tq[3] = 1.0 / 12.0;
+    mem.cv_tq[4] = mem.cv_nlscoef / mem.cv_tq[2];
+    return;
+  }
+
+  const hsum = cvAdamsStart(mem, m);
+  M[0] = cvAltSum(mem.cv_q - 1, m, 1);
+  M[1] = cvAltSum(mem.cv_q - 1, m, 2);
+  cvAdamsFinish(mem, m, M, hsum);
+}
+
+/**
+ * cvAdjustAdams - Adjust history array for Adams order change
+ */
+export function cvAdjustAdams(mem: CvodeMem, deltaq: number): void {
+  const N = mem.cv_N;
+
+  if (deltaq === 1) {
+    // On order increase, set new column of zn to zero
+    mem.cv_zn[mem.cv_L].fill(0.0, 0, N);
+    return;
+  }
+
+  // Order decrease: adjust zn[j] by multiples of zn[q]
+  for (let i = 0; i <= mem.cv_qmax; i++) mem.cv_l[i] = 0.0;
+  mem.cv_l[1] = 1.0;
+  let hsum = 0.0;
+
+  for (let j = 1; j <= mem.cv_q - 2; j++) {
+    hsum += mem.cv_tau[j];
+    const xi = hsum / mem.cv_hscale;
+    for (let i = j + 1; i >= 1; i--)
+      mem.cv_l[i] = mem.cv_l[i] * xi + mem.cv_l[i - 1];
+  }
+
+  for (let j = 1; j <= mem.cv_q - 2; j++)
+    mem.cv_l[j + 1] = mem.cv_q * (mem.cv_l[j] / (j + 1));
+
+
+  if (mem.cv_q > 2) {
+    const znq = mem.cv_zn[mem.cv_q];
+    for (let j = 2; j < mem.cv_q; j++) {
+      const c = -mem.cv_l[j];
+      const znj = mem.cv_zn[j];
+      for (let i = 0; i < N; i++) znj[i] += c * znq[i];
+    }
+  }
+}

package/src/solver-tools/cvode/cvode_bdf.js

@@ -0,0 +1,153 @@
+/**
+ * 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.
+ */
+/**
+ * cvSetBDF - Compute BDF method coefficients l[] and tq[]
+ *
+ * The polynomial l(x) = l_0 + l_1*x + ... + l_q*x^q is given by:
+ *   Lambda(x) = (1 + x/xi*_q) * product_{i=1}^{q-1} (1 + x/xi_i)
+ * where xi_i = (t_n - t_{n-i}) / h
+ *
+ * The tq array holds test quantities for error and convergence tests:
+ *   tq[1] = coefficient for error at order q-1
+ *   tq[2] = coefficient for error at order q
+ *   tq[3] = coefficient for error at order q+1
+ *   tq[4] = convergence test coefficient
+ *   tq[5] = derivative for order q+2 vector
+ */
+export function cvSetBDF(mem) {
+    let alpha0;
+    let alpha0_hat;
+    let xi_inv;
+    let xistar_inv;
+    let hsum;
+    mem.cv_l[0] = 1.0;
+    mem.cv_l[1] = 1.0;
+    xi_inv = 1.0;
+    xistar_inv = 1.0;
+    for (let i = 2; i <= mem.cv_q; i++)
+        mem.cv_l[i] = 0.0;
+    alpha0 = -1.0;
+    alpha0_hat = -1.0;
+    hsum = mem.cv_h;
+    if (mem.cv_q > 1) {
+        for (let j = 2; j < mem.cv_q; j++) {
+            hsum += mem.cv_tau[j - 1];
+            xi_inv = mem.cv_h / hsum;
+            alpha0 -= 1.0 / j;
+            for (let i = j; i >= 1; i--)
+                mem.cv_l[i] += mem.cv_l[i - 1] * xi_inv;
+        }
+        // j = q
+        alpha0 -= 1.0 / mem.cv_q;
+        xistar_inv = -mem.cv_l[1] - alpha0;
+        hsum += mem.cv_tau[mem.cv_q - 1];
+        xi_inv = mem.cv_h / hsum;
+        alpha0_hat = -mem.cv_l[1] - xi_inv;
+        for (let i = mem.cv_q; i >= 1; i--)
+            mem.cv_l[i] += mem.cv_l[i - 1] * xistar_inv;
+    }
+    cvSetTqBDF(mem, hsum, alpha0, alpha0_hat, xi_inv, xistar_inv);
+}
+/**
+ * cvSetTqBDF - Set test quantity array for BDF method
+ */
+function cvSetTqBDF(mem, hsum, alpha0, alpha0_hat, xi_inv, xistar_inv) {
+    const A1 = 1.0 - alpha0_hat + alpha0;
+    const A2 = 1.0 + mem.cv_q * A1;
+    mem.cv_tq[2] = Math.abs(A1 / (alpha0 * A2));
+    mem.cv_tq[5] = Math.abs(A2 * xistar_inv / (mem.cv_l[mem.cv_q] * xi_inv));
+    if (mem.cv_qwait === 1) {
+        if (mem.cv_q > 1) {
+            const C = xistar_inv / mem.cv_l[mem.cv_q];
+            const A3 = alpha0 + 1.0 / mem.cv_q;
+            const A4 = alpha0_hat + xi_inv;
+            const Cpinv = (1.0 - A4 + A3) / A3;
+            mem.cv_tq[1] = Math.abs(C * Cpinv);
+        }
+        else
+            mem.cv_tq[1] = 1.0;
+        hsum += mem.cv_tau[mem.cv_q];
+        const xi_inv2 = mem.cv_h / hsum;
+        const A5 = alpha0 - 1.0 / (mem.cv_q + 1);
+        const A6 = alpha0_hat - xi_inv2;
+        const Cppinv = (1.0 - A6 + A5) / A2;
+        mem.cv_tq[3] = Math.abs(Cppinv / (xi_inv2 * (mem.cv_q + 2) * A5));
+    }
+    mem.cv_tq[4] = mem.cv_nlscoef / mem.cv_tq[2];
+}
+/**
+ * cvIncreaseBDF - Adjust history array for BDF order increase
+ */
+export function cvIncreaseBDF(mem) {
+    const N = mem.cv_N;
+    for (let i = 0; i <= mem.cv_qmax; i++)
+        mem.cv_l[i] = 0.0;
+    mem.cv_l[2] = 1.0;
+    let alpha1 = 1.0;
+    let prod = 1.0;
+    let xiold = 1.0;
+    let alpha0 = -1.0;
+    let hsum = mem.cv_hscale;
+    if (mem.cv_q > 1) {
+        for (let j = 1; j < mem.cv_q; j++) {
+            hsum += mem.cv_tau[j + 1];
+            const xi = hsum / mem.cv_hscale;
+            prod *= xi;
+            alpha0 -= 1.0 / (j + 1);
+            alpha1 += 1.0 / xi;
+            for (let i = j + 2; i >= 2; i--)
+                mem.cv_l[i] = mem.cv_l[i] * xiold + mem.cv_l[i - 1];
+            xiold = xi;
+        }
+    }
+    const A1 = (-alpha0 - alpha1) / prod;
+    const znL = mem.cv_zn[mem.cv_L];
+    const znSaved = mem.cv_zn[mem.cv_indx_acor];
+    for (let i = 0; i < N; i++)
+        znL[i] = A1 * znSaved[i];
+    if (mem.cv_q > 1) {
+        for (let j = 2; j <= mem.cv_q; j++) {
+            const c = mem.cv_l[j];
+            const znj = mem.cv_zn[j];
+            for (let i = 0; i < N; i++)
+                znj[i] += c * znL[i];
+        }
+    }
+}
+/**
+ * cvDecreaseBDF - Adjust history array for BDF order decrease
+ */
+export function cvDecreaseBDF(mem) {
+    const N = mem.cv_N;
+    for (let i = 0; i <= mem.cv_qmax; i++)
+        mem.cv_l[i] = 0.0;
+    mem.cv_l[2] = 1.0;
+    let hsum = 0.0;
+    for (let j = 1; j <= mem.cv_q - 2; j++) {
+        hsum += mem.cv_tau[j];
+        const xi = hsum / mem.cv_hscale;
+        for (let i = j + 2; i >= 2; i--)
+            mem.cv_l[i] = mem.cv_l[i] * xi + mem.cv_l[i - 1];
+    }
+    if (mem.cv_q > 2) {
+        const znq = mem.cv_zn[mem.cv_q];
+        for (let j = 2; j < mem.cv_q; j++) {
+            const c = -mem.cv_l[j];
+            const znj = mem.cv_zn[j];
+            for (let i = 0; i < N; i++)
+                znj[i] += c * znq[i];
+        }
+    }
+}
+//# sourceMappingURL=cvode_bdf.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"cvode_bdf.js","sourceRoot":"","sources":["cvode_bdf.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAIH;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,QAAQ,CAAC,GAAa;IACpC,IAAI,MAAc,CAAC;IAAC,IAAI,UAAkB,CAAC;IAAC,IAAI,MAAc,CAAC;IAAC,IAAI,UAAkB,CAAC;IAAC,IAAI,IAAY,CAAC;IAEzG,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAClB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAClB,MAAM,GAAG,GAAG,CAAC;IACb,UAAU,GAAG,GAAG,CAAC;IACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE;QAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAEtD,MAAM,GAAG,CAAC,GAAG,CAAC;IACd,UAAU,GAAG,CAAC,GAAG,CAAC;IAClB,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;IAEhB,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;QACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,IAAI,IAAI,GAAG,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC1B,MAAM,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;YACzB,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC;YAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;gBACzB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QAC5C,CAAC;QAED,QAAQ;QACR,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,IAAI,CAAC;QACzB,UAAU,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QACnC,IAAI,IAAI,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QACjC,MAAM,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;QACzB,UAAU,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC;QAEnC,KAAK,IAAI,CAAC,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAChC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC;IAChD,CAAC;IAED,UAAU,CAAC,GAAG,EAAE,IAAI,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,CAAC,CAAC;AAChE,CAAC;AAED;;GAEG;AACH,SAAS,UAAU,CAAC,GAAa,EAAE,IAAY,EAAE,MAAc,EAC7D,UAAkB,EAAE,MAAc,EAAE,UAAkB;IACtD,MAAM,EAAE,GAAG,GAAG,GAAG,UAAU,GAAG,MAAM,CAAC;IACrC,MAAM,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,GAAG,EAAE,CAAC;IAC/B,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC;IAC5C,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,UAAU,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;IAEzE,IAAI,GAAG,CAAC,QAAQ,KAAK,CAAC,EAAE,CAAC;QACvB,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;YACjB,MAAM,CAAC,GAAG,UAAU,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;YAC1C,MAAM,EAAE,GAAG,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,IAAI,CAAC;YACnC,MAAM,EAAE,GAAG,UAAU,GAAG,MAAM,CAAC;YAC/B,MAAM,KAAK,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;YACnC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QACrC,CAAC;;YACC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAErB,IAAI,IAAI,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC7B,MAAM,OAAO,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;QAChC,MAAM,EAAE,GAAG,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QACzC,MAAM,EAAE,GAAG,UAAU,GAAG,OAAO,CAAC;QAChC,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;QACpC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACpE,CAAC;IACD,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,UAAU,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;AAC/C,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,aAAa,CAAC,GAAa;IACzC,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IAEnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,OAAO,EAAE,CAAC,EAAE;QAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACzD,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAClB,IAAI,MAAM,GAAG,GAAG,CAAC;IACjB,IAAI,IAAI,GAAG,GAAG,CAAC;IACf,IAAI,KAAK,GAAG,GAAG,CAAC;IAChB,IAAI,MAAM,GAAG,CAAC,GAAG,CAAC;IAClB,IAAI,IAAI,GAAG,GAAG,CAAC,SAAS,CAAC;IAEzB,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;QACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,IAAI,IAAI,GAAG,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC1B,MAAM,EAAE,GAAG,IAAI,GAAG,GAAG,CAAC,SAAS,CAAC;YAChC,IAAI,IAAI,EAAE,CAAC;YACX,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACxB,MAAM,IAAI,GAAG,GAAG,EAAE,CAAC;YACnB,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;gBAC7B,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAEtD,KAAK,GAAG,EAAE,CAAC;QACb,CAAC;IACH,CAAC;IAED,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,IAAI,CAAC;IACrC,MAAM,GAAG,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAChC,MAAM,OAAO,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;IAC5C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;IAErD,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;QACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YACnC,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACtB,MAAM,GAAG,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACnD,CAAC;IACH,CAAC;AACH,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,aAAa,CAAC,GAAa;IACzC,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IAEnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,OAAO,EAAE,CAAC,EAAE;QAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACzD,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAClB,IAAI,IAAI,GAAG,GAAG,CAAC;IAEf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,IAAI,IAAI,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACtB,MAAM,EAAE,GAAG,IAAI,GAAG,GAAG,CAAC,SAAS,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;YAC7B,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IACrD,CAAC;IAED,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;QACjB,MAAM,GAAG,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACvB,MAAM,GAAG,GAAG,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;gBAAE,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACnD,CAAC;IACH,CAAC;AACH,CAAC"}

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

@@ -0,0 +1,164 @@
+/**
+ * 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';
+
+/**
+ * cvSetBDF - Compute BDF method coefficients l[] and tq[]
+ *
+ * The polynomial l(x) = l_0 + l_1*x + ... + l_q*x^q is given by:
+ *   Lambda(x) = (1 + x/xi*_q) * product_{i=1}^{q-1} (1 + x/xi_i)
+ * where xi_i = (t_n - t_{n-i}) / h
+ *
+ * The tq array holds test quantities for error and convergence tests:
+ *   tq[1] = coefficient for error at order q-1
+ *   tq[2] = coefficient for error at order q
+ *   tq[3] = coefficient for error at order q+1
+ *   tq[4] = convergence test coefficient
+ *   tq[5] = derivative for order q+2 vector
+ */
+export function cvSetBDF(mem: CvodeMem): void {
+  let alpha0: number; let alpha0_hat: number; let xi_inv: number; let xistar_inv: number; let hsum: number;
+
+  mem.cv_l[0] = 1.0;
+  mem.cv_l[1] = 1.0;
+  xi_inv = 1.0;
+  xistar_inv = 1.0;
+  for (let i = 2; i <= mem.cv_q; i++) mem.cv_l[i] = 0.0;
+
+  alpha0 = -1.0;
+  alpha0_hat = -1.0;
+  hsum = mem.cv_h;
+
+  if (mem.cv_q > 1) {
+    for (let j = 2; j < mem.cv_q; j++) {
+      hsum += mem.cv_tau[j - 1];
+      xi_inv = mem.cv_h / hsum;
+      alpha0 -= 1.0 / j;
+      for (let i = j; i >= 1; i--)
+        mem.cv_l[i] += mem.cv_l[i - 1] * xi_inv;
+    }
+
+    // j = q
+    alpha0 -= 1.0 / mem.cv_q;
+    xistar_inv = -mem.cv_l[1] - alpha0;
+    hsum += mem.cv_tau[mem.cv_q - 1];
+    xi_inv = mem.cv_h / hsum;
+    alpha0_hat = -mem.cv_l[1] - xi_inv;
+
+    for (let i = mem.cv_q; i >= 1; i--)
+      mem.cv_l[i] += mem.cv_l[i - 1] * xistar_inv;
+  }
+
+  cvSetTqBDF(mem, hsum, alpha0, alpha0_hat, xi_inv, xistar_inv);
+}
+
+/**
+ * cvSetTqBDF - Set test quantity array for BDF method
+ */
+function cvSetTqBDF(mem: CvodeMem, hsum: number, alpha0: number,
+  alpha0_hat: number, xi_inv: number, xistar_inv: number): void {
+  const A1 = 1.0 - alpha0_hat + alpha0;
+  const A2 = 1.0 + mem.cv_q * A1;
+  mem.cv_tq[2] = Math.abs(A1 / (alpha0 * A2));
+  mem.cv_tq[5] = Math.abs(A2 * xistar_inv / (mem.cv_l[mem.cv_q] * xi_inv));
+
+  if (mem.cv_qwait === 1) {
+    if (mem.cv_q > 1) {
+      const C = xistar_inv / mem.cv_l[mem.cv_q];
+      const A3 = alpha0 + 1.0 / mem.cv_q;
+      const A4 = alpha0_hat + xi_inv;
+      const Cpinv = (1.0 - A4 + A3) / A3;
+      mem.cv_tq[1] = Math.abs(C * Cpinv);
+    } else
+      mem.cv_tq[1] = 1.0;
+
+    hsum += mem.cv_tau[mem.cv_q];
+    const xi_inv2 = mem.cv_h / hsum;
+    const A5 = alpha0 - 1.0 / (mem.cv_q + 1);
+    const A6 = alpha0_hat - xi_inv2;
+    const Cppinv = (1.0 - A6 + A5) / A2;
+    mem.cv_tq[3] = Math.abs(Cppinv / (xi_inv2 * (mem.cv_q + 2) * A5));
+  }
+  mem.cv_tq[4] = mem.cv_nlscoef / mem.cv_tq[2];
+}
+
+/**
+ * cvIncreaseBDF - Adjust history array for BDF order increase
+ */
+export function cvIncreaseBDF(mem: CvodeMem): void {
+  const N = mem.cv_N;
+
+  for (let i = 0; i <= mem.cv_qmax; i++) mem.cv_l[i] = 0.0;
+  mem.cv_l[2] = 1.0;
+  let alpha1 = 1.0;
+  let prod = 1.0;
+  let xiold = 1.0;
+  let alpha0 = -1.0;
+  let hsum = mem.cv_hscale;
+
+  if (mem.cv_q > 1) {
+    for (let j = 1; j < mem.cv_q; j++) {
+      hsum += mem.cv_tau[j + 1];
+      const xi = hsum / mem.cv_hscale;
+      prod *= xi;
+      alpha0 -= 1.0 / (j + 1);
+      alpha1 += 1.0 / xi;
+      for (let i = j + 2; i >= 2; i--)
+        mem.cv_l[i] = mem.cv_l[i] * xiold + mem.cv_l[i - 1];
+
+      xiold = xi;
+    }
+  }
+
+  const A1 = (-alpha0 - alpha1) / prod;
+  const znL = mem.cv_zn[mem.cv_L];
+  const znSaved = mem.cv_zn[mem.cv_indx_acor];
+  for (let i = 0; i < N; i++) znL[i] = A1 * znSaved[i];
+
+  if (mem.cv_q > 1) {
+    for (let j = 2; j <= mem.cv_q; j++) {
+      const c = mem.cv_l[j];
+      const znj = mem.cv_zn[j];
+      for (let i = 0; i < N; i++) znj[i] += c * znL[i];
+    }
+  }
+}
+
+/**
+ * cvDecreaseBDF - Adjust history array for BDF order decrease
+ */
+export function cvDecreaseBDF(mem: CvodeMem): void {
+  const N = mem.cv_N;
+
+  for (let i = 0; i <= mem.cv_qmax; i++) mem.cv_l[i] = 0.0;
+  mem.cv_l[2] = 1.0;
+  let hsum = 0.0;
+
+  for (let j = 1; j <= mem.cv_q - 2; j++) {
+    hsum += mem.cv_tau[j];
+    const xi = hsum / mem.cv_hscale;
+    for (let i = j + 2; i >= 2; i--)
+      mem.cv_l[i] = mem.cv_l[i] * xi + mem.cv_l[i - 1];
+  }
+
+  if (mem.cv_q > 2) {
+    const znq = mem.cv_zn[mem.cv_q];
+    for (let j = 2; j < mem.cv_q; j++) {
+      const c = -mem.cv_l[j];
+      const znj = mem.cv_zn[j];
+      for (let i = 0; i < N; i++) znj[i] += c * znq[i];
+    }
+  }
+}