diff-grok 1.2.0 → 1.3.1

package/src/examples/check-methods.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"check-methods.js","sourceRoot":"","sources":["check-methods.ts"],"names":[],"mappings":"AAAA,oCAAoC;AAEpC,OAAO,EAAO,SAAS,EAAe,GAAG,EAAE,OAAO,EAAE,QAAQ,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,KAAK,EAAE,KAAK,EAClG,SAAS,EAAE,SAAS,EAAC,MAAM,aAAa,CAAC;AAE3C,wGAAwG;AACxG,SAAS,QAAQ,CAAC,MAAsC,EAAE,OAAoB;IAC5E,IAAI,KAAK,GAAG,CAAC,CAAC;IAEd,yBAAyB;IACzB,MAAM,cAAc,GAAG,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;IAE5C,MAAM,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC;IAE5B,MAAM,GAAG,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC;IAE9B,MAAM,WAAW,GAAG,GAAG,CAAC,MAAM,CAAC;IAC/B,MAAM,UAAU,GAAG,cAAc,CAAC,MAAM,GAAG,CAAC,CAAC;IAE7C,gBAAgB;IAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,EAAE,CAAC,EAAE,CAAC;QACrC,MAAM,aAAa,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAEpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,EAAE,EAAE,CAAC;YACjC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,aAAa,CAAC,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnF,CAAC;IAED,OAAO,KAAK,CAAC;AACf,CAAC,CAAC,WAAW;AAEb,SAAS,8BAA8B,CAAC,QAAwB,EAAE,QAAsB;IACtF,IAAI,QAAQ,GAAG,CAAC,CAAC;IACjB,MAAM,OAAO,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IACvC,IAAI,GAAG,GAAG,CAAC,CAAC;IAEZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,EAAE,CAAC,EAAE,CAAC;QACzC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;QACnD,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,QAAQ,CAAC,CAAC;IACrC,CAAC;IAED,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC,iCAAiC;AAEnC,MAAM,eAAe,GAAG,IAAI,GAAG,CAAC;IAC9B,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,SAAS,EAAE,OAAO,CAAC;IACpB,CAAC,UAAU,EAAE,QAAQ,CAAC;IACtB,CAAC,OAAO,EAAE,KAAK,CAAC;IAChB,CAAC,OAAO,EAAE,KAAK,CAAC;CACjB,CAAC,CAAC;AAEH,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC;IACtB,GAAG,eAAe;IAClB,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,KAAK,EAAE,GAAG,CAAC;IACZ,CAAC,MAAM,EAAE,IAAI,CAAC;IACd,CAAC,KAAK,EAAE,GAAG,CAAC;CACb,CAAC,CAAC;AAEH,OAAO,CAAC,GAAG,CAAC,eAAe,CAAC,CAAC;AAE7B,eAAe,CAAC,OAAO,CAAC,CAAC,MAAM,EAAE,IAAI,EAAE,EAAE;IACvC,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;IAEvB,OAAO,CAAC,GAAG,CAAC,qCAAqC,CAAC,CAAC;IAEnD,SAAS,CAAC,OAAO,CAAC,CAAC,IAAI,EAAE,GAAG,EAAE,EAAE;QAC9B,IAAI,IAAI,KAAK,OAAO,IAAI,IAAI,CAAC,IAAI,KAAK,IAAI;YACxC,OAAO;QAET,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;QACzB,MAAM,QAAQ,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;QAC9B,MAAM,MAAM,GAAG,8BAA8B,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;QAClG,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC;QAC1B,OAAO,CAAC,GAAG,CAAC,KAAK,IAAI,CAAC,IAAI,QAAQ,MAAM,GAAG,KAAK,QAAQ,MAAM,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IACjF,CAAC,CAAC,CAAC;IAEH,OAAO,CAAC,GAAG,EAAE,CAAC;AAChB,CAAC,CAAC,CAAC;AAEH,OAAO,CAAC,GAAG,CAAC,mDAAmD,CAAC,CAAC;AAEjE,OAAO,CAAC,GAAG,CAAC,mEAAmE,CAAC,CAAC;AAEjF,OAAO,CAAC,OAAO,CAAC,CAAC,MAAM,EAAE,IAAI,EAAE,EAAE;IAC/B,OAAO,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IACxB,OAAO,CAAC,GAAG,CAAC,yBAAyB,CAAC,CAAC;IAEvC,SAAS,CAAC,OAAO,CAAC,CAAC,OAAO,EAAE,EAAE;QAC5B,MAAM,KAAK,GAAG,QAAQ,CAAC,MAAM,EAAE,OAAO,CAAC,CAAC;QACxC,OAAO,CAAC,GAAG,CAAC,QAAQ,OAAO,CAAC,IAAI,CAAC,IAAI,OAAO,KAAK,CAAC,aAAa,CAAC,CAAC,CAAC,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IACrF,CAAC,CAAC,CAAC;IAEH,OAAO,CAAC,GAAG,EAAE,CAAC;AAChB,CAAC,CAAC,CAAC"}

package/src/examples/check-methods.ts

@@ -1,96 +0,0 @@
-// A script for checking performance
-
-import {ODEs, corrProbs, CorrProblem, mrt, ros3prw, ros34prw, rk4, ab5, ab4, rkdp, rk3, lsoda, cvode,
-  perfProbs, refPoints} 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;
-} // getError
-
-function getDeviationFromReferencePoint(solution: Float64Array[], refPoint: Float64Array): number {
-  let absolute = 0;
-  const lastIdx = solution[0].length - 1;
-  let cur = 0;
-
-  for (let k = 0; k < solution.length; ++k) {
-    cur = Math.abs(solution[k][lastIdx] - refPoint[k]);
-    absolute = Math.max(cur, absolute);
-  }
-
-  return absolute;
-} // getDeviationFromReferencePoint
-
-const implicitMethods = new Map([
-  ['MRT', mrt],
-  ['ROS3PRw', ros3prw],
-  ['ROS34PRw', ros34prw],
-  ['LSODA', lsoda],
-  ['CVODE', cvode],
-]);
-
-const methods = new Map([
-  ...implicitMethods,
-  ['RK4', rk4],
-  ['AB5', ab5],
-  ['AB4', ab4],
-  ['RKDP', rkdp],
-  ['RK3', rk3],
-]);
-
-console.log('Performance\n');
-
-implicitMethods.forEach((method, name) => {
-  console.log(' ', name);
-
-  console.log('           PROBLEM   TIME,MS   ERR*');
-
-  perfProbs.forEach((odes, idx) => {
-    if (name === 'CVODE' && odes.name === 'E5')
-      return;
-
-    const start = Date.now();
-    const solution = method(odes);
-    const absErr = getDeviationFromReferencePoint(solution.slice(1), refPoints[idx]).toExponential(2);
-    const finish = Date.now();
-    console.log(`  ${odes.name}     ${finish - start}     ${absErr}`.padStart(37));
-  });
-
-  console.log();
-});
-
-console.log('* approximate vs. reference (also, numerical)\n\n');
-
-console.log('Correctness (maximum absolute deviation): approximate vs. exact\n');
-
-methods.forEach((method, name) => {
-  console.log('  ', name);
-  console.log('          PROBLEM   ERR');
-
-  corrProbs.forEach((problem) => {
-    const error = getError(method, problem);
-    console.log(`     ${problem.odes.name}    ${error.toExponential(2)}`.padStart(28));
-  });
-
-  console.log();
-});

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

@@ -1 +0,0 @@
-{"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;AAUF;;GAEG;AACH,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

@@ -1,155 +0,0 @@
-// Problems with known exact solution (for testing correctness)
-
-import {ODEs} from '../../index';
-
-/** 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
- * @internal
- */
-export type CorrProblem = {
-    odes: ODEs,
-    exact: (t: number) => Float64Array,
-  };
-
-/** Problems for testing correctness
- * @internal
- */
-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

@@ -1,102 +0,0 @@
-/* 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
-
-#meta.solver: {method: 'rkdp'; maxTimeMs: 5000}`;
-/** 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

@@ -1 +0,0 @@
-{"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;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;gDAyCkC,CAAC;AAEjD,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

@@ -1,112 +0,0 @@
-/* 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
-
-#meta.solver: {method: 'rkdp'; maxTimeMs: 5000}`;
-
-/** 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

@@ -1,36 +0,0 @@
-/** 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'],
-};
-export const e5ReferencePoint = new Float64Array([
-    0.1152903278711829E-290,
-    0.8867655517642120E-22,
-    0.8854814626268838E-22,
-    0,
-]);
-//# sourceMappingURL=e5.js.map

package/src/examples/e5.js.map

@@ -1 +0,0 @@
-{"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;AAEF,MAAM,CAAC,MAAM,gBAAgB,GAAG,IAAI,YAAY,CAAC;IAC/C,uBAAuB;IACvB,sBAAsB;IACtB,sBAAsB;IACtB,CAAC;CACF,CAAC,CAAC"}

package/src/examples/e5.ts

@@ -1,36 +0,0 @@
- /** 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'],
-};
-
-export const e5ReferencePoint = new Float64Array([
-  0.1152903278711829E-290,
-  0.8867655517642120E-22,
-  0.8854814626268838E-22,
-  0,
-]);

package/src/examples/hires.js

@@ -1,39 +0,0 @@
-/** 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
-export const hiresReferencePoint = new Float64Array([
-    0.7371312573325668E-3,
-    0.1442485726316185E-3,
-    0.5888729740967575E-4,
-    0.1175651343283149E-2,
-    0.2386356198831331E-2,
-    0.6238968252742796E-2,
-    0.2849998395185769E-2,
-    0.2850001604814231E-2,
-]);
-//# sourceMappingURL=hires.js.map

package/src/examples/hires.js.map

@@ -1 +0,0 @@
-{"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;AAEX,MAAM,CAAC,MAAM,mBAAmB,GAAG,IAAI,YAAY,CAAC;IAClD,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;IACrB,qBAAqB;CACtB,CAAC,CAAC"}

package/src/examples/hires.ts

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

package/src/examples/index.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,SAAS,EAAE,uBAAuB,EAAC,MAAM,aAAa,CAAC;AAC/D,OAAO,EAAC,KAAK,EAAE,mBAAmB,EAAC,MAAM,SAAS,CAAC;AACnD,OAAO,EAAC,KAAK,EAAE,mBAAmB,EAAC,MAAM,SAAS,CAAC;AACnD,OAAO,EAAC,KAAK,EAAE,mBAAmB,EAAC,MAAM,SAAS,CAAC;AACnD,OAAO,EAAC,EAAE,EAAE,gBAAgB,EAAC,MAAM,MAAM,CAAC;AAC1C,OAAO,EAAC,SAAS,EAAE,uBAAuB,EAAC,MAAM,aAAa,CAAC;AAE/D;;GAEG;AACH,MAAM,CAAC,MAAM,SAAS,GAAG;IACvB,SAAS;IACT,KAAK;IACL,KAAK;IACL,KAAK;IACL,EAAE;IACF,SAAS;CACV,CAAC;AAEF;;GAEG;AACH,MAAM,CAAC,MAAM,SAAS,GAAG;IACvB,uBAAuB;IACvB,mBAAmB;IACnB,mBAAmB;IACnB,mBAAmB;IACnB,gBAAgB;IAChB,uBAAuB;CACxB,CAAC;AAEF,OAAO,EAAC,SAAS,EAAc,MAAM,cAAc,CAAC;AAEpD,OAAO,EAAC,cAAc,EAAE,UAAU,EAAE,UAAU,EAAE,OAAO,EAAE,UAAU,EAAE,cAAc,EAAC,MAAM,mBAAmB,CAAC"}