diff-grok 1.1.0 → 1.3.0

package/src/solver-tools/lsoda/methodswitch.ts

@@ -1,107 +0,0 @@
-import {LsodaContext, ETA, RATIO, sm1, cm1, cm2} from './common';
-import {vmnorm} from './blas';
-
-/** Mutable rh output */
-export interface MethodSwitchResult {
-  rh: number;
-}
-
-/**
- * Consider switching between Adams (meth=1) and BDF (meth=2).
- * Modifies ctx.common.meth, miter, nq, icount, pdlast on switch.
- */
-export function methodswitch(
-  ctx: LsodaContext, dsm: number, pnorm: number, out: MethodSwitchResult,
-): void {
-  const c = ctx.common!;
-  const neq = ctx.neq;
-  const mxordn = ctx.opt!.mxordn;
-  const mxords = ctx.opt!.mxords;
-
-  if (c.meth === 1) {
-    // Currently Adams, consider switching to BDF
-    if (c.nq > 5) return;
-
-    let rh2: number;
-    let nqm2: number;
-
-    if (dsm <= (100. * pnorm * ETA) || c.pdest === 0.) {
-      if (c.irflag === 0) return;
-      rh2 = 2.;
-      nqm2 = Math.min(c.nq, mxords);
-    } else {
-      const exsm = 1. / (c.nq + 1);
-      let rh1 = 1. / (1.2 * Math.pow(dsm, exsm) + 0.0000012);
-      let rh1it = 2. * rh1;
-      const pdh = c.pdlast * Math.abs(c.h);
-      if ((pdh * rh1) > 0.00001)
-        rh1it = sm1[c.nq] / pdh;
-      rh1 = Math.min(rh1, rh1it);
-
-      if (c.nq > mxords) {
-        nqm2 = mxords;
-        const lm2 = mxords + 1;
-        const exm2 = 1. / lm2;
-        const lm2p1 = lm2 + 1;
-        const dm2 = vmnorm(neq, c.yh[lm2p1], c.ewt) / cm2[mxords];
-        rh2 = 1. / (1.2 * Math.pow(dm2, exm2) + 0.0000012);
-      } else {
-        const dm2 = dsm * (cm1[c.nq] / cm2[c.nq]);
-        rh2 = 1. / (1.2 * Math.pow(dm2, exsm) + 0.0000012);
-        nqm2 = c.nq;
-      }
-      if (rh2 < RATIO * rh1) return;
-    }
-
-    // Switch to BDF
-    out.rh = rh2;
-    c.icount = 20;
-    c.meth = 2;
-    c.miter = 2;
-    c.pdlast = 0.;
-    c.nq = nqm2;
-    return;
-  }
-
-  // Currently BDF (meth=2), consider switching to Adams
-  const exsm = 1. / (c.nq + 1);
-  let nqm1: number;
-  let exm1: number;
-  let rh1: number;
-  let dm1: number;
-
-  if (mxordn < c.nq) {
-    nqm1 = mxordn;
-    const lm1 = mxordn + 1;
-    exm1 = 1. / lm1;
-    const lm1p1 = lm1 + 1;
-    dm1 = vmnorm(neq, c.yh[lm1p1], c.ewt) / cm1[mxordn];
-    rh1 = 1. / (1.2 * Math.pow(dm1, exm1) + 0.0000012);
-  } else {
-    dm1 = dsm * (cm2[c.nq] / cm1[c.nq]);
-    rh1 = 1. / (1.2 * Math.pow(dm1, exsm) + 0.0000012);
-    nqm1 = c.nq;
-    exm1 = exsm;
-  }
-
-  let rh1it = 2. * rh1;
-  const pdh = c.pdnorm * Math.abs(c.h);
-  if ((pdh * rh1) > 0.00001)
-    rh1it = sm1[nqm1] / pdh;
-  rh1 = Math.min(rh1, rh1it);
-
-  const rh2 = 1. / (1.2 * Math.pow(dsm, exsm) + 0.0000012);
-  if ((rh1 * RATIO) < (5. * rh2)) return;
-
-  const alpha = Math.max(0.001, rh1);
-  dm1 *= Math.pow(alpha, exm1);
-  if (dm1 <= 1000. * ETA * pnorm) return;
-
-  // Switch to Adams
-  out.rh = rh1;
-  c.icount = 20;
-  c.meth = 1;
-  c.miter = 0;
-  c.pdlast = 0.;
-  c.nq = nqm1;
-}

package/src/solver-tools/lsoda/orderswitch.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"orderswitch.js","sourceRoot":"","sources":["orderswitch.ts"],"names":[],"mappings":"AAAA,OAAO,EAAe,GAAG,EAAC,MAAM,UAAU,CAAC;AAC3C,OAAO,EAAC,MAAM,EAAC,MAAM,QAAQ,CAAC;AAO9B;;;GAGG;AACH,MAAM,UAAU,WAAW,CACzB,GAAiB,EAAE,IAAY,EAAE,GAAW,EAC5C,GAAsB,EAAE,KAAa,EAAE,MAAc;IAErD,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IAEpB,MAAM,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IAC7B,IAAI,IAAI,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC;IAExD,IAAI,IAAI,GAAG,EAAE,CAAC;IACd,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;QACf,MAAM,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAClE,MAAM,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;QACvB,IAAI,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC;IACtD,CAAC;IAED,sDAAsD;IACtD,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;QACjB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,QAAQ,CAAC,CAAC;QACzD,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC;YACvB,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QAC7C,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;QACvC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC;YACV,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QAC7C,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC;IACf,CAAC;IAED,IAAI,IAAY,CAAC;IACjB,IAAI,IAAI,IAAI,IAAI,EAAE,CAAC;QACjB,IAAI,IAAI,IAAI,IAAI,EAAE,CAAC;YACjB,IAAI,GAAG,CAAC,CAAC,EAAE,CAAC;YACZ,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC;QAChB,CAAC;aAAM,CAAC;YACN,IAAI,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;YAChB,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC;YACd,IAAI,KAAK,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE;gBAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC;QAC5C,CAAC;IACH,CAAC;SAAM,CAAC;QACN,IAAI,IAAI,IAAI,IAAI,EAAE,CAAC;YACjB,IAAI,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;YAChB,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC;YACd,IAAI,KAAK,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,GAAG,EAAE;gBAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC;QAC5C,CAAC;aAAM,CAAC;YACN,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC;YACd,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,EAAE,CAAC;gBAClB,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;gBACtC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;oBAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACpC,OAAO,CAAC,CAAC;YACX,CAAC;iBAAM,CAAC;gBACN,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;gBACZ,OAAO,CAAC,CAAC;YACX,CAAC;QACH,CAAC;IACH,CAAC;IAED,gEAAgE;IAChE,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;QACjB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,QAAQ,CAAC,CAAC;QACzD,IAAI,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,OAAO,GAAG,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;YACvC,IAAI,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC;gBAChC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;gBACZ,OAAO,CAAC,CAAC;YACX,CAAC;QACH,CAAC;IACH,CAAC;SAAM,CAAC;QACN,IAAI,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC;YAChC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;YACZ,OAAO,CAAC,CAAC;QACX,CAAC;IACH,CAAC;IAED,IAAI,KAAK,IAAI,CAAC,CAAC;QACb,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IAEjC,IAAI,IAAI,KAAK,CAAC,CAAC,EAAE;QACf,OAAO,CAAC,CAAC;IAEX,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC;IACZ,OAAO,CAAC,CAAC;AACX,CAAC"}

package/src/solver-tools/lsoda/orderswitch.ts

@@ -1,95 +0,0 @@
-import {LsodaContext, sm1} from './common';
-import {vmnorm} from './blas';
-
-/** Mutable rh output */
-export interface OrderSwitchResult {
-  rh: number;
-}
-
-/**
- * Decide whether to change the order and/or step size.
- * Returns orderflag: 0 = no change, 1 = change h only, 2 = change both h and nq.
- */
-export function orderswitch(
-  ctx: LsodaContext, rhup: number, dsm: number,
-  out: OrderSwitchResult, kflag: number, maxord: number,
-): number {
-  const c = ctx.common!;
-  const neq = ctx.neq;
-
-  const exsm = 1. / (c.nq + 1);
-  let rhsm = 1. / (1.2 * Math.pow(dsm, exsm) + 0.0000012);
-
-  let rhdn = 0.;
-  if (c.nq !== 1) {
-    const ddn = vmnorm(neq, c.yh[c.nq + 1], c.ewt) / c.tesco[c.nq][1];
-    const exdn = 1. / c.nq;
-    rhdn = 1. / (1.3 * Math.pow(ddn, exdn) + 0.0000013);
-  }
-
-  // If meth = 1, limit rh according to stability region
-  if (c.meth === 1) {
-    const pdh = Math.max(Math.abs(c.h) * c.pdlast, 0.000001);
-    if (c.nq + 1 < maxord + 1)
-      rhup = Math.min(rhup, sm1[c.nq + 1] / pdh);
-    rhsm = Math.min(rhsm, sm1[c.nq] / pdh);
-    if (c.nq > 1)
-      rhdn = Math.min(rhdn, sm1[c.nq - 1] / pdh);
-    c.pdest = 0.;
-  }
-
-  let newq: number;
-  if (rhsm >= rhup) {
-    if (rhsm >= rhdn) {
-      newq = c.nq;
-      out.rh = rhsm;
-    } else {
-      newq = c.nq - 1;
-      out.rh = rhdn;
-      if (kflag < 0 && out.rh > 1.) out.rh = 1.;
-    }
-  } else {
-    if (rhup <= rhdn) {
-      newq = c.nq - 1;
-      out.rh = rhdn;
-      if (kflag < 0 && out.rh > 1.) out.rh = 1.;
-    } else {
-      out.rh = rhup;
-      if (out.rh >= 1.1) {
-        const r = c.el[c.nq + 1] / (c.nq + 1);
-        c.nq = c.nq + 1;
-        for (let i = 1; i <= neq; i++)
-          c.yh[c.nq + 1][i] = c.acor[i] * r;
-        return 2;
-      } else {
-        c.ialth = 3;
-        return 0;
-      }
-    }
-  }
-
-  // If meth = 1 and h is restricted by stability, bypass 10% test
-  if (c.meth === 1) {
-    const pdh = Math.max(Math.abs(c.h) * c.pdlast, 0.000001);
-    if (out.rh * pdh * 1.00001 < sm1[newq]) {
-      if (kflag === 0 && out.rh < 1.1) {
-        c.ialth = 3;
-        return 0;
-      }
-    }
-  } else {
-    if (kflag === 0 && out.rh < 1.1) {
-      c.ialth = 3;
-      return 0;
-    }
-  }
-
-  if (kflag <= -2)
-    out.rh = Math.min(out.rh, 0.2);
-
-  if (newq === c.nq)
-    return 1;
-
-  c.nq = newq;
-  return 2;
-}

package/src/solver-tools/lsoda/prja.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"prja.js","sourceRoot":"","sources":["prja.ts"],"names":[],"mappings":"AAAA,OAAO,EAAe,GAAG,EAAE,OAAO,EAAC,MAAM,UAAU,CAAC;AACpD,OAAO,EAAC,MAAM,EAAE,KAAK,EAAE,KAAK,EAAC,MAAM,QAAQ,CAAC;AAE5C;;;;GAIG;AACH,MAAM,UAAU,IAAI,CAAC,GAAiB,EAAE,CAAe;IACrD,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IAEpB,CAAC,CAAC,GAAG,EAAE,CAAC;IACR,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAE1B,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QAClB,OAAO,CAAC,KAAK,CAAC,mBAAmB,CAAC,CAAC;QACnC,OAAO,CAAC,CAAC;IACX,CAAC;IAED,IAAI,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;IACrC,IAAI,EAAE,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;IACjD,IAAI,EAAE,KAAK,EAAE;QAAE,EAAE,GAAG,EAAE,CAAC;IAEvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1D,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QACV,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;QACf,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;QAC5D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAC7C,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACZ,CAAC;IACD,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC;IAEb,2BAA2B;IAC3B,CAAC,CAAC,MAAM,GAAG,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAEnD,sBAAsB;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;QAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;IAEnB,wBAAwB;IACxB,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;IACrC,IAAI,GAAG,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IAExB,OAAO,CAAC,CAAC;AACX,CAAC"}

package/src/solver-tools/lsoda/prja.ts

@@ -1,49 +0,0 @@
-import {LsodaContext, ETA, SQRTETA} from './common';
-import {vmnorm, fnorm, dgefa} from './blas';
-
-/**
- * Compute and process the matrix P = I - h * el[1] * J,
- * where J is approximated by finite differences.
- * Returns 1 on success, 0 on failure (singular matrix).
- */
-export function prja(ctx: LsodaContext, y: Float64Array): number {
-  const c = ctx.common!;
-  const neq = ctx.neq;
-
-  c.nje++;
-  const hl0 = c.h * c.el[1];
-
-  if (c.miter !== 2) {
-    console.error('[prja] miter != 2');
-    return 0;
-  }
-
-  let fac = vmnorm(neq, c.savf, c.ewt);
-  let r0 = 1000. * Math.abs(c.h) * ETA * neq * fac;
-  if (r0 === 0.) r0 = 1.;
-
-  for (let j = 1; j <= neq; j++) {
-    const yj = y[j];
-    const r = Math.max(SQRTETA * Math.abs(yj), r0 / c.ewt[j]);
-    y[j] += r;
-    fac = -hl0 / r;
-    ctx.func(c.tn, y.subarray(1), c.acor.subarray(1), ctx.data);
-    for (let i = 1; i <= neq; i++)
-      c.wm[i][j] = (c.acor[i] - c.savf[i]) * fac;
-    y[j] = yj;
-  }
-  c.nfe += neq;
-
-  // Compute norm of Jacobian
-  c.pdnorm = fnorm(neq, c.wm, c.ewt) / Math.abs(hl0);
-
-  // Add identity matrix
-  for (let i = 1; i <= neq; i++)
-    c.wm[i][i] += 1.;
-
-  // LU decomposition on P
-  const ier = dgefa(c.wm, neq, c.ipvt);
-  if (ier !== 0) return 0;
-
-  return 1;
-}

package/src/solver-tools/lsoda/scaleh.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"scaleh.js","sourceRoot":"","sources":["scaleh.ts"],"names":[],"mappings":"AAAA,OAAO,EAAe,GAAG,EAAC,MAAM,UAAU,CAAC;AAE3C;;GAEG;AACH,MAAM,UAAU,MAAM,CAAC,GAAiB,EAAE,EAAU;IAClD,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IACpB,MAAM,IAAI,GAAG,GAAG,CAAC,GAAI,CAAC,IAAI,CAAC;IAE3B,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;IAC1B,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC;IAElD,4CAA4C;IAC5C,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;QACjB,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;QACb,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,QAAQ,CAAC,CAAC;QACzD,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,OAAO,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC;YACtC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;YACrB,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;QACf,CAAC;IACH,CAAC;IAED,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACnC,CAAC,IAAI,EAAE,CAAC;QACR,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IACpB,CAAC;IACD,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;IACV,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC;IACX,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;AACrB,CAAC"}

package/src/solver-tools/lsoda/scaleh.ts

@@ -1,33 +0,0 @@
-import {LsodaContext, sm1} from './common';
-
-/**
- * Scale step size. Rescales the yh array when h changes.
- */
-export function scaleh(ctx: LsodaContext, rh: number): void {
-  const c = ctx.common!;
-  const neq = ctx.neq;
-  const hmxi = ctx.opt!.hmxi;
-
-  rh = Math.min(rh, c.rmax);
-  rh = rh / Math.max(1., Math.abs(c.h) * hmxi * rh);
-
-  // If meth = 1, restrict by stability region
-  if (c.meth === 1) {
-    c.irflag = 0;
-    const pdh = Math.max(Math.abs(c.h) * c.pdlast, 0.000001);
-    if ((rh * pdh * 1.00001) >= sm1[c.nq]) {
-      rh = sm1[c.nq] / pdh;
-      c.irflag = 1;
-    }
-  }
-
-  let r = 1.;
-  for (let j = 2; j <= c.nq + 1; j++) {
-    r *= rh;
-    for (let i = 1; i <= neq; i++)
-      c.yh[j][i] *= r;
-  }
-  c.h *= rh;
-  c.rc *= rh;
-  c.ialth = c.nq + 1;
-}

package/src/solver-tools/lsoda/solsy.js

@@ -1,14 +0,0 @@
-import { dgesl } from './blas';
-/**
- * Solve the linear system arising from a chord iteration.
- * Only miter=2 (dense) is supported.
- */
-export function solsy(ctx, y) {
-    const c = ctx.common;
-    const neq = ctx.neq;
-    if (c.miter !== 2)
-        throw new Error('[solsy] miter != 2 not implemented');
-    dgesl(c.wm, neq, c.ipvt, y, 0);
-    return 1;
-}
-//# sourceMappingURL=solsy.js.map

package/src/solver-tools/lsoda/solsy.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"solsy.js","sourceRoot":"","sources":["solsy.ts"],"names":[],"mappings":"AACA,OAAO,EAAC,KAAK,EAAC,MAAM,QAAQ,CAAC;AAE7B;;;GAGG;AACH,MAAM,UAAU,KAAK,CAAC,GAAiB,EAAE,CAAe;IACtD,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IAEpB,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC;QACf,MAAM,IAAI,KAAK,CAAC,oCAAoC,CAAC,CAAC;IAGxD,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC/B,OAAO,CAAC,CAAC;AACX,CAAC"}

package/src/solver-tools/lsoda/solsy.ts

@@ -1,18 +0,0 @@
-import {LsodaContext} from './common';
-import {dgesl} from './blas';
-
-/**
- * Solve the linear system arising from a chord iteration.
- * Only miter=2 (dense) is supported.
- */
-export function solsy(ctx: LsodaContext, y: Float64Array): number {
-  const c = ctx.common!;
-  const neq = ctx.neq;
-
-  if (c.miter !== 2)
-    throw new Error('[solsy] miter != 2 not implemented');
-
-
-  dgesl(c.wm, neq, c.ipvt, y, 0);
-  return 1;
-}

package/src/solver-tools/lsoda/stoda.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"stoda.js","sourceRoot":"","sources":["stoda.ts"],"names":[],"mappings":"AAAA,OAAO,EAAe,KAAK,EAAE,IAAI,EAAC,MAAM,UAAU,CAAC;AACnD,OAAO,EAAC,MAAM,EAAC,MAAM,QAAQ,CAAC;AAC9B,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAC9B,OAAO,EAAC,MAAM,EAAC,MAAM,UAAU,CAAC;AAChC,OAAO,EAAC,UAAU,EAAkB,MAAM,cAAc,CAAC;AACzD,OAAO,EAAC,WAAW,EAAoB,MAAM,eAAe,CAAC;AAC7D,OAAO,EAAC,YAAY,EAAqB,MAAM,gBAAgB,CAAC;AAEhE;;;;;GAKG;AACH,MAAM,UAAU,KAAK,CAAC,GAAiB,EAAE,CAAe,EAAE,MAAc;IACtE,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,IAAI,GAAG,GAAG,CAAC,GAAI,CAAC,IAAI,CAAC;IAC3B,MAAM,MAAM,GAAG,GAAG,CAAC,GAAI,CAAC,MAAM,CAAC;IAC/B,MAAM,MAAM,GAAG,GAAG,CAAC,GAAI,CAAC,MAAM,CAAC;IAC/B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IAEpB,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,MAAM,IAAI,GAAG,CAAC,CAAC,EAAE,CAAC;IAClB,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IACV,IAAI,IAAI,GAAG,EAAE,CAAC;IAEd,IAAI,MAAM,GAAG,MAAM,CAAC;IACpB,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;QAAE,MAAM,GAAG,MAAM,CAAC;IAElC,+CAA+C;IAC/C,SAAS,QAAQ;QACf,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAC3B,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QACjB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;IACf,CAAC;IAED,6DAA6D;IAC7D,SAAS,UAAU;QACjB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE;YAChC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5B,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAC9B,CAAC;IAED,yCAAyC;IAEzC,IAAI,MAAM,KAAK,CAAC,EAAE,CAAC;QACjB,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACT,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;QACZ,CAAC,CAAC,IAAI,GAAG,MAAM,CAAC;QAChB,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC;QACV,CAAC,CAAC,KAAK,GAAG,GAAG,CAAC;QACd,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;QACb,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;QACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;QACjB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,CAAC,CAAC,MAAM,GAAG,EAAE,CAAC;QACd,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;QACb,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC;QACb,CAAC,CAAC,MAAM,GAAG,EAAE,CAAC;QACd,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACd,UAAU,EAAE,CAAC;IACf,CAAC;IAED,IAAI,MAAM,KAAK,CAAC,CAAC,EAAE,CAAC;QAClB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;QACjB,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC;YAAE,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;QAC/B,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,EAAE,CAAC;YACvB,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;YACnB,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;YACnB,UAAU,EAAE,CAAC;QACf,CAAC;QACD,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC;YACnB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACxB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACb,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAClB,CAAC;IACH,CAAC;IAED,IAAI,MAAM,KAAK,CAAC,CAAC,EAAE,CAAC;QAClB,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC;YACnB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACxB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACb,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAClB,CAAC;IACH,CAAC;IAED,gCAAgC;IAEhC,IAAI,GAAG,GAAG,GAAG,CAAC;IACd,MAAM,EAAE,GAAoB,EAAC,GAAG,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAC,CAAC;IACvD,MAAM,QAAQ,GAAsB,EAAC,EAAE,EAAE,CAAC,EAAC,CAAC;IAC5C,MAAM,QAAQ,GAAuB,EAAC,EAAE,EAAE,CAAC,EAAC,CAAC;IAC7C,IAAI,EAAU,CAAC;IAEf,KAAK,EACL,OAAO,IAAI,EAAE,CAAC;QACZ,8BAA8B;QAC9B,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;QAEX,KAAK,EACL,OAAO,IAAI,EAAE,CAAC;YACZ,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,KAAK;gBAC7B,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;YACnB,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,IAAI,GAAG,IAAI;gBACxB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;YAEnB,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;YAEZ,oDAAoD;YACpD,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;oBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnC,CAAC;YACH,CAAC;YAED,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;YAE1C,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;YACX,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC;YACf,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;YACT,MAAM,OAAO,GAAG,UAAU,CAAC,GAAG,EAAE,CAAC,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;YACpD,IAAI,GAAG,EAAE,CAAC,IAAI,CAAC;YAEf,IAAI,OAAO,KAAK,CAAC;gBAAE,MAAM,KAAK,CAAC;YAE/B,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;gBAClB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1C,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;gBAChB,SAAS,KAAK,CAAC;YACjB,CAAC;YAED,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;gBAClB,KAAK,GAAG,CAAC,CAAC,CAAC;gBACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;gBACb,MAAM,GAAG,CAAC,CAAC;gBACX,OAAO,KAAK,CAAC;YACf,CAAC;QACH,CAAC;QAED,mBAAmB;QACnB,IAAI,EAAE,CAAC,CAAC,KAAK,CAAC;YACZ,GAAG,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAClC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC;YACV,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAEtD,IAAI,GAAG,IAAI,EAAE,EAAE,CAAC;YACd,kBAAkB;YAClB,KAAK,GAAG,CAAC,CAAC;YACV,CAAC,CAAC,GAAG,EAAE,CAAC;YACR,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACX,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC;YACb,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,IAAI,CAAC;YAEjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACnC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;oBAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAChC,CAAC;YAED,CAAC,CAAC,MAAM,EAAE,CAAC;YACX,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;gBACjB,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,YAAY,CAAC,GAAG,EAAE,GAAG,EAAE,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,CAAC;gBAC9D,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,EAAE,CAAC;oBACvB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;oBAChB,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC;oBACb,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;YACH,CAAC;YAED,6CAA6C;YAC7C,CAAC,CAAC,KAAK,EAAE,CAAC;YACV,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;gBAClB,IAAI,IAAI,GAAG,EAAE,CAAC;gBACd,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,KAAK,MAAM,GAAG,CAAC,EAAE,CAAC;oBAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC9C,MAAM,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC1D,MAAM,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;oBAC7B,IAAI,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC;gBACtD,CAAC;gBAED,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,MAAM,SAAS,GAAG,WAAW,CAAC,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;gBAEvE,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;oBAChB,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC;oBACb,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,UAAU,EAAE,CAAC;oBACb,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;oBAChB,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC;oBACb,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;YACH,CAAC;YAED,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,KAAK,MAAM,GAAG,CAAC,EAAE,CAAC;gBAC3C,QAAQ,EAAE,CAAC;gBACX,MAAM,KAAK,CAAC;YACd,CAAC;YAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;gBAC3B,CAAC,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAClC,QAAQ,EAAE,CAAC;YACX,MAAM,KAAK,CAAC;QACd,CAAC;aAAM,CAAC;YACN,oBAAoB;YACpB,KAAK,EAAE,CAAC;YACR,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC;YAEZ,mBAAmB;YACnB,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;oBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnC,CAAC;YACH,CAAC;YAED,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC;YACZ,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,OAAO,EAAE,CAAC;gBACpC,KAAK,GAAG,CAAC,CAAC,CAAC;gBACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;gBACb,MAAM,GAAG,CAAC,CAAC;gBACX,MAAM,KAAK,CAAC;YACd,CAAC;YAED,IAAI,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC;gBACf,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,MAAM,SAAS,GAAG,WAAW,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;gBACrE,IAAI,SAAS,KAAK,CAAC,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACvC,IAAI,SAAS,KAAK,CAAC;wBACjB,QAAQ,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;oBAC3C,QAAQ,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC1D,MAAM,CAAC,GAAG,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC;gBAC3B,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,UAAU,EAAE,CAAC;oBACb,QAAQ,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC1D,MAAM,CAAC,GAAG,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC;gBAC3B,CAAC;gBACD,SAAS,KAAK,CAAC;YACjB,CAAC;iBAAM,CAAC;gBACN,qBAAqB;gBACrB,IAAI,KAAK,KAAK,CAAC,EAAE,EAAE,CAAC;oBAClB,KAAK,GAAG,CAAC,CAAC,CAAC;oBACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;oBACb,MAAM,GAAG,CAAC,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;qBAAM,CAAC;oBACN,EAAE,GAAG,GAAG,CAAC;oBACT,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;oBACxC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;oBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;oBAC5D,CAAC,CAAC,GAAG,EAAE,CAAC;oBACR,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;oBAC/B,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;oBACjB,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;oBACZ,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC;wBAAE,SAAS,KAAK,CAAC;oBAC/B,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;oBACT,UAAU,EAAE,CAAC;oBACb,SAAS,KAAK,CAAC;gBACjB,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,KAAK,CAAC;AACf,CAAC"}

package/src/solver-tools/lsoda/stoda.ts

@@ -1,286 +0,0 @@
-import {LsodaContext, CCMAX, MSBP} from './common';
-import {vmnorm} from './blas';
-import {cfode} from './cfode';
-import {scaleh} from './scaleh';
-import {correction, CorrectionState} from './correction';
-import {orderswitch, OrderSwitchResult} from './orderswitch';
-import {methodswitch, MethodSwitchResult} from './methodswitch';
-
-/**
- * Perform one step of the integration.
- *
- * jstart:  0 = first step, >0 = continue, -1 = new h/meth/miter, -2 = new h only
- * Returns kflag: 0 = success, -1 = error test failed, -2 = convergence failed.
- */
-export function stoda(ctx: LsodaContext, y: Float64Array, jstart: number): number {
-  const c = ctx.common!;
-  const hmin = ctx.opt!.hmin;
-  const mxords = ctx.opt!.mxords;
-  const mxordn = ctx.opt!.mxordn;
-  const neq = ctx.neq;
-
-  let kflag = 0;
-  const told = c.tn;
-  c.ncf = 0;
-  let delp = 0.;
-
-  let maxord = mxordn;
-  if (c.meth === 2) maxord = mxords;
-
-  // Helper: endstoda — finalize step, scale acor
-  function endstoda(): void {
-    const r = 1. / c.tesco[c.nqu][2];
-    for (let i = 1; i <= neq; i++)
-      c.acor[i] *= r;
-    c.hold = c.h;
-  }
-
-  // Helper: resetcoeff — reload el from elco for current order
-  function resetcoeff(): void {
-    const el0 = c.el[1];
-    for (let i = 1; i <= c.nq + 1; i++)
-      c.el[i] = c.elco[c.nq][i];
-    c.rc = c.rc * c.el[1] / el0;
-  }
-
-  // --- Initialization based on jstart ---
-
-  if (jstart === 0) {
-    c.nq = 1;
-    c.ialth = 2;
-    c.rmax = 10000.;
-    c.rc = 0.;
-    c.crate = 0.7;
-    c.hold = c.h;
-    c.nslp = 0;
-    c.ipup = c.miter;
-    c.el[1] = 1.0;
-    c.icount = 20;
-    c.irflag = 0;
-    c.pdest = 0.;
-    c.pdlast = 0.;
-    cfode(ctx, 1);
-    resetcoeff();
-  }
-
-  if (jstart === -1) {
-    c.ipup = c.miter;
-    if (c.ialth === 1) c.ialth = 2;
-    if (c.meth !== c.mused) {
-      cfode(ctx, c.meth);
-      c.ialth = c.nq + 1;
-      resetcoeff();
-    }
-    if (c.h !== c.hold) {
-      const rh = c.h / c.hold;
-      c.h = c.hold;
-      scaleh(ctx, rh);
-    }
-  }
-
-  if (jstart === -2) {
-    if (c.h !== c.hold) {
-      const rh = c.h / c.hold;
-      c.h = c.hold;
-      scaleh(ctx, rh);
-    }
-  }
-
-  // --- Main integration loop ---
-
-  let dsm = 0.0;
-  const cs: CorrectionState = {del: 0, delp: delp, m: 0};
-  const osResult: OrderSwitchResult = {rh: 0};
-  const msResult: MethodSwitchResult = {rh: 0};
-  let rh: number;
-
-  outer:
-  while (true) {
-    // Before the corrector starts
-    c.jcur = 0;
-
-    inner:
-    while (true) {
-      if (Math.abs(c.rc - 1.) > CCMAX)
-        c.ipup = c.miter;
-      if (c.nst >= c.nslp + MSBP)
-        c.ipup = c.miter;
-
-      c.tn += c.h;
-
-      // Prediction: multiply yh by Pascal triangle matrix
-      for (let j = c.nq; j >= 1; j--) {
-        for (let i1 = j; i1 <= c.nq; i1++) {
-          for (let i = 1; i <= neq; i++)
-            c.yh[i1][i] += c.yh[i1 + 1][i];
-        }
-      }
-
-      const pnorm = vmnorm(neq, c.yh[1], c.ewt);
-
-      cs.del = 0;
-      cs.delp = delp;
-      cs.m = 0;
-      const corflag = correction(ctx, y, pnorm, cs, told);
-      delp = cs.delp;
-
-      if (corflag === 0) break inner;
-
-      if (corflag === 1) {
-        rh = Math.max(0.25, hmin / Math.abs(c.h));
-        scaleh(ctx, rh);
-        continue inner;
-      }
-
-      if (corflag === 2) {
-        kflag = -2;
-        c.hold = c.h;
-        jstart = 1;
-        return kflag;
-      }
-    }
-
-    // Local error test
-    if (cs.m === 0)
-      dsm = cs.del / c.tesco[c.nq][2];
-    if (cs.m > 0)
-      dsm = vmnorm(neq, c.acor, c.ewt) / c.tesco[c.nq][2];
-
-    if (dsm <= 1.) {
-      // Successful step
-      kflag = 0;
-      c.nst++;
-      c.hu = c.h;
-      c.nqu = c.nq;
-      c.mused = c.meth;
-
-      for (let j = 1; j <= c.nq + 1; j++) {
-        const r = c.el[j];
-        for (let i = 1; i <= neq; i++)
-          c.yh[j][i] += r * c.acor[i];
-      }
-
-      c.icount--;
-      if (c.icount < 0) {
-        msResult.rh = 0;
-        methodswitch(ctx, dsm, vmnorm(neq, c.yh[1], c.ewt), msResult);
-        if (c.meth !== c.mused) {
-          rh = Math.max(msResult.rh, hmin / Math.abs(c.h));
-          scaleh(ctx, rh);
-          c.rmax = 10.;
-          endstoda();
-          break outer;
-        }
-      }
-
-      // No method switch — do step/order selection
-      c.ialth--;
-      if (c.ialth === 0) {
-        let rhup = 0.;
-        if (c.nq + 1 !== maxord + 1) {
-          for (let i = 1; i <= neq; i++)
-            c.savf[i] = c.acor[i] - c.yh[maxord + 1][i];
-          const dup = vmnorm(neq, c.savf, c.ewt) / c.tesco[c.nq][3];
-          const exup = 1. / (c.nq + 2);
-          rhup = 1. / (1.4 * Math.pow(dup, exup) + 0.0000014);
-        }
-
-        osResult.rh = 0;
-        const orderflag = orderswitch(ctx, rhup, dsm, osResult, kflag, maxord);
-
-        if (orderflag === 0) {
-          endstoda();
-          break outer;
-        }
-        if (orderflag === 1) {
-          rh = Math.max(osResult.rh, hmin / Math.abs(c.h));
-          scaleh(ctx, rh);
-          c.rmax = 10.;
-          endstoda();
-          break outer;
-        }
-        if (orderflag === 2) {
-          resetcoeff();
-          rh = Math.max(osResult.rh, hmin / Math.abs(c.h));
-          scaleh(ctx, rh);
-          c.rmax = 10.;
-          endstoda();
-          break outer;
-        }
-      }
-
-      if (c.ialth > 1 || c.nq + 1 === maxord + 1) {
-        endstoda();
-        break outer;
-      }
-
-      for (let i = 1; i <= neq; i++)
-        c.yh[maxord + 1][i] = c.acor[i];
-      endstoda();
-      break outer;
-    } else {
-      // Error test failed
-      kflag--;
-      c.tn = told;
-
-      // Retract yh array
-      for (let j = c.nq; j >= 1; j--) {
-        for (let i1 = j; i1 <= c.nq; i1++) {
-          for (let i = 1; i <= neq; i++)
-            c.yh[i1][i] -= c.yh[i1 + 1][i];
-        }
-      }
-
-      c.rmax = 2.;
-      if (Math.abs(c.h) <= hmin * 1.00001) {
-        kflag = -1;
-        c.hold = c.h;
-        jstart = 1;
-        break outer;
-      }
-
-      if (kflag > -3) {
-        osResult.rh = 0;
-        const orderflag = orderswitch(ctx, 0., dsm, osResult, kflag, maxord);
-        if (orderflag === 1 || orderflag === 0) {
-          if (orderflag === 0)
-            osResult.rh = Math.min(osResult.rh, 0.2);
-          osResult.rh = Math.max(osResult.rh, hmin / Math.abs(c.h));
-          scaleh(ctx, osResult.rh);
-        }
-        if (orderflag === 2) {
-          resetcoeff();
-          osResult.rh = Math.max(osResult.rh, hmin / Math.abs(c.h));
-          scaleh(ctx, osResult.rh);
-        }
-        continue outer;
-      } else {
-        // 3 or more failures
-        if (kflag === -10) {
-          kflag = -1;
-          c.hold = c.h;
-          jstart = 1;
-          break outer;
-        } else {
-          rh = 0.1;
-          rh = Math.max(hmin / Math.abs(c.h), rh);
-          c.h *= rh;
-          for (let i = 1; i <= neq; i++)
-            y[i] = c.yh[1][i];
-          ctx.func(c.tn, y.subarray(1), c.savf.subarray(1), ctx.data);
-          c.nfe++;
-          for (let i = 1; i <= neq; i++)
-            c.yh[2][i] = c.h * c.savf[i];
-          c.ipup = c.miter;
-          c.ialth = 5;
-          if (c.nq === 1) continue outer;
-          c.nq = 1;
-          resetcoeff();
-          continue outer;
-        }
-      }
-    }
-  }
-
-  return kflag;
-}

package/src/solver-tools/lsoda-method.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"lsoda-method.js","sourceRoot":"","sources":["lsoda-method.ts"],"names":[],"mappings":"AAAA;;;;;;;;;EASE;AAEF,OAAO,EAAO,SAAS,EAAC,MAAM,eAAe,CAAC;AAE9C,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAE9B,2FAA2F;AAC3F,SAAS,QAAQ,CAAC,CAAe;IAC/B,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE;QACpB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QACd,OAAO,CAAC,CAAC;IACX,CAAC,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EA0CE;AACF,MAAM,UAAU,KAAK,CAAC,IAAU,EAAE,QAAmB;IACnD,gCAAgC;IAChC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;IAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC;IAC3B,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC;IAC3B,MAAM,SAAS,GAAG,IAAI,CAAC,SAAS,CAAC;IACjC,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,CAAC,MAAM,CAAC;IAChC,MAAM,MAAM,GAAG,KAAK,CAAC;IAErB,mDAAmD;IACnD,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IACnD,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IAEnD,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE;QACjD,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM;KAC1C,CAAC,CAAC;IAEH,qFAAqF;IACrF,iFAAiF;IACjF,+DAA+D;IAC/D,IAAI,QAAQ;QACV,QAAQ,CAAC,gBAAgB,EAAE,CAAC;IAE9B,IAAI,CAAC,GAAsB,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC;IAC7C,IAAI,CAAC,GAAG,EAAE,CAAC;IAEX,qFAAqF;IACrF,MAAM,UAAU,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC;IACnD,IAAI,UAAU,GAAG,EAAE,IAAI,UAAU,GAAG,EAAE,EAAE,CAAC;QACvC,MAAM,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QAC1C,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACT,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACT,IAAI,MAAM,CAAC,KAAK,IAAI,CAAC;YACnB,MAAM,IAAI,KAAK,CAAC,SAAS,CAAC,WAAW,CAAC,CAAC;IAC3C,CAAC;IAED,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;IACpD,MAAM,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC;IAClD,MAAM,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,cAAc,CAAC;IAE1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,cAAc,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;QACzD,MAAM,MAAM,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QACxC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QACb,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QAEb,IAAI,MAAM,CAAC,KAAK,IAAI,CAAC;YACnB,MAAM,IAAI,KAAK,CAAC,SAAS,CAAC,WAAW,CAAC,CAAC;IAC3C,CAAC;IAED,gEAAgE;IAChE,MAAM,KAAK,GAAG,MAAM,CAAC,cAAc,EAAE,CAAC;IACtC,MAAM,IAAI,GAAG,KAAK,CAAC,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;IAE7C,IAAI,QAAQ;QACV,QAAQ,CAAC,uBAAuB,EAAE,CAAC;IAErC,gDAAgD;IAChD,MAAM,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC;IAC/B,MAAM,QAAQ,GAAG,KAAK,CAAe,GAAG,GAAG,CAAC,CAAC,CAAC;IAC9C,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;IACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAE9B,8DAA8D;IAC9D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAEvC,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC,QAAQ"}