diff-grok 1.1.0 → 1.3.0

package/src/tests/test-defs.js

@@ -1,207 +0,0 @@
-/* eslint-disable max-len */
-import { mrt, ros3prw, ros34prw, rk4, ab5, ab4, rkdp, rk3, lsoda } from '../../index';
-export const methods = new Map([
-    ['MRT', mrt],
-    ['ROS3PRw', ros3prw],
-    ['ROS34PRw', ros34prw],
-    ['RK4', rk4],
-    ['AB5', ab5],
-    ['AB4', ab4],
-    ['RKDP', rkdp],
-    ['RK3', rk3],
-    ['LSODA', lsoda],
-]);
-/** Implicit methods suitable for stiff performance benchmarks */
-export const implicitMethods = new Map([
-    ['MRT', mrt],
-    ['ROS3PRw', ros3prw],
-    ['ROS34PRw', ros34prw],
-    ['LSODA', lsoda],
-]);
-export const MAX_MAD = 0.1;
-export const TIMEOUT_MS = 10000;
-const BASIC_MODEL = `#name: M|M|2|2
-#description: Modelling the system with two servers & two waiting places (EXPLAINED)
-#equations:  
-  dp0/dt = -la * p0 + mu * p1
-  dp1/dt = la * p0 - (la + mu) * p1 + 2 * mu * p2
-  dp2/dt = la * p1 - (la + 2 * mu) * p2 + 2 * mu * p3
-  dp3/dt = la * p2 - (la + 2 * mu) * p3 + 2 * mu * p4
-
-#expressions:
-  la = 1 / arrival
-  mu = 1 / service
-  p4 = 1 - p0 - p1 - p2 - p3
-  busy = 1 - p0 - p1
-  queue = p3 + 2 * p4
-
-#argument: t
-  _t0 = 0     {min: 0;      max: 10;   caption: start;   category: Time; units: min}  [Initial time of simulation]
-  _t1 = 60    {min: 20;     max: 100;  caption: finish;  category: Time; units: min}  [Final time of simulation]
-  _h  = 1     {min: 0.01;   max: 0.1;  caption: step;    category: Time; units: min}  [Time step of simulation] 
-
-#inits:
-  p0 = 1   {min: 0; max: 1; category: Initial state; caption: empty}   [Probability that initially there are NO customers]
-  p1 = 0   {min: 0; max: 1; category: Initial state; caption: single}  [Probability that initially there is ONE customer]
-  p2 = 0   {min: 0; max: 1; category: Initial state; caption: two}     [Probability that initially there are TWO customers]
-  p3 = 0   {min: 0; max: 1; category: Initial state; caption: three}   [Probability that initially there are THREE customers]
-
-#output:
-  t     {caption: Time, min}
-  p0    {caption: P(Empty)}
-  busy  {caption: P(Full load)}
-  queue {caption: E(Queue)}
-
-#parameters:
-  arrival = 10   {min: 1; max: 100; category: Means; caption: arrival;  units: min} [Mean arrival time]
-  service = 100  {min: 1; max: 100; category: Means; caption: service; units: min} [Mean service time]`;
-const BASIC_MODEL_INPUTS = {
-    _t0: 0,
-    _t1: 60,
-    _h: 1,
-    p0: 1,
-    p1: 0,
-    p2: 0,
-    p3: 0,
-    arrival: 10,
-    service: 100,
-};
-const BASIC_MODEL_OUTPUTS = 4;
-const CYCLIC_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`;
-const CYCLIC_MODEL_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 CYCLIC_MODEL_OUTPUTS = 5;
-const MULTISTAGE_MODEL = `#name: GA-production
-#tags: model
-#description: Gluconic acid (GA) production by Aspergillus niger modeling
-#equations:
-  dX/dt = rX
-  dS/dt = -gamma * rX - lambda * X
-  dO/dt = Kla * (Cod - O) - delta * rX - phi * X
-  dP/dt = alpha * rX + beta * X
-
-#expressions:
-  mu = muM * S / (Ks + S) * O / (Ko + O)
-  rX = mu * X
-
-#argument: t, 1-st stage
-  _t0 = 0 {units: h; caption: initial; category: Misc} [Start of the process]
-  _t1 = 60 {units: h; caption: 1-st stage; category: Durations; min: 20; max: 80} [Duration of the 1-st stage]
-  step = 0.1 {units: h; caption: step; category: Misc; min: 0.01; max: 1} [Time step of simulation]
-
-#update: 2-nd stage
-  duration = overall - _t1
-  S += 70
-
-#inits:  
-  X = 5 {units: kg/m³; caption: biomass; category: Initial concentrations; min: 1; max: 10} [Aspergillus niger biomass]
-  S = 150 {units: kg/m³; caption: glucose; category: Initial concentrations; min: 50; max: 200} [Glucose]
-  O = 7 {units: kg/m³; caption: oxygen; category: Initial concentrations; min: 1; max: 10} [Dissolved oxygen]
-  P = 0 {units: kg/m³; caption: acid; category: Initial concentrations; min: 0; max: 0.1} [Gluconic acid]
-
-#output:
-  t {caption: time}
-  X {caption: biomass}
-  S {caption: glucose}
-  O {caption: oxygen}
-  P {caption: acid}
-
-#parameters:
-  overall = 100 {units: h; category: Durations; min: 100; max: 140} [Overall duration]
-  muM = 0.668 {units: 1/h; category: Parameters} [Monod type model parameter]
-  alpha = 2.92 {category: Parameters} [Monod type model parameter]
-  beta = 0.131 {units: 1/h; category: Parameters} [Monod type model parameter]
-  gamma = 2.12 {category: Parameters} [Monod type model parameter]
-  lambda = 0.232 {units: 1/h; category: Parameters} [Monod type model parameter]
-  delta = 0.278 {category: Parameters} [Monod type model parameter]
-  phi = 4.87e-3 {units: 1/h; category: Parameters} [Monod type model parameter]
-  Ks = 1.309e2 {units: g/L; category: Parameters} [Monod type model parameter]
-  Ko = 3.63e-4 {units: g/L; category: Parameters} [Monod type model parameter]
-  Kla = 1.7e-2 {units: 1/s; category: Parameters} [Volumetric mass transfer coefficient]
-  Cod = 15 {units: kg/m³; category: Parameters} [Liquid phase dissolved oxygen saturation concentration]
-  
-#tolerance: 1e-9`;
-const MULTISTAGE_MODEL_INPUTS = {
-    _t0: 0,
-    _t1: 60,
-    _h: 1,
-    X: 5,
-    S: 150,
-    O: 7,
-    P: 0,
-    overall: 100,
-    muM: 0.668,
-    alpha: 2.92,
-    beta: 0.131,
-    gamma: 2.12,
-    lambda: 0.232,
-    delta: 0.278,
-    phi: 4.87e-3,
-    Ks: 1.309e2,
-    Ko: 3.63e-4,
-    Kla: 1.7e-2,
-    Cod: 15,
-};
-const MULTISTAGE_MODEL_OUTPUTS = 5;
-export const problems = new Map([
-    ['basic', { model: BASIC_MODEL, inputs: BASIC_MODEL_INPUTS, outputsCount: BASIC_MODEL_OUTPUTS }],
-    ['cyclic', { model: CYCLIC_MODEL, inputs: CYCLIC_MODEL_INPUTS, outputsCount: CYCLIC_MODEL_OUTPUTS }],
-    ['multistage', { model: MULTISTAGE_MODEL, inputs: MULTISTAGE_MODEL_INPUTS, outputsCount: MULTISTAGE_MODEL_OUTPUTS }],
-]);
-//# sourceMappingURL=test-defs.js.map

package/src/tests/test-defs.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"test-defs.js","sourceRoot":"","sources":["test-defs.ts"],"names":[],"mappings":"AAAA,4BAA4B;AAC5B,OAAO,EAAC,GAAG,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,KAAK,EAAC,MAAM,aAAa,CAAC;AAEpF,MAAM,CAAC,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC;IAC7B,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,SAAS,EAAE,OAAO,CAAC;IACpB,CAAC,UAAU,EAAE,QAAQ,CAAC;IACtB,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,MAAM,EAAE,IAAI,CAAC;IACd,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,OAAO,EAAE,KAAK,CAAC;CACjB,CAAC,CAAC;AAEH,iEAAiE;AACjE,MAAM,CAAC,MAAM,eAAe,GAAG,IAAI,GAAG,CAAC;IACrC,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,SAAS,EAAE,OAAO,CAAC;IACpB,CAAC,UAAU,EAAE,QAAQ,CAAC;IACtB,CAAC,OAAO,EAAE,KAAK,CAAC;CACjB,CAAC,CAAC;AAEH,MAAM,CAAC,MAAM,OAAO,GAAG,GAAG,CAAC;AAC3B,MAAM,CAAC,MAAM,UAAU,GAAG,KAAK,CAAC;AAEhC,MAAM,WAAW,GAAG;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;uGAkCmF,CAAC;AAExG,MAAM,kBAAkB,GAAG;IACzB,GAAG,EAAE,CAAC;IACN,GAAG,EAAE,EAAE;IACP,EAAE,EAAE,CAAC;IACL,EAAE,EAAE,CAAC;IACL,EAAE,EAAE,CAAC;IACL,EAAE,EAAE,CAAC;IACL,EAAE,EAAE,CAAC;IACL,OAAO,EAAE,EAAE;IACX,OAAO,EAAE,GAAG;CACb,CAAC;AAEF,MAAM,mBAAmB,GAAG,CAAC,CAAC;AAE9B,MAAM,YAAY,GAAG;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;iBAuCJ,CAAC;AAElB,MAAM,mBAAmB,GAAG;IAC1B,MAAM,EAAE,EAAE;IACV,GAAG,EAAE,CAAC;IACN,GAAG,EAAE,EAAE;IACP,EAAE,EAAE,CAAC;IACL,KAAK,EAAE,CAAC;IACR,KAAK,EAAE,CAAC;IACR,IAAI,EAAE,CAAC;IACP,GAAG,EAAE,GAAG;IACR,IAAI,EAAE,KAAK;IACX,EAAE,EAAE,IAAI;IACR,EAAE,EAAE,GAAG;IACP,EAAE,EAAE,IAAI;IACR,CAAC,EAAE,KAAK;IACR,EAAE,EAAE,IAAI;IACR,IAAI,EAAE,IAAI;IACV,GAAG,EAAE,KAAK;IACV,IAAI,EAAE,KAAK;CACZ,CAAC;AAEF,MAAM,oBAAoB,GAAG,CAAC,CAAC;AAE/B,MAAM,gBAAgB,GAAG;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;iBAiDR,CAAC;AAElB,MAAM,uBAAuB,GAAG;IAC9B,GAAG,EAAE,CAAC;IACN,GAAG,EAAE,EAAE;IACP,EAAE,EAAE,CAAC;IACL,CAAC,EAAE,CAAC;IACJ,CAAC,EAAE,GAAG;IACN,CAAC,EAAE,CAAC;IACJ,CAAC,EAAE,CAAC;IACJ,OAAO,EAAE,GAAG;IACZ,GAAG,EAAE,KAAK;IACV,KAAK,EAAE,IAAI;IACX,IAAI,EAAE,KAAK;IACX,KAAK,EAAE,IAAI;IACX,MAAM,EAAE,KAAK;IACb,KAAK,EAAE,KAAK;IACZ,GAAG,EAAE,OAAO;IACZ,EAAE,EAAE,OAAO;IACX,EAAE,EAAE,OAAO;IACX,GAAG,EAAE,MAAM;IACX,GAAG,EAAE,EAAE;CACR,CAAC;AAEF,MAAM,wBAAwB,GAAG,CAAC,CAAC;AAQnC,MAAM,CAAC,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAsB;IACnD,CAAC,OAAO,EAAE,EAAC,KAAK,EAAE,WAAW,EAAE,MAAM,EAAE,kBAAkB,EAAE,YAAY,EAAE,mBAAmB,EAAC,CAAC;IAC9F,CAAC,QAAQ,EAAE,EAAC,KAAK,EAAE,YAAY,EAAE,MAAM,EAAE,mBAAmB,EAAE,YAAY,EAAE,oBAAoB,EAAC,CAAC;IAClG,CAAC,YAAY,EAAE,EAAC,KAAK,EAAE,gBAAgB,EAAE,MAAM,EAAE,uBAAuB,EAAE,YAAY,EAAE,wBAAwB,EAAC,CAAC;CACnH,CAAC,CAAC"}

package/src/tests/test-defs.ts

@@ -1,225 +0,0 @@
-/* eslint-disable max-len */
-import {mrt, ros3prw, ros34prw, rk4, ab5, ab4, rkdp, rk3, lsoda} from '../../index';
-
-export const methods = new Map([
-  ['MRT', mrt],
-  ['ROS3PRw', ros3prw],
-  ['ROS34PRw', ros34prw],
-  ['RK4', rk4],
-  ['AB5', ab5],
-  ['AB4', ab4],
-  ['RKDP', rkdp],
-  ['RK3', rk3],
-  ['LSODA', lsoda],
-]);
-
-/** Implicit methods suitable for stiff performance benchmarks */
-export const implicitMethods = new Map([
-  ['MRT', mrt],
-  ['ROS3PRw', ros3prw],
-  ['ROS34PRw', ros34prw],
-  ['LSODA', lsoda],
-]);
-
-export const MAX_MAD = 0.1;
-export const TIMEOUT_MS = 10000;
-
-const BASIC_MODEL = `#name: M|M|2|2
-#description: Modelling the system with two servers & two waiting places (EXPLAINED)
-#equations:  
-  dp0/dt = -la * p0 + mu * p1
-  dp1/dt = la * p0 - (la + mu) * p1 + 2 * mu * p2
-  dp2/dt = la * p1 - (la + 2 * mu) * p2 + 2 * mu * p3
-  dp3/dt = la * p2 - (la + 2 * mu) * p3 + 2 * mu * p4
-
-#expressions:
-  la = 1 / arrival
-  mu = 1 / service
-  p4 = 1 - p0 - p1 - p2 - p3
-  busy = 1 - p0 - p1
-  queue = p3 + 2 * p4
-
-#argument: t
-  _t0 = 0     {min: 0;      max: 10;   caption: start;   category: Time; units: min}  [Initial time of simulation]
-  _t1 = 60    {min: 20;     max: 100;  caption: finish;  category: Time; units: min}  [Final time of simulation]
-  _h  = 1     {min: 0.01;   max: 0.1;  caption: step;    category: Time; units: min}  [Time step of simulation] 
-
-#inits:
-  p0 = 1   {min: 0; max: 1; category: Initial state; caption: empty}   [Probability that initially there are NO customers]
-  p1 = 0   {min: 0; max: 1; category: Initial state; caption: single}  [Probability that initially there is ONE customer]
-  p2 = 0   {min: 0; max: 1; category: Initial state; caption: two}     [Probability that initially there are TWO customers]
-  p3 = 0   {min: 0; max: 1; category: Initial state; caption: three}   [Probability that initially there are THREE customers]
-
-#output:
-  t     {caption: Time, min}
-  p0    {caption: P(Empty)}
-  busy  {caption: P(Full load)}
-  queue {caption: E(Queue)}
-
-#parameters:
-  arrival = 10   {min: 1; max: 100; category: Means; caption: arrival;  units: min} [Mean arrival time]
-  service = 100  {min: 1; max: 100; category: Means; caption: service; units: min} [Mean service time]`;
-
-const BASIC_MODEL_INPUTS = {
-  _t0: 0,
-  _t1: 60,
-  _h: 1,
-  p0: 1,
-  p1: 0,
-  p2: 0,
-  p3: 0,
-  arrival: 10,
-  service: 100,
-};
-
-const BASIC_MODEL_OUTPUTS = 4;
-
-const CYCLIC_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`;
-
-const CYCLIC_MODEL_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 CYCLIC_MODEL_OUTPUTS = 5;
-
-const MULTISTAGE_MODEL = `#name: GA-production
-#tags: model
-#description: Gluconic acid (GA) production by Aspergillus niger modeling
-#equations:
-  dX/dt = rX
-  dS/dt = -gamma * rX - lambda * X
-  dO/dt = Kla * (Cod - O) - delta * rX - phi * X
-  dP/dt = alpha * rX + beta * X
-
-#expressions:
-  mu = muM * S / (Ks + S) * O / (Ko + O)
-  rX = mu * X
-
-#argument: t, 1-st stage
-  _t0 = 0 {units: h; caption: initial; category: Misc} [Start of the process]
-  _t1 = 60 {units: h; caption: 1-st stage; category: Durations; min: 20; max: 80} [Duration of the 1-st stage]
-  step = 0.1 {units: h; caption: step; category: Misc; min: 0.01; max: 1} [Time step of simulation]
-
-#update: 2-nd stage
-  duration = overall - _t1
-  S += 70
-
-#inits:  
-  X = 5 {units: kg/m³; caption: biomass; category: Initial concentrations; min: 1; max: 10} [Aspergillus niger biomass]
-  S = 150 {units: kg/m³; caption: glucose; category: Initial concentrations; min: 50; max: 200} [Glucose]
-  O = 7 {units: kg/m³; caption: oxygen; category: Initial concentrations; min: 1; max: 10} [Dissolved oxygen]
-  P = 0 {units: kg/m³; caption: acid; category: Initial concentrations; min: 0; max: 0.1} [Gluconic acid]
-
-#output:
-  t {caption: time}
-  X {caption: biomass}
-  S {caption: glucose}
-  O {caption: oxygen}
-  P {caption: acid}
-
-#parameters:
-  overall = 100 {units: h; category: Durations; min: 100; max: 140} [Overall duration]
-  muM = 0.668 {units: 1/h; category: Parameters} [Monod type model parameter]
-  alpha = 2.92 {category: Parameters} [Monod type model parameter]
-  beta = 0.131 {units: 1/h; category: Parameters} [Monod type model parameter]
-  gamma = 2.12 {category: Parameters} [Monod type model parameter]
-  lambda = 0.232 {units: 1/h; category: Parameters} [Monod type model parameter]
-  delta = 0.278 {category: Parameters} [Monod type model parameter]
-  phi = 4.87e-3 {units: 1/h; category: Parameters} [Monod type model parameter]
-  Ks = 1.309e2 {units: g/L; category: Parameters} [Monod type model parameter]
-  Ko = 3.63e-4 {units: g/L; category: Parameters} [Monod type model parameter]
-  Kla = 1.7e-2 {units: 1/s; category: Parameters} [Volumetric mass transfer coefficient]
-  Cod = 15 {units: kg/m³; category: Parameters} [Liquid phase dissolved oxygen saturation concentration]
-  
-#tolerance: 1e-9`;
-
-const MULTISTAGE_MODEL_INPUTS = {
-  _t0: 0,
-  _t1: 60,
-  _h: 1,
-  X: 5,
-  S: 150,
-  O: 7,
-  P: 0,
-  overall: 100,
-  muM: 0.668,
-  alpha: 2.92,
-  beta: 0.131,
-  gamma: 2.12,
-  lambda: 0.232,
-  delta: 0.278,
-  phi: 4.87e-3,
-  Ks: 1.309e2,
-  Ko: 3.63e-4,
-  Kla: 1.7e-2,
-  Cod: 15,
-};
-
-const MULTISTAGE_MODEL_OUTPUTS = 5;
-
-type TestProblem = {
-  model: string,
-  inputs: Record<string, number>,
-  outputsCount: number,
-};
-
-export const problems = new Map<string, TestProblem>([
-  ['basic', {model: BASIC_MODEL, inputs: BASIC_MODEL_INPUTS, outputsCount: BASIC_MODEL_OUTPUTS}],
-  ['cyclic', {model: CYCLIC_MODEL, inputs: CYCLIC_MODEL_INPUTS, outputsCount: CYCLIC_MODEL_OUTPUTS}],
-  ['multistage', {model: MULTISTAGE_MODEL, inputs: MULTISTAGE_MODEL_INPUTS, outputsCount: MULTISTAGE_MODEL_OUTPUTS}],
-]);

package/src/tests/test-utils.js

@@ -1,28 +0,0 @@
-import { getIVP, getIvp2WebWorker, getInputVector, getPipelineCreator, applyPipeline } from '../../index';
-/** Return numerical solution error: maximum absolute deviation between approximate & exact solutions */
-export function getError(method, corProb) {
-    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
-/** Evaluate model  */
-export function evalModel(model, inputs) {
-    const ivp = getIVP(model);
-    const ivpWW = getIvp2WebWorker(ivp);
-    const inputVector = getInputVector(inputs, ivp);
-    const creator = getPipelineCreator(ivp);
-    const pipeline = creator.getPipeline(inputVector);
-    return applyPipeline(pipeline, ivpWW, inputVector);
-} // evalModel
-//# sourceMappingURL=test-utils.js.map

package/src/tests/test-utils.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"test-utils.js","sourceRoot":"","sources":["test-utils.ts"],"names":[],"mappings":"AAAA,OAAO,EAAoB,MAAM,EAAE,gBAAgB,EAAE,cAAc,EAAE,kBAAkB,EACrF,aAAa,EAAC,MAAM,aAAa,CAAC;AAEpC,wGAAwG;AACxG,MAAM,UAAU,QAAQ,CAAC,MAAsC,EAAE,OAAoB;IACnF,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,sBAAsB;AACtB,MAAM,UAAU,SAAS,CAAC,KAAa,EAAE,MAA8B;IACrE,MAAM,GAAG,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;IAC1B,MAAM,KAAK,GAAG,gBAAgB,CAAC,GAAG,CAAC,CAAC;IACpC,MAAM,WAAW,GAAG,cAAc,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;IAChD,MAAM,OAAO,GAAG,kBAAkB,CAAC,GAAG,CAAC,CAAC;IACxC,MAAM,QAAQ,GAAG,OAAO,CAAC,WAAW,CAAC,WAAW,CAAC,CAAC;IAClD,OAAO,aAAa,CAAC,QAAQ,EAAE,KAAK,EAAE,WAAW,CAAC,CAAC;AACrD,CAAC,CAAC,YAAY"}

package/src/tests/test-utils.ts

@@ -1,37 +0,0 @@
-import {ODEs, CorrProblem, getIVP, getIvp2WebWorker, getInputVector, getPipelineCreator,
-  applyPipeline} from '../../index';
-
-/** Return numerical solution error: maximum absolute deviation between approximate & exact solutions */
-export 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
-
-/** Evaluate model  */
-export function evalModel(model: string, inputs: Record<string, number>): Float64Array[] {
-  const ivp = getIVP(model);
-  const ivpWW = getIvp2WebWorker(ivp);
-  const inputVector = getInputVector(inputs, ivp);
-  const creator = getPipelineCreator(ivp);
-  const pipeline = creator.getPipeline(inputVector);
-  return applyPipeline(pipeline, ivpWW, inputVector);
-} // evalModel

package/src/worker-tools/index.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAgB,gBAAgB,EAAC,MAAM,aAAa,CAAC;AAC5D,OAAO,EAAC,QAAQ,EAAC,MAAM,WAAW,CAAC"}

package/src/worker-tools/index.ts

@@ -1,2 +0,0 @@
-export {IVP2WebWorker, getIvp2WebWorker} from './scripting';
-export {solveIvp} from './solving';

package/src/worker-tools/scripting.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"scripting.js","sourceRoot":"","sources":["scripting.ts"],"names":[],"mappings":"AAAA,4BAA4B;AAC5B,kDAAkD;AAElD,OAAO,EAAM,WAAW,EAAE,OAAO,EAAE,UAAU,EAAC,MAAM,oCAAoC,CAAC;AAoBzF,MAAM,cAAc,GAAG,UAAU,CAAC;AAClC,MAAM,mBAAmB,GAAG,KAAK,CAAC;AAElC,4BAA4B;AAC5B,SAAS,aAAa,CAAC,IAAY;IACjC,IAAI,CAAC;QACH,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI;aACzB,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC;aACpB,OAAO,CAAC,QAAQ,EAAE,UAAU,CAAC;aAC7B,OAAO,CAAC,WAAW,EAAE,aAAa,CAAC;aACnC,OAAO,CAAC,eAAe,EAAE,iBAAiB,CAAC;aAC3C,OAAO,CAAC,GAAG,EAAE,GAAG,CAAC,CACnB,CAAC;QAEF,OAAO,IAAI,CAAC;IACd,CAAC;IAAC,OAAO,CAAC,EAAE,CAAC;QACX,OAAO;YACL,MAAM,EAAE,cAAc;YACtB,SAAS,EAAE,mBAAmB;SAC/B,CAAC;IACJ,CAAC;AACH,CAAC;AAED,4DAA4D;AAC5D,SAAS,eAAe,CAAC,GAAQ;IAC/B,MAAM,KAAK,GAAa,EAAE,CAAC;IAE3B,kBAAkB;IAClB,GAAG,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE;QAChC,IAAI,GAAG,KAAK,OAAO;YACjB,KAAK,CAAC,IAAI,CAAC,SAAS,UAAU,CAAC,GAAG,CAAC,kBAAkB,UAAU,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;;YAE5E,KAAK,CAAC,IAAI,CAAC,uCAAuC,CAAC,CAAC;IACxD,CAAC,CAAC,CAAC;IAEH,mBAAmB;IACnB,GAAG,CAAC,cAAc,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,WAAW,CAAC,GAAG,CAAC,WAAW,WAAW,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAEzG,kBAAkB;IAClB,IAAI,GAAG,CAAC,MAAM,KAAK,IAAI;QACrB,GAAG,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,KAAK,EAAE,IAAI,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,IAAI,MAAM,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IAErF,gCAAgC;IAChC,KAAK,CAAC,IAAI,CAAC,SAAS,GAAG,CAAC,GAAG,CAAC,IAAI,kBAAkB,CAAC,CAAC;IACpD,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,IAAI,mBAAmB,GAAG,IAAI,CAAC,CAAC,CAAC;IAEnG,uBAAuB;IACvB,IAAI,GAAG,CAAC,KAAK,KAAK,IAAI;QACpB,GAAG,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,IAAI,EAAE,IAAI,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,IAAI,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;IAE9E,iBAAiB;IACjB,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,gBAAgB,GAAG,OAAO,GAAG,CAAC,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAErH,OAAO,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AAC1B,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA6FG;AACH,MAAM,UAAU,gBAAgB,CAAC,GAAQ;IACvC,MAAM,UAAU,GAAa,EAAE,CAAC;IAChC,IAAI,GAAG,CAAC,MAAM,KAAK,IAAI;QACrB,GAAG,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;IAEzD,MAAM,KAAK,GAAkB;QAC3B,SAAS,EAAE,MAAM,CAAC,GAAG,CAAC,SAAS,CAAC;QAChC,UAAU,EAAE,aAAa,CAAC,GAAG,CAAC,cAAc,CAAC;QAC7C,YAAY,EAAE,eAAe,CAAC,GAAG,CAAC;QAClC,UAAU,EAAE,UAAU;QACtB,SAAS,EAAE,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM;KACzC,CAAC;IAEF,OAAO,KAAK,CAAC;AACf,CAAC;AAED,uFAAuF;AACvF,MAAM,UAAU,OAAO,CAAC,GAAkB,EAAE,SAAuB;IACjE,MAAM,KAAK,GAAa,EAAE,CAAC;IAE3B,GAAG,CAAC,UAAU,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,IAAI,MAAM,SAAS,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAExF,OAAO,KAAK,CAAC,MAAM,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AACnD,CAAC"}

package/src/worker-tools/solving.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"solving.js","sourceRoot":"","sources":["solving.ts"],"names":[],"mappings":"AAAA,wCAAwC;AAExC,OAAO,EAAe,GAAG,EAAE,OAAO,EAAE,QAAQ,EAAQ,WAAW,EAC7D,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,KAAK,EAAC,MAAM,iBAAiB,CAAC;AAE1D,OAAO,EAAgB,OAAO,EAAC,MAAM,aAAa,CAAC;AAEnD,MAAM,QAAQ,GAAG,CAAC,CAAC;AACnB,MAAM,SAAS,GAAG,CAAC,CAAC;AACpB,MAAM,QAAQ,GAAG,CAAC,CAAC;AACnB,MAAM,aAAa,GAAG,CAAC,CAAC;AAExB,oCAAoC;AACpC,SAAS,SAAS,CAAC,IAAwB;IACzC,QAAQ,IAAI,EAAE,CAAC;QACf,KAAK,KAAK;YACR,OAAO,GAAG,CAAC;QAEb,KAAK,SAAS;YACZ,OAAO,OAAO,CAAC;QAEjB,KAAK,KAAK;YACR,OAAO,GAAG,CAAC;QAEb,KAAK,KAAK;YACR,OAAO,GAAG,CAAC;QAEb,KAAK,KAAK;YACR,OAAO,GAAG,CAAC;QAEb,KAAK,MAAM;YACT,OAAO,IAAI,CAAC;QAEd,KAAK,KAAK;YACR,OAAO,GAAG,CAAC;QAEb,KAAK,OAAO;YACV,OAAO,KAAK,CAAC;QAEf;YACE,OAAO,QAAQ,CAAC;IAClB,CAAC;AACH,CAAC;AAED,gDAAgD;AAChD,SAAS,WAAW,CAAC,GAAkB,EAAE,MAAoB;IAC3D,MAAM,QAAQ,GAAG,aAAa,GAAG,GAAG,CAAC,SAAS,GAAG,GAAG,CAAC,UAAU,CAAC,MAAM,CAAC;IAEvE,IAAI,QAAQ,GAAG,MAAM,CAAC,MAAM;QAC1B,MAAM,IAAI,KAAK,CAAC,qCAAqC,QAAQ,cAAc,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC;AAChG,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,QAAQ,CAAC,GAAkB,EAAE,MAAoB;IAC/D,WAAW,CAAC,GAAG,EAAE,MAAM,CAAC,CAAC;IACzB,MAAM,mBAAmB,GAAG,MAAM,CAAC,MAAM,GAAG,GAAG,CAAC,UAAU,CAAC,MAAM,CAAC;IAElE,MAAM,QAAQ,GAAG,OAAO,CAAC,GAAG,EAAE,MAAM,CAAC,KAAK,CAAC,mBAAmB,CAAC,CAAC,CAAC;IACjE,MAAM,IAAI,GAAG,IAAI,QAAQ,CAAC,QAAQ,CAAS,CAAC;IAE5C,MAAM,IAAI,GAAS;QACjB,IAAI,EAAE,EAAE;QACR,GAAG,EAAE;YACH,IAAI,EAAE,EAAE;YACR,KAAK,EAAE,MAAM,CAAC,QAAQ,CAAC;YACvB,MAAM,EAAE,MAAM,CAAC,SAAS,CAAC;YACzB,IAAI,EAAE,MAAM,CAAC,QAAQ,CAAC;SACvB;QACD,OAAO,EAAE,MAAM,CAAC,KAAK,CAAC,QAAQ,GAAG,CAAC,EAAE,mBAAmB,CAAC;QACxD,IAAI,EAAE,IAAI;QACV,SAAS,EAAE,GAAG,CAAC,SAAS;QACxB,gBAAgB,EAAE,EAAE;KACrB,CAAC;IAEF,MAAM,MAAM,GAAG,SAAS,CAAC,GAAG,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;IAChD,MAAM,QAAQ,GAAG,WAAW,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC;IAE7C,OAAO,MAAM,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;AAChC,CAAC"}

package/src/worker-tools/solving.ts

@@ -1,83 +0,0 @@
-// Tools for solving IVPs in web-workers
-
-import {SolverMethod, mrt, ros3prw, ros34prw, ODEs, getCallback,
-  rk4, ab5, ab4, rkdp, rk3, lsoda} from '../solver-tools';
-import {Func} from '../solver-tools/solver-defs';
-import {IVP2WebWorker, getFunc} from './scripting';
-
-const ARG_INIT = 0;
-const ARG_FINAL = 1;
-const ARG_STEP = 2;
-const ARG_INP_COUNT = 3;
-
-/** Return method for solving IVP */
-function getMethod(name: string | undefined): SolverMethod {
-  switch (name) {
-  case 'mrt':
-    return mrt;
-
-  case 'ros3prw':
-    return ros3prw;
-
-  case 'rk4':
-    return rk4;
-
-  case 'ab5':
-    return ab5;
-
-  case 'ab4':
-    return ab4;
-
-  case 'rkdp':
-    return rkdp;
-
-  case 'rk3':
-    return rk3;
-
-  case 'lsoda':
-    return lsoda;
-
-  default:
-    return ros34prw;
-  }
-}
-
-/** Check the correspondence of inputs to IVP */
-function checkInputs(ivp: IVP2WebWorker, inputs: Float64Array): void {
-  const expected = ARG_INP_COUNT + ivp.deqsCount + ivp.paramNames.length;
-
-  if (expected > inputs.length)
-    throw new Error(`Incorrect inputs count, expected: ${expected}, current: ${inputs.length}`);
-}
-
-/** Return solution of IVP with a specified inputs
- * @param ivp - initial value problem to be solved
- * @param inputs - vector of input values
- * @returns solution of the given initial value problem
- */
-export function solveIvp(ivp: IVP2WebWorker, inputs: Float64Array): Float64Array[] {
-  checkInputs(ivp, inputs);
-  const nonParamInputsCount = inputs.length - ivp.paramNames.length;
-
-  const funcCode = getFunc(ivp, inputs.slice(nonParamInputsCount));
-  const func = new Function(funcCode) as Func;
-
-  const odes: ODEs = {
-    name: '',
-    arg: {
-      name: '',
-      start: inputs[ARG_INIT],
-      finish: inputs[ARG_FINAL],
-      step: inputs[ARG_STEP],
-    },
-    initial: inputs.slice(ARG_STEP + 1, nonParamInputsCount),
-    func: func,
-    tolerance: ivp.tolerance,
-    solutionColNames: [],
-  };
-
-  const method = getMethod(ivp.solverOpts.method);
-  const callback = getCallback(ivp.solverOpts);
-
-  return method(odes, callback);
-}