diff-grok 1.0.8 → 1.0.10

package/src/tests/test-utils.ts

@@ -0,0 +1,37 @@
+import {ODEs, CorrProblem, getIVP, getIvp2WebWorker, getInputVector, getPipelineCreator,
+  applyPipeline} from '../../index';
+
+/** Return numerical solution error: maximum absolute deviation between approximate & exact solutions */
+export function getError(method: (odes: ODEs) => Float64Array[], corProb: CorrProblem): number {
+  let error = 0;
+
+  // Get numerical solution
+  const approxSolution = method(corProb.odes);
+
+  const exact = corProb.exact;
+
+  const arg = approxSolution[0];
+
+  const pointsCount = arg.length;
+  const funcsCount = approxSolution.length - 1;
+
+  // Compute error
+  for (let i = 0; i < pointsCount; ++i) {
+    const exactSolution = exact(arg[i]);
+
+    for (let j = 0; j < funcsCount; ++j)
+      error = Math.max(error, Math.abs(exactSolution[j] - approxSolution[j + 1][i]));
+  }
+
+  return error;
+} // getError
+
+/** Evaluate model  */
+export function evalModel(model: string, inputs: Record<string, number>): Float64Array[] {
+  const ivp = getIVP(model);
+  const ivpWW = getIvp2WebWorker(ivp);
+  const inputVector = getInputVector(inputs, ivp);
+  const creator = getPipelineCreator(ivp);
+  const pipeline = creator.getPipeline(inputVector);
+  return applyPipeline(pipeline, ivpWW, inputVector);
+} // evalModel

package/src/worker-tools/scripting.js

@@ -1,3 +1,4 @@
+/* eslint-disable max-len */
 // Scripting tools for solving IVPs in web-workers
 import { MATH_CONSTS, POW_IDX, MATH_FUNCS } from '../scripting-tools/scripting-tools';
 const DEFAULT_METHOD = 'ros34prw';
@@ -45,7 +46,100 @@ function getFuncMainBody(ivp) {
     ivp.deqs.solutionNames.forEach((name, idx) => lines.push(`arguments[2][${idx}] = ${ivp.deqs.equations.get(name)};`));
     return lines.join('\n');
 }
-/** Return IVP-structure to be passed to web-worker */
+/** Return IVP-structure to be passed to web-worker
+ * @param ivp - initial value problem to be solved
+ * @example
+ *  // Here, we consider modeling queues: https://en.wikipedia.org/wiki/M/M/c_queue
+ * import * as DGL from 'diff-grok';
+ *
+ * // 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!');
+ * }
+ */
 export function getIvp2WebWorker(ivp) {
     const paramNames = [];
     if (ivp.params !== null)

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

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

package/src/worker-tools/scripting.ts

@@ -1,9 +1,18 @@
+/* eslint-disable max-len */
 // Scripting tools for solving IVPs in web-workers
 
 import {IVP, MATH_CONSTS, POW_IDX, MATH_FUNCS} from '../scripting-tools/scripting-tools';
 import {SolverOptions} from '../solver-tools/solver-defs';
 
-/** Initial value problem structure to be passed to web-worker */
+/**
+ * Represents a fully parsed simulation prepared for execution in Web Workers.
+ *
+ * @property deqsCount - The number of differential equations in the system.
+ * @property funcMainBody - The main body of the function implementing the ODE system.
+ * @property paramNames - Names of the model parameters.
+ * @property solverOpts - Configuration options for the numerical solver (partial).
+ * @property tolerance - Tolerance for the numerical solver.
+ */
 export type IVP2WebWorker = {
     tolerance: number,
     solverOpts: Partial<SolverOptions>,
@@ -68,7 +77,100 @@ function getFuncMainBody(ivp: IVP): string {
   return lines.join('\n');
 }
 
-/** Return IVP-structure to be passed to web-worker */
+/** Return IVP-structure to be passed to web-worker
+ * @param ivp - initial value problem to be solved
+ * @example
+ *  // Here, we consider modeling queues: https://en.wikipedia.org/wiki/M/M/c_queue
+ * import * as DGL from 'diff-grok';
+ *
+ * // 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!');
+ * }
+ */
 export function getIvp2WebWorker(ivp: IVP): IVP2WebWorker {
   const paramNames: string[] = [];
   if (ivp.params !== null)

package/src/worker-tools/solving.js

@@ -22,7 +22,11 @@ function checkInputs(ivp, inputs) {
     if (expected > inputs.length)
         throw new Error(`Incorrect inputs count, expected: ${expected}, current: ${inputs.length}`);
 }
-/** Return solution of IVP with a specified inputs */
+/** Return solution of IVP with a specified inputs
+ * @param ivp - initial value problem to be solved
+ * @param inputs - vector of input values
+ * @returns solution of the given initial value problem
+ */
 export function solveIvp(ivp, inputs) {
     checkInputs(ivp, inputs);
     const nonParamInputsCount = inputs.length - ivp.paramNames.length;

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

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

package/src/worker-tools/solving.ts

@@ -31,7 +31,11 @@ function checkInputs(ivp: IVP2WebWorker, inputs: Float64Array): void {
     throw new Error(`Incorrect inputs count, expected: ${expected}, current: ${inputs.length}`);
 }
 
-/** Return solution of IVP with a specified inputs */
+/** Return solution of IVP with a specified inputs
+ * @param ivp - initial value problem to be solved
+ * @param inputs - vector of input values
+ * @returns solution of the given initial value problem
+ */
 export function solveIvp(ivp: IVP2WebWorker, inputs: Float64Array): Float64Array[] {
   checkInputs(ivp, inputs);
   const nonParamInputsCount = inputs.length - ivp.paramNames.length;

package/tsconfig.json

@@ -48,7 +48,7 @@
     // "paths": {},                                 /* A series of entries which re-map imports to lookup locations relative to the 'baseUrl'. */
     // "rootDirs": [],                              /* List of root folders whose combined content represents the structure of the project at runtime. */
     // "typeRoots": [],                             /* List of folders to include type definitions from. */    
-    // "types": [],                                 /* Type declaration files to be included in compilation. */
+    "types": ["jest", "node"],                                 /* Type declaration files to be included in compilation. */
     // "allowSyntheticDefaultImports": true,        /* Allow default imports from modules with no default export. This does not affect code emit, just typechecking. */
     "esModuleInterop": true,                        /* Enables emit interoperability between CommonJS and ES Modules via creation of namespace objects for all imports. Implies 'allowSyntheticDefaultImports'. */
     // "preserveSymlinks": true,                    /* Do not resolve the real path of symlinks. */

package/typedoc.json

@@ -0,0 +1,8 @@
+{
+  "entryPoints": ["index.ts"],
+  "out": "docs",
+  "includeVersion": true,
+  "excludePrivate": true,
+  "excludeProtected": true,
+  "excludeInternal": true
+}