diff-grok 1.0.8

package/src/examples/check-methods.ts

@@ -0,0 +1,62 @@
+// A script for checking performance
+
+import {ODEs, corrProbs, CorrProblem, mrt, ros3prw, ros34prw, perfProbs} from '../../index';
+
+/** Return numerical solution error: maximum absolute deviation between approximate & exact solutions */
+function getError(method: (odes: ODEs) => Float64Array[], corProb: CorrProblem): number {
+  let error = 0;
+
+  // Get numerical solution
+  const approxSolution = method(corProb.odes);
+
+  const exact = corProb.exact;
+
+  const arg = approxSolution[0];
+
+  const pointsCount = arg.length;
+  const funcsCount = approxSolution.length - 1;
+
+  // Compute error
+  for (let i = 0; i < pointsCount; ++i) {
+    const exactSolution = exact(arg[i]);
+
+    for (let j = 0; j < funcsCount; ++j)
+      error = Math.max(error, Math.abs(exactSolution[j] - approxSolution[j + 1][i]));
+  }
+
+  return error;
+}
+
+const methods = new Map([
+  ['MRT', mrt],
+  ['ROS3PRw', ros3prw],
+  ['ROS34PRw', ros34prw],
+]);
+
+console.log('Performance:\n');
+
+methods.forEach((method, name) => {
+  console.log('  ', name);
+
+  perfProbs.forEach((odes) => {
+    const start = Date.now();
+    method(odes);
+    const finish = Date.now();
+    console.log(`     ${odes.name}: ${finish - start} ms.`);
+  });
+
+  console.log();
+});
+
+console.log('Correctness (maximum absolute deviation):\n');
+
+methods.forEach((method, name) => {
+  console.log('  ', name);
+
+  corrProbs.forEach((problem) => {
+    const error = getError(method, problem);
+    console.log(`     ${problem.odes.name}: ${error}`);
+  });
+
+  console.log();
+});

package/src/examples/corr-probs.js

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

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

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

package/src/examples/corr-probs.ts

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

package/src/examples/cyclic-model.js

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

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

@@ -0,0 +1 @@
+{"version":3,"file":"cyclic-model.js","sourceRoot":"","sources":["cyclic-model.ts"],"names":[],"mappings":"AAAA,4BAA4B;AAC5B;;;;EAIE;AAEF,OAAO,KAAK,GAAG,MAAM,aAAa,CAAC;AAEnC,6BAA6B;AAC7B,MAAM,KAAK,GAAG;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;iBAuCG,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,MAAM,EAAE,EAAE;QACV,GAAG,EAAE,CAAC;QACN,GAAG,EAAE,EAAE;QACP,EAAE,EAAE,CAAC;QACL,KAAK,EAAE,CAAC;QACR,KAAK,EAAE,CAAC;QACR,IAAI,EAAE,CAAC;QACP,GAAG,EAAE,GAAG;QACR,IAAI,EAAE,KAAK;QACX,EAAE,EAAE,IAAI;QACR,EAAE,EAAE,GAAG;QACP,EAAE,EAAE,IAAI;QACR,CAAC,EAAE,KAAK;QACR,EAAE,EAAE,IAAI;QACR,IAAI,EAAE,IAAI;QACV,GAAG,EAAE,KAAK;QACV,IAAI,EAAE,KAAK;KACZ,CAAC;IACF,MAAM,WAAW,GAAG,GAAG,CAAC,cAAc,CAAC,MAAM,EAAE,GAAG,CAAC,CAAC;IAEpD,yBAAyB;IACzB,MAAM,OAAO,GAAG,GAAG,CAAC,kBAAkB,CAAC,GAAG,CAAC,CAAC;IAC5C,MAAM,QAAQ,GAAG,OAAO,CAAC,WAAW,CAAC,WAAW,CAAC,CAAC;IAElD,8CAA8C;IAC9C,MAAM,QAAQ,GAAG,GAAG,CAAC,aAAa,CAAC,QAAQ,EAAE,KAAK,EAAE,WAAW,CAAC,CAAC;IAEjE,qBAAqB;IAErB,sBAAsB;IACtB,IAAI,IAAI,GAAG,WAAW,CAAC;IACvB,WAAW,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,IAAI,IAAI,IAAI,GAAG,aAAa,CAAC,CAAC;IAC5D,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAElB,6BAA6B;IAC7B,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;QAChC,IAAI,GAAG,EAAE,CAAC;QAEV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,EAAE,CAAC;YAC9B,IAAI,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;QAE9C,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACpB,CAAC;AACH,CAAC;AAAC,OAAO,GAAG,EAAE,CAAC;IACb,OAAO,CAAC,GAAG,CAAC,qBAAqB,EAAE,GAAG,YAAY,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,kBAAkB,CAAC,CAAC;AAC9F,CAAC"}

package/src/examples/cyclic-model.ts

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

package/src/examples/e5.js

@@ -0,0 +1,30 @@
+/** Kintetic constants for the E5 model */
+var E5;
+(function (E5) {
+    E5[E5["K1"] = 7.89e-10] = "K1";
+    E5[E5["K2"] = 1130000000] = "K2";
+    E5[E5["K3"] = 11000000] = "K3";
+    E5[E5["K4"] = 1130] = "K4";
+})(E5 || (E5 = {}));
+;
+/** The E5 model (chemical pyrolysis: https://archimede.uniba.it/~testset/report/e5.pdf) */
+export const e5 = {
+    name: 'E5',
+    arg: { name: 't', start: 0, finish: 1e13, step: 2.5e8 },
+    initial: [0.00176, 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];
+        // compute output
+        output[0] = -E5.K1 * y1 - E5.K3 * y1 * y3;
+        output[1] = E5.K1 * y1 - E5.K2 * y2 * y3;
+        output[2] = E5.K1 * y1 - E5.K2 * y2 * y3 - E5.K3 * y1 * y3 + E5.K4 * y4;
+        output[3] = E5.K3 * y1 * y3 - E5.K4 * y4;
+    },
+    tolerance: 1e-6,
+    solutionColNames: ['y1', 'y2', 'y3', 'y4'],
+};
+//# sourceMappingURL=e5.js.map

package/src/examples/e5.js.map

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

package/src/examples/e5.ts

@@ -0,0 +1,29 @@
+ /** Kintetic constants for the E5 model */
+ enum E5 {
+    K1 = 7.89e-10,
+    K2 = 1.13e9,
+    K3 = 1.1e7,
+    K4 = 1.13e3,
+  };
+
+/** The E5 model (chemical pyrolysis: https://archimede.uniba.it/~testset/report/e5.pdf) */
+export const e5 = {
+  name: 'E5',
+  arg: {name: 't', start: 0, finish: 1e13, step: 2.5e8},
+  initial: [0.00176, 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];
+
+    // compute output
+    output[0] = -E5.K1 * y1 - E5.K3 * y1 * y3;
+    output[1] = E5.K1 * y1 - E5.K2 * y2 * y3;
+    output[2] = E5.K1 * y1 - E5.K2 * y2 * y3 - E5.K3 * y1 * y3 + E5.K4 * y4;
+    output[3] = E5.K3 * y1 * y3 - E5.K4 * y4;
+  },
+  tolerance: 1e-6,
+  solutionColNames: ['y1', 'y2', 'y3', 'y4'],
+};

package/src/examples/hires.js

@@ -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, 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];
+        // 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
+//# sourceMappingURL=hires.js.map

package/src/examples/hires.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"hires.js","sourceRoot":"","sources":["hires.ts"],"names":[],"mappings":"AAAA,mHAAmH;AACnH,MAAM,CAAC,MAAM,KAAK,GAAG;IACnB,IAAI,EAAE,OAAO;IACb,GAAG,EAAE,EAAC,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,QAAQ,EAAE,IAAI,EAAE,IAAI,EAAC;IACxD,OAAO,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC;IACtC,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;QAEhB,iBAAiB;QACjB,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,MAAM,CAAC;QACxD,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;QAClC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,KAAK,GAAG,EAAE,CAAC;QACjD,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;QAC9C,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;QAChD,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;QAC3E,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;QACtC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,CAAC;IACzC,CAAC;IACD,SAAS,EAAE,KAAK;IAChB,gBAAgB,EAAE,CAAC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,CAAC;CACnE,CAAC,CAAC,QAAQ"}