diff-grok 1.1.0 → 1.3.0

package/src/examples/bioreactor.js

@@ -1,125 +0,0 @@
-/** Example. Solve the following initial value problem
-
-      dx / dt = x + y - t
-      dy / dt = x * y + t
-      x(0) = 1
-      y(0) = -1
-
-    on [0, 2] with the step 0.01.
- */
-import { lsoda } from '../../index';
-// constants
-const VLinit = 7.2;
-const VTV = 10;
-const speed = 400;
-const diam = 6;
-const power = 2.1;
-const pH = 7.4;
-const k1red = 0.05604;
-const k1ox = 0.0108;
-const k2Fd = 1.35;
-const k2Fa = 110400000;
-const k2Kd = 0.04038;
-const k2Ka = 120000000;
-const k3FKa = 181200000;
-const k3FKd = 0.01188;
-const k4ox = 0.0108;
-const k4red = 0.05604;
-const ktox = 0.005;
-const krcyst = 0;
-const pO2sat = 100;
-const pKa2MEA = 8.18;
-const H = 1.072069378;
-const R = 0.082;
-const _t0 = 0;
-const _t1 = 1000;
-const _h = 1;
-const FFox = 0.2;
-const KKox = 0.2;
-const FFred = 0.1;
-const KKred = 0.1;
-const Ffree = 0;
-const Kfree = 0;
-const FKred = 0;
-const FKox = 0;
-const MEAthiol = 15;
-const CO2 = 0.12;
-const yO2P = 0.209;
-const CYST = 0;
-const VL = 7.2;
-const qin = 1;
-const yO2in = 0.21;
-const T = 300;
-const P = 1;
-const switchTime = 135;
-// the problem definition
-const task = {
-    name: 'Bioreactor',
-    arg: { name: 't', start: _t0, finish: _t1, step: _h },
-    initial: [FFox, KKox, FFred, KKred, Ffree, Kfree, FKred, FKox, MEAthiol, CO2, yO2P, CYST, VL],
-    func: (t, _y, _output) => {
-        // extract function values
-        const FFox = _y[0];
-        const KKox = _y[1];
-        const FFred = _y[2];
-        const KKred = _y[3];
-        const Ffree = _y[4];
-        const Kfree = _y[5];
-        const FKred = _y[6];
-        const FKox = _y[7];
-        const MEAthiol = _y[8];
-        const CO2 = _y[9];
-        const yO2P = _y[10];
-        const CYST = _y[11];
-        const VL = _y[12];
-        // evaluate expressions
-        const KF = pow(VL, -0.65) * 0.065 * pow(Math.pow(speed, 3) * Math.pow(diam, 5) * power / 2.16e12, 0.361);
-        const MA = MEAthiol * pow(10, pH - pKa2MEA);
-        const qout = qin - KF * (yO2P * H - CO2) * VL * R * T / (P * 1000);
-        const OTR = KF * (yO2P * H - CO2);
-        const Fin = t < switchTime ? 0 : 0.025;
-        const Fper = t < switchTime ? 0.025 : Fin;
-        const E0 = VLinit / VL;
-        const E1 = Math.pow((MA * E0), 2);
-        const E11 = k1red * FFox * E0 * E1;
-        const E12 = k1ox * FFred * E0;
-        const E31 = k2Fd * FFred * E0;
-        const E32 = k2Fa * Math.pow((Ffree * E0), 2) * E1;
-        const E21 = k1red * KKox * E0 * E1;
-        const E22 = k1ox * KKred * E0;
-        const E41 = k2Kd * KKred * E0;
-        const E42 = k2Ka * Math.pow((Kfree * E0), 2) * E1;
-        const E51 = k3FKa * Ffree * E0 * Kfree * E0;
-        const E52 = k3FKd * FKred * E0;
-        const E61 = k4ox * FKred * E0 * Math.pow((CYST * E0), 2);
-        const E62 = k4red * FKox * E0 * E1;
-        const E70 = E0 * CO2;
-        const E71 = Math.pow((MEAthiol * E0), 2);
-        const E72 = (E70 >= 0) ? sqrt(E70) : 0;
-        // compute output
-        _output[0] = -E11 + E12;
-        _output[1] = -E21 + E22;
-        _output[2] = E11 - E12 - E31 + E32;
-        _output[3] = E21 - E22 - E41 + E42;
-        _output[4] = 2 * (E31 - E32) - E51 + E52;
-        _output[5] = 2 * (E41 - E42) - E51 + E52;
-        _output[6] = E51 - E52 - E61 + E62;
-        _output[7] = E61 - E62;
-        _output[8] = 2 * (-E11 + E12 - E21 + E22 + E31 + E41 - E32 - E42 - E62 - ktox * E71 * E72) -
-            (MEAthiol + MA) * (Fin + Fper) / VL;
-        _output[9] = (Fin * pO2sat * 0.0022 - 2 * Fper * CO2) / VL + OTR - 0.5 * ktox * E71 * E72;
-        _output[10] = -OTR * (VL / (VTV - VL)) * R * T * P + yO2in * qin - yO2P * qout;
-        _output[11] = ktox * E71 * E72 - krcyst * CYST * E0 - (Fin + Fper) * CYST / VL;
-        _output[12] = Fin - Fper;
-    },
-    tolerance: 0.00005,
-    solutionColNames: ['FFox', 'KKox', 'FFred', 'KKred', 'Ffree', 'Kfree', 'FKred',
-        'FKox', 'MEAthiol', 'CO2', 'yO2P', 'CYST', 'VL'],
-};
-// used Math-functions
-const pow = (x, y) => Math.pow(x, y);
-const tan = (x) => Math.tan(x);
-const sqrt = (x) => Math.sqrt(x);
-const exp = (x) => Math.exp(x);
-const solution = lsoda(task);
-//# sourceMappingURL=bioreactor.js.map

package/src/examples/bioreactor.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"bioreactor.js","sourceRoot":"","sources":["bioreactor.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,OAAO,EAAO,KAAK,EAAM,MAAM,aAAa,CAAC;AAE7C,YAAY;AACZ,MAAM,MAAM,GAAG,GAAG,CAAC;AACnB,MAAM,GAAG,GAAG,EAAE,CAAC;AACf,MAAM,KAAK,GAAG,GAAG,CAAC;AAClB,MAAM,IAAI,GAAG,CAAC,CAAC;AACf,MAAM,KAAK,GAAG,GAAG,CAAC;AAClB,MAAM,EAAE,GAAG,GAAG,CAAC;AACf,MAAM,KAAK,GAAG,OAAO,CAAC;AACtB,MAAM,IAAI,GAAG,MAAM,CAAC;AACpB,MAAM,IAAI,GAAG,IAAI,CAAC;AAClB,MAAM,IAAI,GAAG,SAAS,CAAC;AACvB,MAAM,IAAI,GAAG,OAAO,CAAC;AACrB,MAAM,IAAI,GAAG,SAAS,CAAC;AACvB,MAAM,KAAK,GAAG,SAAS,CAAC;AACxB,MAAM,KAAK,GAAG,OAAO,CAAC;AACtB,MAAM,IAAI,GAAG,MAAM,CAAC;AACpB,MAAM,KAAK,GAAG,OAAO,CAAC;AACtB,MAAM,IAAI,GAAG,KAAK,CAAC;AACnB,MAAM,MAAM,GAAG,CAAC,CAAC;AACjB,MAAM,MAAM,GAAG,GAAG,CAAC;AACnB,MAAM,OAAO,GAAG,IAAI,CAAC;AACrB,MAAM,CAAC,GAAG,WAAW,CAAC;AACtB,MAAM,CAAC,GAAG,KAAK,CAAC;AAEhB,MAAM,GAAG,GAAG,CAAC,CAAC;AACd,MAAM,GAAG,GAAG,IAAI,CAAC;AACjB,MAAM,EAAE,GAAG,CAAC,CAAC;AACb,MAAM,IAAI,GAAG,GAAG,CAAC;AACjB,MAAM,IAAI,GAAG,GAAG,CAAC;AACjB,MAAM,KAAK,GAAG,GAAG,CAAC;AAClB,MAAM,KAAK,GAAG,GAAG,CAAC;AAClB,MAAM,KAAK,GAAG,CAAC,CAAC;AAChB,MAAM,KAAK,GAAG,CAAC,CAAC;AAChB,MAAM,KAAK,GAAG,CAAC,CAAC;AAChB,MAAM,IAAI,GAAG,CAAC,CAAC;AACf,MAAM,QAAQ,GAAG,EAAE,CAAC;AACpB,MAAM,GAAG,GAAG,IAAI,CAAC;AACjB,MAAM,IAAI,GAAG,KAAK,CAAC;AACnB,MAAM,IAAI,GAAG,CAAC,CAAC;AACf,MAAM,EAAE,GAAG,GAAG,CAAC;AACf,MAAM,GAAG,GAAG,CAAC,CAAC;AACd,MAAM,KAAK,GAAG,IAAI,CAAC;AACnB,MAAM,CAAC,GAAG,GAAG,CAAC;AACd,MAAM,CAAC,GAAG,CAAC,CAAC;AACZ,MAAM,UAAU,GAAG,GAAG,CAAC;AAEvB,yBAAyB;AACzB,MAAM,IAAI,GAAG;IACX,IAAI,EAAE,YAAY;IAClB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,EAAC;IACnD,OAAO,EAAE,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,IAAI,EAAE,QAAQ,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,CAAC;IAC7F,IAAI,EAAE,CAAC,CAAS,EAAE,EAAgB,EAAE,OAAqB,EAAE,EAAE;QAC3D,0BAA0B;QAC1B,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACpB,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACpB,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACpB,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACpB,MAAM,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACpB,MAAM,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,QAAQ,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAClB,MAAM,IAAI,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;QACpB,MAAM,IAAI,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;QACpB,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC;QAElB,uBAAuB;QACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,KAAK,GAAG,GAAG,CAAC,SAAA,KAAK,EAAE,CAAC,CAAA,GAAG,SAAA,IAAI,EAAE,CAAC,CAAA,GAAG,KAAK,GAAG,OAAO,EAAE,KAAK,CAAC,CAAC;QACrF,MAAM,EAAE,GAAG,QAAQ,GAAG,GAAG,CAAC,EAAE,EAAE,EAAE,GAAG,OAAO,CAAC,CAAC;QAC5C,MAAM,IAAI,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;QACnE,MAAM,GAAG,GAAG,EAAE,GAAG,CAAC,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;QAClC,MAAM,GAAG,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC;QACvC,MAAM,IAAI,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC;QAC1C,MAAM,EAAE,GAAG,MAAM,GAAG,EAAE,CAAC;QACvB,MAAM,EAAE,GAAG,SAAA,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,CAAA,CAAC;QACxB,MAAM,GAAG,GAAG,KAAK,GAAG,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;QACnC,MAAM,GAAG,GAAG,IAAI,GAAG,KAAK,GAAG,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAG,IAAI,GAAG,KAAK,GAAG,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAG,IAAI,GAAG,SAAA,CAAC,KAAK,GAAG,EAAE,CAAC,EAAE,CAAC,CAAA,GAAG,EAAE,CAAC;QACxC,MAAM,GAAG,GAAG,KAAK,GAAG,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;QACnC,MAAM,GAAG,GAAG,IAAI,GAAG,KAAK,GAAG,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAG,IAAI,GAAG,KAAK,GAAG,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAG,IAAI,GAAG,SAAA,CAAC,KAAK,GAAG,EAAE,CAAC,EAAE,CAAC,CAAA,GAAG,EAAE,CAAC;QACxC,MAAM,GAAG,GAAG,KAAK,GAAG,KAAK,GAAG,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC;QAC5C,MAAM,GAAG,GAAG,KAAK,GAAG,KAAK,GAAG,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,IAAI,GAAG,KAAK,GAAG,EAAE,GAAG,SAAA,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,CAAC,CAAA,CAAC;QAC/C,MAAM,GAAG,GAAG,KAAK,GAAG,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;QACnC,MAAM,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;QACrB,MAAM,GAAG,GAAG,SAAA,CAAC,QAAQ,GAAG,EAAE,CAAC,EAAE,CAAC,CAAA,CAAC;QAC/B,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAEvC,iBAAiB;QACjB,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QACxB,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QACxB,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACnC,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACnC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;QACzC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;QACzC,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACnC,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;QACvB,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC;YACxF,CAAC,QAAQ,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;QACtC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,MAAM,GAAG,MAAM,GAAG,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC;QAC1F,OAAO,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,GAAG,GAAG,IAAI,GAAG,IAAI,CAAC;QAC/E,OAAO,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,MAAM,GAAG,IAAI,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC;QAC/E,OAAO,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC;IAC3B,CAAC;IACD,SAAS,EAAE,OAAO;IAClB,gBAAgB,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO;QAC5E,MAAM,EAAE,UAAU,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,CAAC;CACnD,CAAC;AAEF,sBAAsB;AACtB,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;AACrD,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACvC,MAAM,IAAI,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;AACzC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AAEvC,MAAM,QAAQ,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC"}

package/src/examples/bioreactor.ts

@@ -1,132 +0,0 @@
-/** Example. Solve the following initial value problem
-
-      dx / dt = x + y - t
-      dy / dt = x * y + t
-      x(0) = 1
-      y(0) = -1
-
-    on [0, 2] with the step 0.01.
- */
-
-import {ODEs, lsoda, mrt} from '../../index';
-
-// constants
-const VLinit = 7.2;
-const VTV = 10;
-const speed = 400;
-const diam = 6;
-const power = 2.1;
-const pH = 7.4;
-const k1red = 0.05604;
-const k1ox = 0.0108;
-const k2Fd = 1.35;
-const k2Fa = 110400000;
-const k2Kd = 0.04038;
-const k2Ka = 120000000;
-const k3FKa = 181200000;
-const k3FKd = 0.01188;
-const k4ox = 0.0108;
-const k4red = 0.05604;
-const ktox = 0.005;
-const krcyst = 0;
-const pO2sat = 100;
-const pKa2MEA = 8.18;
-const H = 1.072069378;
-const R = 0.082;
-
-const _t0 = 0;
-const _t1 = 1000;
-const _h = 1;
-const FFox = 0.2;
-const KKox = 0.2;
-const FFred = 0.1;
-const KKred = 0.1;
-const Ffree = 0;
-const Kfree = 0;
-const FKred = 0;
-const FKox = 0;
-const MEAthiol = 15;
-const CO2 = 0.12;
-const yO2P = 0.209;
-const CYST = 0;
-const VL = 7.2;
-const qin = 1;
-const yO2in = 0.21;
-const T = 300;
-const P = 1;
-const switchTime = 135;
-
-// the problem definition
-const task = {
-  name: 'Bioreactor',
-  arg: {name: 't', start: _t0, finish: _t1, step: _h},
-  initial: [FFox, KKox, FFred, KKred, Ffree, Kfree, FKred, FKox, MEAthiol, CO2, yO2P, CYST, VL],
-  func: (t: number, _y: Float64Array, _output: Float64Array) => {
-    // extract function values
-    const FFox = _y[0];
-    const KKox = _y[1];
-    const FFred = _y[2];
-    const KKred = _y[3];
-    const Ffree = _y[4];
-    const Kfree = _y[5];
-    const FKred = _y[6];
-    const FKox = _y[7];
-    const MEAthiol = _y[8];
-    const CO2 = _y[9];
-    const yO2P = _y[10];
-    const CYST = _y[11];
-    const VL = _y[12];
-
-    // evaluate expressions
-    const KF = pow(VL, -0.65) * 0.065 * pow(speed**3 * diam**5 * power / 2.16e12, 0.361);
-    const MA = MEAthiol * pow(10, pH - pKa2MEA);
-    const qout = qin - KF * (yO2P * H - CO2) * VL * R * T / (P * 1000);
-    const OTR = KF * (yO2P * H - CO2);
-    const Fin = t < switchTime ? 0 : 0.025;
-    const Fper = t < switchTime ? 0.025 : Fin;
-    const E0 = VLinit / VL;
-    const E1 = (MA * E0)**2;
-    const E11 = k1red * FFox * E0 * E1;
-    const E12 = k1ox * FFred * E0;
-    const E31 = k2Fd * FFred * E0;
-    const E32 = k2Fa * (Ffree * E0)**2 * E1;
-    const E21 = k1red * KKox * E0 * E1;
-    const E22 = k1ox * KKred * E0;
-    const E41 = k2Kd * KKred * E0;
-    const E42 = k2Ka * (Kfree * E0)**2 * E1;
-    const E51 = k3FKa * Ffree * E0 * Kfree * E0;
-    const E52 = k3FKd * FKred * E0;
-    const E61 = k4ox * FKred * E0 * (CYST * E0)**2;
-    const E62 = k4red * FKox * E0 * E1;
-    const E70 = E0 * CO2;
-    const E71 = (MEAthiol * E0)**2;
-    const E72 = (E70 >= 0) ? sqrt(E70) : 0;
-
-    // compute output
-    _output[0] = -E11 + E12;
-    _output[1] = -E21 + E22;
-    _output[2] = E11 - E12 - E31 + E32;
-    _output[3] = E21 - E22 - E41 + E42;
-    _output[4] = 2 * (E31 - E32) - E51 + E52;
-    _output[5] = 2 * (E41 - E42) - E51 + E52;
-    _output[6] = E51 - E52 - E61 + E62;
-    _output[7] = E61 - E62;
-    _output[8] = 2 * (-E11 + E12 - E21 + E22 + E31 + E41 - E32 - E42 - E62 - ktox * E71 * E72) -
-      (MEAthiol + MA) * (Fin + Fper) / VL;
-    _output[9] = (Fin * pO2sat * 0.0022 - 2 * Fper * CO2) / VL + OTR - 0.5 * ktox * E71 * E72;
-    _output[10] = -OTR * (VL / (VTV - VL)) * R * T * P + yO2in * qin - yO2P * qout;
-    _output[11] = ktox * E71 * E72 - krcyst * CYST * E0 - (Fin + Fper) * CYST / VL;
-    _output[12] = Fin - Fper;
-  },
-  tolerance: 0.00005,
-  solutionColNames: ['FFox', 'KKox', 'FFred', 'KKred', 'Ffree', 'Kfree', 'FKred',
-    'FKox', 'MEAthiol', 'CO2', 'yO2P', 'CYST', 'VL'],
-};
-
-// used Math-functions
-const pow = (x: number, y: number) => Math.pow(x, y);
-const tan = (x: number) => Math.tan(x);
-const sqrt = (x: number) => Math.sqrt(x);
-const exp = (x: number) => Math.exp(x);
-
-const solution = lsoda(task);

package/src/examples/check-methods.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"check-methods.js","sourceRoot":"","sources":["check-methods.ts"],"names":[],"mappings":"AAAA,oCAAoC;AAEpC,OAAO,EAAO,SAAS,EAAe,GAAG,EAAE,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,SAAS,EAAC,MAAM,aAAa,CAAC;AAEvG,wGAAwG;AACxG,SAAS,QAAQ,CAAC,MAAsC,EAAE,OAAoB;IAC5E,IAAI,KAAK,GAAG,CAAC,CAAC;IAEd,yBAAyB;IACzB,MAAM,cAAc,GAAG,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;IAE5C,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC;IAE5B,MAAM,GAAG,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC;IAE9B,MAAM,WAAW,GAAG,GAAG,CAAC,MAAM,CAAC;IAC/B,MAAM,UAAU,GAAG,cAAc,CAAC,MAAM,GAAG,CAAC,CAAC;IAE7C,gBAAgB;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,EAAE,CAAC,EAAE,CAAC;QACrC,MAAM,aAAa,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAEpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,EAAE,CAAC;YACjC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnF,CAAC;IAED,OAAO,KAAK,CAAC;AACf,CAAC,CAAC,WAAW;AAEb,SAAS,8BAA8B,CAAC,QAAwB,EAAE,QAAsB;IACtF,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,MAAM,OAAO,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IACvC,IAAI,GAAG,GAAG,CAAC,CAAC;IAEZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;QACzC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;QACnD,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,QAAQ,CAAC,CAAC;IACrC,CAAC;IAED,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC,iCAAiC;AAEnC,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC;IACtB,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,SAAS,EAAE,OAAO,CAAC;IACpB,CAAC,UAAU,EAAE,QAAQ,CAAC;CACvB,CAAC,CAAC;AAEH,OAAO,CAAC,GAAG,CAAC,eAAe,CAAC,CAAC;AAE7B,OAAO,CAAC,OAAO,CAAC,CAAC,MAAM,EAAE,IAAI,EAAE,EAAE;IAC/B,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IAEvB,OAAO,CAAC,GAAG,CAAC,qCAAqC,CAAC,CAAC;IAEnD,SAAS,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,EAAE;QAC9B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;QACzB,MAAM,QAAQ,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;QAC9B,MAAM,MAAM,GAAG,8BAA8B,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;QAClG,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;QAC1B,OAAO,CAAC,GAAG,CAAC,KAAK,IAAI,CAAC,IAAI,QAAQ,MAAM,GAAG,KAAK,QAAQ,MAAM,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IACjF,CAAC,CAAC,CAAC;IAEH,OAAO,CAAC,GAAG,EAAE,CAAC;AAChB,CAAC,CAAC,CAAC;AAEH,OAAO,CAAC,GAAG,CAAC,mDAAmD,CAAC,CAAC;AAEjE,OAAO,CAAC,GAAG,CAAC,mEAAmE,CAAC,CAAC;AAEjF,OAAO,CAAC,OAAO,CAAC,CAAC,MAAM,EAAE,IAAI,EAAE,EAAE;IAC/B,OAAO,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IACxB,OAAO,CAAC,GAAG,CAAC,yBAAyB,CAAC,CAAC;IAEvC,SAAS,CAAC,OAAO,CAAC,CAAC,OAAO,EAAE,EAAE;QAC5B,MAAM,KAAK,GAAG,QAAQ,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;QACxC,OAAO,CAAC,GAAG,CAAC,QAAQ,OAAO,CAAC,IAAI,CAAC,IAAI,OAAO,KAAK,CAAC,aAAa,CAAC,CAAC,CAAC,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IACrF,CAAC,CAAC,CAAC;IAEH,OAAO,CAAC,GAAG,EAAE,CAAC;AAChB,CAAC,CAAC,CAAC"}

package/src/examples/check-methods.ts

@@ -1,81 +0,0 @@
-// A script for checking performance
-
-import {ODEs, corrProbs, CorrProblem, mrt, ros3prw, ros34prw, perfProbs, refPoints} from '../../index';
-
-/** Return numerical solution error: maximum absolute deviation between approximate & exact solutions */
-function getError(method: (odes: ODEs) => Float64Array[], corProb: CorrProblem): number {
-  let error = 0;
-
-  // Get numerical solution
-  const approxSolution = method(corProb.odes);
-
-  const exact = corProb.exact;
-
-  const arg = approxSolution[0];
-
-  const pointsCount = arg.length;
-  const funcsCount = approxSolution.length - 1;
-
-  // Compute error
-  for (let i = 0; i < pointsCount; ++i) {
-    const exactSolution = exact(arg[i]);
-
-    for (let j = 0; j < funcsCount; ++j)
-      error = Math.max(error, Math.abs(exactSolution[j] - approxSolution[j + 1][i]));
-  }
-
-  return error;
-} // getError
-
-function getDeviationFromReferencePoint(solution: Float64Array[], refPoint: Float64Array): number {
-  let absolute = 0;
-  const lastIdx = solution[0].length - 1;
-  let cur = 0;
-
-  for (let k = 0; k < solution.length; ++k) {
-    cur = Math.abs(solution[k][lastIdx] - refPoint[k]);
-    absolute = Math.max(cur, absolute);
-  }
-
-  return absolute;
-} // getDeviationFromReferencePoint
-
-const methods = new Map([
-  ['MRT', mrt],
-  ['ROS3PRw', ros3prw],
-  ['ROS34PRw', ros34prw],
-]);
-
-console.log('Performance\n');
-
-methods.forEach((method, name) => {
-  console.log(' ', name);
-
-  console.log('           PROBLEM   TIME,MS   ERR*');
-
-  perfProbs.forEach((odes, idx) => {
-    const start = Date.now();
-    const solution = method(odes);
-    const absErr = getDeviationFromReferencePoint(solution.slice(1), refPoints[idx]).toExponential(2);
-    const finish = Date.now();
-    console.log(`  ${odes.name}     ${finish - start}     ${absErr}`.padStart(37));
-  });
-
-  console.log();
-});
-
-console.log('* approximate vs. reference (also, numerical)\n\n');
-
-console.log('Correctness (maximum absolute deviation): approximate vs. exact\n');
-
-methods.forEach((method, name) => {
-  console.log('  ', name);
-  console.log('          PROBLEM   ERR');
-
-  corrProbs.forEach((problem) => {
-    const error = getError(method, problem);
-    console.log(`     ${problem.odes.name}    ${error.toExponential(2)}`.padStart(28));
-  });
-
-  console.log();
-});

package/src/examples/corr-probs.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"corr-probs.js","sourceRoot":"","sources":["corr-probs.ts"],"names":[],"mappings":"AAAA,+DAA+D;AAI/D,4DAA4D;AAC5D,MAAM,UAAU,GAAS;IACvB,IAAI,EAAE,cAAc;IACpB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,EAAC;IACjD,OAAO,EAAE,CAAC,CAAC,CAAC;IACZ,IAAI,EAAE,CAAC,CAAS,EAAE,CAAe,EAAE,MAAoB,EAAE,EAAE;QACzD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,CAAC;IACD,SAAS,EAAE,OAAO;IAClB,gBAAgB,EAAE,CAAC,GAAG,CAAC;CACxB,CAAC;AAEF,kEAAkE;AAClE,MAAM,eAAe,GAAG,CAAC,CAAS,EAAE,EAAE;IACpC,OAAO,IAAI,YAAY,CAAC;QACtB,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;KAC5E,CAAC,CAAC;AACL,CAAC,CAAC;AAEF,0CAA0C;AAC1C,MAAM,UAAU,GAAS;IACvB,IAAI,EAAE,cAAc;IACpB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,EAAC;IACjD,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IACf,IAAI,EAAE,CAAC,CAAS,EAAE,CAAe,EAAE,MAAoB,EAAE,EAAE;QACzD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1B,CAAC;IACD,SAAS,EAAE,WAAW;IACtB,gBAAgB,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC;CAC7B,CAAC;AAEF,gDAAgD;AAChD,MAAM,eAAe,GAAG,CAAC,CAAS,EAAE,EAAE;IACpC,OAAO,IAAI,YAAY,CAAC;QACtB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACzC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;KAC1C,CAAC,CAAC;AACL,CAAC,CAAC;AAEF,2CAA2C;AAC3C,MAAM,UAAU,GAAS;IACvB,IAAI,EAAE,cAAc;IACpB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,KAAK,EAAC;IAClD,OAAO,EAAE,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;IACvB,IAAI,EAAE,CAAC,CAAS,EAAE,CAAe,EAAE,MAAoB,EAAE,EAAE;QACzD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC9C,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/C,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,SAAA,CAAC,EAAE,CAAC,CAAA,GAAG,CAAC,GAAG,SAAA,CAAC,EAAE,CAAC,CAAA,GAAG,CAAC,GAAG,CAAC,CAAC;IAC1C,CAAC;IACD,SAAS,EAAE,UAAU;IACrB,gBAAgB,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC;CAClC,CAAC;AAEF,gDAAgD;AAChD,MAAM,eAAe,GAAG,CAAC,CAAS,EAAE,EAAE;IACpC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACtB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAEtB,OAAO,IAAI,YAAY,CAAC;QACtB,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE;QACpC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE;QAC7C,SAAA,CAAC,EAAE,CAAC,CAAA,GAAG,SAAA,CAAC,EAAE,CAAC,CAAA,GAAG,SAAA,CAAC,EAAE,CAAC,CAAA;KACnB,CAAC,CAAC;AACL,CAAC,CAAC;AAEF,wDAAwD;AACxD,MAAM,OAAO,GAAS;IACpB,IAAI,EAAE,UAAU;IAChB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,EAAC;IACjD,OAAO,EAAE,CAAC,CAAC,CAAC;IACZ,IAAI,EAAE,CAAC,CAAS,EAAE,CAAe,EAAE,MAAoB,EAAE,EAAE;QACzD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACxD,CAAC;IACD,SAAS,EAAE,SAAS;IACpB,gBAAgB,EAAE,CAAC,GAAG,CAAC;CACxB,CAAC;AAEF,8DAA8D;AAC9D,MAAM,YAAY,GAAG,CAAC,CAAS,EAAE,EAAE;IACjC,OAAO,IAAI,YAAY,CAAC;QACtB,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;KACvD,CAAC,CAAC;AACL,CAAC,CAAC;AAEF,yDAAyD;AACzD,MAAM,OAAO,GAAS;IACpB,IAAI,EAAE,UAAU;IAChB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,EAAC;IACjD,OAAO,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC;IACvB,IAAI,EAAE,CAAC,CAAS,EAAE,CAAe,EAAE,MAAoB,EAAE,EAAE;QACzD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACtC,CAAC;IACD,SAAS,EAAE,SAAS;IACpB,gBAAgB,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC;CAC7B,CAAC;AAEF,8DAA8D;AAC9D,MAAM,YAAY,GAAG,CAAC,CAAS,EAAE,EAAE;IACjC,OAAO,IAAI,YAAY,CAAC;QACtB,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC;QAC9D,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC;KAChE,CAAC,CAAC;AACL,CAAC,CAAC;AAEF,uCAAuC;AACvC,MAAM,OAAO,GAAS;IACpB,IAAI,EAAE,UAAU;IAChB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,EAAC;IACjD,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAClB,IAAI,EAAE,CAAC,CAAS,EAAE,CAAe,EAAE,MAAoB,EAAE,EAAE;QACzD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5C,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACpB,CAAC;IACD,SAAS,EAAE,UAAU;IACrB,gBAAgB,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC;CAClC,CAAC;AAEF,4CAA4C;AAC5C,MAAM,YAAY,GAAG,CAAC,CAAS,EAAE,EAAE;IACjC,OAAO,IAAI,YAAY,CAAC;QACtB,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QAChC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QACX,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;KACb,CAAC,CAAC;AACL,CAAC,CAAC;AAUF;;GAEG;AACH,MAAM,CAAC,MAAM,SAAS,GAAG;IACvB,EAAC,IAAI,EAAE,UAAU,EAAE,KAAK,EAAE,eAAe,EAAC;IAC1C,EAAC,IAAI,EAAE,UAAU,EAAE,KAAK,EAAE,eAAe,EAAC;IAC1C,EAAC,IAAI,EAAE,UAAU,EAAE,KAAK,EAAE,eAAe,EAAC;IAC1C,EAAC,IAAI,EAAE,OAAO,EAAE,KAAK,EAAE,YAAY,EAAC;IACpC,EAAC,IAAI,EAAE,OAAO,EAAE,KAAK,EAAE,YAAY,EAAC;IACpC,EAAC,IAAI,EAAE,OAAO,EAAE,KAAK,EAAE,YAAY,EAAC;CACrC,CAAC"}

package/src/examples/corr-probs.ts

@@ -1,155 +0,0 @@
-// Problems with known exact solution (for testing correctness)
-
-import {ODEs} from '../../index';
-
-/** Non-stiff 1D test problem: equation (see [1], p. 736) */
-const nonStiff1D: ODEs = {
-  name: 'non-stiff 1D',
-  arg: {name: 't', start: 0, finish: 4, step: 0.01},
-  initial: [2],
-  func: (t: number, y: Float64Array, output: Float64Array) => {
-    output[0] = 4 * Math.exp(0.8 * t) - 0.5 * y[0];
-  },
-  tolerance: 0.00001,
-  solutionColNames: ['y'],
-};
-
-/** Non-stiff 1D test problem: exact solution (see [1], p. 736) */
-const exactNonStiff1D = (t: number) => {
-  return new Float64Array([
-    (Math.exp(0.8 * t) - Math.exp(-0.5 * t)) * 4 / 1.3 + 2 * Math.exp(-0.5 * t),
-  ]);
-};
-
-/** Non-stiff 2D test problem: equation */
-const nonStiff2D: ODEs = {
-  name: 'non-stiff 2D',
-  arg: {name: 't', start: 0, finish: 4, step: 0.01},
-  initial: [1, 1],
-  func: (t: number, y: Float64Array, output: Float64Array) => {
-    output[0] = y[0] + y[1];
-    output[1] = y[1] - y[0];
-  },
-  tolerance: 0.000000001,
-  solutionColNames: ['x', 'y'],
-};
-
-/** Non-stiff 2D test problem: exact solution */
-const exactNonStiff2D = (t: number) => {
-  return new Float64Array([
-    Math.exp(t) * (Math.cos(t) + Math.sin(t)),
-    Math.exp(t) * (Math.cos(t) - Math.sin(t)),
-  ]);
-};
-
-/** Non-stiff 3D test problem: equations */
-const nonStiff3D: ODEs = {
-  name: 'non-stiff 3D',
-  arg: {name: 't', start: 0, finish: 2, step: 0.001},
-  initial: [0.3, -0.8, 0],
-  func: (t: number, y: Float64Array, output: Float64Array) => {
-    output[0] = 5 * y[0] + 2 * y[1] + Math.sin(t);
-    output[1] = -4 * y[0] - y[1] + Math.exp(2 * t);
-    output[2] = 5 * t**4 - 3 * t**2 + 2 * t;
-  },
-  tolerance: 0.00000001,
-  solutionColNames: ['x', 'y', 'z'],
-};
-
-/** Non-stiff 3D test problem: exact solution */
-const exactNonStiff3D = (t: number) => {
-  const e1 = Math.exp(t);
-  const e2 = Math.exp(2 * t);
-  const e3 = Math.exp(3 * t);
-  const c = Math.cos(t);
-  const s = Math.sin(t);
-
-  return new Float64Array([
-    e1 + e3 + 0.1 * (3 * c - s) - 2 * e2,
-    -2 * e1 - e3 - 0.2 * (2 * s + 4 * c) + 3 * e2,
-    t**5 - t**3 + t**2,
-  ]);
-};
-
-/** Stiff 1D test problem: equation (see [1], p. 767) */
-const stiff1D: ODEs = {
-  name: 'stiff 1D',
-  arg: {name: 't', start: 0, finish: 4, step: 0.01},
-  initial: [0],
-  func: (t: number, y: Float64Array, output: Float64Array) => {
-    output[0] = -1000 * y[0] + 3000 - 2000 * Math.exp(-t);
-  },
-  tolerance: 0.0000005,
-  solutionColNames: ['y'],
-};
-
-/** Stiff 1D test problem: exact solution (see [1], p. 767) */
-const exactStiff1D = (t: number) => {
-  return new Float64Array([
-    3 - 0.998 * Math.exp(-1000 * t) - 2.002 * Math.exp(-t),
-  ]);
-};
-
-/** Stiff 2D test problem: equations (see [1], p. 770) */
-const stiff2D: ODEs = {
-  name: 'stiff 2D',
-  arg: {name: 't', start: 0, finish: 4, step: 0.01},
-  initial: [52.29, 83.82],
-  func: (t: number, y: Float64Array, output: Float64Array) => {
-    output[0] = -5 * y[0] + 3 * y[1];
-    output[1] = 100 * y[0] - 301 * y[1];
-  },
-  tolerance: 0.0000005,
-  solutionColNames: ['x', 'y'],
-};
-
-/** Stiff 2D test problem: exact solution (see [1], p. 770) */
-const exactStiff2D = (t: number) => {
-  return new Float64Array([
-    52.96 * Math.exp(-3.9899 * t) - 0.67 * Math.exp(-302.0101 * t),
-    17.83 * Math.exp(-3.9899 * t) + 65.99 * Math.exp(-302.0101 * t),
-  ]);
-};
-
-/** Stiff 3D test problem: equations */
-const stiff3D: ODEs = {
-  name: 'stiff 3D',
-  arg: {name: 't', start: 0, finish: 4, step: 0.01},
-  initial: [2, 1, 0],
-  func: (t: number, y: Float64Array, output: Float64Array) => {
-    output[0] = -100 * y[0] + 100 * y[1] + y[2];
-    output[1] = y[2];
-    output[2] = -y[1];
-  },
-  tolerance: 0.00000001,
-  solutionColNames: ['x', 'y', 'z'],
-};
-
-/** Stiff 3D test problem: exact solution */
-const exactStiff3D = (t: number) => {
-  return new Float64Array([
-    Math.exp(-100 * t) + Math.cos(t),
-    Math.cos(t),
-    -Math.sin(t),
-  ]);
-};
-
-/** Correctness problem
- * @internal
- */
-export type CorrProblem = {
-    odes: ODEs,
-    exact: (t: number) => Float64Array,
-  };
-
-/** Problems for testing correctness
- * @internal
- */
-export const corrProbs = [
-  {odes: nonStiff1D, exact: exactNonStiff1D},
-  {odes: nonStiff2D, exact: exactNonStiff2D},
-  {odes: nonStiff3D, exact: exactNonStiff3D},
-  {odes: stiff1D, exact: exactStiff1D},
-  {odes: stiff2D, exact: exactStiff2D},
-  {odes: stiff3D, exact: exactStiff3D},
-];

package/src/examples/cyclic-model.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"cyclic-model.js","sourceRoot":"","sources":["cyclic-model.ts"],"names":[],"mappings":"AAAA,4BAA4B;AAC5B;;;;EAIE;AAEF,OAAO,KAAK,GAAG,MAAM,aAAa,CAAC;AAEnC,6BAA6B;AAC7B,MAAM,KAAK,GAAG;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;gDAyCkC,CAAC;AAEjD,oFAAoF;AACpF,MAAM,GAAG,GAAG,GAAG,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC;AAC9B,MAAM,KAAK,GAAG,GAAG,CAAC,gBAAgB,CAAC,GAAG,CAAC,CAAC;AAExC,8BAA8B;AAC9B,IAAI,CAAC;IACH,gCAAgC;IAChC,MAAM,WAAW,GAAG,GAAG,CAAC,cAAc,CAAC,GAAG,CAAC,CAAC;IAC5C,MAAM,OAAO,GAAG,WAAW,CAAC,MAAM,CAAC;IAEnC,wBAAwB;IACxB,MAAM,MAAM,GAAG;QACb,MAAM,EAAE,EAAE;QACV,GAAG,EAAE,CAAC;QACN,GAAG,EAAE,EAAE;QACP,EAAE,EAAE,CAAC;QACL,KAAK,EAAE,CAAC;QACR,KAAK,EAAE,CAAC;QACR,IAAI,EAAE,CAAC;QACP,GAAG,EAAE,GAAG;QACR,IAAI,EAAE,KAAK;QACX,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,GAAG;QACP,EAAE,EAAE,IAAI;QACR,CAAC,EAAE,KAAK;QACR,EAAE,EAAE,IAAI;QACR,IAAI,EAAE,IAAI;QACV,GAAG,EAAE,KAAK;QACV,IAAI,EAAE,KAAK;KACZ,CAAC;IACF,MAAM,WAAW,GAAG,GAAG,CAAC,cAAc,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;IAEpD,yBAAyB;IACzB,MAAM,OAAO,GAAG,GAAG,CAAC,kBAAkB,CAAC,GAAG,CAAC,CAAC;IAC5C,MAAM,QAAQ,GAAG,OAAO,CAAC,WAAW,CAAC,WAAW,CAAC,CAAC;IAElD,8CAA8C;IAC9C,MAAM,QAAQ,GAAG,GAAG,CAAC,aAAa,CAAC,QAAQ,EAAE,KAAK,EAAE,WAAW,CAAC,CAAC;IAEjE,qBAAqB;IAErB,sBAAsB;IACtB,IAAI,IAAI,GAAG,WAAW,CAAC;IACvB,WAAW,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,IAAI,IAAI,GAAG,aAAa,CAAC,CAAC;IAC5D,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAElB,6BAA6B;IAC7B,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;QAChC,IAAI,GAAG,EAAE,CAAC;QAEV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,EAAE,CAAC;YAC9B,IAAI,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;QAE9C,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACpB,CAAC;AACH,CAAC;AAAC,OAAO,GAAG,EAAE,CAAC;IACb,OAAO,CAAC,GAAG,CAAC,qBAAqB,EAAE,GAAG,YAAY,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,kBAAkB,CAAC,CAAC;AAC9F,CAAC"}

package/src/examples/cyclic-model.ts

@@ -1,112 +0,0 @@
-/* eslint-disable max-len */
-/** This example shows how to apply pipelines and cyclic models.
-    This approach can be used for in-webworkers analysis of models.
-
-    Here, we consider pharmacokinetic-pharmacodynamic (PK-PD) simulation: two-compartment model.
-*/
-
-import * as DGL from '../../index';
-
-/** 1. Model specification */
-const model = `#name: PK-PD
-#tags: model
-#description: Pharmacokinetic-pharmacodynamic (PK-PD) simulation: two-compartment model
-#equations:
-  d(depot)/dt = -KA * depot
-  d(centr)/dt = KA * depot - CL * C2 - Q * C2 + Q * C3
-  d(peri)/dt  = Q * C2 - Q * C3
-  d(eff)/dt  = Kin - Kout * (1 - C2/(EC50 + C2)) * eff
-
-#expressions:
-  C2 = centr / V2
-  C3 = peri / V3
-
-#loop:
-  _count = 10 {caption: count; category: Dosing; min: 1; max: 20} [Number of doses]
-  depot += dose
-
-#argument: t
-  _t0 = 0 {units: h; caption: begin; category: Dosing; min: 0; max: 1} [Begin of dosing interval]
-  _t1 = 12 {units: h; caption: end; category: Dosing; min: 5; max: 15} [End of dosing interval]
-  _h = 1 {units: h; caption: step; category: Dosing; min: 0.01; max: 0.1} [Time step of simulation]  
-
-#inits:  
-  depot = 0 {category: Initial values}
-  centr = 0 {category: Initial values} [Central]
-  peri = 0 {category: Initial values} [Peripheral]
-  eff = 0.2 {category: Initial values} [Effective compartment rate]
-
-#parameters:  
-  dose = 1e4 {category: Dosing; min: 1e3; max: 2e4; step: 1e3} [Dosage]
-  KA = 0.3 {caption: rate constant; category: Parameters; min: 0.1; max: 1}
-  CL = 2 {caption: clearance; category: Parameters; min: 1; max: 5}
-  V2 = 4 {caption: central volume; category: Parameters; min: 1; max: 10} [Central compartment volume]
-  Q = 1 {caption: inter rate; category: Parameters; min: 0.1; max: 1} [Intercompartmental rate]
-  V3 = 30 {caption: peri volume; category: Parameters; min: 20; max: 40} [Peripheral compartment volume]
-  EC50 = 8 {caption: effect; category: Parameters; min: 1; max: 10}
-  Kin = 0.2 {caption: Kin; category: Parameters; min: 0.1; max: 0.5} [The first-order production constant]
-  Kout = 0.2 {caption: Kout; category: Parameters; min: 0.1; max: 0.5} [The first-order dissipation rate constant]
-  
-#tolerance: 1e-9
-
-#meta.solver: {method: 'rkdp'; maxTimeMs: 5000}`;
-
-/** 2. Generate IVP-objects: for the main thread & for computations in webworkers */
-const ivp = DGL.getIVP(model);
-const ivpWW = DGL.getIvp2WebWorker(ivp);
-
-/** 3. Perform computations */
-try {
-  // 3.1) Extract names of outputs
-  const outputNames = DGL.getOutputNames(ivp);
-  const outSize = outputNames.length;
-
-  // 3.2) Set model inputs
-  const inputs = {
-    _count: 10,
-    _t0: 0,
-    _t1: 16,
-    _h: 1,
-    depot: 0,
-    centr: 0,
-    peri: 0,
-    eff: 0.2,
-    dose: 10000,
-    KA: 0.63,
-    CL: 3.2,
-    V2: 6.58,
-    Q: 0.622,
-    V3: 35.6,
-    EC50: 5.41,
-    Kin: 0.272,
-    Kout: 0.276,
-  };
-  const inputVector = DGL.getInputVector(inputs, ivp);
-
-  // 3.3) Create a pipeline
-  const creator = DGL.getPipelineCreator(ivp);
-  const pipeline = creator.getPipeline(inputVector);
-
-  // 3.4) Apply pipeline to perform computations
-  const solution = DGL.applyPipeline(pipeline, ivpWW, inputVector);
-
-  // 3.5) Print results
-
-  // 3.5.1) Table header
-  let line = '         ';
-  outputNames.forEach((name) => line += name + '           ');
-  console.log(line);
-
-  // 3.5.2) Table with solution
-  const length = solution[0].length;
-  for (let i = 0; i < length; ++i) {
-    line = '';
-
-    for (let j = 0; j < outSize; ++j)
-      line += solution[j][i].toFixed(8) + '     ';
-
-    console.log(line);
-  }
-} catch (err) {
-  console.log('Simulation failed: ', err instanceof Error ? err.message : 'Unknown problem!');
-}

package/src/examples/e5.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"e5.js","sourceRoot":"","sources":["e5.ts"],"names":[],"mappings":"AAAC,0CAA0C;AAC1C,IAAK,EAKH;AALF,WAAK,EAAE;IACJ,8BAAa,CAAA;IACb,gCAAW,CAAA;IACX,8BAAU,CAAA;IACV,0BAAW,CAAA;AACb,CAAC,EALG,EAAE,KAAF,EAAE,QAKL;AAAA,CAAC;AAEJ,2FAA2F;AAC3F,MAAM,CAAC,MAAM,EAAE,GAAG;IAChB,IAAI,EAAE,IAAI;IACV,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,KAAK,EAAC;IACrD,OAAO,EAAE,CAAC,OAAO,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAC3B,IAAI,EAAE,CAAC,CAAS,EAAE,CAAe,EAAE,MAAoB,EAAE,EAAE;QACzD,0BAA0B;QAC1B,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAEhB,iBAAiB;QACjB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC1C,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QACzC,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC;QACxE,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC;IAC3C,CAAC;IACD,SAAS,EAAE,IAAI;IACf,gBAAgB,EAAE,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC;CAC3C,CAAC;AAEF,MAAM,CAAC,MAAM,gBAAgB,GAAG,IAAI,YAAY,CAAC;IAC/C,uBAAuB;IACvB,sBAAsB;IACtB,sBAAsB;IACtB,CAAC;CACF,CAAC,CAAC"}