diff-grok 1.0.8

package/src/examples/hires.ts

@@ -0,0 +1,29 @@
+/** High Irradiance Responses of photomorphogenesis (see https://archimede.uniba.it/~testset/problems/hires.php) */
+export const hires = {
+  name: 'HIRES',
+  arg: {name: 't', start: 0, finish: 321.8122, step: 0.01},
+  initial: [1, 0, 0, 0, 0, 0, 0, 0.0057],
+  func: (t: number, y: Float64Array, output: Float64Array) => {
+    // extract function values
+    const y1 = y[0];
+    const y2 = y[1];
+    const y3 = y[2];
+    const y4 = y[3];
+    const y5 = y[4];
+    const y6 = y[5];
+    const y7 = y[6];
+    const y8 = y[7];
+
+    // compute output
+    output[0] = -1.71 * y1 + 0.43 * y2 + 8.32 * y3 + 0.0007;
+    output[1] = 1.71 * y1 - 8.75 * y2;
+    output[2] = -10.03 * y3 + 0.43 * y4 + 0.035 * y5;
+    output[3] = 8.32 * y2 + 1.71 * y3 - 1.12 * y4;
+    output[4] = -1.745 * y5 + 0.43 * y6 + 0.43 * y7;
+    output[5] = -280 * y6 * y8 + 0.69 * y4 + 1.71 * y5 - 0.43 * y6 + 0.69 * y7;
+    output[6] = 280 * y6 * y8 - 1.81 * y7;
+    output[7] = -280 * y6 * y8 + 1.81 * y7;
+  },
+  tolerance: 1e-10,
+  solutionColNames: ['y1', 'y2', 'y3', 'y4', 'y5', 'y6', 'y7', 'y8'],
+}; // hires

package/src/examples/index.js

@@ -0,0 +1,17 @@
+import { robertson } from './robertson';
+import { hires } from './hires';
+import { vdpol } from './vdpol';
+import { orego } from './orego';
+import { e5 } from './e5';
+import { pollution } from './pollution';
+/** Problems for performance check */
+export const perfProbs = [
+    robertson,
+    hires,
+    vdpol,
+    orego,
+    e5,
+    pollution,
+];
+export { corrProbs } from './corr-probs';
+//# sourceMappingURL=index.js.map

package/src/examples/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,SAAS,EAAC,MAAM,aAAa,CAAC;AACtC,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAC9B,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAC9B,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAC9B,OAAO,EAAC,EAAE,EAAC,MAAM,MAAM,CAAC;AACxB,OAAO,EAAC,SAAS,EAAC,MAAM,aAAa,CAAC;AAEtC,qCAAqC;AACrC,MAAM,CAAC,MAAM,SAAS,GAAG;IACvB,SAAS;IACT,KAAK;IACL,KAAK;IACL,KAAK;IACL,EAAE;IACF,SAAS;CACV,CAAC;AAEF,OAAO,EAAC,SAAS,EAAc,MAAM,cAAc,CAAC"}

package/src/examples/index.ts

@@ -0,0 +1,18 @@
+import {robertson} from './robertson';
+import {hires} from './hires';
+import {vdpol} from './vdpol';
+import {orego} from './orego';
+import {e5} from './e5';
+import {pollution} from './pollution';
+
+/** Problems for performance check */
+export const perfProbs = [
+  robertson,
+  hires,
+  vdpol,
+  orego,
+  e5,
+  pollution,
+];
+
+export {corrProbs, CorrProblem} from './corr-probs';

package/src/examples/model-updates.js

@@ -0,0 +1,112 @@
+/* eslint-disable max-len */
+/** This example shows how to apply pipelines and models with updates.
+    This approach can be used for in-webworkers analysis of models.
+
+    Here, we consider gluconic acid (GA) production by Aspergillus niger modeling.
+*/
+import * as DGL from '../../index';
+/** 1. Model specification */
+const 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`;
+/** 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 = {
+        _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 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!');
+}
+//# sourceMappingURL=model-updates.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"model-updates.js","sourceRoot":"","sources":["model-updates.ts"],"names":[],"mappings":"AAAA,4BAA4B;AAC5B;;;;EAIE;AAEF,OAAO,KAAK,GAAG,MAAM,aAAa,CAAC;AAEnC,6BAA6B;AAC7B,MAAM,KAAK,GAAG;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;iBAiDG,CAAC;AAElB,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,GAAG,EAAE,CAAC;QACN,GAAG,EAAE,EAAE;QACP,EAAE,EAAE,CAAC;QACL,CAAC,EAAE,CAAC;QACJ,CAAC,EAAE,GAAG;QACN,CAAC,EAAE,CAAC;QACJ,CAAC,EAAE,CAAC;QACJ,OAAO,EAAE,GAAG;QACZ,GAAG,EAAE,KAAK;QACV,KAAK,EAAE,IAAI;QACX,IAAI,EAAE,KAAK;QACX,KAAK,EAAE,IAAI;QACX,MAAM,EAAE,KAAK;QACb,KAAK,EAAE,KAAK;QACZ,GAAG,EAAE,OAAO;QACZ,EAAE,EAAE,OAAO;QACX,EAAE,EAAE,OAAO;QACX,GAAG,EAAE,MAAM;QACX,GAAG,EAAE,EAAE;KACR,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,OAAO,CAAC;IACnB,WAAW,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,IAAI,IAAI,GAAG,WAAW,CAAC,CAAC;IAC1D,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/model-updates.ts

@@ -0,0 +1,122 @@
+/* eslint-disable max-len */
+/** This example shows how to apply pipelines and models with updates.
+    This approach can be used for in-webworkers analysis of models.
+
+    Here, we consider gluconic acid (GA) production by Aspergillus niger modeling.
+*/
+
+import * as DGL from '../../index';
+
+/** 1. Model specification */
+const 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`;
+
+/** 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 = {
+    _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 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/orego.js

@@ -0,0 +1,19 @@
+/** The OREGO model (see https://archimede.uniba.it/~testset/report/orego.pdf) */
+export const orego = {
+    name: 'OREGO',
+    arg: { name: 't', start: 0, finish: 360, step: 0.01 },
+    initial: [1, 2, 3],
+    func: (t, y, output) => {
+        // extract function values
+        const y1 = y[0];
+        const y2 = y[1];
+        const y3 = y[2];
+        // compute output
+        output[0] = 77.27 * (y2 - y1 * y2 + y1 - 0.000008375 * y1 * y1);
+        output[1] = 1 / 77.27 * (-y2 - y1 * y2 + y3);
+        output[2] = 0.161 * (y1 - y3);
+    },
+    tolerance: 1e-8,
+    solutionColNames: ['y1', 'y2', 'y3'],
+};
+//# sourceMappingURL=orego.js.map

package/src/examples/orego.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"orego.js","sourceRoot":"","sources":["orego.ts"],"names":[],"mappings":"AAAA,iFAAiF;AACjF,MAAM,CAAC,MAAM,KAAK,GAAG;IACnB,IAAI,EAAE,OAAO;IACb,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,EAAC;IACnD,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IAClB,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;QAEhB,iBAAiB;QACjB,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,WAAW,GAAG,EAAE,GAAG,EAAE,CAAE,CAAC;QACjE,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QAC7C,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;IAChC,CAAC;IACD,SAAS,EAAE,IAAI;IACf,gBAAgB,EAAE,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC;CACrC,CAAC"}

package/src/examples/orego.ts

@@ -0,0 +1,19 @@
+/** The OREGO model (see https://archimede.uniba.it/~testset/report/orego.pdf) */
+export const orego = {
+  name: 'OREGO',
+  arg: {name: 't', start: 0, finish: 360, step: 0.01},
+  initial: [1, 2, 3],
+  func: (t: number, y: Float64Array, output: Float64Array) => {
+    // extract function values
+    const y1 = y[0];
+    const y2 = y[1];
+    const y3 = y[2];
+
+    // compute output
+    output[0] = 77.27 * (y2 - y1 * y2 + y1 - 0.000008375 * y1 * y1 );
+    output[1] = 1 / 77.27 * (-y2 - y1 * y2 + y3);
+    output[2] = 0.161 * (y1 - y3);
+  },
+  tolerance: 1e-8,
+  solutionColNames: ['y1', 'y2', 'y3'],
+};

package/src/examples/pipeline-use.js

@@ -0,0 +1,87 @@
+/* eslint-disable max-len */
+/** This example shows how to apply pipelines and models with customized outputs.
+    This approach can be used for in-webworkers analysis of models.
+
+    Here, we consider modeling queues: https://en.wikipedia.org/wiki/M/M/c_queue
+*/
+import * as DGL from '../../index';
+/** 1. Model specification: the M|M|2|2 model  */
+const 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]`;
+/** 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 = {
+        _t0: 0, // Initial time of simulation
+        _t1: 60, // Final time of simulation
+        _h: 1, // Time step of simulation
+        p0: 1, // Probability that initially there are NO customers
+        p1: 0, // Probability that initially there is ONE customer
+        p2: 0, // Probability that initially there are TWO customers
+        p3: 0, // Probability that initially there are THREE customers
+        arrival: 10, // Mean arrival time
+        service: 100, // Mean service time
+    };
+    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!');
+}
+//# sourceMappingURL=pipeline-use.js.map

package/src/examples/pipeline-use.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"pipeline-use.js","sourceRoot":"","sources":["pipeline-use.ts"],"names":[],"mappings":"AAAA,4BAA4B;AAC5B;;;;EAIE;AAEF,OAAO,KAAK,GAAG,MAAM,aAAa,CAAC;AAEnC,iDAAiD;AACjD,MAAM,KAAK,GAAG;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;uGAkCyF,CAAC;AAExG,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,GAAG,EAAE,CAAC,EAAE,6BAA6B;QACrC,GAAG,EAAE,EAAE,EAAE,2BAA2B;QACpC,EAAE,EAAE,CAAC,EAAE,0BAA0B;QACjC,EAAE,EAAE,CAAC,EAAE,oDAAoD;QAC3D,EAAE,EAAE,CAAC,EAAE,mDAAmD;QAC1D,EAAE,EAAE,CAAC,EAAE,qDAAqD;QAC5D,EAAE,EAAE,CAAC,EAAE,uDAAuD;QAC9D,OAAO,EAAE,EAAE,EAAE,oBAAoB;QACjC,OAAO,EAAE,GAAG,EAAE,oBAAoB;KACnC,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,EAAE,CAAC;IACd,WAAW,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,IAAI,IAAI,GAAG,QAAQ,CAAC,CAAC;IACvD,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/pipeline-use.ts

@@ -0,0 +1,98 @@
+/* eslint-disable max-len */
+/** This example shows how to apply pipelines and models with customized outputs.
+    This approach can be used for in-webworkers analysis of models.
+
+    Here, we consider modeling queues: https://en.wikipedia.org/wiki/M/M/c_queue
+*/
+
+import * as DGL from '../../index';
+
+/** 1. Model specification: the M|M|2|2 model  */
+const 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]`;
+
+/** 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 = {
+    _t0: 0, // Initial time of simulation
+    _t1: 60, // Final time of simulation
+    _h: 1, // Time step of simulation
+    p0: 1, // Probability that initially there are NO customers
+    p1: 0, // Probability that initially there is ONE customer
+    p2: 0, // Probability that initially there are TWO customers
+    p3: 0, // Probability that initially there are THREE customers
+    arrival: 10, // Mean arrival time
+    service: 100, // Mean service time
+  };
+  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/pollution.js

@@ -0,0 +1,110 @@
+/** Kintetic constants for the Pollution model */
+var POL;
+(function (POL) {
+    POL[POL["K1"] = 0.35] = "K1";
+    POL[POL["K2"] = 26.6] = "K2";
+    POL[POL["K3"] = 12300] = "K3";
+    POL[POL["K4"] = 0.00086] = "K4";
+    POL[POL["K5"] = 0.00082] = "K5";
+    POL[POL["K6"] = 15000] = "K6";
+    POL[POL["K7"] = 0.00013] = "K7";
+    POL[POL["K8"] = 24000] = "K8";
+    POL[POL["K9"] = 16500] = "K9";
+    POL[POL["K10"] = 9000] = "K10";
+    POL[POL["K11"] = 0.022] = "K11";
+    POL[POL["K12"] = 12000] = "K12";
+    POL[POL["K13"] = 1.88] = "K13";
+    POL[POL["K14"] = 16300] = "K14";
+    POL[POL["K15"] = 4800000] = "K15";
+    POL[POL["K16"] = 0.00035] = "K16";
+    POL[POL["K17"] = 0.0175] = "K17";
+    POL[POL["K18"] = 100000000] = "K18";
+    POL[POL["K19"] = 444000000000] = "K19";
+    POL[POL["K20"] = 1240] = "K20";
+    POL[POL["K21"] = 2.1] = "K21";
+    POL[POL["K22"] = 5.78] = "K22";
+    POL[POL["K23"] = 0.0474] = "K23";
+    POL[POL["K24"] = 1780] = "K24";
+    POL[POL["K25"] = 3.12] = "K25";
+})(POL || (POL = {}));
+; // POL
+/** The chemical reaction part of the air pollution model (https://archimede.uniba.it/~testset/report/pollu.pdf) */
+export const pollution = {
+    name: 'Pollution',
+    arg: { name: 't', start: 0, finish: 60, step: 0.002 },
+    initial: [0, 0.2, 0, 0.04, 0, 0, 0.1, 0.3, 0.01, 0, 0, 0, 0, 0, 0, 0, 0.007, 0, 0, 0],
+    func: (t, y, output) => {
+        // extract function values
+        const y1 = y[0];
+        const y2 = y[1];
+        const y3 = y[2];
+        const y4 = y[3];
+        const y5 = y[4];
+        const y6 = y[5];
+        const y7 = y[6];
+        const y8 = y[7];
+        const y9 = y[8];
+        const y10 = y[9];
+        const y11 = y[10];
+        const y12 = y[11];
+        const y13 = y[12];
+        const y14 = y[13];
+        const y15 = y[14];
+        const y16 = y[15];
+        const y17 = y[16];
+        const y18 = y[17];
+        const y19 = y[18];
+        const y20 = y[19];
+        // evaluate expressions
+        const r1 = POL.K1 * y1;
+        const r2 = POL.K2 * y2 * y4;
+        const r3 = POL.K3 * y5 * y2;
+        const r4 = POL.K4 * y7;
+        const r5 = POL.K5 * y7;
+        const r6 = POL.K6 * y7 * y6;
+        const r7 = POL.K7 * y9;
+        const r8 = POL.K8 * y9 * y6;
+        const r9 = POL.K9 * y11 * y2;
+        const r10 = POL.K10 * y11 * y1;
+        const r11 = POL.K11 * y13;
+        const r12 = POL.K12 * y10 * y2;
+        const r13 = POL.K13 * y14;
+        const r14 = POL.K14 * y1 * y6;
+        const r15 = POL.K15 * y3;
+        const r16 = POL.K16 * y4;
+        const r17 = POL.K17 * y4;
+        const r18 = POL.K18 * y16;
+        const r19 = POL.K19 * y16;
+        const r20 = POL.K20 * y17 * y6;
+        const r21 = POL.K21 * y19;
+        const r22 = POL.K22 * y19;
+        const r23 = POL.K23 * y1 * y4;
+        const r24 = POL.K24 * y19 * y1;
+        const r25 = POL.K25 * y20;
+        // compute output
+        output[0] = -(r1 + r10 + r14 + r23 + r24) + (r2 + r3 + r9 + r11 + r12 + r22 + r25);
+        output[1] = -r2 - r3 - r9 - r12 + r1 + r21;
+        output[2] = -r15 + r1 + r17 + r19 + r22;
+        output[3] = -r2 - r16 - r17 - r23 + r15;
+        output[4] = -r3 + 2 * r4 + r6 + r7 + r13 + r20;
+        output[5] = -r6 - r8 - r14 - r20 + r3 + 2 * r18;
+        output[6] = -r4 - r5 - r6 + r13;
+        output[7] = r4 + r5 + r6 + r7;
+        output[8] = -r7 - r8;
+        output[9] = -r12 + r7 + r9;
+        output[10] = -r9 - r10 + r8 + r11;
+        output[11] = r9;
+        output[12] = -r11 + r10;
+        output[13] = -r13 + r12;
+        output[14] = r14;
+        output[15] = -r18 - r19 + r16;
+        output[16] = -r20;
+        output[17] = r20;
+        output[18] = -r21 - r22 - r24 + r23 + r25;
+        output[19] = -r25 + r24;
+    },
+    tolerance: 1e-6,
+    solutionColNames: ['y1', 'y2', 'y3', 'y4', 'y5', 'y6', 'y7', 'y8', 'y9', 'y10', 'y11',
+        'y12', 'y13', 'y14', 'y15', 'y16', 'y17', 'y18', 'y19', 'y20'],
+}; // pollution
+//# sourceMappingURL=pollution.js.map

package/src/examples/pollution.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"pollution.js","sourceRoot":"","sources":["pollution.ts"],"names":[],"mappings":"AAAE,iDAAiD;AACjD,IAAK,GA0BJ;AA1BD,WAAK,GAAG;IACN,4BAAS,CAAA;IACT,4BAAS,CAAA;IACT,6BAAW,CAAA;IACX,+BAAW,CAAA;IACX,+BAAW,CAAA;IACX,6BAAU,CAAA;IACV,+BAAW,CAAA;IACX,6BAAU,CAAA;IACV,6BAAW,CAAA;IACX,8BAAS,CAAA;IACT,+BAAW,CAAA;IACX,+BAAW,CAAA;IACX,8BAAU,CAAA;IACV,+BAAY,CAAA;IACZ,iCAAW,CAAA;IACX,iCAAY,CAAA;IACZ,gCAAY,CAAA;IACZ,mCAAS,CAAA;IACT,sCAAa,CAAA;IACb,8BAAU,CAAA;IACV,6BAAS,CAAA;IACT,8BAAU,CAAA;IACV,gCAAY,CAAA;IACZ,8BAAU,CAAA;IACV,8BAAU,CAAA;AACZ,CAAC,EA1BI,GAAG,KAAH,GAAG,QA0BP;AAAA,CAAC,CAAC,MAAM;AAEX,mHAAmH;AACnH,MAAM,CAAC,MAAM,SAAS,GAAG;IACvB,IAAI,EAAE,WAAW;IACjB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAC;IACnD,OAAO,EAAE,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACrF,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;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;QAChB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAClB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAElB,uBAAuB;QACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC;QACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC5B,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC5B,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC;QACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC;QACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC5B,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC;QACvB,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC5B,MAAM,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;QAC7B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QAC1B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QAC1B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC;QACzB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC;QACzB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC;QACzB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QAC1B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QAC1B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QAC1B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QAC1B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QAC9B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QAE1B,iBAAiB;QACjB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC;QACnF,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;QAC3C,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACxC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QACxC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAE,EAAE,GAAE,GAAG,GAAG,GAAG,CAAC;QAC5C,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,CAAC;QAChD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC;QAChC,MAAM,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAC9B,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC;QACrB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAE,EAAE,GAAG,EAAE,CAAC;QAC1B,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;QAClC,MAAM,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;QAChB,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QACxB,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;QACxB,MAAM,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;QACjB,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAC9B,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC;QAClB,MAAM,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;QACjB,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAC1C,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;IAC1B,CAAC;IACD,SAAS,EAAE,IAAI;IACf,gBAAgB,EAAE,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK;QACnF,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC;CACjE,CAAC,CAAC,YAAY"}