diff-grok 1.1.0 → 1.2.0

package/docs/media/orego.ts

@@ -0,0 +1,25 @@
+/** 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'],
+};
+
+export const oregoReferencePoint = new Float64Array([
+  1.000814870318523,
+  1228.178521549917,
+  132.0554942846706,
+]);

package/docs/media/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/docs/media/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/docs/media/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/docs/media/pollution.js

@@ -0,0 +1,132 @@
+/** 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
+export const pollutionReferencePoint = new Float64Array([
+    0.05646255480022769,
+    0.1342484130422339,
+    0.4139734331099427E-8,
+    0.5523140207484359E-2,
+    0.2018977262302196E-6,
+    0.1464541863493966E-6,
+    0.7784249118997964E-1,
+    0.3245075353396018,
+    0.7494013383880406E-2,
+    0.1622293157301561E-7,
+    0.1135863833257075E-7,
+    0.2230505975721359E-2,
+    0.2087162882798630E-3,
+    0.1396921016840158E-4,
+    0.8964884856898295E-2,
+    0.4352846369330103E-17,
+    0.6899219696263405E-2,
+    0.1007803037365946E-3,
+    0.1772146513969984E-5,
+    0.5682943292316392E-4,
+]);
+//# sourceMappingURL=pollution.js.map

package/docs/media/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;AAEf,MAAM,CAAC,MAAM,uBAAuB,GAAG,IAAI,YAAY,CAAC;IACtD,mBAAmB;IACnB,kBAAkB;IAClB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,kBAAkB;IAClB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,sBAAsB;IACtB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;CACtB,CAAC,CAAC"}

package/docs/media/pollution.ts

@@ -0,0 +1,133 @@
+  /** Kintetic constants for the Pollution model */
+  enum POL {
+    K1 = 0.35,
+    K2 = 26.6,
+    K3 = 1.23e4,
+    K4 = 8.6e-4,
+    K5 = 8.2e-4,
+    K6 = 1.5e4,
+    K7 = 1.3e-4,
+    K8 = 2.4e4,
+    K9 = 1.65e4,
+    K10 = 9e3,
+    K11 = 0.022,
+    K12 = 1.2e4,
+    K13 = 1.88,
+    K14 = 1.63e4,
+    K15 = 4.8e6,
+    K16 = 3.5e-4,
+    K17 = 0.0175,
+    K18 = 1e8,
+    K19 = 4.44e11,
+    K20 = 1240,
+    K21 = 2.1,
+    K22 = 5.78,
+    K23 = 0.0474,
+    K24 = 1780,
+    K25 = 3.12,
+  }; // 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: 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];
+    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
+
+export const pollutionReferencePoint = new Float64Array([
+  0.05646255480022769,
+  0.1342484130422339,
+  0.4139734331099427E-8,
+  0.5523140207484359E-2,
+  0.2018977262302196E-6,
+  0.1464541863493966E-6,
+  0.7784249118997964E-1,
+  0.3245075353396018,
+  0.7494013383880406E-2,
+  0.1622293157301561E-7,
+  0.1135863833257075E-7,
+  0.2230505975721359E-2,
+  0.2087162882798630E-3,
+  0.1396921016840158E-4,
+  0.8964884856898295E-2,
+  0.4352846369330103E-17,
+  0.6899219696263405E-2,
+  0.1007803037365946E-3,
+  0.1772146513969984E-5,
+  0.5682943292316392E-4,
+]);