diff-grok 1.0.8

package/src/examples/pollution.ts

@@ -0,0 +1,110 @@
+  /** 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

package/src/examples/robertson.js

@@ -0,0 +1,14 @@
+/** Robertson chemical reaction, updated version (see https://archimede.uniba.it/~testset/problems/rober.php) */
+export const robertson = {
+    name: 'Robertson',
+    arg: { name: 't', start: 0, finish: 10e11, step: 2.5e7 },
+    initial: [1, 0, 0],
+    func: (t, y, output) => {
+        output[0] = -0.04 * y[0] + 1e4 * y[1] * y[2];
+        output[1] = 0.04 * y[0] - 1e4 * y[1] * y[2] - 3e7 * Math.pow(y[1], 2);
+        output[2] = 3e7 * Math.pow(y[1], 2);
+    },
+    tolerance: 1e-7,
+    solutionColNames: ['A', 'B', 'C'],
+};
+//# sourceMappingURL=robertson.js.map

package/src/examples/robertson.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"robertson.js","sourceRoot":"","sources":["robertson.ts"],"names":[],"mappings":"AAAA,gHAAgH;AAChH,MAAM,CAAC,MAAM,SAAS,GAAG;IACvB,IAAI,EAAE,WAAW;IACjB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAC;IACtD,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,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7C,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,SAAA,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAA,CAAC;QAC5D,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,SAAA,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAA,CAAC;IAC5B,CAAC;IACD,SAAS,EAAE,IAAI;IACf,gBAAgB,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC;CAClC,CAAC"}

package/src/examples/robertson.ts

@@ -0,0 +1,13 @@
+/** Robertson chemical reaction, updated version (see https://archimede.uniba.it/~testset/problems/rober.php) */
+export const robertson = {
+  name: 'Robertson',
+  arg: {name: 't', start: 0, finish: 10e11, step: 2.5e7},
+  initial: [1, 0, 0],
+  func: (t: number, y: Float64Array, output: Float64Array) => {
+    output[0] = -0.04 * y[0] + 1e4 * y[1] * y[2];
+    output[1] = 0.04 * y[0] - 1e4 * y[1] * y[2] - 3e7 * y[1]**2;
+    output[2] = 3e7 * y[1]**2;
+  },
+  tolerance: 1e-7,
+  solutionColNames: ['A', 'B', 'C'],
+};

package/src/examples/scripting.js

@@ -0,0 +1,31 @@
+import { getIVP, getJScode } from '../../index';
+/** Example. Solve the following initial value problem
+     
+      dx / dt = x + y - cos(t)
+      dy / dt = x - y + sin(t)
+      x(0) = 1
+      y(0) = 0
+
+    on [0, 1] with the step 0.01.
+ */
+/** Diff Studio mode specifying the given problem */
+const model = `#name: Example
+#equations:
+  dx/dt = x + y - cos(t)
+  dy/dt = x - y + sin(t)
+
+#inits:
+  x = 1
+  y = 0
+
+#argument: t
+  start = 0
+  finish = 1
+  step = 0.01`;
+// Get parsed formulas
+const ivp = getIVP(model);
+console.log(ivp);
+// Get lines of JS-code specifying ODEs object
+const lines = getJScode(ivp);
+console.log(lines);
+//# sourceMappingURL=scripting.js.map

package/src/examples/scripting.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"scripting.js","sourceRoot":"","sources":["scripting.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,MAAM,EAAE,SAAS,EAAC,MAAM,aAAa,CAAC;AAE9C;;;;;;;;GAQG;AAEH,oDAAoD;AACpD,MAAM,KAAK,GAAG;;;;;;;;;;;;cAYA,CAAC;AAEf,sBAAsB;AACtB,MAAM,GAAG,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC;AAE1B,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;AAEjB,8CAA8C;AAC9C,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;AAE7B,OAAO,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC"}

package/src/examples/scripting.ts

@@ -0,0 +1,36 @@
+import {getIVP, getJScode} from '../../index';
+
+/** Example. Solve the following initial value problem
+     
+      dx / dt = x + y - cos(t)
+      dy / dt = x - y + sin(t)
+      x(0) = 1
+      y(0) = 0
+
+    on [0, 1] with the step 0.01.
+ */
+
+/** Diff Studio mode specifying the given problem */
+const model = `#name: Example
+#equations:
+  dx/dt = x + y - cos(t)
+  dy/dt = x - y + sin(t)
+
+#inits:
+  x = 1
+  y = 0
+
+#argument: t
+  start = 0
+  finish = 1
+  step = 0.01`;
+
+// Get parsed formulas
+const ivp = getIVP(model);
+
+console.log(ivp);
+
+// Get lines of JS-code specifying ODEs object
+const lines = getJScode(ivp);
+
+console.log(lines);

package/src/examples/vdpol.js

@@ -0,0 +1,13 @@
+/** Van der Pol oscillator (see https://archimede.uniba.it/~testset/problems/vdpol.php) */
+export const vdpol = {
+    name: 'van der Pol',
+    arg: { name: 't', start: 0, finish: 2000, step: 0.1 },
+    initial: [-1, 1],
+    func: (t, y, output) => {
+        output[0] = y[1];
+        output[1] = -y[0] + 1000 * (1 - y[0] * y[0]) * y[1];
+    },
+    tolerance: 1e-12,
+    solutionColNames: ['x1', 'x2'],
+};
+//# sourceMappingURL=vdpol.js.map

package/src/examples/vdpol.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"vdpol.js","sourceRoot":"","sources":["vdpol.ts"],"names":[],"mappings":"AAAA,0FAA0F;AAC1F,MAAM,CAAC,MAAM,KAAK,GAAG;IACnB,IAAI,EAAE,aAAa;IACnB,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAC;IACnD,OAAO,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAChB,IAAI,EAAE,CAAC,CAAS,EAAE,CAAe,EAAE,MAAoB,EAAE,EAAE;QACzD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACtD,CAAC;IACD,SAAS,EAAE,KAAK;IAChB,gBAAgB,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC;CAC/B,CAAC"}

package/src/examples/vdpol.ts

@@ -0,0 +1,12 @@
+/** Van der Pol oscillator (see https://archimede.uniba.it/~testset/problems/vdpol.php) */
+export const vdpol = {
+  name: 'van der Pol',
+  arg: {name: 't', start: 0, finish: 2000, step: 0.1},
+  initial: [-1, 1],
+  func: (t: number, y: Float64Array, output: Float64Array) => {
+    output[0] = y[1];
+    output[1] = -y[0] + 1000 * (1 - y[0] * y[0]) * y[1];
+  },
+  tolerance: 1e-12,
+  solutionColNames: ['x1', 'x2'],
+};

package/src/pipeline/basic-pipeline-creator.js

@@ -0,0 +1,21 @@
+import { getOutputCode } from './output';
+import { PipelineCreator } from './pipeline-creator';
+/** Pipeline creator for model without loops & updates */
+export class BasicModelPipelineCreator extends PipelineCreator {
+    constructor(ivp) {
+        super(ivp);
+    }
+    getPipeline(inputs) {
+        return {
+            wrappers: [
+                {
+                    preproc: null,
+                    out: null,
+                    postproc: null,
+                },
+            ],
+            out: getOutputCode(this.ivp),
+        };
+    }
+} // BasicModelPipelineCreator
+//# sourceMappingURL=basic-pipeline-creator.js.map

package/src/pipeline/basic-pipeline-creator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"basic-pipeline-creator.js","sourceRoot":"","sources":["basic-pipeline-creator.ts"],"names":[],"mappings":"AACA,OAAO,EAAC,aAAa,EAAC,MAAM,UAAU,CAAC;AAEvC,OAAO,EAAC,eAAe,EAAC,MAAM,oBAAoB,CAAC;AAEnD,yDAAyD;AACzD,MAAM,OAAO,yBAA0B,SAAQ,eAAe;IAC5D,YAAY,GAAQ;QAClB,KAAK,CAAC,GAAG,CAAC,CAAC;IACb,CAAC;IAEM,WAAW,CAAC,MAAoB;QACrC,OAAO;YACL,QAAQ,EAAE;gBACR;oBACE,OAAO,EAAE,IAAI;oBACb,GAAG,EAAE,IAAI;oBACT,QAAQ,EAAE,IAAI;iBACf;aACF;YACD,GAAG,EAAE,aAAa,CAAC,IAAI,CAAC,GAAG,CAAC;SAC7B,CAAC;IACJ,CAAC;CACF,CAAC,4BAA4B"}

package/src/pipeline/basic-pipeline-creator.ts

@@ -0,0 +1,24 @@
+import {IVP} from '../scripting-tools';
+import {getOutputCode} from './output';
+import {Pipeline} from './pipeline';
+import {PipelineCreator} from './pipeline-creator';
+
+/** Pipeline creator for model without loops & updates */
+export class BasicModelPipelineCreator extends PipelineCreator {
+  constructor(ivp: IVP) {
+    super(ivp);
+  }
+
+  public getPipeline(inputs: Float64Array): Pipeline {
+    return {
+      wrappers: [
+        {
+          preproc: null,
+          out: null,
+          postproc: null,
+        },
+      ],
+      out: getOutputCode(this.ivp),
+    };
+  }
+} // BasicModelPipelineCreator

package/src/pipeline/constants.js

@@ -0,0 +1,21 @@
+export const ARG_INP_COUNT = 3;
+export const ARG_COL_IDX = 0;
+export const CORRECTION_FACTOR = 0.5;
+export const DEFAULT_CORRECTION = 0.0000001;
+export const LOOP_PARAMS_CONT = 1;
+export const LOOP_MIN_COUNT = 1;
+/** Default argument names*/
+export var ARG;
+(function (ARG) {
+    ARG["START"] = "_t0";
+    ARG["FINISH"] = "_t1";
+    ARG["STEP"] = "_h";
+})(ARG || (ARG = {}));
+;
+export const LOOP_COUNT_NAME = '_count';
+export const argName2IdxMap = new Map([
+    [ARG.START, 0],
+    [ARG.FINISH, 1],
+    [ARG.STEP, 2],
+]);
+//# sourceMappingURL=constants.js.map

package/src/pipeline/constants.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"constants.js","sourceRoot":"","sources":["constants.ts"],"names":[],"mappings":"AAAA,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,CAAC;AAC/B,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAC;AAC7B,MAAM,CAAC,MAAM,iBAAiB,GAAG,GAAG,CAAC;AACrC,MAAM,CAAC,MAAM,kBAAkB,GAAG,SAAS,CAAC;AAC5C,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAAC,CAAC;AAClC,MAAM,CAAC,MAAM,cAAc,GAAG,CAAC,CAAC;AAEhC,4BAA4B;AAC5B,MAAM,CAAN,IAAY,GAIX;AAJD,WAAY,GAAG;IACb,oBAAa,CAAA;IACb,qBAAc,CAAA;IACd,kBAAW,CAAA;AACb,CAAC,EAJW,GAAG,KAAH,GAAG,QAId;AAAA,CAAC;AAEF,MAAM,CAAC,MAAM,eAAe,GAAG,QAAQ,CAAC;AAExC,MAAM,CAAC,MAAM,cAAc,GAAG,IAAI,GAAG,CAAiB;IACpD,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC;IACd,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC;IACf,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC;CACd,CAAC,CAAC"}

package/src/pipeline/constants.ts

@@ -0,0 +1,21 @@
+export const ARG_INP_COUNT = 3;
+export const ARG_COL_IDX = 0;
+export const CORRECTION_FACTOR = 0.5;
+export const DEFAULT_CORRECTION = 0.0000001;
+export const LOOP_PARAMS_CONT = 1;
+export const LOOP_MIN_COUNT = 1;
+
+/** Default argument names*/
+export enum ARG {
+  START = '_t0',
+  FINISH = '_t1',
+  STEP = '_h',
+};
+
+export const LOOP_COUNT_NAME = '_count';
+
+export const argName2IdxMap = new Map<string, number>([
+  [ARG.START, 0],
+  [ARG.FINISH, 1],
+  [ARG.STEP, 2],
+]);

package/src/pipeline/index.js

@@ -0,0 +1,8 @@
+export { applyPipeline, getPipelineCreator } from './pipeline';
+export { getOutputCode, getOutputNames } from './output';
+export { PipelineCreator } from './pipeline-creator';
+export { BasicModelPipelineCreator } from './basic-pipeline-creator';
+export { CyclicModelPipelineCreator } from './loops-pipeline-creator';
+export { UpdatesModelPipelineCreator } from './updates-pipeline-creator';
+export { getInputVector } from './input';
+//# sourceMappingURL=index.js.map

package/src/pipeline/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAoB,aAAa,EAAE,kBAAkB,EAAC,MAAM,YAAY,CAAC;AAChF,OAAO,EAAC,aAAa,EAAE,cAAc,EAAC,MAAM,UAAU,CAAC;AACvD,OAAO,EAAC,eAAe,EAAC,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAC,yBAAyB,EAAC,MAAM,0BAA0B,CAAC;AACnE,OAAO,EAAC,0BAA0B,EAAC,MAAM,0BAA0B,CAAC;AACpE,OAAO,EAAC,2BAA2B,EAAC,MAAM,4BAA4B,CAAC;AACvE,OAAO,EAAC,cAAc,EAAC,MAAM,SAAS,CAAC"}

package/src/pipeline/index.ts

@@ -0,0 +1,7 @@
+export {Pipeline, Wrapper, applyPipeline, getPipelineCreator} from './pipeline';
+export {getOutputCode, getOutputNames} from './output';
+export {PipelineCreator} from './pipeline-creator';
+export {BasicModelPipelineCreator} from './basic-pipeline-creator';
+export {CyclicModelPipelineCreator} from './loops-pipeline-creator';
+export {UpdatesModelPipelineCreator} from './updates-pipeline-creator';
+export {getInputVector} from './input';

package/src/pipeline/input.js

@@ -0,0 +1,45 @@
+import { ARG_INP_COUNT, argName2IdxMap, LOOP_PARAMS_CONT, LOOP_COUNT_NAME } from './constants';
+/** Return inputs in a form of a vector */
+export function getInputVector(inputs, ivp) {
+    const eqsCount = ivp.deqs.solutionNames.length;
+    const paramsCount = (ivp.params !== null) ? ivp.params.size : 0;
+    const loopParamsCount = (ivp.loop !== null) ? LOOP_PARAMS_CONT : 0;
+    const size = ARG_INP_COUNT + eqsCount + paramsCount + loopParamsCount;
+    const inputVector = new Float64Array(size);
+    // Argument
+    argName2IdxMap.forEach((index, key) => {
+        if (key in inputs)
+            inputVector[index] = inputs[key];
+        else
+            throw new Error(`Inconsistent inputs: "${key}" is missing`);
+    });
+    let idx = ARG_INP_COUNT;
+    // Init values
+    ivp.deqs.solutionNames.forEach((name) => {
+        if (name in inputs)
+            inputVector[idx] = inputs[name];
+        else
+            throw new Error(`Inconsistent inputs: "${name}" is missing`);
+        ++idx;
+    });
+    // Parameters
+    if (ivp.params !== null) {
+        ivp.params.forEach((_, name) => {
+            if (name in inputs)
+                inputVector[idx] = inputs[name];
+            else
+                throw new Error(`Inconsistent inputs: "${name}" is missing`);
+            ++idx;
+        });
+    }
+    // Loop params
+    if (ivp.loop !== null) {
+        if (LOOP_COUNT_NAME in inputs)
+            inputVector[idx] = inputs[LOOP_COUNT_NAME];
+        else
+            throw new Error(`Inconsistent inputs, loop repetitions count is not defined: "${LOOP_COUNT_NAME}" is missing`);
+        ++idx;
+    }
+    return inputVector;
+} // getInputVector
+//# sourceMappingURL=input.js.map

package/src/pipeline/input.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"input.js","sourceRoot":"","sources":["input.ts"],"names":[],"mappings":"AACA,OAAO,EAAC,aAAa,EAAE,cAAc,EAAE,gBAAgB,EAAE,eAAe,EAAC,MAAM,aAAa,CAAC;AAE7F,0CAA0C;AAC1C,MAAM,UAAU,cAAc,CAAC,MAA8B,EAAE,GAAQ;IACrE,MAAM,QAAQ,GAAG,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC;IAC/C,MAAM,WAAW,GAAG,CAAC,GAAG,CAAC,MAAM,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAChE,MAAM,eAAe,GAAG,CAAC,GAAG,CAAC,IAAI,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,gBAAgB,CAAC,CAAC,CAAC,CAAC,CAAC;IACnE,MAAM,IAAI,GAAG,aAAa,GAAG,QAAQ,GAAG,WAAW,GAAG,eAAe,CAAC;IACtE,MAAM,WAAW,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,CAAC;IAE3C,WAAW;IACX,cAAc,CAAC,OAAO,CAAC,CAAC,KAAK,EAAE,GAAG,EAAE,EAAE;QACpC,IAAI,GAAG,IAAI,MAAM;YACf,WAAW,CAAC,KAAK,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;;YAEjC,MAAM,IAAI,KAAK,CAAC,yBAAyB,GAAG,cAAc,CAAC,CAAC;IAChE,CAAC,CAAC,CAAC;IAEH,IAAI,GAAG,GAAG,aAAa,CAAC;IAExB,cAAc;IACd,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE;QACtC,IAAI,IAAI,IAAI,MAAM;YAChB,WAAW,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;;YAEhC,MAAM,IAAI,KAAK,CAAC,yBAAyB,IAAI,cAAc,CAAC,CAAC;QAE/D,EAAE,GAAG,CAAC;IACR,CAAC,CAAC,CAAC;IAEH,aAAa;IACb,IAAI,GAAG,CAAC,MAAM,KAAK,IAAI,EAAE,CAAC;QACxB,GAAG,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE;YAC7B,IAAI,IAAI,IAAI,MAAM;gBAChB,WAAW,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;;gBAEhC,MAAM,IAAI,KAAK,CAAC,yBAAyB,IAAI,cAAc,CAAC,CAAC;YAE/D,EAAE,GAAG,CAAC;QACR,CAAC,CAAC,CAAC;IACL,CAAC;IAED,cAAc;IACd,IAAI,GAAG,CAAC,IAAI,KAAK,IAAI,EAAE,CAAC;QACtB,IAAI,eAAe,IAAI,MAAM;YAC3B,WAAW,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,eAAe,CAAC,CAAC;;YAE3C,MAAM,IAAI,KAAK,CAAC,gEAAgE,eAAe,cAAc,CAAC,CAAC;QAEjH,EAAE,GAAG,CAAC;IACR,CAAC;IAED,OAAO,WAAW,CAAC;AACrB,CAAC,CAAC,iBAAiB"}

package/src/pipeline/input.ts

@@ -0,0 +1,55 @@
+import {IVP} from '../scripting-tools';
+import {ARG_INP_COUNT, argName2IdxMap, LOOP_PARAMS_CONT, LOOP_COUNT_NAME} from './constants';
+
+/** Return inputs in a form of a vector */
+export function getInputVector(inputs: Record<string, number>, ivp: IVP): Float64Array {
+  const eqsCount = ivp.deqs.solutionNames.length;
+  const paramsCount = (ivp.params !== null) ? ivp.params.size : 0;
+  const loopParamsCount = (ivp.loop !== null) ? LOOP_PARAMS_CONT : 0;
+  const size = ARG_INP_COUNT + eqsCount + paramsCount + loopParamsCount;
+  const inputVector = new Float64Array(size);
+
+  // Argument
+  argName2IdxMap.forEach((index, key) => {
+    if (key in inputs)
+      inputVector[index] = inputs[key];
+    else
+      throw new Error(`Inconsistent inputs: "${key}" is missing`);
+  });
+
+  let idx = ARG_INP_COUNT;
+
+  // Init values
+  ivp.deqs.solutionNames.forEach((name) => {
+    if (name in inputs)
+      inputVector[idx] = inputs[name];
+    else
+      throw new Error(`Inconsistent inputs: "${name}" is missing`);
+
+    ++idx;
+  });
+
+  // Parameters
+  if (ivp.params !== null) {
+    ivp.params.forEach((_, name) => {
+      if (name in inputs)
+        inputVector[idx] = inputs[name];
+      else
+        throw new Error(`Inconsistent inputs: "${name}" is missing`);
+
+      ++idx;
+    });
+  }
+
+  // Loop params
+  if (ivp.loop !== null) {
+    if (LOOP_COUNT_NAME in inputs)
+      inputVector[idx] = inputs[LOOP_COUNT_NAME];
+    else
+      throw new Error(`Inconsistent inputs, loop repetitions count is not defined: "${LOOP_COUNT_NAME}" is missing`);
+
+    ++idx;
+  }
+
+  return inputVector;
+} // getInputVector

package/src/pipeline/loops-pipeline-creator.js

@@ -0,0 +1,107 @@
+import { MATH_CONSTS, MATH_FUNCS, POW_IDX } from '../scripting-tools/scripting-tools';
+import { getOutputCode } from './output';
+import { PipelineCreator } from './pipeline-creator';
+import { argName2IdxMap, ARG, ARG_COL_IDX, CORRECTION_FACTOR, DEFAULT_CORRECTION, LOOP_MIN_COUNT } from './constants';
+/** Pipeline creator for cyclic model */
+export class CyclicModelPipelineCreator extends PipelineCreator {
+    constructor(ivp) {
+        super(ivp);
+        if (ivp.loop === null)
+            throw new Error('Incorrect use of CyclicModelPipelineCreator: model has no loop');
+    }
+    /** Return computations pipeline */
+    getPipeline(inputs) {
+        if (this.ivp.loop === null)
+            throw new Error('Incorrect use of CyclicModelPipelineCreator: model has no loop');
+        const argSize = argName2IdxMap.size;
+        const restricted = inputs.length - 1;
+        const repetitions = inputs[restricted];
+        if (repetitions < LOOP_MIN_COUNT)
+            throw new Error(`Incorrect repetitions count: ${repetitions}, expected: > 0`);
+        const preProcLines = [];
+        // Used math funcs
+        this.ivp.usedMathFuncs.forEach((idx) => {
+            if (idx !== POW_IDX)
+                preProcLines.push(`const ${MATH_FUNCS[idx]} = (x) => Math.${MATH_FUNCS[idx]}(x);`);
+            else
+                preProcLines.push(`const pow = (x, y) => Math.pow(x, y);`);
+        });
+        // Used math consts
+        this.ivp.usedMathConsts.forEach((idx) => {
+            preProcLines.push(`const ${MATH_CONSTS[idx]} = Math.${MATH_CONSTS[idx]};`);
+        });
+        // Model constants
+        if (this.ivp.consts !== null)
+            this.ivp.consts.forEach((input, name) => preProcLines.push(`const ${name} = ${input.value};`));
+        preProcLines.push(`const inputs = arguments[1].slice(0, ${restricted});`);
+        // Extract inputs
+        argName2IdxMap.forEach((idx, name) => preProcLines.push(`let ${name} = inputs[${idx}];`));
+        this.ivp.deqs.solutionNames.forEach((name, idx) => {
+            preProcLines.push(`let ${name} = inputs[${idx + argSize}];`);
+        });
+        if (this.ivp.params !== null) {
+            const shift = argSize + this.ivp.deqs.solutionNames.length;
+            let idx = 0;
+            this.ivp.params.forEach((_, name) => {
+                preProcLines.push(`let ${name} = inputs[${idx + shift}];`);
+                ++idx;
+            });
+        }
+        // Updates
+        this.ivp.loop.updates.forEach((expr) => preProcLines.push(`${expr};`));
+        // Input vector update
+        argName2IdxMap.forEach((idx, name) => preProcLines.push(`inputs[${idx}] = ${name};`));
+        this.ivp.deqs.solutionNames.forEach((name, idx) => {
+            preProcLines.push(`inputs[${idx + argSize}] = ${name};`);
+        });
+        if (this.ivp.params !== null) {
+            const shift = argSize + this.ivp.deqs.solutionNames.length;
+            let idx = 0;
+            this.ivp.params.forEach((_, name) => {
+                preProcLines.push(`inputs[${idx + shift}] = ${name};`);
+                ++idx;
+            });
+        }
+        preProcLines.push(`return inputs;`);
+        const postProcLines = [
+            'const solution = arguments[0];',
+            'const lastIdx = solution[0].length - 1;',
+            'const inputs = new Float64Array(arguments[1]);',
+            `const duration = inputs[${argName2IdxMap.get(ARG.FINISH)}] - inputs[${argName2IdxMap.get(ARG.START)}];`,
+            `inputs[${argName2IdxMap.get(ARG.START)}] = inputs[${argName2IdxMap.get(ARG.FINISH)}];`,
+            `inputs[${argName2IdxMap.get(ARG.FINISH)}] += duration;`,
+        ];
+        this.ivp.deqs.solutionNames.forEach((_, idx) => {
+            postProcLines.push(`inputs[${idx + argSize}] = solution[${idx + 1}][lastIdx];`);
+        });
+        postProcLines.push('return inputs;');
+        const preCode = preProcLines.join('\n');
+        const postCode = postProcLines.join('\n');
+        const wrappers = new Array(repetitions);
+        for (let i = 0; i < repetitions; ++i) {
+            wrappers[i] = {
+                preproc: preCode,
+                out: null,
+                postproc: postCode,
+            };
+        }
+        wrappers[repetitions - 1].postproc = this.getArgumentCorrectionCode();
+        return {
+            wrappers: wrappers,
+            out: getOutputCode(this.ivp),
+        };
+    } // getPipeline
+    /** Return code for correction of an argument */
+    getArgumentCorrectionCode() {
+        return `const arg = arguments[0][${ARG_COL_IDX}];
+    const length = arg.length;
+    const step = arg[1] - arg[0];
+    const correction = Math.min(step * ${CORRECTION_FACTOR}, ${DEFAULT_CORRECTION});
+    for (let i = 0; i < length - 1; ++i) {
+      if (arg[i + 1] - arg[i] < step)
+        arg[i] -= correction;
+    }
+    return arguments[0];`;
+    } // getArgumentCorrectionCode
+} // CyclicModelPipelineCreator
+//# sourceMappingURL=loops-pipeline-creator.js.map

package/src/pipeline/loops-pipeline-creator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"loops-pipeline-creator.js","sourceRoot":"","sources":["loops-pipeline-creator.ts"],"names":[],"mappings":"AACA,OAAO,EAAC,WAAW,EAAE,UAAU,EAAE,OAAO,EAAC,MAAM,oCAAoC,CAAC;AACpF,OAAO,EAAC,aAAa,EAAC,MAAM,UAAU,CAAC;AAEvC,OAAO,EAAC,eAAe,EAAC,MAAM,oBAAoB,CAAC;AAEnD,OAAO,EAAC,cAAc,EAAE,GAAG,EAAE,WAAW,EAAE,iBAAiB,EAAE,kBAAkB,EAAE,cAAc,EAAC,MAAM,aAAa,CAAC;AAEpH,wCAAwC;AACxC,MAAM,OAAO,0BAA2B,SAAQ,eAAe;IAC7D,YAAY,GAAQ;QAClB,KAAK,CAAC,GAAG,CAAC,CAAC;QAEX,IAAI,GAAG,CAAC,IAAI,KAAK,IAAI;YACnB,MAAM,IAAI,KAAK,CAAC,gEAAgE,CAAC,CAAC;IACtF,CAAC;IAED,mCAAmC;IAC5B,WAAW,CAAC,MAAoB;QACrC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,KAAK,IAAI;YACxB,MAAM,IAAI,KAAK,CAAC,gEAAgE,CAAC,CAAC;QAEpF,MAAM,OAAO,GAAG,cAAc,CAAC,IAAI,CAAC;QACpC,MAAM,UAAU,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC;QACrC,MAAM,WAAW,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC;QAEvC,IAAI,WAAW,GAAG,cAAc;YAC9B,MAAM,IAAI,KAAK,CAAC,gCAAgC,WAAW,iBAAiB,CAAC,CAAC;QAEhF,MAAM,YAAY,GAAa,EAAE,CAAC;QAElC,kBAAkB;QAClB,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE;YACrC,IAAI,GAAG,KAAK,OAAO;gBACjB,YAAY,CAAC,IAAI,CAAC,SAAS,UAAU,CAAC,GAAG,CAAC,kBAAkB,UAAU,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;;gBAEnF,YAAY,CAAC,IAAI,CAAC,uCAAuC,CAAC,CAAC;QAC/D,CAAC,CAAC,CAAC;QAEH,mBAAmB;QACnB,IAAI,CAAC,GAAG,CAAC,cAAc,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,EAAE;YACtC,YAAY,CAAC,IAAI,CAAC,SAAS,WAAW,CAAC,GAAG,CAAC,WAAW,WAAW,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC7E,CAAC,CAAC,CAAC;QAEH,kBAAkB;QAClB,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,KAAK,IAAI;YAC1B,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,KAAK,EAAE,IAAI,EAAE,EAAE,CAAC,YAAY,CAAC,IAAI,CAAC,SAAS,IAAI,MAAM,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAEjG,YAAY,CAAC,IAAI,CAAC,wCAAwC,UAAU,IAAI,CAAC,CAAC;QAE1E,iBAAiB;QACjB,cAAc,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC,YAAY,CAAC,IAAI,CAAC,OAAO,IAAI,aAAa,GAAG,IAAI,CAAC,CAAC,CAAC;QAE1F,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,EAAE;YAChD,YAAY,CAAC,IAAI,CAAC,OAAO,IAAI,aAAa,GAAG,GAAG,OAAO,IAAI,CAAC,CAAC;QAC/D,CAAC,CAAC,CAAC;QAEH,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,KAAK,IAAI,EAAE,CAAC;YAC7B,MAAM,KAAK,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC;YAC3D,IAAI,GAAG,GAAG,CAAC,CAAC;YAEZ,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE;gBAClC,YAAY,CAAC,IAAI,CAAC,OAAO,IAAI,aAAa,GAAG,GAAG,KAAK,IAAI,CAAC,CAAC;gBAC3D,EAAE,GAAG,CAAC;YACR,CAAC,CAAC,CAAC;QACL,CAAC;QAED,UAAU;QACV,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,YAAY,CAAC,IAAI,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC;QAEvE,sBAAsB;QACtB,cAAc,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC,YAAY,CAAC,IAAI,CAAC,UAAU,GAAG,OAAO,IAAI,GAAG,CAAC,CAAC,CAAC;QAEtF,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,EAAE;YAChD,YAAY,CAAC,IAAI,CAAC,UAAU,GAAG,GAAG,OAAO,OAAO,IAAI,GAAG,CAAC,CAAC;QAC3D,CAAC,CAAC,CAAC;QAEH,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,KAAK,IAAI,EAAE,CAAC;YAC7B,MAAM,KAAK,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC;YAC3D,IAAI,GAAG,GAAG,CAAC,CAAC;YAEZ,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE;gBAClC,YAAY,CAAC,IAAI,CAAC,UAAU,GAAG,GAAG,KAAK,OAAO,IAAI,GAAG,CAAC,CAAC;gBACvD,EAAE,GAAG,CAAC;YACR,CAAC,CAAC,CAAC;QACL,CAAC;QAED,YAAY,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;QAEpC,MAAM,aAAa,GAAG;YACpB,gCAAgC;YAChC,yCAAyC;YACzC,gDAAgD;YAChD,2BAA2B,cAAc,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,cAAc,cAAc,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI;YACxG,UAAU,cAAc,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,cAAc,cAAc,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI;YACvF,UAAU,cAAc,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,gBAAgB;SACzD,CAAC;QAEF,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,aAAa,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,EAAE;YAC7C,aAAa,CAAC,IAAI,CAAC,UAAU,GAAG,GAAG,OAAO,gBAAgB,GAAG,GAAG,CAAC,aAAa,CAAC,CAAC;QAClF,CAAC,CAAC,CAAC;QAEH,aAAa,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;QAErC,MAAM,OAAO,GAAG,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACxC,MAAM,QAAQ,GAAG,aAAa,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAE1C,MAAM,QAAQ,GAAG,IAAI,KAAK,CAAU,WAAW,CAAC,CAAC;QAEjD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,EAAE,CAAC,EAAE,CAAC;YACrC,QAAQ,CAAC,CAAC,CAAC,GAAG;gBACZ,OAAO,EAAE,OAAO;gBAChB,GAAG,EAAE,IAAI;gBACT,QAAQ,EAAE,QAAQ;aACnB,CAAC;QACJ,CAAC;QAED,QAAQ,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,yBAAyB,EAAE,CAAC;QAEtE,OAAO;YACL,QAAQ,EAAE,QAAQ;YAClB,GAAG,EAAE,aAAa,CAAC,IAAI,CAAC,GAAG,CAAC;SAC7B,CAAC;IACJ,CAAC,CAAC,cAAc;IAEhB,gDAAgD;IACxC,yBAAyB;QAC/B,OAAO,4BAA4B,WAAW;;;yCAGT,iBAAiB,KAAK,kBAAkB;;;;;yBAKxD,CAAC;IACxB,CAAC,CAAC,4BAA4B;CAC/B,CAAC,6BAA6B"}