diff --git a/DESCRIPTION b/DESCRIPTION
index 325ba1de..bf2195fd 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: odin
 Title: ODE Generation and Integration
-Version: 1.4.1
+Version: 1.4.2
 Authors@R: c(person("Rich", "FitzJohn", role = c("aut", "cre"),
                     email = "rich.fitzjohn@gmail.com"),
              person("Thibaut", "Jombart", role = "ctb"),
diff --git a/inst/js/dust.js b/inst/js/dust.js
index cdccc0bd..907a7a1a 100644
--- a/inst/js/dust.js
+++ b/inst/js/dust.js
@@ -1,2 +1,2 @@
 /*! For license information please see dust.js.LICENSE.txt */
-var dust;(()=>{var t={8006:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.interpolateSearch=e.InterpolatorBase=void 0;var r=function(){function t(t,e){(function(t){return Array.isArray(t[0])})(e)||(e=[e]),this._i=0,this._x=t,this._y=e,this.nX=this._x.length,this.nY=this._y.length;for(var r=0;r<this.nY;++r)if(this._y[r].length!==this.nX)throw Error("Invalid length for 'y', expected ".concat(this.nX))}return t.prototype.evalAll=function(t){for(var e=Array(this.nY),r=0;r<this.nY;++r)e[r]=this.eval(t,r);return e},t.prototype.search=function(t,e){var r=n(t,this._x,this._i);if(r<0||!e&&r===this.nX)throw Error("Interpolation failed as 'x' is out of range");return this._i=r,r},t}();function n(t,e,r){var n=r,i=r,o=1,a=e.length;if(e[n]<=t){if(n>=a-1)return a;for(i=n+o;e[i]<t;)if(n=i,(i+=o*=2)>=a){if(e[i=a-1]<t)return a;break}}else{if(0===n)return-1;for(n-=o;e[n]>t;){if(i=n,n<(o*=2)){if(e[n=0]>t)return-1;break}n-=o}}for(;i-n>1;){var s=Math.floor((i+n)/2);e[s]<t?n=s:i=s}return n}e.InterpolatorBase=r,e.interpolateSearch=n},4731:function(t,e,r){"use strict";var n,i=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorConstant=void 0;var o=function(t){function e(){return null!==t&&t.apply(this,arguments)||this}return i(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=Math.min(this.search(t,!0),this.nX-1);return r!==this.nX-1&&this._x[r+1]===t&&++r,this._y[e][r]},e}(r(8006).InterpolatorBase);e.InterpolatorConstant=o},9673:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorSpline=e.InterpolatorLinear=e.InterpolatorConstant=e.InterpolatorBase=void 0;var n=r(8006);Object.defineProperty(e,"InterpolatorBase",{enumerable:!0,get:function(){return n.InterpolatorBase}});var i=r(4731);Object.defineProperty(e,"InterpolatorConstant",{enumerable:!0,get:function(){return i.InterpolatorConstant}});var o=r(9477);Object.defineProperty(e,"InterpolatorLinear",{enumerable:!0,get:function(){return o.InterpolatorLinear}});var a=r(7869);Object.defineProperty(e,"InterpolatorSpline",{enumerable:!0,get:function(){return a.InterpolatorSpline}})},9477:function(t,e,r){"use strict";var n,i=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorLinear=void 0;var o=function(t){function e(){return null!==t&&t.apply(this,arguments)||this}return i(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=this.search(t,!1),n=this._x[r],i=(t-n)/(this._x[r+1]-n),o=this._y[e],a=o[r];return a+(o[r+1]-a)*i},e}(r(8006).InterpolatorBase);e.InterpolatorLinear=o},7869:function(t,e,r){"use strict";var n,i=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)}),o=this&&this.__importDefault||function(t){return t&&t.__esModule?t:{default:t}};Object.defineProperty(e,"__esModule",{value:!0}),e.solve=e.InterpolatorSpline=void 0;var a=o(r(8863)),s=function(t){function e(e,r){var n=t.call(this,e,r)||this,i=function(t){var e=t.length,r=Array(e),n=Array(e),i=Array(e);r[0]=0,n[0]=2/(t[1]-t[0]),i[0]=1/(t[1]-t[0]);for(var o=1;o<e-1;++o){var a=t[o-1],s=t[o],u=t[o+1];r[o]=1/(s-a),n[o]=2*(1/(s-a)+1/(u-s)),i[o]=1/(u-s)}return r[e-1]=1/(t[e-1]-t[e-2]),n[e-1]=2/(t[e-1]-t[e-2]),i[e-1]=0,[r,n,i]}(n._x),o=function(t,e){for(var r=t.length,n=e.length,i=r-1,o=[],a=0;a<n;++a){var s=Array(r),u=e[a];s[0]=3*(u[1]-u[0])/Math.pow(t[1]-t[0],2);for(var c=1;c<i;++c){var l=t[c-1],f=t[c],h=t[c+1],p=u[c-1],d=u[c],v=u[c+1];s[c]=3*((d-p)/Math.pow(f-l,2)+(v-d)/Math.pow(h-f,2))}s[i]=3*(u[i]-u[i-1])/Math.pow(t[i]-t[i-1],2),o.push(s)}return o}(n._x,n._y);return n._k=function(t,e){for(var r=t[0],n=t[1],i=t[2],o=r.length,a=0,s=e;a<s.length;a++){var c=s[a];u(o,r,n.slice(),i,c)}return e}(i,o),n}return i(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=this.search(t,!1);return Array(this.nY),function(t,e,r,n,i){var o=(e-r[t])/(r[t+1]-r[t]),a=i[t]*(r[t+1]-r[t])-(n[t+1]-n[t]),s=-i[t+1]*(r[t+1]-r[t])+(n[t+1]-n[t]);return(1-o)*n[t]+o*n[t+1]+o*(1-o)*(a*(1-o)+s*o)}(r,t,this._x,this._y[e],this._k[e])},e}(r(8006).InterpolatorBase);function u(t,e,r,n,i){if(!(0,a.default)(t,e,r,n,i))throw Error("solve failed: singular matrix?")}e.InterpolatorSpline=s,e.solve=u},8251:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.base=void 0;const n=r(3350),i=r(320),o=r(9928),a=r(2630);e.base={delay:n.delay,interpolate:i.interpolate,maths:o.maths,user:a.user}},7601:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.updatePars=e.batchParsDisplace=e.batchParsRange=e.batchRun=e.Batch=void 0;const n=r(46),i=r(1729),o=r(7466);class a{constructor(t,e,r,n,i){this.pars=e,this.tStart=r,this.tEnd=n,this.solutions=e.values.map((a=>{const s=u(e.base,e.name,a);return(0,o.wodinRun)(t,s,r,n,i)}))}valueAtTime(t){const e=this.solutions.map((e=>e({mode:n.TimeMode.Given,times:[t]})));return{names:e[0].names,x:this.pars.values,y:e[0].names.map(((t,r)=>e.map((t=>t.y[r][0]))))}}extreme(t){return this.findExtremes()[t]}findExtremes(){if(void 0===this._extremes){const t=501,e={mode:n.TimeMode.Grid,nPoints:t,tEnd:this.tEnd,tStart:this.tStart},r=this.solutions.map((t=>t(e))),o=r[0].names,a=this.pars.values,s=(0,i.loop)(o.length,(t=>r.map((e=>function(t,e){const r=e.x,n=e.y[t],o=(0,i.whichMin)(n),a=(0,i.whichMax)(n),s=r[o],u=r[a],c=n[o];return{tMax:u,tMin:s,yMax:n[a],yMin:c}}(t,e)))));this._extremes={tMax:l("tMax",o,a,s),tMin:l("tMin",o,a,s),yMax:l("yMax",o,a,s),yMin:l("yMin",o,a,s)}}return this._extremes}}function s(t,e,r,n,o,a){const s=c(t,e);if(o>s)throw Error(`Expected lower bound to be no greater than ${s}`);if(a<s)throw Error(`Expected upper bound to be no less than ${s}`);if(o>=a)throw Error("Expected upper bound to be greater than lower bound");if(r<2)throw Error("Must include at least 2 traces in the batch");if(n&&o<=0)throw Error("Lower bound must be greater than 0 for logarithmic scale");return{base:t,name:e,values:n?(0,i.gridLog)(o,a,r):(0,i.grid)(o,a,r)}}function u(t,e,r){const n=Object.assign({},t);return n[e]=r,n}function c(t,e){const r=t[e];if(void 0===r)throw Error(`Expected a value for '${e}'`);if("number"!=typeof r)throw Error(`Expected a number for '${e}'`);return r}function l(t,e,r,n){return{names:e,x:r,y:(0,i.loop)(e.length,(e=>n[e].map((e=>e[t]))))}}e.Batch=a,e.batchRun=function(t,e,r,n,i){return new a(t,e,r,n,i)},e.batchParsRange=s,e.batchParsDisplace=function(t,e,r,n,i){const o=c(t,e),a=i/100;return s(t,e,r,n,o*(1-a),o*(1+a))},e.updatePars=u},3350:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.delay=void 0,e.delay=function(t,e,r,n){const i=t(e);for(let t=0;t<r.length;++t)n[t]=i[r[t]]}},7771:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.sumOfSquares=e.updatePars=e.fitTarget=void 0;const n=r(8251),i=r(1267),o=r(46);function a(t,e){const r=Object.assign({},t.base);for(let n=0;n<t.vary.length;++n)r[t.vary[n]]=e[n];return r}function s(t,e){let r=0;for(let n=0;n<t.length;++n){const i=t[n];isNaN(i)||(r+=Math.pow(i-e[n],2))}return r}e.fitTarget=function(t,e,r,u,c){if(e.time[0]<0)throw Error("Expected the first time to be at least 0");const l=e.time[e.time.length-1];return f=>{const h=a(r,f),p=new t(n.base,h,"error"),d=(0,i.runModel)(p,null,0,l,c).solution,v=p.names(),y=v.indexOf(u),m=d(e.time).map((t=>t[y]));return{data:{endTime:l,names:v,pars:h,solution:(0,o.interpolatedSolution)(d,v,0,l),startTime:0},value:s(e.value,m)}}},e.updatePars=a,e.sumOfSquares=s},1688:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.versions=e.batchRun=e.batchParsRange=e.batchParsDisplace=e.Batch=e.base=e.PkgWrapper=e.wodinRun=e.wodinFitValue=e.wodinFit=void 0;var n=r(7466);Object.defineProperty(e,"wodinFit",{enumerable:!0,get:function(){return n.wodinFit}}),Object.defineProperty(e,"wodinFitValue",{enumerable:!0,get:function(){return n.wodinFitValue}}),Object.defineProperty(e,"wodinRun",{enumerable:!0,get:function(){return n.wodinRun}});var i=r(8714);Object.defineProperty(e,"PkgWrapper",{enumerable:!0,get:function(){return i.PkgWrapper}});var o=r(8251);Object.defineProperty(e,"base",{enumerable:!0,get:function(){return o.base}});var a=r(7601);Object.defineProperty(e,"Batch",{enumerable:!0,get:function(){return a.Batch}}),Object.defineProperty(e,"batchParsDisplace",{enumerable:!0,get:function(){return a.batchParsDisplace}}),Object.defineProperty(e,"batchParsRange",{enumerable:!0,get:function(){return a.batchParsRange}}),Object.defineProperty(e,"batchRun",{enumerable:!0,get:function(){return a.batchRun}});var s=r(7497);Object.defineProperty(e,"versions",{enumerable:!0,get:function(){return s.versions}})},320:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.interpolate=e.interpolateTimes=e.interpolateCheckT=e.interpolateCheckY=e.interpolateAlloc=void 0;const n=r(9673);function i(t,e,r){const i=function(t,e){if(t.length===e)return[t];const r=[],n=t.length/e;for(let i=0;i<n;++i)r.push(t.slice(e*i,e*(i+1)));return r}(r,e.length);if("constant"===t)return new n.InterpolatorConstant(e,i);if("linear"===t)return new n.InterpolatorLinear(e,i);if("spline"===t)return new n.InterpolatorSpline(e,i);throw Error(`Invalid interpolation method '${t}'`)}function o(t,e,r,n){if(1===e.length){if(t[0]!==e[0])throw Error(`Expected ${r} to have length ${t[0]} (for ${n})`)}else for(let i=0;i<e.length;++i)if(t[i]!==e[i])throw Error(`Expected dimension ${i+1} of ${r} to have size ${t[i]} (for ${n})`)}function a(t,e){const r=Math.max(...t);return{max:Math.min(...e),min:r}}e.interpolateAlloc=i,e.interpolateCheckY=o,e.interpolateCheckT=function(t,e,r,n){if(void 0===r)return 1/0;if(t<r.min)throw Error(`Integration times do not span interpolation range; min: ${r.min}`);if(e>r.max)throw Error(`Integration times do not span interpolation range; max: ${r.max}`);return void 0===n&&(n=r.max),n},e.interpolateTimes=a,e.interpolate={alloc:i,checkY:o,times:a}},9928:(t,e)=>{"use strict";function r(t){return.5===i(t,1)?2*Math.round(t/2):Math.round(t)}function n(t,e){if(void 0===e||0===e)return r(t);{const n=Math.pow(10,e);return r(t*n)/n}}function i(t,e){let r=t%e;return r*e<0&&(r+=e),r}function o(t,e){return Math.floor(t/e)}function a(t,e,r){let n=0;for(let i=e;i<r;++i)n+=t[i];return n}function s(t,e,r,n,i,o){let a=0;for(let s=n;s<i;++s){const n=s*o;for(let i=e;i<r;++i)a+=t[i+n]}return a}function u(t,e,r,n,i,o,a,s,u){let c=0;for(let l=o;l<a;++l){const o=l*u;for(let a=n;a<i;++a){const n=a*s+o;for(let i=e;i<r;++i)c+=t[i+n]}}return c}Object.defineProperty(e,"__esModule",{value:!0}),e.maths=e.odinSum3=e.odinSum2=e.odinSum1=e.intdivr=e.modr=e.round2=void 0,e.round2=n,e.modr=i,e.intdivr=o,e.odinSum1=a,e.odinSum2=s,e.odinSum3=u,e.maths={intdivr:o,modr:i,odinSum1:a,odinSum2:s,odinSum3:u,round2:n}},1267:function(t,e,r){"use strict";var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var i=Object.getOwnPropertyDescriptor(e,r);i&&!("get"in i?!e.__esModule:i.writable||i.configurable)||(i={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,i)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),i=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),o=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return i(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.runModel=e.isODEModel=e.isDDEModel=void 0;const a=o(r(6245)),s=r(320);function u(t){return 4===t.rhs.length}e.isDDEModel=u,e.isODEModel=function(t){return 3===t.rhs.length},e.runModel=function(t,e,r,n,i){const o=t.getMetadata().interpolateTimes;return i.tcrit=(0,s.interpolateCheckT)(r,n,o,i.tcrit),u(t)?function(t,e,r,n,i){let o=null;"function"==typeof t.output&&(o=(e,r,n)=>t.output(e,r,n)),null===e&&(e=t.initial(r)),t.getInternal().initial_t=r;const s=new a.DDE((function(e,r,n,i){t.rhs(e,r,n,i)}),e.length,i,o);return s.initialise(r,e),{solution:s.run(n),statistics:s.statistics()}}(t,e,r,n,i):function(t,e,r,n,i){let o=null;"function"==typeof t.output&&(o=(e,r)=>t.output(e,r)),null===e&&(e=t.initial(r));const s=new a.Dopri((function(e,r,n){t.rhs(e,r,n)}),e.length,i,o);return s.initialise(r,e),{solution:s.run(n),statistics:s.statistics()}}(t,e,r,n,i)}},8714:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.PkgWrapper=void 0;const n=r(8251),i=r(1267);e.PkgWrapper=class{constructor(t,e,r){this.model=new t(n.base,e,r)}initial(t){return this.model.initial(t)}rhs(t,e){const r=new Array(e.length).fill(0);let n=null;if(!(0,i.isODEModel)(this.model))throw Error("Can't use rhs() with delay models");return this.model,this.model.rhs(t,e,r),this.model.output&&(n=this.model.output(t,e)),{output:n,state:r}}getMetadata(){return this.model.getMetadata()}getInternal(){return this.model.getInternal()}setUser(t,e){this.model.setUser(t,e)}run(t,e,r){const n=t[0],o=t[t.length-1],a=(0,i.runModel)(this.model,e,n,o,r);return{names:this.model.names(),statistics:a.statistics,y:a.solution(t)}}}},46:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.interpolatedSolution=e.TimeMode=void 0;const n=r(1729);var i;!function(t){t.Grid="grid",t.Given="given"}(i=e.TimeMode||(e.TimeMode={})),e.interpolatedSolution=function(t,e,r,o){return a=>{const s=function(t,e,r){switch(t.mode){case i.Grid:return(0,n.grid)(Math.max(e,t.tStart),Math.min(r,t.tEnd),t.nPoints);case i.Given:return t.times}}(a,r,o),u=t(s),c=u[0].map(((t,e)=>u.map((t=>t[e]))));return{names:e,x:s,y:c}}}},2630:(t,e)=>{"use strict";function r(t,e,r){if("ignore"===r)return;const n=[];for(const r of Object.keys(t))e.includes(r)||n.push(r);if(n.length>0){const t="Unknown user parameters: "+n.join(", ");if("message"===r)console.log(t);else{if("warning"!==r)throw"stop"===r?Error(t):Error(t+" (and invalid value for unusedUserAction)");console.warn(t)}}}function n(t,e,r,n,i,o,a){const s=t[e];if(void 0===s){if(void 0!==r[e])return;if(null===n)throw Error(`Expected a value for '${e}'`);r[e]=n}else{if("number"!=typeof s)throw Error(`Expected a number for '${e}'`);c(s,i,o,a,e),r[e]=s}}function i(t,e,r,n,i,o,c){let l=t[e];if(void 0===l){if(void 0!==r[e])return;throw Error(`Expected a value for '${e}'`)}{const t=n.length-1;l=a(l),s(t,l,e),function(t,e,r){const n=t.length-1;for(let i=0;i<n;++i){const o=t[i+1];if(e.dim[i]!==o)throw 1===n?Error(`Expected length ${o} value for '${r}'`):Error(`Incorrect size of dimension ${i+1} of '${r}' (expected ${o})`)}}(n,l,e),u(l,i,o,c,e),r[e]=l.data.slice()}}function o(t,e,r,n,i,o,c){let l=t[e];if(void 0===l){if(void 0!==r[e])return;throw Error("Expected a value for '"+e+"'")}{const t=n.length-1;l=a(l),s(t,l,e),u(l,i,o,c,e),n[0]=l.data.length;for(let e=0;e<t;++e)n[e+1]=l.dim[e];r[e]=l.data.slice()}}function a(t,e){return Array.isArray(t)?t={data:t,dim:[t.length]}:"number"==typeof t&&(t={data:[t],dim:[1]}),t}function s(t,e,r){if(e.dim.length!==t)throw 1===t?Error(`Expected a numeric vector for '${r}'`):2===t?Error(`Expected a numeric matrix for '${r}'`):Error(`Expected a numeric array of rank ${t} for '${r}'`)}function u(t,e,r,n,i){for(const o of t.data){if(null===o)throw Error(`'${i}' must not contain any NA values`);c(o,e,r,n,i)}}function c(t,e,r,n,i){if("number"!=typeof t)throw Error(`Expected a number for '${i}'`);if(t<e)throw Error(`Expected '${i}' to be at least ${e}`);if(t>r)throw Error(`Expected '${i}' to be at most ${r}`);if(n&&!Number.isInteger(t))throw Error(`Expected '${i}' to be integer-like`)}Object.defineProperty(e,"__esModule",{value:!0}),e.user=e.setUserArrayVariable=e.setUserArrayFixed=e.setUserScalar=e.checkUser=void 0,e.checkUser=r,e.setUserScalar=n,e.setUserArrayFixed=i,e.setUserArrayVariable=o,e.user={checkUser:r,setUserArrayFixed:i,setUserArrayVariable:o,setUserScalar:n}},1729:(t,e)=>{"use strict";function r(t,e,r){const n=(e-t)/(r-1),i=[];for(let e=0;e<r-1;++e)i.push(t+e*n);return i.push(e),i}Object.defineProperty(e,"__esModule",{value:!0}),e.loop=e.whichMax=e.whichMin=e.gridLog=e.grid=void 0,e.grid=r,e.gridLog=function(t,e,n){return r(Math.log(t),Math.log(e),n).map(Math.exp)},e.whichMin=function(t){let e=-1,r=1/0;for(let n=0;n<t.length;++n)t[n]<r&&(e=n,r=t[n]);return e},e.whichMax=function(t){let e=-1,r=-1/0;for(let n=0;n<t.length;++n)t[n]>r&&(e=n,r=t[n]);return e},e.loop=function(t,e){const r=[];for(let n=0;n<t;++n)r.push(e(n));return r}},7497:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.versions=void 0,e.versions=function(){return{dfoptim:"0.0.5",dopri:"0.0.12",odinjs:"0.0.15"}}},7466:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.wodinFitValue=e.wodinFit=e.wodinRun=void 0;const n=r(4153),i=r(8251),o=r(7771),a=r(1267),s=r(46);e.wodinRun=function(t,e,r,n,o){const u=new t(i.base,e,"error"),c=(0,a.runModel)(u,null,r,n,o).solution,l=u.names();return(0,s.interpolatedSolution)(c,l,r,n)},e.wodinFit=function(t,e,r,i,a,s){const u=(0,o.fitTarget)(t,e,r,i,a),c=r.vary.map((t=>r.base[t]));return new n.Simplex(u,c,s)},e.wodinFitValue=function(t,e,r){const{names:n,y:i}=t({mode:s.TimeMode.Given,times:e.time}),a=n.indexOf(r);return(0,o.sumOfSquares)(e.value,i[a])}},7689:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.binomial=void 0,e.binomial=function(t,e,r){if(0===e||0===r)return 0;if(1===r)return e;const o=r<=.5?r:1-r,a=e*o>=10?function(t,e,r){const i=Math.sqrt(e*r*(1-r)),o=1.15+2.53*i,a=.0248*o-.0873+.01*r,s=e*r+.5,u=.92-4.2/o,c=r/(1-r),l=(2.83+5.1/o)*i,f=Math.floor((e+1)*r);let h;for(;;){const r=t.random()-.5;let i=t.random();const p=.5-Math.abs(r),d=Math.floor((2*a/p+o)*r+s);if(p>=.07&&i<=u){h=d;break}if(!(d<0||d>e)&&(i=Math.log(i*l/(a/(p*p)+o)),i<=(f+.5)*Math.log((f+1)/(c*(e-f+1)))+(e+1)*Math.log((e-f+1)/(e-d+1))+(d+.5)*Math.log(c*(e-d+1)/(d+1))+n(f)+n(e-f)-n(d)-n(e-d))){h=d;break}}return h}(t,e,o):function(t,e,r){let n=-1;do{n=i(t.random(),e,r)}while(n<0);return n}(t,e,o);return r<=.5?a:e-a};const r=[.08106146679532733,.04134069595540946,.02767792568499805,.02079067210376495,.016644691189821703,.013876128823071987,.011896709945892425,.010411265261971892,.00925546218270945,.008330563433359472];function n(t){return t<r.length?r[t]:function(t){const e=(t+1)*(t+1);return(1/12-(1/360-1/1260/e)/e)/(t+1)}(t)}function i(t,e,r){const n=1-r,i=r/n,o=i*(e+1);let a=Math.pow(n,e),s=0,u=a;for(;t>=a;)t-=a,s++,a*=o/s-i,u=a;return s}},9386:(t,e)=>{"use strict";function r(t){return-Math.log(t.random())}Object.defineProperty(e,"__esModule",{value:!0}),e.exponential=e.randomExponential=void 0,e.randomExponential=r,e.exponential=function(t,e){return r(t)/e}},6532:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.uniform=e.randomUniform=e.randomNormal=e.normal=e.randomExponential=e.exponential=e.binomial=e.RngStateReplay=e.RngStateObserved=e.RngStateBuiltin=e.RngState=e.Random=void 0;var n=r(1929);Object.defineProperty(e,"Random",{enumerable:!0,get:function(){return n.Random}});var i=r(3726);Object.defineProperty(e,"RngState",{enumerable:!0,get:function(){return i.RngState}});var o=r(2035);Object.defineProperty(e,"RngStateBuiltin",{enumerable:!0,get:function(){return o.RngStateBuiltin}});var a=r(2301);Object.defineProperty(e,"RngStateObserved",{enumerable:!0,get:function(){return a.RngStateObserved}});var s=r(5843);Object.defineProperty(e,"RngStateReplay",{enumerable:!0,get:function(){return s.RngStateReplay}});var u=r(7689);Object.defineProperty(e,"binomial",{enumerable:!0,get:function(){return u.binomial}});var c=r(9386);Object.defineProperty(e,"exponential",{enumerable:!0,get:function(){return c.exponential}}),Object.defineProperty(e,"randomExponential",{enumerable:!0,get:function(){return c.randomExponential}});var l=r(2638);Object.defineProperty(e,"normal",{enumerable:!0,get:function(){return l.normal}}),Object.defineProperty(e,"randomNormal",{enumerable:!0,get:function(){return l.randomNormal}});var f=r(6390);Object.defineProperty(e,"randomUniform",{enumerable:!0,get:function(){return f.randomUniform}}),Object.defineProperty(e,"uniform",{enumerable:!0,get:function(){return f.uniform}})},2638:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.normal=e.randomNormal=void 0;const r=2*Math.PI;function n(t){let e=0,n=0;do{e=t.random(),n=t.random()}while(e<Number.EPSILON);return Math.sqrt(-2*Math.log(e))*Math.cos(r*n)}e.randomNormal=n,e.normal=function(t,e,r){return e+n(t)*r}},5789:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.poisson=void 0;const n=r(7788);e.poisson=function(t,e){return e<10?function(t,e){let r=0;const n=Math.exp(-e);let o=1;for(;o*=t.random(),!(o<=n&&r<=i);)r++;return r}(t,e):function(t,e){let r=0;const i=Math.log(e),o=.931+2.53*Math.sqrt(e),a=.02483*o-.059,s=1.1239+1.1328/(o-3.4);for(;;){let u=t.random();u-=.5;const c=t.random(),l=.5-Math.abs(u),f=Math.floor((2*a/l+o)*u+e+.43);if(l>=.07&&c<=.9277-3.6224/(o-2)){r=f;break}if(!(f<0||l<.013&&c>l)&&Math.log(c*s/(a/(l*l)+o))<=f*i-e-(0,n.lngamma)(f+1)){r=f;break}}return r}(t,e)};const i=2**31-1},1929:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.Random=void 0;const n=r(7689),i=r(9386),o=r(2638),a=r(5789),s=r(6390);e.Random=class{state;constructor(t){this.state=t}randomExponential(){return(0,i.randomExponential)(this.state)}randomNormal(){return(0,o.randomNormal)(this.state)}randomUniform(){return(0,s.randomUniform)(this.state)}binomial(t,e){return(0,n.binomial)(this.state,t,e)}exponential(t){return(0,i.exponential)(this.state,t)}normal(t,e){return(0,o.normal)(this.state,t,e)}poisson(t){return(0,a.poisson)(this.state,t)}uniform(t,e){return(0,s.uniform)(this.state,t,e)}}},2035:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.RngStateBuiltin=void 0;const n=r(3726);class i extends n.RngState{constructor(){super()}random(){return Math.random()}}e.RngStateBuiltin=i},2301:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.RngStateObserved=void 0;const n=r(3726),i=r(2035),o=r(5843);class a extends n.RngState{_state;_values;constructor(t){super(),void 0===t&&(t=new i.RngStateBuiltin),this._state=t,this._values=[]}random(){const t=this._state.random();return this._values.push(t),t}replay(){return new o.RngStateReplay(this._values)}length(){return this._values.length}values(){return this._values}}e.RngStateObserved=a},5843:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.RngStateReplay=void 0;const n=r(3726);class i extends n.RngState{_i;_values;constructor(t){super(),this._i=0,this._values=t}random(){const t=this._values[this._i++];if(this._i>this._values.length)throw Error("Exhausted stream");return t}length(){return this._values.length-this._i}}e.RngStateReplay=i},3726:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.RngState=void 0,e.RngState=class{setSeed(t){throw Error("This generator cannot be seeded (really!)")}}},6390:(t,e)=>{"use strict";function r(t){return t.random()}Object.defineProperty(e,"__esModule",{value:!0}),e.uniform=e.randomUniform=void 0,e.randomUniform=r,e.uniform=function(t,e,n){return r(t)*(n-e)+e}},4394:(t,e)=>{"use strict";function r(t,e){return void 0===t?e:t}Object.defineProperty(e,"__esModule",{value:!0}),e.simplexControl=void 0,e.simplexControl=function(t){void 0===t&&(t={});var e={deltaNonZero:r(t.deltaNonZero,.05),deltaZero:r(t.deltaZero,.001),errorOnFailure:r(t.errorOnFailure,!1),tolerance:r(t.tolerance,1e-5)};if(e.tolerance<=0)throw new Error("Invalid control parameter: 'tolerance' must be strictly positive");return e}},4153:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0});var n=r(2169);Object.defineProperty(e,"Simplex",{enumerable:!0,get:function(){return n.Simplex}}),Object.defineProperty(e,"fitSimplex",{enumerable:!0,get:function(){return n.fitSimplex}})},2169:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.Simplex=e.fitSimplex=void 0;var n=r(4394),i=r(2003);function o(t,e,r){for(var n=e.length,i=new Array(n),o=1+t,a=-t,s=0;s<i.length;++s)i[s]=o*e[s]+a*r[s];return i}e.fitSimplex=function(t,e,r,n){return void 0===r&&(r={}),void 0===n&&(n=200),new a(t,e,r).run(n*e.length)};var a=function(){function t(t,e,r){void 0===r&&(r={}),this.rho=1,this.chi=2,this.psi=-.5,this.sigma=.5,this._iterations=0,this._id=0,this._converged=!1,this._target=t,this._control=n.simplexControl(r),this._n=e.length,this._simplex=[],this._simplex.push(this._point(e.slice())),r.errorOnFailure||(this._target=function(t){return function(e){try{return t(e)}catch(t){return 1/0}}}(t));for(var i=0;i<this._n;++i){var o=e.slice();o[i]?o[i]*=1+this._control.deltaNonZero:o[i]=this._control.deltaZero,this._simplex.push(this._point(o))}this._sort()}return t.prototype.step=function(){if(this._iterations++,this._isConverged())return this._converged=!0,!0;var t=this._n,e=this._simplex[0],r=this._simplex[t],n=this._centroid(),i=this._reflect(r,n);if(i.value<e.value){var o=this._expand(r,n);this._update(o.value<i.value?o:i)}else if(i.value>=this._simplex[t-1].value){var a=i.value>r.value?this._contractInside(r,n):this._contractOutside(r,n);if(a.value<r.value)this._update(a);else{for(var s=1;s<=t;++s)this._simplex[s]=this._shrink(this._simplex[s],e.location);this._sort()}}else this._update(i);return!1},t.prototype.run=function(t){if(!this._converged)for(var e=0;e<t&&!this.step();++e);return this.result()},t.prototype.result=function(){var t=this._simplex[0];return{converged:this._converged,data:t.data,evaluations:this._id,iterations:this._iterations,location:t.location,value:t.value}},t.prototype.simplex=function(){return this._simplex.map((function(t){return{location:t.location,value:t.value}}))},t.prototype._point=function(t){var e=i.checkResult(this._target(t));return{data:e.data,id:this._id++,location:t,value:e.value}},t.prototype._sort=function(){this._simplex.sort((function(t,e){return t.value-e.value}))},t.prototype._centroid=function(){for(var t=this._n,e=new Array(t),r=0;r<t;++r){e[r]=0;for(var n=0;n<t;++n)e[r]+=this._simplex[n].location[r];e[r]/=t}return e},t.prototype._update=function(t){this._simplex[this._n]=t,this._sort()},t.prototype._reflect=function(t,e){return this._point(o(this.rho,e,t.location))},t.prototype._expand=function(t,e){return this._point(o(this.chi,e,t.location))},t.prototype._contractInside=function(t,e){return this._point(o(this.psi,e,t.location))},t.prototype._contractOutside=function(t,e){return this._point(o(-this.psi*this.rho,e,t.location))},t.prototype._shrink=function(t,e){return this._point(o(-this.sigma,e,t.location))},t.prototype._isConverged=function(){for(var t=this._control.tolerance,e=this._simplex[0],r=this._simplex[this._n],n=r.value-e.value<t||1-e.value/r.value<t,i=0,o=0;o<this._n;++o)i=Math.max(i,Math.abs(e.location[o]-r.location[o]));return i<t&&n},t}();e.Simplex=a},2003:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.checkResult=void 0,e.checkResult=function(t){return"number"==typeof t?{data:null,value:t}:t}},2886:(t,e)=>{"use strict";function r(t,e){return new Error("Invalid control parameter: '".concat(t,"' ").concat(e))}function n(t,e){return void 0===t?e:t}Object.defineProperty(e,"__esModule",{value:!0}),e.dopriControl=void 0,e.dopriControl=function(t){void 0===t&&(t={});var e=1e-6,i=1e4,o=1e-6,a=1/0,s=1e-8,u=!1,c=0,l=1/0,f={atol:n(t.atol,e),maxSteps:n(t.maxSteps,i),rtol:n(t.rtol,o),stepSizeMax:n(t.stepSizeMax,a),stepSizeMin:n(t.stepSizeMin,s),stepSizeMinAllow:n(t.stepSizeMinAllow,u),stiffCheck:n(t.stiffCheck,c),tcrit:n(t.tcrit,l)};if(f.maxSteps<1)throw r("maxSteps","must be at least 1");if(f.atol<=0)throw r("atol","must be strictly positive");if(f.rtol<=0)throw r("rtol","must be strictly positive");return f}},8742:function(t,e,r){"use strict";var n,i=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.DDE=e.integrateDDE=void 0;var o=r(2223),a=r(4660);e.integrateDDE=function(t,e,r,n,i,o){void 0===i&&(i={}),void 0===o&&(o=null);var a=new s(t,e.length,i,o);return a.initialise(r,e),a.run(n)};var s=function(t){function e(e,r,n,i){void 0===n&&(n={}),void 0===i&&(i=null);var o=this,a=function(t){return o._interpolate(t)},s=null===i?null:function(t,e){return i(t,e,a)};return(o=t.call(this,(function(t,r,n){return e(t,r,n,a)}),r,n,s)||this)._y0=new Array(r),o}return i(e,t),e.prototype.initialise=function(e,r){return this._y0=r,t.prototype.initialise.call(this,e,r),this},e.prototype._interpolate=function(t){var e=this._findHistory(t);return e<0?this._y0.slice(0):this._stepper.interpolate(t,this._history[e])},e.prototype._findHistory=function(t){return(0,a.search)(this._history,(function(e){return e.t>t}))},e}(o.Dopri);e.DDE=s},2223:function(t,e,r){"use strict";var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var i=Object.getOwnPropertyDescriptor(e,r);i&&!("get"in i?!e.__esModule:i.writable||i.configurable)||(i={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,i)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),i=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),o=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return i(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri=e.integrateDopri=void 0;var a=r(2886),s=o(r(2922)),u=r(1216),c=o(r(4660)),l=Math.pow(2,-52);function f(t,e){return new Error("Integration failure: ".concat(t," at ").concat(e))}e.integrateDopri=function(t,e,r,n,i,o){void 0===i&&(i={}),void 0===o&&(o=null);var a=new h(t,e.length,i,o);return a.initialise(r,e),a.run(n)};var h=function(){function t(t,e,r,n){void 0===r&&(r={}),void 0===n&&(n=null),this._history=[],this._t=0,this._h=0,this._nSteps=0,this._nStepsAccepted=0,this._nStepsRejected=0,this._stiffNStiff=0,this._stiffNNonstiff=0,this._lastError=0,this._stepper=new s.Dopri5(t,e),this._control=(0,a.dopriControl)(r),this._output=n}return t.prototype.initialise=function(t,e){var r=this._stepper.n;if(e.length!==r)throw Error("Invalid size 'y' - expected a length ".concat(r," array"));return this._stepper.reset(t,e),this._reset(),this._h=function(t,e,r,n,i,o){var a=new Array(t.n),s=new Array(t.n),u=new Array(t.n);t.rhs(e,r,a),t.nEval++;var l=0,f=0,h=0;for(h=0;h<t.n;++h){var p=n+i*Math.abs(r[h]);l+=c.square(a[h]/p),f+=c.square(r[h]/p)}var d=l<=1e-10||f<=1e-10?1e-6:.01*Math.sqrt(f/l);for(d=Math.min(d,o),h=0;h<t.n;++h)u[h]=r[h]+d*a[h];t.rhs(e+d,u,s),t.nEval++;var v=0;for(h=0;h<t.n;++h)p=n+i*Math.abs(r[h]),v+=c.square((s[h]-a[h])/p);v=Math.sqrt(v)/d;var y=Math.max(Math.abs(v),Math.sqrt(l)),m=y<=1e-15?Math.max(1e-6,.001*Math.abs(d)):Math.pow(.01/y,1/t.order);return Math.min(Math.min(100*Math.abs(d),m),o)}(this._stepper,t,e,this._control.atol,this._control.rtol,this._control.stepSizeMax),this._t=t,this._history=[],this},t.prototype.run=function(t){for(;this._t<t;)this._step(),this._history.push(this._stepper.history.clone());return(0,u.interpolator)(this._history.slice(0),this._stepper,this._output)},t.prototype.statistics=function(){return{lastError:this._lastError,nEval:this._stepper.nEval,nSteps:this._nSteps,nStepsAccepted:this._nStepsAccepted,nStepsRejected:this._nStepsRejected,stiffNNonstiff:this._stiffNNonstiff,stiffNStiff:this._stiffNStiff}},t.prototype._reset=function(){this._nSteps=0,this._nStepsAccepted=0,this._nStepsRejected=0,this._stiffNStiff=0,this._stiffNNonstiff=0,this._lastError=0},t.prototype._step=function(){for(var t=this._t,e=this._h,r=!1,n=!1,i=Math.max(this._lastError,1e-4),o=this._stepper.stepControl,a=this._control;!r;){if(this._nSteps>a.maxSteps)throw f("too many steps",t);if(e<a.stepSizeMin){if(!a.stepSizeMinAllow)throw f("step too small",t);e=a.stepSizeMin}if(e<=Math.abs(t)*l)throw f("step size vanished",t);t+e>a.tcrit&&(e=a.tcrit-t),this._stepper.step(t,e),this._nSteps++;var s=this._stepper.error(a.atol,a.rtol),u=Math.pow(s,o.constant),h=1/o.factorMin,p=1/o.factorMax;if(s<=1){if(r=!0,this._nStepsAccepted++,this._isStiff(e))throw f("problem became stiff",t);this._stepper.stepComplete(t,e);var d=u/Math.pow(i,o.beta),v=e/(d=c.constrain(d/o.factorSafe,p,h));this._t+=e,this._h=n?Math.min(v,e):Math.min(v,a.stepSizeMax),this._lastError=s}else n=!0,this._nStepsAccepted>=1&&this._nStepsRejected++,e/=Math.min(h,u/o.factorSafe)}return this._t},t.prototype._isStiff=function(t){if(this._stiffNStiff>0||this._nStepsAccepted%this._control.stiffCheck==0)if(this._stepper.isStiff(t)){if(this._stiffNNonstiff=0,this._stiffNStiff++>=15)return!0}else this._stiffNStiff>0&&this._stiffNNonstiff++>=6&&(this._stiffNStiff=0,this._stiffNNonstiff=0);return!1},t}();e.Dopri=h},7134:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri5StepControl=void 0;e.Dopri5StepControl=function(){this.factorSafe=.9,this.factorMin=.2,this.factorMax=10,this.beta=.04,this.constant=.2-.75*this.beta}},2922:function(t,e,r){"use strict";var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var i=Object.getOwnPropertyDescriptor(e,r);i&&!("get"in i?!e.__esModule:i.writable||i.configurable)||(i={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,i)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),i=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),o=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return i(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri5=void 0;var a=o(r(4660)),s=o(r(7134)),u=r(8430),c=function(){function t(t,e){this.order=5,this.stepControl=new s.Dopri5StepControl,this.nEval=0,this.rhs=t,this.n=e,this.y=new Array(e),this.yNext=new Array(e),this.yStiff=new Array(e),this.k1=new Array(e),this.k2=new Array(e),this.k3=new Array(e),this.k4=new Array(e),this.k5=new Array(e),this.k6=new Array(e),this.history=new u.HistoryElement(this.order*e)}return t.prototype.step=function(t,e){var r=this.n,n=this.y,i=this.yNext,o=this.k1,a=this.k2,s=this.k3,u=this.k4,c=this.k5,l=this.k6,f=this.history.data,h=0;for(h=0;h<r;++h)i[h]=n[h]+.2*e*o[h];for(this.rhs(t+.2*e,i,a),h=0;h<r;++h)i[h]=n[h]+e*(.075*o[h]+.225*a[h]);for(this.rhs(t+.3*e,i,s),h=0;h<r;++h)i[h]=n[h]+e*(.9777777777777777*o[h]+-3.7333333333333334*a[h]+3.5555555555555554*s[h]);for(this.rhs(t+.8*e,i,u),h=0;h<r;++h)i[h]=n[h]+e*(2.9525986892242035*o[h]+-11.595793324188385*a[h]+9.822892851699436*s[h]+-.2908093278463649*u[h]);for(this.rhs(t+.8888888888888888*e,i,c),h=0;h<r;++h)this.yStiff[h]=n[h]+e*(2.8462752525252526*o[h]+-10.757575757575758*a[h]+8.906422717743473*s[h]+.2784090909090909*u[h]+-.2735313036020583*c[h]);var p=t+e;for(this.rhs(p,this.yStiff,l),h=0;h<r;++h)i[h]=n[h]+e*(.09114583333333333*o[h]+.44923629829290207*s[h]+.6510416666666666*u[h]+-.322376179245283*c[h]+.13095238095238096*l[h]);this.rhs(p,i,a);var d=4*r;for(h=0;h<r;++h)f[d++]=e*(-1.1270175653862835*o[h]+2.675424484351598*s[h]+-5.685526961588504*u[h]+3.5219323679207912*c[h]+-1.7672812570757455*l[h]+2.382468931778144*a[h]);for(h=0;h<r;++h)u[h]=e*(.0012326388888888888*o[h]+-.0042527702905061394*s[h]+.03697916666666667*u[h]+-.05086379716981132*c[h]+.0419047619047619*l[h]+-.025*a[h]);this.nEval+=6},t.prototype.stepComplete=function(t,e){this.saveHistory(t,e),a.copyArray(this.k1,this.k2),a.copyArray(this.y,this.yNext)},t.prototype.error=function(t,e){var r=0,n=0;for(n=0;n<this.n;++n){var i=t+e*Math.max(Math.abs(this.y[n]),Math.abs(this.yNext[n]));r+=a.square(this.k4[n]/i)}return Math.sqrt(r/this.n)},t.prototype.interpolate=function(t,e){for(var r=e.data,n=(t-e.t)/e.h,i=1-n,o=this.n,a=new Array(o),s=0;s<o;++s)a[s]=r[s]+n*(r[o+s]+i*(r[2*o+s]+n*(r[3*o+s]+i*r[4*o+s])));return a},t.prototype.isStiff=function(t){for(var e=0,r=0,n=0;n<this.n;++n)e+=a.square(this.k2[n]-this.k6[n]),r+=a.square(this.yNext[n]-this.yStiff[n]);return r>0&&Math.abs(t)*Math.sqrt(e/r)>3.25},t.prototype.reset=function(t,e){for(var r=0;r<this.n;++r)this.y[r]=e[r];this.rhs(t,e,this.k1),this.nEval=1},t.prototype.saveHistory=function(t,e){for(var r=this.history,n=this.n,i=0;i<n;++i){var o=this.yNext[i]-this.y[i],a=e*this.k1[i]-o;r.data[i]=this.y[i],r.data[n+i]=o,r.data[2*n+i]=a,r.data[3*n+i]=-e*this.k2[i]+o-a}r.t=t,r.h=e},t}();e.Dopri5=c},1606:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.HistoryElement=void 0;var r=function(){function t(t){this.t=0,this.h=0,this.data=new Array(t)}return t.prototype.clone=function(){var e=new t(this.data.length);return e.t=this.t,e.h=this.h,e.data=this.data.slice(0),e},t}();e.HistoryElement=r},6245:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.integrate=e.integrateDDE=e.DDE=e.integrateDopri=e.Dopri=void 0;var n=r(2223);Object.defineProperty(e,"Dopri",{enumerable:!0,get:function(){return n.Dopri}}),Object.defineProperty(e,"integrateDopri",{enumerable:!0,get:function(){return n.integrateDopri}});var i=r(8742);Object.defineProperty(e,"DDE",{enumerable:!0,get:function(){return i.DDE}}),Object.defineProperty(e,"integrateDDE",{enumerable:!0,get:function(){return i.integrateDDE}});var o=r(5961);Object.defineProperty(e,"integrate",{enumerable:!0,get:function(){return o.integrate}})},5961:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.integrate=void 0;var n=r(8742),i=r(2223);e.integrate=function(t,e,r,o,a,s){if(void 0===a&&(a={}),void 0===s&&(s=null),function(t){return 3===t.length}(t)){if(!function(t){return null===t||2===t.length}(s))throw new Error("Can't used delayed output with non-delayed rhs");return(0,i.integrateDopri)(t,e,r,o,a,s)}if(!function(t){return null===t||3===t.length}(s))throw new Error("Can't used non-delayed output with delayed rhs");return(0,n.integrateDDE)(t,e,r,o,a,s)}},1216:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.interpolator=void 0,e.interpolator=function(t,e,r){return function(n){return function(t,e,r,n){for(var i=[],o=e,a=0,s=0,u=t;s<u.length;s++){for(var c=u[s];o[a].t+o[a].h<c;)a++;var l=r.interpolate(c,o[a]);null!==n&&(l=l.concat(n(c,l))),i.push(l)}return i}(n,t,e,r)}}},8430:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.HistoryElement=void 0;var n=r(1606);Object.defineProperty(e,"HistoryElement",{enumerable:!0,get:function(){return n.HistoryElement}})},4660:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.search=e.seqLen=e.approxEqualArray=e.approxEqual=e.zeros=e.copyArray=e.constrain=e.square=void 0;var r=Math.pow(2,-26);e.square=function(t){return t*t},e.constrain=function(t,e,r){return Math.max(Math.min(t,r),e)},e.copyArray=function(t,e){for(var r=t.length,n=0;n<r;n++)t[n]=e[n]},e.zeros=function(t){for(var e=Array(t),r=0;r<t;++r)e[r]=0;return e},e.approxEqual=function(t,e,n){void 0===n&&(n=r);var i=Math.abs(t-e),o=Math.abs(t);return o>n&&(i/=o),i<n},e.approxEqualArray=function(t,e,n){if(void 0===n&&(n=r),e.length!==t.length)throw Error("Incompatible arrays");for(var i=0,o=0,a=0,s=0;s<t.length;++s)t[s]!==e[s]&&(i+=Math.abs(t[s]),o+=Math.abs(t[s]-e[s]),a++);return 0===a||(o/=a,(i/=a)>n&&(o/=i),o<n)},e.seqLen=function(t,e,r){for(var n=[],i=0;i<r;++i){var o=i/(r-1);n.push((1-o)*t+o*e)}return n},e.search=function(t,e){var r=0,n=t.length-1;if(0===t.length||e(t[r]))return-1;if(!e(t[n]))return n;for(;n-r>1;){var i=Math.floor((r+n)/2);e(t[i])?n=i:r=i}return r}},7788:(t,e)=>{var r,n=function(){"use strict";var t=new DataView(new ArrayBuffer(8)),e=function(e,r,n){var i,o,a,s=-.16666666666666632,u=0;return t.setFloat64(0,e),u=t.getInt32(0),(u&=2147483647)<1044381696&&0==(0|e)?e:(a=(i=e*e)*e,o=.00833333333332249+i*(i*(27557313707070068e-22+i*(1.58969099521155e-10*i-2.5050760253406863e-8))-.0001984126982985795),0===n?e+a*(s+i*o):e-(i*(.5*r-a*o)-r-a*s))},r=function(e,r){var n,i,o,a=0;return t.setFloat64(0,e),a=t.getInt32(0),(a&=2147483647)<1044381696&&0==(0|e)?1:(i=(n=e*e)*(.0416666666666666+n*(n*(2480158728947673e-20+n*(n*(2.087572321298175e-9+-11359647557788195e-27*n)-2.7557314351390663e-7))-.001388888888887411)),a<1070805811?1-(.5*n-(n*i-e*r)):(a>1072234496?o=.28125:(t.setFloat64(0,o),t.setUint32(0,a-2097152),t.setUint32(4,0),o=t.getFloat64(0)),1-o-(.5*n-o-(n*i-e*r))))},n=function(n){var i,o,a,s,u,c,l,f,h=1.4616321449683622,p=0,d=0;if(t.setFloat64(0,n),c=t.getInt32(0),l=t.getInt32(4),(f=2147483647&c)>=2146435072)return n*n;if(0==(f|l))return 1/0;if(f<999292928)return c<0?-Math.log(-n):-Math.log(n);if(c<0){if(f>=1127219200)return 1/0;if(i=function(n){var i,o=0,a=0;if(t.setFloat64(0,n),a=t.getInt32(0),(a&=2147483647)<1070596096)return e(Math.PI*n,0,0);switch(i=-n,Math.floor(i)!=i?o=4*(i=2*((i*=.5)-Math.floor(i)))|0:a>=1128267776?(i=0,o=0):(o=t.getInt32(4),i=o&=1,o<<=2),o){case 0:i=e(Math.PI*i,0,0);break;case 1:case 2:i=r(Math.PI*(.5-i),0);break;case 3:case 4:i=e(Math.PI*(1-i),0,0);break;case 5:case 6:i=-r(Math.PI*(i-1.5),0);break;default:i=e(Math.PI*(i-2),0,0)}return-i}(n),0==i)return 1/0;p=Math.log(3.141592653589793/Math.abs(i*n)),n=-n}if(0==(f-1072693248|l)||0==(f-1073741824|l))s=0;else if(f<1073741824)switch(f<=1072483532?(s=-Math.log(n),f>=1072130372?(o=1-n,d=0):f>=1070442081?(o=n-(h-1),d=1):(o=n,d=2)):(s=0,f>=1073460419?(o=2-n,d=0):f>=1072936132?(o=n-h,d=1):(o=n-1,d=2)),d){case 0:s+=o*(.07721566490153287+(a=o*o)*(.06735230105312927+a*(.007385550860814029+a*(.0011927076318336207+a*(.00022086279071390839+25214456545125733e-21*a)))))+a*(.3224670334241136+a*(.020580808432516733+a*(.0028905138367341563+a*(.0005100697921535113+a*(.00010801156724758394+44864094961891516e-21*a)))))-.5*o;break;case 1:s+=(a=o*o)*(.48383612272381005+(u=a*o)*(u*(.006100538702462913+u*(.00031563207090362595*u-.0014034646998923284))-.032788541075985965))-(-3638676997039505e-33-u*(u*(.01797067508118204+u*(u*(.000881081882437654+-.00031275416837512086*u)-.0036845201678113826))-.1475877229945939+o*(.06462494023913339+u*(u*(.0022596478090061247+u*(.0003355291926355191*u-.0005385953053567405))-.010314224129834144))))-.12148629053584961;break;case 2:s+=-.5*o+o*(o*(.6328270640250934+o*(1.4549225013723477+o*(.9777175279633727+o*(.22896372806469245+.013381091853678766*o))))-.07721566490153287)/(1+o*(2.4559779371304113+o*(2.128489763798934+o*(.7692851504566728+o*(.10422264559336913+.003217092422824239*o)))))}else if(f<1075838976)switch(i=0,s=.5*(o=n-(d=0|n))+o*(o*(.21498241596060885+o*(.325778796408931+o*(.14635047265246445+o*(.02664227030336386+o*(.0018402845140733772+3194753265841009e-20*o)))))-.07721566490153287)/(1+o*(1.3920053346762105+o*(.7219355475671381+o*(.17193386563280308+o*(.01864591917156529+o*(.0007779424963818936+7326684307446256e-21*o)))))),a=1,d){case 7:a*=o+6;case 6:a*=o+5;case 5:a*=o+4;case 4:a*=o+3;case 3:a*=o+2,s+=Math.log(a)}else s=f<1133510656?(n-.5)*((i=Math.log(n))-1)+(u=.4189385332046727+(a=1/n)*(.08333333333333297+(o=a*a)*(o*(.0007936505586430196+o*(o*(.0008363399189962821+-.0016309293409657527*o)-.00059518755745034))-.0027777777772877554))):n*(Math.log(n)-1);return c<0&&(s=p-s),s};return{lngamma:function(t){return n(t)},gamma:function(t){return Math.exp(n(t))}}}();void 0!==t.exports?t.exports=n:void 0===(r=function(){return n}.apply(e,[]))||(t.exports=r)},6907:t=>{"use strict";t.exports=function(t){for(var e=new Array(t),r=0;r<t;++r)e[r]=r;return e}},8738:t=>{function e(t){return!!t.constructor&&"function"==typeof t.constructor.isBuffer&&t.constructor.isBuffer(t)}t.exports=function(t){return null!=t&&(e(t)||function(t){return"function"==typeof t.readFloatLE&&"function"==typeof t.slice&&e(t.slice(0,0))}(t)||!!t._isBuffer)}},2861:(t,e,r)=>{var n=r(6907),i=r(8738),o="undefined"!=typeof Float64Array;function a(t,e){return t[0]-e[0]}function s(){var t,e=this.stride,r=new Array(e.length);for(t=0;t<r.length;++t)r[t]=[Math.abs(e[t]),t];r.sort(a);var n=new Array(r.length);for(t=0;t<n.length;++t)n[t]=r[t][1];return n}function u(t,e){var r=["View",e,"d",t].join("");e<0&&(r="View_Nil"+t);var i="generic"===t;if(-1===e){var o="function "+r+"(a){this.data=a;};var proto="+r+".prototype;proto.dtype='"+t+"';proto.index=function(){return -1};proto.size=0;proto.dimension=-1;proto.shape=proto.stride=proto.order=[];proto.lo=proto.hi=proto.transpose=proto.step=function(){return new "+r+"(this.data);};proto.get=proto.set=function(){};proto.pick=function(){return null};return function construct_"+r+"(a){return new "+r+"(a);}";return new Function(o)()}if(0===e)return o="function "+r+"(a,d) {this.data = a;this.offset = d};var proto="+r+".prototype;proto.dtype='"+t+"';proto.index=function(){return this.offset};proto.dimension=0;proto.size=1;proto.shape=proto.stride=proto.order=[];proto.lo=proto.hi=proto.transpose=proto.step=function "+r+"_copy() {return new "+r+"(this.data,this.offset)};proto.pick=function "+r+"_pick(){return TrivialArray(this.data);};proto.valueOf=proto.get=function "+r+"_get(){return "+(i?"this.data.get(this.offset)":"this.data[this.offset]")+"};proto.set=function "+r+"_set(v){return "+(i?"this.data.set(this.offset,v)":"this.data[this.offset]=v")+"};return function construct_"+r+"(a,b,c,d){return new "+r+"(a,d)}",new Function("TrivialArray",o)(c[t][0]);o=["'use strict'"];var a=n(e),u=a.map((function(t){return"i"+t})),l="this.offset+"+a.map((function(t){return"this.stride["+t+"]*i"+t})).join("+"),f=a.map((function(t){return"b"+t})).join(","),h=a.map((function(t){return"c"+t})).join(",");o.push("function "+r+"(a,"+f+","+h+",d){this.data=a","this.shape=["+f+"]","this.stride=["+h+"]","this.offset=d|0}","var proto="+r+".prototype","proto.dtype='"+t+"'","proto.dimension="+e),o.push("Object.defineProperty(proto,'size',{get:function "+r+"_size(){return "+a.map((function(t){return"this.shape["+t+"]"})).join("*"),"}})"),1===e?o.push("proto.order=[0]"):(o.push("Object.defineProperty(proto,'order',{get:"),e<4?(o.push("function "+r+"_order(){"),2===e?o.push("return (Math.abs(this.stride[0])>Math.abs(this.stride[1]))?[1,0]:[0,1]}})"):3===e&&o.push("var s0=Math.abs(this.stride[0]),s1=Math.abs(this.stride[1]),s2=Math.abs(this.stride[2]);if(s0>s1){if(s1>s2){return [2,1,0];}else if(s0>s2){return [1,2,0];}else{return [1,0,2];}}else if(s0>s2){return [2,0,1];}else if(s2>s1){return [0,1,2];}else{return [0,2,1];}}})")):o.push("ORDER})")),o.push("proto.set=function "+r+"_set("+u.join(",")+",v){"),i?o.push("return this.data.set("+l+",v)}"):o.push("return this.data["+l+"]=v}"),o.push("proto.get=function "+r+"_get("+u.join(",")+"){"),i?o.push("return this.data.get("+l+")}"):o.push("return this.data["+l+"]}"),o.push("proto.index=function "+r+"_index(",u.join(),"){return "+l+"}"),o.push("proto.hi=function "+r+"_hi("+u.join(",")+"){return new "+r+"(this.data,"+a.map((function(t){return["(typeof i",t,"!=='number'||i",t,"<0)?this.shape[",t,"]:i",t,"|0"].join("")})).join(",")+","+a.map((function(t){return"this.stride["+t+"]"})).join(",")+",this.offset)}");var p=a.map((function(t){return"a"+t+"=this.shape["+t+"]"})),d=a.map((function(t){return"c"+t+"=this.stride["+t+"]"}));o.push("proto.lo=function "+r+"_lo("+u.join(",")+"){var b=this.offset,d=0,"+p.join(",")+","+d.join(","));for(var v=0;v<e;++v)o.push("if(typeof i"+v+"==='number'&&i"+v+">=0){d=i"+v+"|0;b+=c"+v+"*d;a"+v+"-=d}");for(o.push("return new "+r+"(this.data,"+a.map((function(t){return"a"+t})).join(",")+","+a.map((function(t){return"c"+t})).join(",")+",b)}"),o.push("proto.step=function "+r+"_step("+u.join(",")+"){var "+a.map((function(t){return"a"+t+"=this.shape["+t+"]"})).join(",")+","+a.map((function(t){return"b"+t+"=this.stride["+t+"]"})).join(",")+",c=this.offset,d=0,ceil=Math.ceil"),v=0;v<e;++v)o.push("if(typeof i"+v+"==='number'){d=i"+v+"|0;if(d<0){c+=b"+v+"*(a"+v+"-1);a"+v+"=ceil(-a"+v+"/d)}else{a"+v+"=ceil(a"+v+"/d)}b"+v+"*=d}");o.push("return new "+r+"(this.data,"+a.map((function(t){return"a"+t})).join(",")+","+a.map((function(t){return"b"+t})).join(",")+",c)}");var y=new Array(e),m=new Array(e);for(v=0;v<e;++v)y[v]="a[i"+v+"]",m[v]="b[i"+v+"]";for(o.push("proto.transpose=function "+r+"_transpose("+u+"){"+u.map((function(t,e){return t+"=("+t+"===undefined?"+e+":"+t+"|0)"})).join(";"),"var a=this.shape,b=this.stride;return new "+r+"(this.data,"+y.join(",")+","+m.join(",")+",this.offset)}"),o.push("proto.pick=function "+r+"_pick("+u+"){var a=[],b=[],c=this.offset"),v=0;v<e;++v)o.push("if(typeof i"+v+"==='number'&&i"+v+">=0){c=(c+this.stride["+v+"]*i"+v+")|0}else{a.push(this.shape["+v+"]);b.push(this.stride["+v+"])}");return o.push("var ctor=CTOR_LIST[a.length+1];return ctor(this.data,a,b,c)}"),o.push("return function construct_"+r+"(data,shape,stride,offset){return new "+r+"(data,"+a.map((function(t){return"shape["+t+"]"})).join(",")+","+a.map((function(t){return"stride["+t+"]"})).join(",")+",offset)}"),new Function("CTOR_LIST","ORDER",o.join("\n"))(c[t],s)}var c={float32:[],float64:[],int8:[],int16:[],int32:[],uint8:[],uint16:[],uint32:[],array:[],uint8_clamped:[],bigint64:[],biguint64:[],buffer:[],generic:[]};t.exports=function(t,e,r,n){if(void 0===t)return(0,c.array[0])([]);"number"==typeof t&&(t=[t]),void 0===e&&(e=[t.length]);var a=e.length;if(void 0===r){r=new Array(a);for(var s=a-1,l=1;s>=0;--s)r[s]=l,l*=e[s]}if(void 0===n)for(n=0,s=0;s<a;++s)r[s]<0&&(n-=(e[s]-1)*r[s]);for(var f=function(t){if(i(t))return"buffer";if(o)switch(Object.prototype.toString.call(t)){case"[object Float64Array]":return"float64";case"[object Float32Array]":return"float32";case"[object Int8Array]":return"int8";case"[object Int16Array]":return"int16";case"[object Int32Array]":return"int32";case"[object Uint8Array]":return"uint8";case"[object Uint16Array]":return"uint16";case"[object Uint32Array]":return"uint32";case"[object Uint8ClampedArray]":return"uint8_clamped";case"[object BigInt64Array]":return"bigint64";case"[object BigUint64Array]":return"biguint64"}return Array.isArray(t)?"array":"generic"}(t),h=c[f];h.length<=a+1;)h.push(u(f,h.length-1));return(0,h[a+1])(t,e,r,n)}},8863:t=>{"use strict";t.exports=function(t,e,r,n,i){var o,a;for(o=1;o<t;o++){if(0===r[o-1])return console.log("tridiagonalSolve: failed due to lack of diagonal dominance"),!1;a=e[o]/r[o-1],r[o]-=a*n[o-1],i[o]-=a*i[o-1]}if(0===r[t-1])return console.log("tridiagonalSolve: failed due to singular matrix"),!1;for(i[t-1]/=r[t-1],o=t-2;o>=0;o--){if(0===r[o])return console.log("tridiagonalSolve: failed due to singular matrix"),!1;i[o]=(i[o]-n[o]*i[o+1])/r[o]}return!0}},5653:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.base=void 0;var n=r(1688),i=n.base.maths,o=n.base.user;e.base={maths:i,user:o}},9133:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.Dust=void 0;var n=r(6532),i=r(5653),o=r(217),a=r(2367),s=r(8177),u=function(){function t(t,e,r,a,s){if(r<=0)throw Error("Expected at least one particle");this._Model=t,this._particles=[],this._random=s||new n.Random(new n.RngStateBuiltin);for(var u=new t(i.base,e,"ignore"),c=0;c<r;++c)this._particles.push(new o.Particle(u,a))}return t.prototype.nState=function(){return this._particles[0].size},t.prototype.nParticles=function(){return this._particles.length},t.prototype.step=function(){return this._particles[0].step},t.prototype.info=function(){return this._particles[0].info()},t.prototype.setPars=function(t,e){var r=this,n=this.step(),a=(this.nState(),new this._Model(i.base,t,"ignore"));if(a.size()!==this.nState())throw Error("Particle produced unexpected state size");this.forEachParticle((function(t,i){var s=e?void 0:t.state();r._particles[i]=new o.Particle(a,n,s)}))},t.prototype.setStep=function(t){this._particles.forEach((function(e){e.step=t}))},t.prototype.setState=function(t){this.checkState(t),this.forEachParticle((function(e,r){e.setState(t[r])}))},t.prototype.run=function(t){var e=this;this._particles.forEach((function(r){return r.run(t,e._random)}))},t.prototype.simulate=function(t,e){for(var r=null===e?this.nState():e.length,n=this.nParticles(),i=t.length,o=(0,s.dustStateTime)(r,n,i),a=function(r){u.run(t[r]),u.forEachParticle((function(t,n){t.copyState(o.viewParticle(n,r),e)}))},u=this,c=0;c<i;++c)a(c);return o},t.prototype.state=function(t){var e=null===t?this.nState():t.length,r=this.nParticles(),n=(0,a.dustState)(e,r);return this.forEachParticle((function(e,r){e.copyState(n.viewParticle(r),t)})),n},t.prototype.reorder=function(t){var e=this,r=this.nParticles();if(t.length!==r)throw Error("Expected index to have length ".concat(r)+" but given ".concat(t.length));this.forEachParticle((function(n,i){var o=t[i];if(o<0||o>=r)throw Error("Expected index to be in [0, ".concat(r-1,"]")+" but given ".concat(o));n.setState(e._particles[o].state(),!0)})),this.forEachParticle((function(t,e){t.swap()}))},t.prototype.forEachParticle=function(t){this._particles.forEach(t)},t.prototype.checkState=function(t){if(t.length!==this.nParticles())throw Error("Invalid length state, expected ".concat(this.nParticles())+" but given ".concat(t.length));for(var e=0;e<t.length;++e)if(t[e].length!==this.nState())throw Error("Invalid length state for particle ".concat(e,",")+" expected ".concat(this.nState())+" but given ".concat(t[e].length))},t}();e.Dust=u},217:function(t,e){"use strict";var r=this&&this.__spreadArray||function(t,e,r){if(r||2===arguments.length)for(var n,i=0,o=e.length;i<o;i++)!n&&i in e||(n||(n=Array.prototype.slice.call(e,0,i)),n[i]=e[i]);return t.concat(n||Array.prototype.slice.call(e))};Object.defineProperty(e,"__esModule",{value:!0}),e.Particle=void 0;var n=function(){function t(t,e,n){this._model=t,this.step=e,this.size=t.size(),this._y=void 0===n?t.initial(this.step):r([],n,!0),this._yNext=Array(this.size).fill(0)}return t.prototype.info=function(){return this._model.info()},t.prototype.run=function(t,e){for(;this.step<t;)this._model.update(this.step++,this._y,this._yNext,e),this.swap()},t.prototype.state=function(){return this._y},t.prototype.copyState=function(t,e){var r=this._y;if(null===e)for(var n=0;n<r.length;++n)t.set(n,r[n]);else for(n=0;n<e.length;++n)t.set(n,r[e[n]])},t.prototype.setState=function(t,e){void 0===e&&(e=!1);for(var r=e?this._yNext:this._y,n=0;n<t.length;++n)r[n]=t[n]},t.prototype.swap=function(){var t;t=[this._yNext,this._y],this._y=t[0],this._yNext=t[1]},t}();e.Particle=n},586:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.variableNames=e.PkgWrapper=void 0;var n=r(6532),i=r(5653),o=r(9133),a=r(9882),s=function(){function t(t,e,r,o){this.Model=t,this.pars=e,this.random=new n.Random(o||new n.RngStateBuiltin),this.model=new this.Model(i.base,e,r)}return t.random=function(t){return new n.Random(t)},t.prototype.setUser=function(t,e){this.pars=t,this.model=new this.Model(i.base,t,e)},t.prototype.initial=function(t){return this.model.initial(t)},t.prototype.update=function(t,e){var r=Array(e.length).fill(0);return this.model.update(t,e,r,this.random),r},t.prototype.getInternal=function(){return this.model.getInternal()},t.prototype.getMetadata=function(){var t=this.model.info();return{info:t,names:u(t),size:this.model.size()}},t.prototype.run=function(t,e){var r=t[0],n=new o.Dust(this.Model,this.pars,1,r,this.random);null!==e&&n.setState([e]);var i=n.simulate(t,null);return{size:this.model.size(),y:Array.from(i.state.data)}},t}();function u(t){for(var e=[],r=function(t){var r=t.dim,n=t.name;0===r.length?e.push(n):e.push.apply(e,(0,a.combinations)(r).map((function(t){return"".concat(n,"[").concat(t.join(","),"]")})))},n=0,i=t;n<i.length;n++)r(i[n]);return e}e.PkgWrapper=s,e.variableNames=u},8177:function(t,e,r){"use strict";var n=this&&this.__importDefault||function(t){return t&&t.__esModule?t:{default:t}};Object.defineProperty(e,"__esModule",{value:!0}),e.DustStateTime=e.dustStateTime=void 0;var i=n(r(2861)),o=r(2367);e.dustStateTime=function(t,e,r){var n=(0,i.default)(new Float64Array(t*e*r),[r,e,t]);return new a(n,t,e,r)};var a=function(){function t(t,e,r,n){this.nState=e,this.nParticles=r,this.nTime=n,this.state=t}return t.prototype.viewTime=function(t){var e=this.state.pick(t,null,null);return new o.DustState(e,this.nState,this.nParticles)},t.prototype.viewParticle=function(t,e){return this.state.pick(e,t,null)},t.prototype.viewState=function(t,e){return this.state.pick(e,null,t)},t.prototype.viewTrace=function(t,e){return this.state.pick(null,e,t)},t.prototype.getParticle=function(t,e){return(0,o.copyVector)(this.viewParticle(t,e))},t.prototype.getState=function(t,e){return(0,o.copyVector)(this.viewState(t,e))},t.prototype.getTrace=function(t,e){return(0,o.copyVector)(this.viewTrace(t,e))},t}();e.DustStateTime=a},2367:function(t,e,r){"use strict";var n=this&&this.__importDefault||function(t){return t&&t.__esModule?t:{default:t}};Object.defineProperty(e,"__esModule",{value:!0}),e.copyVector=e.DustState=e.dustState=void 0;var i=n(r(2861));e.dustState=function(t,e){var r=(0,i.default)(new Float64Array(t*e),[e,t]);return new o(r,t,e)};var o=function(){function t(t,e,r){this.nState=e,this.nParticles=r,this.state=t}return t.prototype.viewParticle=function(t){return this.state.pick(t,null)},t.prototype.viewState=function(t){return this.state.pick(null,t)},t.prototype.getParticle=function(t){return a(this.viewParticle(t))},t.prototype.getState=function(t){return a(this.viewState(t))},t.prototype.asMatrix=function(){for(var t=[],e=0;e<this.nParticles;++e)t.push(this.getParticle(e));return t},t}();function a(t){for(var e=t.shape[0],r=Array(e),n=0;n<e;++n)r[n]=t.get(n);return r}e.DustState=o,e.copyVector=a},9882:function(t,e){"use strict";var r=this&&this.__spreadArray||function(t,e,r){if(r||2===arguments.length)for(var n,i=0,o=e.length;i<o;i++)!n&&i in e||(n||(n=Array.prototype.slice.call(e,0,i)),n[i]=e[i]);return t.concat(n||Array.prototype.slice.call(e))};function n(t,e){for(var r=[],n=t;n<=e;++n)r.push(n);return r}function i(t,e){return t[0].map((function(r,n){return e(t.map((function(t){return t[n]})))}))}function o(t){return t.reduce((function(t,e){return t+e}),0)/t.length}Object.defineProperty(e,"__esModule",{value:!0}),e.meanArray=e.mean=e.applyArray=e.isEqualArray=e.combinations=e.seq=void 0,e.seq=n,e.combinations=function t(e){var i=n(1,e[0]);if(e.length>1){for(var o=[],a=0,s=t(e.slice(1));a<s.length;a++)for(var u=s[a],c=0,l=i;c<l.length;c++){var f=l[c];o.push(r([f],u,!0))}return o}return i.map((function(t){return[t]}))},e.isEqualArray=function(t,e){return t.every((function(t,r){return t===e[r]}))},e.applyArray=i,e.mean=o,e.meanArray=function(t){return i(t,o)}},9556:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.versions=void 0,e.versions=function(){return{random:"0.0.3",odinjs:"0.0.15",dust:"0.0.1"}}},7915:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.tidyDiscreteSolutionVariable=e.tidyDiscreteSolution=e.wodinRunDiscrete=e.runModelDiscrete=e.DiscreteSeriesMode=void 0;var n,i=r(9133),o=r(9882);function a(t,e,r,n,a,s){var u=Math.floor(r/a),c=Math.ceil(n/a),l=new i.Dust(t,e,s,u),f=l.info(),h=(0,o.seq)(u,c),p=h.map((function(t){return t*a}));return{info:f,state:l.simulate(h,null),times:p}}function s(t){var e=t.times,r=[];return t.info.forEach((function(e){var n=u(e.name,t);r.push.apply(r,n)})),{x:e,values:r}}function u(t,e){for(var r=e.info.map((function(t){return t.name})).indexOf(t),i=e.state,a=i.getTrace(r,0),s=!1,u=[],c=0;c<i.nParticles;++c){var l=i.getTrace(r,c);u.push(l),!s&&c>1&&!(0,o.isEqualArray)(l,a)&&(s=!0)}if(s){var f=u.map((function(e){return{mode:n.Individual,name:t,y:e}}));return f.push({mode:n.Mean,name:t,y:(0,o.meanArray)(u)}),f}return[{mode:n.Deterministic,name:t,y:a}]}!function(t){t.Individual="Individual",t.Mean="Mean",t.Deterministic="Deterministic"}(n=e.DiscreteSeriesMode||(e.DiscreteSeriesMode={})),e.runModelDiscrete=a,e.wodinRunDiscrete=function(t,e,r,n,i,o){return s(a(t,e,r,n,i,o))},e.tidyDiscreteSolution=s,e.tidyDiscreteSolutionVariable=u}},e={};function r(n){var i=e[n];if(void 0!==i)return i.exports;var o=e[n]={exports:{}};return t[n].call(o.exports,o,o.exports,r),o.exports}var n={};(()=>{"use strict";var t=n;Object.defineProperty(t,"__esModule",{value:!0}),t.wodinRunDiscrete=t.DiscreteSeriesMode=t.versions=t.DustStateTime=t.dustStateTime=t.DustState=t.dustState=t.PkgWrapper=t.Particle=t.Dust=t.base=void 0;var e=r(5653);Object.defineProperty(t,"base",{enumerable:!0,get:function(){return e.base}});var i=r(9133);Object.defineProperty(t,"Dust",{enumerable:!0,get:function(){return i.Dust}});var o=r(217);Object.defineProperty(t,"Particle",{enumerable:!0,get:function(){return o.Particle}});var a=r(586);Object.defineProperty(t,"PkgWrapper",{enumerable:!0,get:function(){return a.PkgWrapper}});var s=r(2367);Object.defineProperty(t,"dustState",{enumerable:!0,get:function(){return s.dustState}}),Object.defineProperty(t,"DustState",{enumerable:!0,get:function(){return s.DustState}});var u=r(8177);Object.defineProperty(t,"dustStateTime",{enumerable:!0,get:function(){return u.dustStateTime}}),Object.defineProperty(t,"DustStateTime",{enumerable:!0,get:function(){return u.DustStateTime}});var c=r(9556);Object.defineProperty(t,"versions",{enumerable:!0,get:function(){return c.versions}});var l=r(7915);Object.defineProperty(t,"DiscreteSeriesMode",{enumerable:!0,get:function(){return l.DiscreteSeriesMode}}),Object.defineProperty(t,"wodinRunDiscrete",{enumerable:!0,get:function(){return l.wodinRunDiscrete}})})(),dust=n})();
+var dust;(()=>{var t={8006:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.interpolateSearch=e.InterpolatorBase=void 0;var r=function(){function t(t,e){(function(t){return Array.isArray(t[0])})(e)||(e=[e]),this._i=0,this._x=t,this._y=e,this.nX=this._x.length,this.nY=this._y.length;for(var r=0;r<this.nY;++r)if(this._y[r].length!==this.nX)throw Error("Invalid length for 'y', expected ".concat(this.nX))}return t.prototype.evalAll=function(t){for(var e=Array(this.nY),r=0;r<this.nY;++r)e[r]=this.eval(t,r);return e},t.prototype.search=function(t,e){var r=n(t,this._x,this._i);if(r<0||!e&&r===this.nX)throw Error("Interpolation failed as 'x' is out of range");return this._i=r,r},t}();function n(t,e,r){var n=r,i=r,o=1,a=e.length;if(e[n]<=t){if(n>=a-1)return a;for(i=n+o;e[i]<t;)if(n=i,(i+=o*=2)>=a){if(e[i=a-1]<t)return a;break}}else{if(0===n)return-1;for(n-=o;e[n]>t;){if(i=n,n<(o*=2)){if(e[n=0]>t)return-1;break}n-=o}}for(;i-n>1;){var s=Math.floor((i+n)/2);e[s]<t?n=s:i=s}return n}e.InterpolatorBase=r,e.interpolateSearch=n},4731:function(t,e,r){"use strict";var n,i=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorConstant=void 0;var o=function(t){function e(){return null!==t&&t.apply(this,arguments)||this}return i(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=Math.min(this.search(t,!0),this.nX-1);return r!==this.nX-1&&this._x[r+1]===t&&++r,this._y[e][r]},e}(r(8006).InterpolatorBase);e.InterpolatorConstant=o},9673:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorSpline=e.InterpolatorLinear=e.InterpolatorConstant=e.InterpolatorBase=void 0;var n=r(8006);Object.defineProperty(e,"InterpolatorBase",{enumerable:!0,get:function(){return n.InterpolatorBase}});var i=r(4731);Object.defineProperty(e,"InterpolatorConstant",{enumerable:!0,get:function(){return i.InterpolatorConstant}});var o=r(9477);Object.defineProperty(e,"InterpolatorLinear",{enumerable:!0,get:function(){return o.InterpolatorLinear}});var a=r(7869);Object.defineProperty(e,"InterpolatorSpline",{enumerable:!0,get:function(){return a.InterpolatorSpline}})},9477:function(t,e,r){"use strict";var n,i=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorLinear=void 0;var o=function(t){function e(){return null!==t&&t.apply(this,arguments)||this}return i(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=this.search(t,!1),n=this._x[r],i=(t-n)/(this._x[r+1]-n),o=this._y[e],a=o[r];return a+(o[r+1]-a)*i},e}(r(8006).InterpolatorBase);e.InterpolatorLinear=o},7869:function(t,e,r){"use strict";var n,i=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)}),o=this&&this.__importDefault||function(t){return t&&t.__esModule?t:{default:t}};Object.defineProperty(e,"__esModule",{value:!0}),e.solve=e.InterpolatorSpline=void 0;var a=o(r(8863)),s=function(t){function e(e,r){var n=t.call(this,e,r)||this,i=function(t){var e=t.length,r=Array(e),n=Array(e),i=Array(e);r[0]=0,n[0]=2/(t[1]-t[0]),i[0]=1/(t[1]-t[0]);for(var o=1;o<e-1;++o){var a=t[o-1],s=t[o],u=t[o+1];r[o]=1/(s-a),n[o]=2*(1/(s-a)+1/(u-s)),i[o]=1/(u-s)}return r[e-1]=1/(t[e-1]-t[e-2]),n[e-1]=2/(t[e-1]-t[e-2]),i[e-1]=0,[r,n,i]}(n._x),o=function(t,e){for(var r=t.length,n=e.length,i=r-1,o=[],a=0;a<n;++a){var s=Array(r),u=e[a];s[0]=3*(u[1]-u[0])/Math.pow(t[1]-t[0],2);for(var c=1;c<i;++c){var l=t[c-1],f=t[c],h=t[c+1],p=u[c-1],d=u[c],v=u[c+1];s[c]=3*((d-p)/Math.pow(f-l,2)+(v-d)/Math.pow(h-f,2))}s[i]=3*(u[i]-u[i-1])/Math.pow(t[i]-t[i-1],2),o.push(s)}return o}(n._x,n._y);return n._k=function(t,e){for(var r=t[0],n=t[1],i=t[2],o=r.length,a=0,s=e;a<s.length;a++){var c=s[a];u(o,r,n.slice(),i,c)}return e}(i,o),n}return i(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=this.search(t,!1);return Array(this.nY),function(t,e,r,n,i){var o=(e-r[t])/(r[t+1]-r[t]),a=i[t]*(r[t+1]-r[t])-(n[t+1]-n[t]),s=-i[t+1]*(r[t+1]-r[t])+(n[t+1]-n[t]);return(1-o)*n[t]+o*n[t+1]+o*(1-o)*(a*(1-o)+s*o)}(r,t,this._x,this._y[e],this._k[e])},e}(r(8006).InterpolatorBase);function u(t,e,r,n,i){if(!(0,a.default)(t,e,r,n,i))throw Error("solve failed: singular matrix?")}e.InterpolatorSpline=s,e.solve=u},8251:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.base=void 0;const n=r(3350),i=r(320),o=r(9928),a=r(2630);e.base={delay:n.delay,interpolate:i.interpolate,maths:o.maths,user:a.user}},7601:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.updatePars=e.batchParsDisplace=e.batchParsRange=e.batchRun=e.Batch=void 0;const n=r(46),i=r(1729),o=r(7466);class a{constructor(t,e,r,n,i){this.pars=e,this.tStart=r,this.tEnd=n;const a=[],s=[],c=[];if(e.values.forEach((l=>{const f=u(e.base,e.name,l);try{a.push((0,o.wodinRun)(t,f,r,n,i)),c.push(l)}catch(t){s.push({value:l,error:t.message})}})),0===a.length)throw Error(`All solutions failed; first error: ${s[0].error}`);this.pars=Object.assign(Object.assign({},e),{values:c}),this.solutions=a,this.errors=s}valueAtTime(t){const e=this.solutions.map((e=>e({mode:n.TimeMode.Given,times:[t]})));return{names:e[0].names,x:this.pars.values,y:e[0].names.map(((t,r)=>e.map((t=>t.y[r][0]))))}}extreme(t){return this.findExtremes()[t]}findExtremes(){if(void 0===this._extremes){const t=501,e={mode:n.TimeMode.Grid,nPoints:t,tEnd:this.tEnd,tStart:this.tStart},r=this.solutions.map((t=>t(e))),o=r[0].names,a=this.pars.values,s=(0,i.loop)(o.length,(t=>r.map((e=>function(t,e){const r=e.x,n=e.y[t],o=(0,i.whichMin)(n),a=(0,i.whichMax)(n),s=r[o],u=r[a],c=n[o];return{tMax:u,tMin:s,yMax:n[a],yMin:c}}(t,e)))));this._extremes={tMax:l("tMax",o,a,s),tMin:l("tMin",o,a,s),yMax:l("yMax",o,a,s),yMin:l("yMin",o,a,s)}}return this._extremes}}function s(t,e,r,n,o,a){const s=c(t,e);if(o>s)throw Error(`Expected lower bound to be no greater than ${s}`);if(a<s)throw Error(`Expected upper bound to be no less than ${s}`);if(o>=a)throw Error("Expected upper bound to be greater than lower bound");if(r<2)throw Error("Must include at least 2 traces in the batch");if(n&&o<=0)throw Error("Lower bound must be greater than 0 for logarithmic scale");return{base:t,name:e,values:n?(0,i.gridLog)(o,a,r):(0,i.grid)(o,a,r)}}function u(t,e,r){const n=Object.assign({},t);return n[e]=r,n}function c(t,e){const r=t[e];if(void 0===r)throw Error(`Expected a value for '${e}'`);if("number"!=typeof r)throw Error(`Expected a number for '${e}'`);return r}function l(t,e,r,n){return{names:e,x:r,y:(0,i.loop)(e.length,(e=>n[e].map((e=>e[t]))))}}e.Batch=a,e.batchRun=function(t,e,r,n,i){return new a(t,e,r,n,i)},e.batchParsRange=s,e.batchParsDisplace=function(t,e,r,n,i){const o=c(t,e),a=i/100;return s(t,e,r,n,o*(1-a),o*(1+a))},e.updatePars=u},3350:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.delay=void 0,e.delay=function(t,e,r,n){const i=t(e);for(let t=0;t<r.length;++t)n[t]=i[r[t]]}},7771:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.sumOfSquares=e.updatePars=e.fitTarget=void 0;const n=r(8251),i=r(1267),o=r(46);function a(t,e){const r=Object.assign({},t.base);for(let n=0;n<t.vary.length;++n)r[t.vary[n]]=e[n];return r}function s(t,e){let r=0;for(let n=0;n<t.length;++n){const i=t[n];isNaN(i)||(r+=Math.pow(i-e[n],2))}return r}e.fitTarget=function(t,e,r,u,c){if(e.time[0]<0)throw Error("Expected the first time to be at least 0");const l=e.time[e.time.length-1];return f=>{const h=a(r,f),p=new t(n.base,h,"error"),d=(0,i.runModel)(p,null,0,l,c).solution,v=p.names(),y=v.indexOf(u),m=d(e.time).map((t=>t[y]));return{data:{endTime:l,names:v,pars:h,solution:(0,o.interpolatedSolution)(d,v,0,l),startTime:0},value:s(e.value,m)}}},e.updatePars=a,e.sumOfSquares=s},1688:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.versions=e.batchRun=e.batchParsRange=e.batchParsDisplace=e.Batch=e.TimeMode=e.base=e.PkgWrapper=e.wodinRun=e.wodinFitValue=e.wodinFit=void 0;var n=r(7466);Object.defineProperty(e,"wodinFit",{enumerable:!0,get:function(){return n.wodinFit}}),Object.defineProperty(e,"wodinFitValue",{enumerable:!0,get:function(){return n.wodinFitValue}}),Object.defineProperty(e,"wodinRun",{enumerable:!0,get:function(){return n.wodinRun}});var i=r(8714);Object.defineProperty(e,"PkgWrapper",{enumerable:!0,get:function(){return i.PkgWrapper}});var o=r(8251);Object.defineProperty(e,"base",{enumerable:!0,get:function(){return o.base}});var a=r(46);Object.defineProperty(e,"TimeMode",{enumerable:!0,get:function(){return a.TimeMode}});var s=r(7601);Object.defineProperty(e,"Batch",{enumerable:!0,get:function(){return s.Batch}}),Object.defineProperty(e,"batchParsDisplace",{enumerable:!0,get:function(){return s.batchParsDisplace}}),Object.defineProperty(e,"batchParsRange",{enumerable:!0,get:function(){return s.batchParsRange}}),Object.defineProperty(e,"batchRun",{enumerable:!0,get:function(){return s.batchRun}});var u=r(7497);Object.defineProperty(e,"versions",{enumerable:!0,get:function(){return u.versions}})},320:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.interpolate=e.interpolateTimes=e.interpolateCheckT=e.interpolateCheckY=e.interpolateAlloc=void 0;const n=r(9673);function i(t,e,r){const i=function(t,e){if(t.length===e)return[t];const r=[],n=t.length/e;for(let i=0;i<n;++i)r.push(t.slice(e*i,e*(i+1)));return r}(r,e.length);if("constant"===t)return new n.InterpolatorConstant(e,i);if("linear"===t)return new n.InterpolatorLinear(e,i);if("spline"===t)return new n.InterpolatorSpline(e,i);throw Error(`Invalid interpolation method '${t}'`)}function o(t,e,r,n){if(1===e.length){if(t[0]!==e[0])throw Error(`Expected ${r} to have length ${t[0]} (for ${n})`)}else for(let i=0;i<e.length;++i)if(t[i]!==e[i])throw Error(`Expected dimension ${i+1} of ${r} to have size ${t[i]} (for ${n})`)}function a(t,e){const r=Math.max(...t);return{max:Math.min(...e),min:r}}e.interpolateAlloc=i,e.interpolateCheckY=o,e.interpolateCheckT=function(t,e,r,n){if(void 0===r)return 1/0;if(t<r.min)throw Error(`Integration times do not span interpolation range; min: ${r.min}`);if(e>r.max)throw Error(`Integration times do not span interpolation range; max: ${r.max}`);return void 0===n&&(n=r.max),n},e.interpolateTimes=a,e.interpolate={alloc:i,checkY:o,times:a}},9928:(t,e)=>{"use strict";function r(t){return.5===i(t,1)?2*Math.round(t/2):Math.round(t)}function n(t,e){if(void 0===e||0===e)return r(t);{const n=Math.pow(10,e);return r(t*n)/n}}function i(t,e){let r=t%e;return r*e<0&&(r+=e),r}function o(t,e){return Math.floor(t/e)}function a(t,e,r){let n=0;for(let i=e;i<r;++i)n+=t[i];return n}function s(t,e,r,n,i,o){let a=0;for(let s=n;s<i;++s){const n=s*o;for(let i=e;i<r;++i)a+=t[i+n]}return a}function u(t,e,r,n,i,o,a,s,u){let c=0;for(let l=o;l<a;++l){const o=l*u;for(let a=n;a<i;++a){const n=a*s+o;for(let i=e;i<r;++i)c+=t[i+n]}}return c}Object.defineProperty(e,"__esModule",{value:!0}),e.maths=e.odinSum3=e.odinSum2=e.odinSum1=e.intdivr=e.modr=e.round2=void 0,e.round2=n,e.modr=i,e.intdivr=o,e.odinSum1=a,e.odinSum2=s,e.odinSum3=u,e.maths={intdivr:o,modr:i,odinSum1:a,odinSum2:s,odinSum3:u,round2:n}},1267:function(t,e,r){"use strict";var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var i=Object.getOwnPropertyDescriptor(e,r);i&&!("get"in i?!e.__esModule:i.writable||i.configurable)||(i={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,i)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),i=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),o=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return i(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.runModel=e.isODEModel=e.isDDEModel=void 0;const a=o(r(6245)),s=r(320);function u(t){return 4===t.rhs.length}e.isDDEModel=u,e.isODEModel=function(t){return 3===t.rhs.length},e.runModel=function(t,e,r,n,i){const o=t.getMetadata().interpolateTimes;return i.tcrit=(0,s.interpolateCheckT)(r,n,o,i.tcrit),u(t)?function(t,e,r,n,i){let o=null;"function"==typeof t.output&&(o=(e,r,n)=>t.output(e,r,n)),null===e&&(e=t.initial(r)),t.getInternal().initial_t=r;const s=new a.DDE((function(e,r,n,i){t.rhs(e,r,n,i)}),e.length,i,o);return s.initialise(r,e),{solution:s.run(n),statistics:s.statistics()}}(t,e,r,n,i):function(t,e,r,n,i){let o=null;"function"==typeof t.output&&(o=(e,r)=>t.output(e,r)),null===e&&(e=t.initial(r));const s=new a.Dopri((function(e,r,n){t.rhs(e,r,n)}),e.length,i,o);return s.initialise(r,e),{solution:s.run(n),statistics:s.statistics()}}(t,e,r,n,i)}},8714:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.PkgWrapper=void 0;const n=r(8251),i=r(1267);e.PkgWrapper=class{constructor(t,e,r){this.model=new t(n.base,e,r)}initial(t){return this.model.initial(t)}rhs(t,e){const r=new Array(e.length).fill(0);let n=null;if(!(0,i.isODEModel)(this.model))throw Error("Can't use rhs() with delay models");return this.model,this.model.rhs(t,e,r),this.model.output&&(n=this.model.output(t,e)),{output:n,state:r}}getMetadata(){return this.model.getMetadata()}getInternal(){return this.model.getInternal()}setUser(t,e){this.model.setUser(t,e)}run(t,e,r){const n=t[0],o=t[t.length-1],a=(0,i.runModel)(this.model,e,n,o,r);return{names:this.model.names(),statistics:a.statistics,y:a.solution(t)}}}},46:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.interpolatedSolution=e.TimeMode=void 0;const n=r(1729);var i;!function(t){t.Grid="grid",t.Given="given"}(i=e.TimeMode||(e.TimeMode={})),e.interpolatedSolution=function(t,e,r,o){return a=>{const s=function(t,e,r){switch(t.mode){case i.Grid:return(0,n.grid)(Math.max(e,t.tStart),Math.min(r,t.tEnd),t.nPoints);case i.Given:return t.times}}(a,r,o),u=t(s),c=u[0].map(((t,e)=>u.map((t=>t[e]))));return{names:e,x:s,y:c}}}},2630:(t,e)=>{"use strict";function r(t,e,r){if("ignore"===r)return;const n=[];for(const r of Object.keys(t))e.includes(r)||n.push(r);if(n.length>0){const t="Unknown user parameters: "+n.join(", ");if("message"===r)console.log(t);else{if("warning"!==r)throw"stop"===r?Error(t):Error(t+" (and invalid value for unusedUserAction)");console.warn(t)}}}function n(t,e,r,n,i,o,a){const s=t[e];if(void 0===s){if(void 0!==r[e])return;if(null===n)throw Error(`Expected a value for '${e}'`);r[e]=n}else{if("number"!=typeof s)throw Error(`Expected a number for '${e}'`);c(s,i,o,a,e),r[e]=s}}function i(t,e,r,n,i,o,c){let l=t[e];if(void 0===l){if(void 0!==r[e])return;throw Error(`Expected a value for '${e}'`)}{const t=n.length-1;l=a(l),s(t,l,e),function(t,e,r){const n=t.length-1;for(let i=0;i<n;++i){const o=t[i+1];if(e.dim[i]!==o)throw 1===n?Error(`Expected length ${o} value for '${r}'`):Error(`Incorrect size of dimension ${i+1} of '${r}' (expected ${o})`)}}(n,l,e),u(l,i,o,c,e),r[e]=l.data.slice()}}function o(t,e,r,n,i,o,c){let l=t[e];if(void 0===l){if(void 0!==r[e])return;throw Error("Expected a value for '"+e+"'")}{const t=n.length-1;l=a(l),s(t,l,e),u(l,i,o,c,e),n[0]=l.data.length;for(let e=0;e<t;++e)n[e+1]=l.dim[e];r[e]=l.data.slice()}}function a(t,e){return Array.isArray(t)?t={data:t,dim:[t.length]}:"number"==typeof t&&(t={data:[t],dim:[1]}),t}function s(t,e,r){if(e.dim.length!==t)throw 1===t?Error(`Expected a numeric vector for '${r}'`):2===t?Error(`Expected a numeric matrix for '${r}'`):Error(`Expected a numeric array of rank ${t} for '${r}'`)}function u(t,e,r,n,i){for(const o of t.data){if(null===o)throw Error(`'${i}' must not contain any NA values`);c(o,e,r,n,i)}}function c(t,e,r,n,i){if("number"!=typeof t)throw Error(`Expected a number for '${i}'`);if(t<e)throw Error(`Expected '${i}' to be at least ${e}`);if(t>r)throw Error(`Expected '${i}' to be at most ${r}`);if(n&&!Number.isInteger(t))throw Error(`Expected '${i}' to be integer-like`)}Object.defineProperty(e,"__esModule",{value:!0}),e.user=e.setUserArrayVariable=e.setUserArrayFixed=e.setUserScalar=e.checkUser=void 0,e.checkUser=r,e.setUserScalar=n,e.setUserArrayFixed=i,e.setUserArrayVariable=o,e.user={checkUser:r,setUserArrayFixed:i,setUserArrayVariable:o,setUserScalar:n}},1729:(t,e)=>{"use strict";function r(t,e,r){const n=(e-t)/(r-1),i=[];for(let e=0;e<r-1;++e)i.push(t+e*n);return i.push(e),i}Object.defineProperty(e,"__esModule",{value:!0}),e.loop=e.whichMax=e.whichMin=e.gridLog=e.grid=void 0,e.grid=r,e.gridLog=function(t,e,n){return r(Math.log(t),Math.log(e),n).map(Math.exp)},e.whichMin=function(t){let e=-1,r=1/0;for(let n=0;n<t.length;++n)t[n]<r&&(e=n,r=t[n]);return e},e.whichMax=function(t){let e=-1,r=-1/0;for(let n=0;n<t.length;++n)t[n]>r&&(e=n,r=t[n]);return e},e.loop=function(t,e){const r=[];for(let n=0;n<t;++n)r.push(e(n));return r}},7497:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.versions=void 0,e.versions=function(){return{dfoptim:"0.0.5",dopri:"0.0.12",odinjs:"0.0.18"}}},7466:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.wodinFitValue=e.wodinFit=e.wodinRun=void 0;const n=r(4153),i=r(8251),o=r(7771),a=r(1267),s=r(46);e.wodinRun=function(t,e,r,n,o){const u=new t(i.base,e,"error"),c=(0,a.runModel)(u,null,r,n,o).solution,l=u.names();return(0,s.interpolatedSolution)(c,l,r,n)},e.wodinFit=function(t,e,r,i,a,s){const u=(0,o.fitTarget)(t,e,r,i,a),c=r.vary.map((t=>r.base[t]));return new n.Simplex(u,c,s)},e.wodinFitValue=function(t,e,r){const{names:n,y:i}=t({mode:s.TimeMode.Given,times:e.time}),a=n.indexOf(r);return(0,o.sumOfSquares)(e.value,i[a])}},7689:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.binomial=void 0,e.binomial=function(t,e,r){if(0===e||0===r)return 0;if(1===r)return e;const o=r<=.5?r:1-r,a=e*o>=10?function(t,e,r){const i=Math.sqrt(e*r*(1-r)),o=1.15+2.53*i,a=.0248*o-.0873+.01*r,s=e*r+.5,u=.92-4.2/o,c=r/(1-r),l=(2.83+5.1/o)*i,f=Math.floor((e+1)*r);let h;for(;;){const r=t.random()-.5;let i=t.random();const p=.5-Math.abs(r),d=Math.floor((2*a/p+o)*r+s);if(p>=.07&&i<=u){h=d;break}if(!(d<0||d>e)&&(i=Math.log(i*l/(a/(p*p)+o)),i<=(f+.5)*Math.log((f+1)/(c*(e-f+1)))+(e+1)*Math.log((e-f+1)/(e-d+1))+(d+.5)*Math.log(c*(e-d+1)/(d+1))+n(f)+n(e-f)-n(d)-n(e-d))){h=d;break}}return h}(t,e,o):function(t,e,r){let n=-1;do{n=i(t.random(),e,r)}while(n<0);return n}(t,e,o);return r<=.5?a:e-a};const r=[.08106146679532733,.04134069595540946,.02767792568499805,.02079067210376495,.016644691189821703,.013876128823071987,.011896709945892425,.010411265261971892,.00925546218270945,.008330563433359472];function n(t){return t<r.length?r[t]:function(t){const e=(t+1)*(t+1);return(1/12-(1/360-1/1260/e)/e)/(t+1)}(t)}function i(t,e,r){const n=1-r,i=r/n,o=i*(e+1);let a=Math.pow(n,e),s=0,u=a;for(;t>=a;)t-=a,s++,a*=o/s-i,u=a;return s}},9386:(t,e)=>{"use strict";function r(t){return-Math.log(t.random())}Object.defineProperty(e,"__esModule",{value:!0}),e.exponential=e.randomExponential=void 0,e.randomExponential=r,e.exponential=function(t,e){return r(t)/e}},6532:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.uniform=e.randomUniform=e.randomNormal=e.normal=e.randomExponential=e.exponential=e.binomial=e.RngStateReplay=e.RngStateObserved=e.RngStateBuiltin=e.RngState=e.Random=void 0;var n=r(1929);Object.defineProperty(e,"Random",{enumerable:!0,get:function(){return n.Random}});var i=r(3726);Object.defineProperty(e,"RngState",{enumerable:!0,get:function(){return i.RngState}});var o=r(2035);Object.defineProperty(e,"RngStateBuiltin",{enumerable:!0,get:function(){return o.RngStateBuiltin}});var a=r(2301);Object.defineProperty(e,"RngStateObserved",{enumerable:!0,get:function(){return a.RngStateObserved}});var s=r(5843);Object.defineProperty(e,"RngStateReplay",{enumerable:!0,get:function(){return s.RngStateReplay}});var u=r(7689);Object.defineProperty(e,"binomial",{enumerable:!0,get:function(){return u.binomial}});var c=r(9386);Object.defineProperty(e,"exponential",{enumerable:!0,get:function(){return c.exponential}}),Object.defineProperty(e,"randomExponential",{enumerable:!0,get:function(){return c.randomExponential}});var l=r(2638);Object.defineProperty(e,"normal",{enumerable:!0,get:function(){return l.normal}}),Object.defineProperty(e,"randomNormal",{enumerable:!0,get:function(){return l.randomNormal}});var f=r(6390);Object.defineProperty(e,"randomUniform",{enumerable:!0,get:function(){return f.randomUniform}}),Object.defineProperty(e,"uniform",{enumerable:!0,get:function(){return f.uniform}})},2638:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.normal=e.randomNormal=void 0;const r=2*Math.PI;function n(t){let e=0,n=0;do{e=t.random(),n=t.random()}while(e<Number.EPSILON);return Math.sqrt(-2*Math.log(e))*Math.cos(r*n)}e.randomNormal=n,e.normal=function(t,e,r){return e+n(t)*r}},5789:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.poisson=void 0;const n=r(7788);e.poisson=function(t,e){return e<10?function(t,e){let r=0;const n=Math.exp(-e);let o=1;for(;o*=t.random(),!(o<=n&&r<=i);)r++;return r}(t,e):function(t,e){let r=0;const i=Math.log(e),o=.931+2.53*Math.sqrt(e),a=.02483*o-.059,s=1.1239+1.1328/(o-3.4);for(;;){let u=t.random();u-=.5;const c=t.random(),l=.5-Math.abs(u),f=Math.floor((2*a/l+o)*u+e+.43);if(l>=.07&&c<=.9277-3.6224/(o-2)){r=f;break}if(!(f<0||l<.013&&c>l)&&Math.log(c*s/(a/(l*l)+o))<=f*i-e-(0,n.lngamma)(f+1)){r=f;break}}return r}(t,e)};const i=2**31-1},1929:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.Random=void 0;const n=r(7689),i=r(9386),o=r(2638),a=r(5789),s=r(6390);e.Random=class{state;constructor(t){this.state=t}randomExponential(){return(0,i.randomExponential)(this.state)}randomNormal(){return(0,o.randomNormal)(this.state)}randomUniform(){return(0,s.randomUniform)(this.state)}binomial(t,e){return(0,n.binomial)(this.state,t,e)}exponential(t){return(0,i.exponential)(this.state,t)}normal(t,e){return(0,o.normal)(this.state,t,e)}poisson(t){return(0,a.poisson)(this.state,t)}uniform(t,e){return(0,s.uniform)(this.state,t,e)}}},2035:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.RngStateBuiltin=void 0;const n=r(3726);class i extends n.RngState{constructor(){super()}random(){return Math.random()}}e.RngStateBuiltin=i},2301:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.RngStateObserved=void 0;const n=r(3726),i=r(2035),o=r(5843);class a extends n.RngState{_state;_values;constructor(t){super(),void 0===t&&(t=new i.RngStateBuiltin),this._state=t,this._values=[]}random(){const t=this._state.random();return this._values.push(t),t}replay(){return new o.RngStateReplay(this._values)}length(){return this._values.length}values(){return this._values}}e.RngStateObserved=a},5843:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.RngStateReplay=void 0;const n=r(3726);class i extends n.RngState{_i;_values;constructor(t){super(),this._i=0,this._values=t}random(){const t=this._values[this._i++];if(this._i>this._values.length)throw Error("Exhausted stream");return t}length(){return this._values.length-this._i}}e.RngStateReplay=i},3726:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.RngState=void 0,e.RngState=class{setSeed(t){throw Error("This generator cannot be seeded (really!)")}}},6390:(t,e)=>{"use strict";function r(t){return t.random()}Object.defineProperty(e,"__esModule",{value:!0}),e.uniform=e.randomUniform=void 0,e.randomUniform=r,e.uniform=function(t,e,n){return r(t)*(n-e)+e}},4394:(t,e)=>{"use strict";function r(t,e){return void 0===t?e:t}Object.defineProperty(e,"__esModule",{value:!0}),e.simplexControl=void 0,e.simplexControl=function(t){void 0===t&&(t={});var e={deltaNonZero:r(t.deltaNonZero,.05),deltaZero:r(t.deltaZero,.001),errorOnFailure:r(t.errorOnFailure,!1),tolerance:r(t.tolerance,1e-5)};if(e.tolerance<=0)throw new Error("Invalid control parameter: 'tolerance' must be strictly positive");return e}},4153:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0});var n=r(2169);Object.defineProperty(e,"Simplex",{enumerable:!0,get:function(){return n.Simplex}}),Object.defineProperty(e,"fitSimplex",{enumerable:!0,get:function(){return n.fitSimplex}})},2169:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.Simplex=e.fitSimplex=void 0;var n=r(4394),i=r(2003);function o(t,e,r){for(var n=e.length,i=new Array(n),o=1+t,a=-t,s=0;s<i.length;++s)i[s]=o*e[s]+a*r[s];return i}e.fitSimplex=function(t,e,r,n){return void 0===r&&(r={}),void 0===n&&(n=200),new a(t,e,r).run(n*e.length)};var a=function(){function t(t,e,r){void 0===r&&(r={}),this.rho=1,this.chi=2,this.psi=-.5,this.sigma=.5,this._iterations=0,this._id=0,this._converged=!1,this._target=t,this._control=n.simplexControl(r),this._n=e.length,this._simplex=[],this._simplex.push(this._point(e.slice())),r.errorOnFailure||(this._target=function(t){return function(e){try{return t(e)}catch(t){return 1/0}}}(t));for(var i=0;i<this._n;++i){var o=e.slice();o[i]?o[i]*=1+this._control.deltaNonZero:o[i]=this._control.deltaZero,this._simplex.push(this._point(o))}this._sort()}return t.prototype.step=function(){if(this._iterations++,this._isConverged())return this._converged=!0,!0;var t=this._n,e=this._simplex[0],r=this._simplex[t],n=this._centroid(),i=this._reflect(r,n);if(i.value<e.value){var o=this._expand(r,n);this._update(o.value<i.value?o:i)}else if(i.value>=this._simplex[t-1].value){var a=i.value>r.value?this._contractInside(r,n):this._contractOutside(r,n);if(a.value<r.value)this._update(a);else{for(var s=1;s<=t;++s)this._simplex[s]=this._shrink(this._simplex[s],e.location);this._sort()}}else this._update(i);return!1},t.prototype.run=function(t){if(!this._converged)for(var e=0;e<t&&!this.step();++e);return this.result()},t.prototype.result=function(){var t=this._simplex[0];return{converged:this._converged,data:t.data,evaluations:this._id,iterations:this._iterations,location:t.location,value:t.value}},t.prototype.simplex=function(){return this._simplex.map((function(t){return{location:t.location,value:t.value}}))},t.prototype._point=function(t){var e=i.checkResult(this._target(t));return{data:e.data,id:this._id++,location:t,value:e.value}},t.prototype._sort=function(){this._simplex.sort((function(t,e){return t.value-e.value}))},t.prototype._centroid=function(){for(var t=this._n,e=new Array(t),r=0;r<t;++r){e[r]=0;for(var n=0;n<t;++n)e[r]+=this._simplex[n].location[r];e[r]/=t}return e},t.prototype._update=function(t){this._simplex[this._n]=t,this._sort()},t.prototype._reflect=function(t,e){return this._point(o(this.rho,e,t.location))},t.prototype._expand=function(t,e){return this._point(o(this.chi,e,t.location))},t.prototype._contractInside=function(t,e){return this._point(o(this.psi,e,t.location))},t.prototype._contractOutside=function(t,e){return this._point(o(-this.psi*this.rho,e,t.location))},t.prototype._shrink=function(t,e){return this._point(o(-this.sigma,e,t.location))},t.prototype._isConverged=function(){for(var t=this._control.tolerance,e=this._simplex[0],r=this._simplex[this._n],n=r.value-e.value<t||1-e.value/r.value<t,i=0,o=0;o<this._n;++o)i=Math.max(i,Math.abs(e.location[o]-r.location[o]));return i<t&&n},t}();e.Simplex=a},2003:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.checkResult=void 0,e.checkResult=function(t){return"number"==typeof t?{data:null,value:t}:t}},2886:(t,e)=>{"use strict";function r(t,e){return new Error("Invalid control parameter: '".concat(t,"' ").concat(e))}function n(t,e){return void 0===t?e:t}Object.defineProperty(e,"__esModule",{value:!0}),e.dopriControl=void 0,e.dopriControl=function(t){void 0===t&&(t={});var e=1e-6,i=1e4,o=1e-6,a=1/0,s=1e-8,u=!1,c=0,l=1/0,f={atol:n(t.atol,e),maxSteps:n(t.maxSteps,i),rtol:n(t.rtol,o),stepSizeMax:n(t.stepSizeMax,a),stepSizeMin:n(t.stepSizeMin,s),stepSizeMinAllow:n(t.stepSizeMinAllow,u),stiffCheck:n(t.stiffCheck,c),tcrit:n(t.tcrit,l)};if(f.maxSteps<1)throw r("maxSteps","must be at least 1");if(f.atol<=0)throw r("atol","must be strictly positive");if(f.rtol<=0)throw r("rtol","must be strictly positive");return f}},8742:function(t,e,r){"use strict";var n,i=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.DDE=e.integrateDDE=void 0;var o=r(2223),a=r(4660);e.integrateDDE=function(t,e,r,n,i,o){void 0===i&&(i={}),void 0===o&&(o=null);var a=new s(t,e.length,i,o);return a.initialise(r,e),a.run(n)};var s=function(t){function e(e,r,n,i){void 0===n&&(n={}),void 0===i&&(i=null);var o=this,a=function(t){return o._interpolate(t)},s=null===i?null:function(t,e){return i(t,e,a)};return(o=t.call(this,(function(t,r,n){return e(t,r,n,a)}),r,n,s)||this)._y0=new Array(r),o}return i(e,t),e.prototype.initialise=function(e,r){return this._y0=r,t.prototype.initialise.call(this,e,r),this},e.prototype._interpolate=function(t){var e=this._findHistory(t);return e<0?this._y0.slice(0):this._stepper.interpolate(t,this._history[e])},e.prototype._findHistory=function(t){return(0,a.search)(this._history,(function(e){return e.t>t}))},e}(o.Dopri);e.DDE=s},2223:function(t,e,r){"use strict";var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var i=Object.getOwnPropertyDescriptor(e,r);i&&!("get"in i?!e.__esModule:i.writable||i.configurable)||(i={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,i)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),i=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),o=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return i(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri=e.integrateDopri=void 0;var a=r(2886),s=o(r(2922)),u=r(1216),c=o(r(4660)),l=Math.pow(2,-52);function f(t,e){return new Error("Integration failure: ".concat(t," at ").concat(e))}e.integrateDopri=function(t,e,r,n,i,o){void 0===i&&(i={}),void 0===o&&(o=null);var a=new h(t,e.length,i,o);return a.initialise(r,e),a.run(n)};var h=function(){function t(t,e,r,n){void 0===r&&(r={}),void 0===n&&(n=null),this._history=[],this._t=0,this._h=0,this._nSteps=0,this._nStepsAccepted=0,this._nStepsRejected=0,this._stiffNStiff=0,this._stiffNNonstiff=0,this._lastError=0,this._stepper=new s.Dopri5(t,e),this._control=(0,a.dopriControl)(r),this._output=n}return t.prototype.initialise=function(t,e){var r=this._stepper.n;if(e.length!==r)throw Error("Invalid size 'y' - expected a length ".concat(r," array"));return this._stepper.reset(t,e),this._reset(),this._h=function(t,e,r,n,i,o){var a=new Array(t.n),s=new Array(t.n),u=new Array(t.n);t.rhs(e,r,a),t.nEval++;var l=0,f=0,h=0;for(h=0;h<t.n;++h){var p=n+i*Math.abs(r[h]);l+=c.square(a[h]/p),f+=c.square(r[h]/p)}var d=l<=1e-10||f<=1e-10?1e-6:.01*Math.sqrt(f/l);for(d=Math.min(d,o),h=0;h<t.n;++h)u[h]=r[h]+d*a[h];t.rhs(e+d,u,s),t.nEval++;var v=0;for(h=0;h<t.n;++h)p=n+i*Math.abs(r[h]),v+=c.square((s[h]-a[h])/p);v=Math.sqrt(v)/d;var y=Math.max(Math.abs(v),Math.sqrt(l)),m=y<=1e-15?Math.max(1e-6,.001*Math.abs(d)):Math.pow(.01/y,1/t.order);return Math.min(Math.min(100*Math.abs(d),m),o)}(this._stepper,t,e,this._control.atol,this._control.rtol,this._control.stepSizeMax),this._t=t,this._history=[],this},t.prototype.run=function(t){for(;this._t<t;)this._step(),this._history.push(this._stepper.history.clone());return(0,u.interpolator)(this._history.slice(0),this._stepper,this._output)},t.prototype.statistics=function(){return{lastError:this._lastError,nEval:this._stepper.nEval,nSteps:this._nSteps,nStepsAccepted:this._nStepsAccepted,nStepsRejected:this._nStepsRejected,stiffNNonstiff:this._stiffNNonstiff,stiffNStiff:this._stiffNStiff}},t.prototype._reset=function(){this._nSteps=0,this._nStepsAccepted=0,this._nStepsRejected=0,this._stiffNStiff=0,this._stiffNNonstiff=0,this._lastError=0},t.prototype._step=function(){for(var t=this._t,e=this._h,r=!1,n=!1,i=Math.max(this._lastError,1e-4),o=this._stepper.stepControl,a=this._control;!r;){if(this._nSteps>a.maxSteps)throw f("too many steps",t);if(e<a.stepSizeMin){if(!a.stepSizeMinAllow)throw f("step too small",t);e=a.stepSizeMin}if(e<=Math.abs(t)*l)throw f("step size vanished",t);t+e>a.tcrit&&(e=a.tcrit-t),this._stepper.step(t,e),this._nSteps++;var s=this._stepper.error(a.atol,a.rtol),u=Math.pow(s,o.constant),h=1/o.factorMin,p=1/o.factorMax;if(s<=1){if(r=!0,this._nStepsAccepted++,this._isStiff(e))throw f("problem became stiff",t);this._stepper.stepComplete(t,e);var d=u/Math.pow(i,o.beta),v=e/(d=c.constrain(d/o.factorSafe,p,h));this._t+=e,this._h=n?Math.min(v,e):Math.min(v,a.stepSizeMax),this._lastError=s}else n=!0,this._nStepsAccepted>=1&&this._nStepsRejected++,e/=Math.min(h,u/o.factorSafe)}return this._t},t.prototype._isStiff=function(t){if(this._stiffNStiff>0||this._nStepsAccepted%this._control.stiffCheck==0)if(this._stepper.isStiff(t)){if(this._stiffNNonstiff=0,this._stiffNStiff++>=15)return!0}else this._stiffNStiff>0&&this._stiffNNonstiff++>=6&&(this._stiffNStiff=0,this._stiffNNonstiff=0);return!1},t}();e.Dopri=h},7134:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri5StepControl=void 0;e.Dopri5StepControl=function(){this.factorSafe=.9,this.factorMin=.2,this.factorMax=10,this.beta=.04,this.constant=.2-.75*this.beta}},2922:function(t,e,r){"use strict";var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var i=Object.getOwnPropertyDescriptor(e,r);i&&!("get"in i?!e.__esModule:i.writable||i.configurable)||(i={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,i)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),i=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),o=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return i(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri5=void 0;var a=o(r(4660)),s=o(r(7134)),u=r(8430),c=function(){function t(t,e){this.order=5,this.stepControl=new s.Dopri5StepControl,this.nEval=0,this.rhs=t,this.n=e,this.y=new Array(e),this.yNext=new Array(e),this.yStiff=new Array(e),this.k1=new Array(e),this.k2=new Array(e),this.k3=new Array(e),this.k4=new Array(e),this.k5=new Array(e),this.k6=new Array(e),this.history=new u.HistoryElement(this.order*e)}return t.prototype.step=function(t,e){var r=this.n,n=this.y,i=this.yNext,o=this.k1,a=this.k2,s=this.k3,u=this.k4,c=this.k5,l=this.k6,f=this.history.data,h=0;for(h=0;h<r;++h)i[h]=n[h]+.2*e*o[h];for(this.rhs(t+.2*e,i,a),h=0;h<r;++h)i[h]=n[h]+e*(.075*o[h]+.225*a[h]);for(this.rhs(t+.3*e,i,s),h=0;h<r;++h)i[h]=n[h]+e*(.9777777777777777*o[h]+-3.7333333333333334*a[h]+3.5555555555555554*s[h]);for(this.rhs(t+.8*e,i,u),h=0;h<r;++h)i[h]=n[h]+e*(2.9525986892242035*o[h]+-11.595793324188385*a[h]+9.822892851699436*s[h]+-.2908093278463649*u[h]);for(this.rhs(t+.8888888888888888*e,i,c),h=0;h<r;++h)this.yStiff[h]=n[h]+e*(2.8462752525252526*o[h]+-10.757575757575758*a[h]+8.906422717743473*s[h]+.2784090909090909*u[h]+-.2735313036020583*c[h]);var p=t+e;for(this.rhs(p,this.yStiff,l),h=0;h<r;++h)i[h]=n[h]+e*(.09114583333333333*o[h]+.44923629829290207*s[h]+.6510416666666666*u[h]+-.322376179245283*c[h]+.13095238095238096*l[h]);this.rhs(p,i,a);var d=4*r;for(h=0;h<r;++h)f[d++]=e*(-1.1270175653862835*o[h]+2.675424484351598*s[h]+-5.685526961588504*u[h]+3.5219323679207912*c[h]+-1.7672812570757455*l[h]+2.382468931778144*a[h]);for(h=0;h<r;++h)u[h]=e*(.0012326388888888888*o[h]+-.0042527702905061394*s[h]+.03697916666666667*u[h]+-.05086379716981132*c[h]+.0419047619047619*l[h]+-.025*a[h]);this.nEval+=6},t.prototype.stepComplete=function(t,e){this.saveHistory(t,e),a.copyArray(this.k1,this.k2),a.copyArray(this.y,this.yNext)},t.prototype.error=function(t,e){var r=0,n=0;for(n=0;n<this.n;++n){var i=t+e*Math.max(Math.abs(this.y[n]),Math.abs(this.yNext[n]));r+=a.square(this.k4[n]/i)}return Math.sqrt(r/this.n)},t.prototype.interpolate=function(t,e){for(var r=e.data,n=(t-e.t)/e.h,i=1-n,o=this.n,a=new Array(o),s=0;s<o;++s)a[s]=r[s]+n*(r[o+s]+i*(r[2*o+s]+n*(r[3*o+s]+i*r[4*o+s])));return a},t.prototype.isStiff=function(t){for(var e=0,r=0,n=0;n<this.n;++n)e+=a.square(this.k2[n]-this.k6[n]),r+=a.square(this.yNext[n]-this.yStiff[n]);return r>0&&Math.abs(t)*Math.sqrt(e/r)>3.25},t.prototype.reset=function(t,e){for(var r=0;r<this.n;++r)this.y[r]=e[r];this.rhs(t,e,this.k1),this.nEval=1},t.prototype.saveHistory=function(t,e){for(var r=this.history,n=this.n,i=0;i<n;++i){var o=this.yNext[i]-this.y[i],a=e*this.k1[i]-o;r.data[i]=this.y[i],r.data[n+i]=o,r.data[2*n+i]=a,r.data[3*n+i]=-e*this.k2[i]+o-a}r.t=t,r.h=e},t}();e.Dopri5=c},1606:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.HistoryElement=void 0;var r=function(){function t(t){this.t=0,this.h=0,this.data=new Array(t)}return t.prototype.clone=function(){var e=new t(this.data.length);return e.t=this.t,e.h=this.h,e.data=this.data.slice(0),e},t}();e.HistoryElement=r},6245:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.integrate=e.integrateDDE=e.DDE=e.integrateDopri=e.Dopri=void 0;var n=r(2223);Object.defineProperty(e,"Dopri",{enumerable:!0,get:function(){return n.Dopri}}),Object.defineProperty(e,"integrateDopri",{enumerable:!0,get:function(){return n.integrateDopri}});var i=r(8742);Object.defineProperty(e,"DDE",{enumerable:!0,get:function(){return i.DDE}}),Object.defineProperty(e,"integrateDDE",{enumerable:!0,get:function(){return i.integrateDDE}});var o=r(5961);Object.defineProperty(e,"integrate",{enumerable:!0,get:function(){return o.integrate}})},5961:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.integrate=void 0;var n=r(8742),i=r(2223);e.integrate=function(t,e,r,o,a,s){if(void 0===a&&(a={}),void 0===s&&(s=null),function(t){return 3===t.length}(t)){if(!function(t){return null===t||2===t.length}(s))throw new Error("Can't used delayed output with non-delayed rhs");return(0,i.integrateDopri)(t,e,r,o,a,s)}if(!function(t){return null===t||3===t.length}(s))throw new Error("Can't used non-delayed output with delayed rhs");return(0,n.integrateDDE)(t,e,r,o,a,s)}},1216:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.interpolator=void 0,e.interpolator=function(t,e,r){return function(n){return function(t,e,r,n){for(var i=[],o=e,a=0,s=0,u=t;s<u.length;s++){for(var c=u[s];o[a].t+o[a].h<c;)a++;var l=r.interpolate(c,o[a]);null!==n&&(l=l.concat(n(c,l))),i.push(l)}return i}(n,t,e,r)}}},8430:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.HistoryElement=void 0;var n=r(1606);Object.defineProperty(e,"HistoryElement",{enumerable:!0,get:function(){return n.HistoryElement}})},4660:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.search=e.seqLen=e.approxEqualArray=e.approxEqual=e.zeros=e.copyArray=e.constrain=e.square=void 0;var r=Math.pow(2,-26);e.square=function(t){return t*t},e.constrain=function(t,e,r){return Math.max(Math.min(t,r),e)},e.copyArray=function(t,e){for(var r=t.length,n=0;n<r;n++)t[n]=e[n]},e.zeros=function(t){for(var e=Array(t),r=0;r<t;++r)e[r]=0;return e},e.approxEqual=function(t,e,n){void 0===n&&(n=r);var i=Math.abs(t-e),o=Math.abs(t);return o>n&&(i/=o),i<n},e.approxEqualArray=function(t,e,n){if(void 0===n&&(n=r),e.length!==t.length)throw Error("Incompatible arrays");for(var i=0,o=0,a=0,s=0;s<t.length;++s)t[s]!==e[s]&&(i+=Math.abs(t[s]),o+=Math.abs(t[s]-e[s]),a++);return 0===a||(o/=a,(i/=a)>n&&(o/=i),o<n)},e.seqLen=function(t,e,r){for(var n=[],i=0;i<r;++i){var o=i/(r-1);n.push((1-o)*t+o*e)}return n},e.search=function(t,e){var r=0,n=t.length-1;if(0===t.length||e(t[r]))return-1;if(!e(t[n]))return n;for(;n-r>1;){var i=Math.floor((r+n)/2);e(t[i])?n=i:r=i}return r}},7788:(t,e)=>{var r,n=function(){"use strict";var t=new DataView(new ArrayBuffer(8)),e=function(e,r,n){var i,o,a,s=-.16666666666666632,u=0;return t.setFloat64(0,e),u=t.getInt32(0),(u&=2147483647)<1044381696&&0==(0|e)?e:(a=(i=e*e)*e,o=.00833333333332249+i*(i*(27557313707070068e-22+i*(1.58969099521155e-10*i-2.5050760253406863e-8))-.0001984126982985795),0===n?e+a*(s+i*o):e-(i*(.5*r-a*o)-r-a*s))},r=function(e,r){var n,i,o,a=0;return t.setFloat64(0,e),a=t.getInt32(0),(a&=2147483647)<1044381696&&0==(0|e)?1:(i=(n=e*e)*(.0416666666666666+n*(n*(2480158728947673e-20+n*(n*(2.087572321298175e-9+-11359647557788195e-27*n)-2.7557314351390663e-7))-.001388888888887411)),a<1070805811?1-(.5*n-(n*i-e*r)):(a>1072234496?o=.28125:(t.setFloat64(0,o),t.setUint32(0,a-2097152),t.setUint32(4,0),o=t.getFloat64(0)),1-o-(.5*n-o-(n*i-e*r))))},n=function(n){var i,o,a,s,u,c,l,f,h=1.4616321449683622,p=0,d=0;if(t.setFloat64(0,n),c=t.getInt32(0),l=t.getInt32(4),(f=2147483647&c)>=2146435072)return n*n;if(0==(f|l))return 1/0;if(f<999292928)return c<0?-Math.log(-n):-Math.log(n);if(c<0){if(f>=1127219200)return 1/0;if(i=function(n){var i,o=0,a=0;if(t.setFloat64(0,n),a=t.getInt32(0),(a&=2147483647)<1070596096)return e(Math.PI*n,0,0);switch(i=-n,Math.floor(i)!=i?o=4*(i=2*((i*=.5)-Math.floor(i)))|0:a>=1128267776?(i=0,o=0):(o=t.getInt32(4),i=o&=1,o<<=2),o){case 0:i=e(Math.PI*i,0,0);break;case 1:case 2:i=r(Math.PI*(.5-i),0);break;case 3:case 4:i=e(Math.PI*(1-i),0,0);break;case 5:case 6:i=-r(Math.PI*(i-1.5),0);break;default:i=e(Math.PI*(i-2),0,0)}return-i}(n),0==i)return 1/0;p=Math.log(3.141592653589793/Math.abs(i*n)),n=-n}if(0==(f-1072693248|l)||0==(f-1073741824|l))s=0;else if(f<1073741824)switch(f<=1072483532?(s=-Math.log(n),f>=1072130372?(o=1-n,d=0):f>=1070442081?(o=n-(h-1),d=1):(o=n,d=2)):(s=0,f>=1073460419?(o=2-n,d=0):f>=1072936132?(o=n-h,d=1):(o=n-1,d=2)),d){case 0:s+=o*(.07721566490153287+(a=o*o)*(.06735230105312927+a*(.007385550860814029+a*(.0011927076318336207+a*(.00022086279071390839+25214456545125733e-21*a)))))+a*(.3224670334241136+a*(.020580808432516733+a*(.0028905138367341563+a*(.0005100697921535113+a*(.00010801156724758394+44864094961891516e-21*a)))))-.5*o;break;case 1:s+=(a=o*o)*(.48383612272381005+(u=a*o)*(u*(.006100538702462913+u*(.00031563207090362595*u-.0014034646998923284))-.032788541075985965))-(-3638676997039505e-33-u*(u*(.01797067508118204+u*(u*(.000881081882437654+-.00031275416837512086*u)-.0036845201678113826))-.1475877229945939+o*(.06462494023913339+u*(u*(.0022596478090061247+u*(.0003355291926355191*u-.0005385953053567405))-.010314224129834144))))-.12148629053584961;break;case 2:s+=-.5*o+o*(o*(.6328270640250934+o*(1.4549225013723477+o*(.9777175279633727+o*(.22896372806469245+.013381091853678766*o))))-.07721566490153287)/(1+o*(2.4559779371304113+o*(2.128489763798934+o*(.7692851504566728+o*(.10422264559336913+.003217092422824239*o)))))}else if(f<1075838976)switch(i=0,s=.5*(o=n-(d=0|n))+o*(o*(.21498241596060885+o*(.325778796408931+o*(.14635047265246445+o*(.02664227030336386+o*(.0018402845140733772+3194753265841009e-20*o)))))-.07721566490153287)/(1+o*(1.3920053346762105+o*(.7219355475671381+o*(.17193386563280308+o*(.01864591917156529+o*(.0007779424963818936+7326684307446256e-21*o)))))),a=1,d){case 7:a*=o+6;case 6:a*=o+5;case 5:a*=o+4;case 4:a*=o+3;case 3:a*=o+2,s+=Math.log(a)}else s=f<1133510656?(n-.5)*((i=Math.log(n))-1)+(u=.4189385332046727+(a=1/n)*(.08333333333333297+(o=a*a)*(o*(.0007936505586430196+o*(o*(.0008363399189962821+-.0016309293409657527*o)-.00059518755745034))-.0027777777772877554))):n*(Math.log(n)-1);return c<0&&(s=p-s),s};return{lngamma:function(t){return n(t)},gamma:function(t){return Math.exp(n(t))}}}();void 0!==t.exports?t.exports=n:void 0===(r=function(){return n}.apply(e,[]))||(t.exports=r)},6907:t=>{"use strict";t.exports=function(t){for(var e=new Array(t),r=0;r<t;++r)e[r]=r;return e}},8738:t=>{function e(t){return!!t.constructor&&"function"==typeof t.constructor.isBuffer&&t.constructor.isBuffer(t)}t.exports=function(t){return null!=t&&(e(t)||function(t){return"function"==typeof t.readFloatLE&&"function"==typeof t.slice&&e(t.slice(0,0))}(t)||!!t._isBuffer)}},2861:(t,e,r)=>{var n=r(6907),i=r(8738),o="undefined"!=typeof Float64Array;function a(t,e){return t[0]-e[0]}function s(){var t,e=this.stride,r=new Array(e.length);for(t=0;t<r.length;++t)r[t]=[Math.abs(e[t]),t];r.sort(a);var n=new Array(r.length);for(t=0;t<n.length;++t)n[t]=r[t][1];return n}function u(t,e){var r=["View",e,"d",t].join("");e<0&&(r="View_Nil"+t);var i="generic"===t;if(-1===e){var o="function "+r+"(a){this.data=a;};var proto="+r+".prototype;proto.dtype='"+t+"';proto.index=function(){return -1};proto.size=0;proto.dimension=-1;proto.shape=proto.stride=proto.order=[];proto.lo=proto.hi=proto.transpose=proto.step=function(){return new "+r+"(this.data);};proto.get=proto.set=function(){};proto.pick=function(){return null};return function construct_"+r+"(a){return new "+r+"(a);}";return new Function(o)()}if(0===e)return o="function "+r+"(a,d) {this.data = a;this.offset = d};var proto="+r+".prototype;proto.dtype='"+t+"';proto.index=function(){return this.offset};proto.dimension=0;proto.size=1;proto.shape=proto.stride=proto.order=[];proto.lo=proto.hi=proto.transpose=proto.step=function "+r+"_copy() {return new "+r+"(this.data,this.offset)};proto.pick=function "+r+"_pick(){return TrivialArray(this.data);};proto.valueOf=proto.get=function "+r+"_get(){return "+(i?"this.data.get(this.offset)":"this.data[this.offset]")+"};proto.set=function "+r+"_set(v){return "+(i?"this.data.set(this.offset,v)":"this.data[this.offset]=v")+"};return function construct_"+r+"(a,b,c,d){return new "+r+"(a,d)}",new Function("TrivialArray",o)(c[t][0]);o=["'use strict'"];var a=n(e),u=a.map((function(t){return"i"+t})),l="this.offset+"+a.map((function(t){return"this.stride["+t+"]*i"+t})).join("+"),f=a.map((function(t){return"b"+t})).join(","),h=a.map((function(t){return"c"+t})).join(",");o.push("function "+r+"(a,"+f+","+h+",d){this.data=a","this.shape=["+f+"]","this.stride=["+h+"]","this.offset=d|0}","var proto="+r+".prototype","proto.dtype='"+t+"'","proto.dimension="+e),o.push("Object.defineProperty(proto,'size',{get:function "+r+"_size(){return "+a.map((function(t){return"this.shape["+t+"]"})).join("*"),"}})"),1===e?o.push("proto.order=[0]"):(o.push("Object.defineProperty(proto,'order',{get:"),e<4?(o.push("function "+r+"_order(){"),2===e?o.push("return (Math.abs(this.stride[0])>Math.abs(this.stride[1]))?[1,0]:[0,1]}})"):3===e&&o.push("var s0=Math.abs(this.stride[0]),s1=Math.abs(this.stride[1]),s2=Math.abs(this.stride[2]);if(s0>s1){if(s1>s2){return [2,1,0];}else if(s0>s2){return [1,2,0];}else{return [1,0,2];}}else if(s0>s2){return [2,0,1];}else if(s2>s1){return [0,1,2];}else{return [0,2,1];}}})")):o.push("ORDER})")),o.push("proto.set=function "+r+"_set("+u.join(",")+",v){"),i?o.push("return this.data.set("+l+",v)}"):o.push("return this.data["+l+"]=v}"),o.push("proto.get=function "+r+"_get("+u.join(",")+"){"),i?o.push("return this.data.get("+l+")}"):o.push("return this.data["+l+"]}"),o.push("proto.index=function "+r+"_index(",u.join(),"){return "+l+"}"),o.push("proto.hi=function "+r+"_hi("+u.join(",")+"){return new "+r+"(this.data,"+a.map((function(t){return["(typeof i",t,"!=='number'||i",t,"<0)?this.shape[",t,"]:i",t,"|0"].join("")})).join(",")+","+a.map((function(t){return"this.stride["+t+"]"})).join(",")+",this.offset)}");var p=a.map((function(t){return"a"+t+"=this.shape["+t+"]"})),d=a.map((function(t){return"c"+t+"=this.stride["+t+"]"}));o.push("proto.lo=function "+r+"_lo("+u.join(",")+"){var b=this.offset,d=0,"+p.join(",")+","+d.join(","));for(var v=0;v<e;++v)o.push("if(typeof i"+v+"==='number'&&i"+v+">=0){d=i"+v+"|0;b+=c"+v+"*d;a"+v+"-=d}");for(o.push("return new "+r+"(this.data,"+a.map((function(t){return"a"+t})).join(",")+","+a.map((function(t){return"c"+t})).join(",")+",b)}"),o.push("proto.step=function "+r+"_step("+u.join(",")+"){var "+a.map((function(t){return"a"+t+"=this.shape["+t+"]"})).join(",")+","+a.map((function(t){return"b"+t+"=this.stride["+t+"]"})).join(",")+",c=this.offset,d=0,ceil=Math.ceil"),v=0;v<e;++v)o.push("if(typeof i"+v+"==='number'){d=i"+v+"|0;if(d<0){c+=b"+v+"*(a"+v+"-1);a"+v+"=ceil(-a"+v+"/d)}else{a"+v+"=ceil(a"+v+"/d)}b"+v+"*=d}");o.push("return new "+r+"(this.data,"+a.map((function(t){return"a"+t})).join(",")+","+a.map((function(t){return"b"+t})).join(",")+",c)}");var y=new Array(e),m=new Array(e);for(v=0;v<e;++v)y[v]="a[i"+v+"]",m[v]="b[i"+v+"]";for(o.push("proto.transpose=function "+r+"_transpose("+u+"){"+u.map((function(t,e){return t+"=("+t+"===undefined?"+e+":"+t+"|0)"})).join(";"),"var a=this.shape,b=this.stride;return new "+r+"(this.data,"+y.join(",")+","+m.join(",")+",this.offset)}"),o.push("proto.pick=function "+r+"_pick("+u+"){var a=[],b=[],c=this.offset"),v=0;v<e;++v)o.push("if(typeof i"+v+"==='number'&&i"+v+">=0){c=(c+this.stride["+v+"]*i"+v+")|0}else{a.push(this.shape["+v+"]);b.push(this.stride["+v+"])}");return o.push("var ctor=CTOR_LIST[a.length+1];return ctor(this.data,a,b,c)}"),o.push("return function construct_"+r+"(data,shape,stride,offset){return new "+r+"(data,"+a.map((function(t){return"shape["+t+"]"})).join(",")+","+a.map((function(t){return"stride["+t+"]"})).join(",")+",offset)}"),new Function("CTOR_LIST","ORDER",o.join("\n"))(c[t],s)}var c={float32:[],float64:[],int8:[],int16:[],int32:[],uint8:[],uint16:[],uint32:[],array:[],uint8_clamped:[],bigint64:[],biguint64:[],buffer:[],generic:[]};t.exports=function(t,e,r,n){if(void 0===t)return(0,c.array[0])([]);"number"==typeof t&&(t=[t]),void 0===e&&(e=[t.length]);var a=e.length;if(void 0===r){r=new Array(a);for(var s=a-1,l=1;s>=0;--s)r[s]=l,l*=e[s]}if(void 0===n)for(n=0,s=0;s<a;++s)r[s]<0&&(n-=(e[s]-1)*r[s]);for(var f=function(t){if(i(t))return"buffer";if(o)switch(Object.prototype.toString.call(t)){case"[object Float64Array]":return"float64";case"[object Float32Array]":return"float32";case"[object Int8Array]":return"int8";case"[object Int16Array]":return"int16";case"[object Int32Array]":return"int32";case"[object Uint8Array]":return"uint8";case"[object Uint16Array]":return"uint16";case"[object Uint32Array]":return"uint32";case"[object Uint8ClampedArray]":return"uint8_clamped";case"[object BigInt64Array]":return"bigint64";case"[object BigUint64Array]":return"biguint64"}return Array.isArray(t)?"array":"generic"}(t),h=c[f];h.length<=a+1;)h.push(u(f,h.length-1));return(0,h[a+1])(t,e,r,n)}},8863:t=>{"use strict";t.exports=function(t,e,r,n,i){var o,a;for(o=1;o<t;o++){if(0===r[o-1])return console.log("tridiagonalSolve: failed due to lack of diagonal dominance"),!1;a=e[o]/r[o-1],r[o]-=a*n[o-1],i[o]-=a*i[o-1]}if(0===r[t-1])return console.log("tridiagonalSolve: failed due to singular matrix"),!1;for(i[t-1]/=r[t-1],o=t-2;o>=0;o--){if(0===r[o])return console.log("tridiagonalSolve: failed due to singular matrix"),!1;i[o]=(i[o]-n[o]*i[o+1])/r[o]}return!0}},5653:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.base=void 0;var n=r(1688),i=n.base.maths,o=n.base.user;e.base={maths:i,user:o}},9133:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.Dust=void 0;var n=r(6532),i=r(5653),o=r(217),a=r(2367),s=r(8177),u=function(){function t(t,e,r,a,s){if(r<=0)throw Error("Expected at least one particle");this._Model=t,this._particles=[],this._random=s||new n.Random(new n.RngStateBuiltin);for(var u=new t(i.base,e,"ignore"),c=0;c<r;++c)this._particles.push(new o.Particle(u,a))}return t.prototype.nState=function(){return this._particles[0].size},t.prototype.nParticles=function(){return this._particles.length},t.prototype.step=function(){return this._particles[0].step},t.prototype.info=function(){return this._particles[0].info()},t.prototype.setPars=function(t,e){var r=this,n=this.step(),a=(this.nState(),new this._Model(i.base,t,"ignore"));if(a.size()!==this.nState())throw Error("Particle produced unexpected state size");this.forEachParticle((function(t,i){var s=e?void 0:t.state();r._particles[i]=new o.Particle(a,n,s)}))},t.prototype.setStep=function(t){this._particles.forEach((function(e){e.step=t}))},t.prototype.setState=function(t){this.checkState(t),this.forEachParticle((function(e,r){e.setState(t[r])}))},t.prototype.run=function(t){var e=this;this._particles.forEach((function(r){return r.run(t,e._random)}))},t.prototype.simulate=function(t,e){for(var r=null===e?this.nState():e.length,n=this.nParticles(),i=t.length,o=(0,s.dustStateTime)(r,n,i),a=function(r){u.run(t[r]),u.forEachParticle((function(t,n){t.copyState(o.viewParticle(n,r),e)}))},u=this,c=0;c<i;++c)a(c);return o},t.prototype.state=function(t){var e=null===t?this.nState():t.length,r=this.nParticles(),n=(0,a.dustState)(e,r);return this.forEachParticle((function(e,r){e.copyState(n.viewParticle(r),t)})),n},t.prototype.reorder=function(t){var e=this,r=this.nParticles();if(t.length!==r)throw Error("Expected index to have length ".concat(r)+" but given ".concat(t.length));this.forEachParticle((function(n,i){var o=t[i];if(o<0||o>=r)throw Error("Expected index to be in [0, ".concat(r-1,"]")+" but given ".concat(o));n.setState(e._particles[o].state(),!0)})),this.forEachParticle((function(t,e){t.swap()}))},t.prototype.forEachParticle=function(t){this._particles.forEach(t)},t.prototype.checkState=function(t){if(t.length!==this.nParticles())throw Error("Invalid length state, expected ".concat(this.nParticles())+" but given ".concat(t.length));for(var e=0;e<t.length;++e)if(t[e].length!==this.nState())throw Error("Invalid length state for particle ".concat(e,",")+" expected ".concat(this.nState())+" but given ".concat(t[e].length))},t}();e.Dust=u},217:function(t,e){"use strict";var r=this&&this.__spreadArray||function(t,e,r){if(r||2===arguments.length)for(var n,i=0,o=e.length;i<o;i++)!n&&i in e||(n||(n=Array.prototype.slice.call(e,0,i)),n[i]=e[i]);return t.concat(n||Array.prototype.slice.call(e))};Object.defineProperty(e,"__esModule",{value:!0}),e.Particle=void 0;var n=function(){function t(t,e,n){this._model=t,this.step=e,this.size=t.size(),this._y=void 0===n?t.initial(this.step):r([],n,!0),this._yNext=Array(this.size).fill(0)}return t.prototype.info=function(){return this._model.info()},t.prototype.run=function(t,e){for(;this.step<t;)this._model.update(this.step++,this._y,this._yNext,e),this.swap()},t.prototype.state=function(){return this._y},t.prototype.copyState=function(t,e){var r=this._y;if(null===e)for(var n=0;n<r.length;++n)t.set(n,r[n]);else for(n=0;n<e.length;++n)t.set(n,r[e[n]])},t.prototype.setState=function(t,e){void 0===e&&(e=!1);for(var r=e?this._yNext:this._y,n=0;n<t.length;++n)r[n]=t[n]},t.prototype.swap=function(){var t;t=[this._yNext,this._y],this._y=t[0],this._yNext=t[1]},t}();e.Particle=n},586:(t,e,r)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.variableNames=e.PkgWrapper=void 0;var n=r(6532),i=r(5653),o=r(9133),a=r(9882),s=function(){function t(t,e,r,o){this.Model=t,this.pars=e,this.random=new n.Random(o||new n.RngStateBuiltin),this.model=new this.Model(i.base,e,r)}return t.random=function(t){return new n.Random(t)},t.prototype.setUser=function(t,e){this.pars=t,this.model=new this.Model(i.base,t,e)},t.prototype.initial=function(t){return this.model.initial(t)},t.prototype.update=function(t,e){var r=Array(e.length).fill(0);return this.model.update(t,e,r,this.random),r},t.prototype.getInternal=function(){return this.model.getInternal()},t.prototype.getMetadata=function(){var t=this.model.info();return{info:t,names:u(t),size:this.model.size()}},t.prototype.run=function(t,e){var r=t[0],n=new o.Dust(this.Model,this.pars,1,r,this.random);null!==e&&n.setState([e]);var i=n.simulate(t,null);return{size:this.model.size(),y:Array.from(i.state.data)}},t}();function u(t){for(var e=[],r=function(t){var r=t.dim,n=t.name;0===r.length?e.push(n):e.push.apply(e,(0,a.combinations)(r).map((function(t){return"".concat(n,"[").concat(t.join(","),"]")})))},n=0,i=t;n<i.length;n++)r(i[n]);return e}e.PkgWrapper=s,e.variableNames=u},8177:function(t,e,r){"use strict";var n=this&&this.__importDefault||function(t){return t&&t.__esModule?t:{default:t}};Object.defineProperty(e,"__esModule",{value:!0}),e.DustStateTime=e.dustStateTime=void 0;var i=n(r(2861)),o=r(2367);e.dustStateTime=function(t,e,r){var n=(0,i.default)(new Float64Array(t*e*r),[r,e,t]);return new a(n,t,e,r)};var a=function(){function t(t,e,r,n){this.nState=e,this.nParticles=r,this.nTime=n,this.state=t}return t.prototype.viewTime=function(t){var e=this.state.pick(t,null,null);return new o.DustState(e,this.nState,this.nParticles)},t.prototype.viewParticle=function(t,e){return this.state.pick(e,t,null)},t.prototype.viewState=function(t,e){return this.state.pick(e,null,t)},t.prototype.viewTrace=function(t,e){return this.state.pick(null,e,t)},t.prototype.getParticle=function(t,e){return(0,o.copyVector)(this.viewParticle(t,e))},t.prototype.getState=function(t,e){return(0,o.copyVector)(this.viewState(t,e))},t.prototype.getTrace=function(t,e){return(0,o.copyVector)(this.viewTrace(t,e))},t}();e.DustStateTime=a},2367:function(t,e,r){"use strict";var n=this&&this.__importDefault||function(t){return t&&t.__esModule?t:{default:t}};Object.defineProperty(e,"__esModule",{value:!0}),e.copyVector=e.DustState=e.dustState=void 0;var i=n(r(2861));e.dustState=function(t,e){var r=(0,i.default)(new Float64Array(t*e),[e,t]);return new o(r,t,e)};var o=function(){function t(t,e,r){this.nState=e,this.nParticles=r,this.state=t}return t.prototype.viewParticle=function(t){return this.state.pick(t,null)},t.prototype.viewState=function(t){return this.state.pick(null,t)},t.prototype.getParticle=function(t){return a(this.viewParticle(t))},t.prototype.getState=function(t){return a(this.viewState(t))},t.prototype.asMatrix=function(){for(var t=[],e=0;e<this.nParticles;++e)t.push(this.getParticle(e));return t},t}();function a(t){for(var e=t.shape[0],r=Array(e),n=0;n<e;++n)r[n]=t.get(n);return r}e.DustState=o,e.copyVector=a},9882:function(t,e){"use strict";var r=this&&this.__spreadArray||function(t,e,r){if(r||2===arguments.length)for(var n,i=0,o=e.length;i<o;i++)!n&&i in e||(n||(n=Array.prototype.slice.call(e,0,i)),n[i]=e[i]);return t.concat(n||Array.prototype.slice.call(e))};function n(t,e){for(var r=[],n=t;n<=e;++n)r.push(n);return r}function i(t,e){return t[0].map((function(r,n){return e(t.map((function(t){return t[n]})))}))}function o(t){return t.reduce((function(t,e){return t+e}),0)/t.length}Object.defineProperty(e,"__esModule",{value:!0}),e.findClosest=e.meanArray=e.mean=e.applyArray=e.isEqualArray=e.combinations=e.seqBy=e.seq=void 0,e.seq=n,e.seqBy=function(t,e,r){for(var n=[],i=t;i<=e;i+=r)n.push(i);return n},e.combinations=function t(e){var i=n(1,e[0]);if(e.length>1){for(var o=[],a=0,s=t(e.slice(1));a<s.length;a++)for(var u=s[a],c=0,l=i;c<l.length;c++){var f=l[c];o.push(r([f],u,!0))}return o}return i.map((function(t){return[t]}))},e.isEqualArray=function(t,e){return t.every((function(t,r){return t===e[r]}))},e.applyArray=i,e.mean=o,e.meanArray=function(t){return i(t,o)},e.findClosest=function(t,e){var r=e.findIndex((function(e){return e>=t}));if(e[r]===t||0===r)return r;var n=e[r-1],i=e[r];return t-n<i-t?r-1:r}},9556:(t,e)=>{"use strict";Object.defineProperty(e,"__esModule",{value:!0}),e.versions=void 0,e.versions=function(){return{random:"0.0.3",odinjs:"0.0.18",dust:"0.0.1"}}},7915:function(t,e,r){"use strict";var n=this&&this.__assign||function(){return n=Object.assign||function(t){for(var e,r=1,n=arguments.length;r<n;r++)for(var i in e=arguments[r])Object.prototype.hasOwnProperty.call(e,i)&&(t[i]=e[i]);return t},n.apply(this,arguments)};Object.defineProperty(e,"__esModule",{value:!0}),e.tidyDiscreteSolutionVariable=e.filterSolution=e.filterIndex=e.tidyDiscreteSolution=e.wodinRunDiscrete=e.runModelDiscrete=e.DiscreteSeriesMode=void 0;var i,o=r(1688),a=r(9133),s=r(9882);function u(t,e,r,n,i,o){var u=Math.floor(r/i),c=Math.ceil(n/i),l=new a.Dust(t,e,o,u),f=l.info(),h=(0,s.seq)(u,c),p=h.map((function(t){return t*i}));return{info:f,state:l.simulate(h,null),times:p}}function c(t){var e=t.times,r=[];return t.info.forEach((function(e){var n=h(e.name,t);r.push.apply(r,n)})),{x:e,values:r}}function l(t,e){if(e.mode===o.TimeMode.Grid){var r=(0,s.findClosest)(e.tStart,t),n=(0,s.findClosest)(e.tEnd,t),i=Math.max(1,Math.ceil((n-r)/e.nPoints));return(0,s.seqBy)(r,n,i)}return e.times.map((function(e){return(0,s.findClosest)(e,t)}))}function f(t){return function(e){var r=l(t.x,e),i=function(t){return r.map((function(e){return t[e]}))};return{x:i(t.x),values:t.values.map((function(t){return n(n({},t),{y:i(t.y)})}))}}}function h(t,e){for(var r=e.info.map((function(t){return t.name})).indexOf(t),n=e.state,o=n.getTrace(r,0),a=!1,u=[],c=0;c<n.nParticles;++c){var l=n.getTrace(r,c);u.push(l),!a&&c>1&&!(0,s.isEqualArray)(l,o)&&(a=!0)}if(a){var f=u.map((function(e){return{mode:i.Individual,name:t,y:e}}));return f.push({mode:i.Mean,name:t,y:(0,s.meanArray)(u)}),f}return[{mode:i.Deterministic,name:t,y:o}]}!function(t){t.Individual="Individual",t.Mean="Mean",t.Deterministic="Deterministic"}(i=e.DiscreteSeriesMode||(e.DiscreteSeriesMode={})),e.runModelDiscrete=u,e.wodinRunDiscrete=function(t,e,r,n,i,o){return f(c(u(t,e,r,n,i,o)))},e.tidyDiscreteSolution=c,e.filterIndex=l,e.filterSolution=f,e.tidyDiscreteSolutionVariable=h}},e={};function r(n){var i=e[n];if(void 0!==i)return i.exports;var o=e[n]={exports:{}};return t[n].call(o.exports,o,o.exports,r),o.exports}var n={};(()=>{"use strict";var t=n;Object.defineProperty(t,"__esModule",{value:!0}),t.wodinRunDiscrete=t.DiscreteSeriesMode=t.versions=t.DustStateTime=t.dustStateTime=t.DustState=t.dustState=t.PkgWrapper=t.Particle=t.Dust=t.base=void 0;var e=r(5653);Object.defineProperty(t,"base",{enumerable:!0,get:function(){return e.base}});var i=r(9133);Object.defineProperty(t,"Dust",{enumerable:!0,get:function(){return i.Dust}});var o=r(217);Object.defineProperty(t,"Particle",{enumerable:!0,get:function(){return o.Particle}});var a=r(586);Object.defineProperty(t,"PkgWrapper",{enumerable:!0,get:function(){return a.PkgWrapper}});var s=r(2367);Object.defineProperty(t,"dustState",{enumerable:!0,get:function(){return s.dustState}}),Object.defineProperty(t,"DustState",{enumerable:!0,get:function(){return s.DustState}});var u=r(8177);Object.defineProperty(t,"dustStateTime",{enumerable:!0,get:function(){return u.dustStateTime}}),Object.defineProperty(t,"DustStateTime",{enumerable:!0,get:function(){return u.DustStateTime}});var c=r(9556);Object.defineProperty(t,"versions",{enumerable:!0,get:function(){return c.versions}});var l=r(7915);Object.defineProperty(t,"DiscreteSeriesMode",{enumerable:!0,get:function(){return l.DiscreteSeriesMode}}),Object.defineProperty(t,"wodinRunDiscrete",{enumerable:!0,get:function(){return l.wodinRunDiscrete}})})(),dust=n})();
diff --git a/inst/js/odin.js b/inst/js/odin.js
index 60a9a50b..097afe8d 100644
--- a/inst/js/odin.js
+++ b/inst/js/odin.js
@@ -1 +1 @@
-var odinjs;(()=>{"use strict";var t={6:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.interpolateSearch=e.InterpolatorBase=void 0;var r=function(){function t(t,e){(function(t){return Array.isArray(t[0])})(e)||(e=[e]),this._i=0,this._x=t,this._y=e,this.nX=this._x.length,this.nY=this._y.length;for(var r=0;r<this.nY;++r)if(this._y[r].length!==this.nX)throw Error("Invalid length for 'y', expected ".concat(this.nX))}return t.prototype.evalAll=function(t){for(var e=Array(this.nY),r=0;r<this.nY;++r)e[r]=this.eval(t,r);return e},t.prototype.search=function(t,e){var r=n(t,this._x,this._i);if(r<0||!e&&r===this.nX)throw Error("Interpolation failed as 'x' is out of range");return this._i=r,r},t}();function n(t,e,r){var n=r,o=r,i=1,a=e.length;if(e[n]<=t){if(n>=a-1)return a;for(o=n+i;e[o]<t;)if(n=o,(o+=i*=2)>=a){if(e[o=a-1]<t)return a;break}}else{if(0===n)return-1;for(n-=i;e[n]>t;){if(o=n,n<(i*=2)){if(e[n=0]>t)return-1;break}n-=i}}for(;o-n>1;){var s=Math.floor((o+n)/2);e[s]<t?n=s:o=s}return n}e.InterpolatorBase=r,e.interpolateSearch=n},731:function(t,e,r){var n,o=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorConstant=void 0;var i=function(t){function e(){return null!==t&&t.apply(this,arguments)||this}return o(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=Math.min(this.search(t,!0),this.nX-1);return r!==this.nX-1&&this._x[r+1]===t&&++r,this._y[e][r]},e}(r(6).InterpolatorBase);e.InterpolatorConstant=i},673:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorSpline=e.InterpolatorLinear=e.InterpolatorConstant=e.InterpolatorBase=void 0;var n=r(6);Object.defineProperty(e,"InterpolatorBase",{enumerable:!0,get:function(){return n.InterpolatorBase}});var o=r(731);Object.defineProperty(e,"InterpolatorConstant",{enumerable:!0,get:function(){return o.InterpolatorConstant}});var i=r(477);Object.defineProperty(e,"InterpolatorLinear",{enumerable:!0,get:function(){return i.InterpolatorLinear}});var a=r(869);Object.defineProperty(e,"InterpolatorSpline",{enumerable:!0,get:function(){return a.InterpolatorSpline}})},477:function(t,e,r){var n,o=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorLinear=void 0;var i=function(t){function e(){return null!==t&&t.apply(this,arguments)||this}return o(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=this.search(t,!1),n=this._x[r],o=(t-n)/(this._x[r+1]-n),i=this._y[e],a=i[r];return a+(i[r+1]-a)*o},e}(r(6).InterpolatorBase);e.InterpolatorLinear=i},869:function(t,e,r){var n,o=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)}),i=this&&this.__importDefault||function(t){return t&&t.__esModule?t:{default:t}};Object.defineProperty(e,"__esModule",{value:!0}),e.solve=e.InterpolatorSpline=void 0;var a=i(r(863)),s=function(t){function e(e,r){var n=t.call(this,e,r)||this,o=function(t){var e=t.length,r=Array(e),n=Array(e),o=Array(e);r[0]=0,n[0]=2/(t[1]-t[0]),o[0]=1/(t[1]-t[0]);for(var i=1;i<e-1;++i){var a=t[i-1],s=t[i],u=t[i+1];r[i]=1/(s-a),n[i]=2*(1/(s-a)+1/(u-s)),o[i]=1/(u-s)}return r[e-1]=1/(t[e-1]-t[e-2]),n[e-1]=2/(t[e-1]-t[e-2]),o[e-1]=0,[r,n,o]}(n._x),i=function(t,e){for(var r=t.length,n=e.length,o=r-1,i=[],a=0;a<n;++a){var s=Array(r),u=e[a];s[0]=3*(u[1]-u[0])/Math.pow(t[1]-t[0],2);for(var l=1;l<o;++l){var c=t[l-1],f=t[l],h=t[l+1],p=u[l-1],d=u[l],v=u[l+1];s[l]=3*((d-p)/Math.pow(f-c,2)+(v-d)/Math.pow(h-f,2))}s[o]=3*(u[o]-u[o-1])/Math.pow(t[o]-t[o-1],2),i.push(s)}return i}(n._x,n._y);return n._k=function(t,e){for(var r=t[0],n=t[1],o=t[2],i=r.length,a=0,s=e;a<s.length;a++){var l=s[a];u(i,r,n.slice(),o,l)}return e}(o,i),n}return o(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=this.search(t,!1);return Array(this.nY),function(t,e,r,n,o){var i=(e-r[t])/(r[t+1]-r[t]),a=o[t]*(r[t+1]-r[t])-(n[t+1]-n[t]),s=-o[t+1]*(r[t+1]-r[t])+(n[t+1]-n[t]);return(1-i)*n[t]+i*n[t+1]+i*(1-i)*(a*(1-i)+s*i)}(r,t,this._x,this._y[e],this._k[e])},e}(r(6).InterpolatorBase);function u(t,e,r,n,o){if(!(0,a.default)(t,e,r,n,o))throw Error("solve failed: singular matrix?")}e.InterpolatorSpline=s,e.solve=u},394:(t,e)=>{function r(t,e){return void 0===t?e:t}Object.defineProperty(e,"__esModule",{value:!0}),e.simplexControl=void 0,e.simplexControl=function(t){void 0===t&&(t={});var e={deltaNonZero:r(t.deltaNonZero,.05),deltaZero:r(t.deltaZero,.001),errorOnFailure:r(t.errorOnFailure,!1),tolerance:r(t.tolerance,1e-5)};if(e.tolerance<=0)throw new Error("Invalid control parameter: 'tolerance' must be strictly positive");return e}},153:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0});var n=r(169);Object.defineProperty(e,"Simplex",{enumerable:!0,get:function(){return n.Simplex}}),Object.defineProperty(e,"fitSimplex",{enumerable:!0,get:function(){return n.fitSimplex}})},169:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.Simplex=e.fitSimplex=void 0;var n=r(394),o=r(3);function i(t,e,r){for(var n=e.length,o=new Array(n),i=1+t,a=-t,s=0;s<o.length;++s)o[s]=i*e[s]+a*r[s];return o}e.fitSimplex=function(t,e,r,n){return void 0===r&&(r={}),void 0===n&&(n=200),new a(t,e,r).run(n*e.length)};var a=function(){function t(t,e,r){void 0===r&&(r={}),this.rho=1,this.chi=2,this.psi=-.5,this.sigma=.5,this._iterations=0,this._id=0,this._converged=!1,this._target=t,this._control=n.simplexControl(r),this._n=e.length,this._simplex=[],this._simplex.push(this._point(e.slice())),r.errorOnFailure||(this._target=function(t){return function(e){try{return t(e)}catch(t){return 1/0}}}(t));for(var o=0;o<this._n;++o){var i=e.slice();i[o]?i[o]*=1+this._control.deltaNonZero:i[o]=this._control.deltaZero,this._simplex.push(this._point(i))}this._sort()}return t.prototype.step=function(){if(this._iterations++,this._isConverged())return this._converged=!0,!0;var t=this._n,e=this._simplex[0],r=this._simplex[t],n=this._centroid(),o=this._reflect(r,n);if(o.value<e.value){var i=this._expand(r,n);this._update(i.value<o.value?i:o)}else if(o.value>=this._simplex[t-1].value){var a=o.value>r.value?this._contractInside(r,n):this._contractOutside(r,n);if(a.value<r.value)this._update(a);else{for(var s=1;s<=t;++s)this._simplex[s]=this._shrink(this._simplex[s],e.location);this._sort()}}else this._update(o);return!1},t.prototype.run=function(t){if(!this._converged)for(var e=0;e<t&&!this.step();++e);return this.result()},t.prototype.result=function(){var t=this._simplex[0];return{converged:this._converged,data:t.data,evaluations:this._id,iterations:this._iterations,location:t.location,value:t.value}},t.prototype.simplex=function(){return this._simplex.map((function(t){return{location:t.location,value:t.value}}))},t.prototype._point=function(t){var e=o.checkResult(this._target(t));return{data:e.data,id:this._id++,location:t,value:e.value}},t.prototype._sort=function(){this._simplex.sort((function(t,e){return t.value-e.value}))},t.prototype._centroid=function(){for(var t=this._n,e=new Array(t),r=0;r<t;++r){e[r]=0;for(var n=0;n<t;++n)e[r]+=this._simplex[n].location[r];e[r]/=t}return e},t.prototype._update=function(t){this._simplex[this._n]=t,this._sort()},t.prototype._reflect=function(t,e){return this._point(i(this.rho,e,t.location))},t.prototype._expand=function(t,e){return this._point(i(this.chi,e,t.location))},t.prototype._contractInside=function(t,e){return this._point(i(this.psi,e,t.location))},t.prototype._contractOutside=function(t,e){return this._point(i(-this.psi*this.rho,e,t.location))},t.prototype._shrink=function(t,e){return this._point(i(-this.sigma,e,t.location))},t.prototype._isConverged=function(){for(var t=this._control.tolerance,e=this._simplex[0],r=this._simplex[this._n],n=r.value-e.value<t||1-e.value/r.value<t,o=0,i=0;i<this._n;++i)o=Math.max(o,Math.abs(e.location[i]-r.location[i]));return o<t&&n},t}();e.Simplex=a},3:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.checkResult=void 0,e.checkResult=function(t){return"number"==typeof t?{data:null,value:t}:t}},886:(t,e)=>{function r(t,e){return new Error("Invalid control parameter: '".concat(t,"' ").concat(e))}function n(t,e){return void 0===t?e:t}Object.defineProperty(e,"__esModule",{value:!0}),e.dopriControl=void 0,e.dopriControl=function(t){void 0===t&&(t={});var e=1e-6,o=1e4,i=1e-6,a=1/0,s=1e-8,u=!1,l=0,c=1/0,f={atol:n(t.atol,e),maxSteps:n(t.maxSteps,o),rtol:n(t.rtol,i),stepSizeMax:n(t.stepSizeMax,a),stepSizeMin:n(t.stepSizeMin,s),stepSizeMinAllow:n(t.stepSizeMinAllow,u),stiffCheck:n(t.stiffCheck,l),tcrit:n(t.tcrit,c)};if(f.maxSteps<1)throw r("maxSteps","must be at least 1");if(f.atol<=0)throw r("atol","must be strictly positive");if(f.rtol<=0)throw r("rtol","must be strictly positive");return f}},742:function(t,e,r){var n,o=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.DDE=e.integrateDDE=void 0;var i=r(223),a=r(660);e.integrateDDE=function(t,e,r,n,o,i){void 0===o&&(o={}),void 0===i&&(i=null);var a=new s(t,e.length,o,i);return a.initialise(r,e),a.run(n)};var s=function(t){function e(e,r,n,o){void 0===n&&(n={}),void 0===o&&(o=null);var i=this,a=function(t){return i._interpolate(t)},s=null===o?null:function(t,e){return o(t,e,a)};return(i=t.call(this,(function(t,r,n){return e(t,r,n,a)}),r,n,s)||this)._y0=new Array(r),i}return o(e,t),e.prototype.initialise=function(e,r){return this._y0=r,t.prototype.initialise.call(this,e,r),this},e.prototype._interpolate=function(t){var e=this._findHistory(t);return e<0?this._y0.slice(0):this._stepper.interpolate(t,this._history[e])},e.prototype._findHistory=function(t){return(0,a.search)(this._history,(function(e){return e.t>t}))},e}(i.Dopri);e.DDE=s},223:function(t,e,r){var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var o=Object.getOwnPropertyDescriptor(e,r);o&&!("get"in o?!e.__esModule:o.writable||o.configurable)||(o={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,o)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),o=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),i=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return o(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri=e.integrateDopri=void 0;var a=r(886),s=i(r(922)),u=r(216),l=i(r(660)),c=Math.pow(2,-52);function f(t,e){return new Error("Integration failure: ".concat(t," at ").concat(e))}e.integrateDopri=function(t,e,r,n,o,i){void 0===o&&(o={}),void 0===i&&(i=null);var a=new h(t,e.length,o,i);return a.initialise(r,e),a.run(n)};var h=function(){function t(t,e,r,n){void 0===r&&(r={}),void 0===n&&(n=null),this._history=[],this._t=0,this._h=0,this._nSteps=0,this._nStepsAccepted=0,this._nStepsRejected=0,this._stiffNStiff=0,this._stiffNNonstiff=0,this._lastError=0,this._stepper=new s.Dopri5(t,e),this._control=(0,a.dopriControl)(r),this._output=n}return t.prototype.initialise=function(t,e){var r=this._stepper.n;if(e.length!==r)throw Error("Invalid size 'y' - expected a length ".concat(r," array"));return this._stepper.reset(t,e),this._reset(),this._h=function(t,e,r,n,o,i){var a=new Array(t.n),s=new Array(t.n),u=new Array(t.n);t.rhs(e,r,a),t.nEval++;var c=0,f=0,h=0;for(h=0;h<t.n;++h){var p=n+o*Math.abs(r[h]);c+=l.square(a[h]/p),f+=l.square(r[h]/p)}var d=c<=1e-10||f<=1e-10?1e-6:.01*Math.sqrt(f/c);for(d=Math.min(d,i),h=0;h<t.n;++h)u[h]=r[h]+d*a[h];t.rhs(e+d,u,s),t.nEval++;var v=0;for(h=0;h<t.n;++h)p=n+o*Math.abs(r[h]),v+=l.square((s[h]-a[h])/p);v=Math.sqrt(v)/d;var _=Math.max(Math.abs(v),Math.sqrt(c)),y=_<=1e-15?Math.max(1e-6,.001*Math.abs(d)):Math.pow(.01/_,1/t.order);return Math.min(Math.min(100*Math.abs(d),y),i)}(this._stepper,t,e,this._control.atol,this._control.rtol,this._control.stepSizeMax),this._t=t,this._history=[],this},t.prototype.run=function(t){for(;this._t<t;)this._step(),this._history.push(this._stepper.history.clone());return(0,u.interpolator)(this._history.slice(0),this._stepper,this._output)},t.prototype.statistics=function(){return{lastError:this._lastError,nEval:this._stepper.nEval,nSteps:this._nSteps,nStepsAccepted:this._nStepsAccepted,nStepsRejected:this._nStepsRejected,stiffNNonstiff:this._stiffNNonstiff,stiffNStiff:this._stiffNStiff}},t.prototype._reset=function(){this._nSteps=0,this._nStepsAccepted=0,this._nStepsRejected=0,this._stiffNStiff=0,this._stiffNNonstiff=0,this._lastError=0},t.prototype._step=function(){for(var t=this._t,e=this._h,r=!1,n=!1,o=Math.max(this._lastError,1e-4),i=this._stepper.stepControl,a=this._control;!r;){if(this._nSteps>a.maxSteps)throw f("too many steps",t);if(e<a.stepSizeMin){if(!a.stepSizeMinAllow)throw f("step too small",t);e=a.stepSizeMin}if(e<=Math.abs(t)*c)throw f("step size vanished",t);t+e>a.tcrit&&(e=a.tcrit-t),this._stepper.step(t,e),this._nSteps++;var s=this._stepper.error(a.atol,a.rtol),u=Math.pow(s,i.constant),h=1/i.factorMin,p=1/i.factorMax;if(s<=1){if(r=!0,this._nStepsAccepted++,this._isStiff(e))throw f("problem became stiff",t);this._stepper.stepComplete(t,e);var d=u/Math.pow(o,i.beta),v=e/(d=l.constrain(d/i.factorSafe,p,h));this._t+=e,this._h=n?Math.min(v,e):Math.min(v,a.stepSizeMax),this._lastError=s}else n=!0,this._nStepsAccepted>=1&&this._nStepsRejected++,e/=Math.min(h,u/i.factorSafe)}return this._t},t.prototype._isStiff=function(t){if(this._stiffNStiff>0||this._nStepsAccepted%this._control.stiffCheck==0)if(this._stepper.isStiff(t)){if(this._stiffNNonstiff=0,this._stiffNStiff++>=15)return!0}else this._stiffNStiff>0&&this._stiffNNonstiff++>=6&&(this._stiffNStiff=0,this._stiffNNonstiff=0);return!1},t}();e.Dopri=h},167:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri5StepControl=void 0;e.Dopri5StepControl=function(){this.factorSafe=.9,this.factorMin=.2,this.factorMax=10,this.beta=.04,this.constant=.2-.75*this.beta}},922:function(t,e,r){var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var o=Object.getOwnPropertyDescriptor(e,r);o&&!("get"in o?!e.__esModule:o.writable||o.configurable)||(o={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,o)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),o=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),i=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return o(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri5=void 0;var a=i(r(660)),s=i(r(167)),u=r(430),l=function(){function t(t,e){this.order=5,this.stepControl=new s.Dopri5StepControl,this.nEval=0,this.rhs=t,this.n=e,this.y=new Array(e),this.yNext=new Array(e),this.yStiff=new Array(e),this.k1=new Array(e),this.k2=new Array(e),this.k3=new Array(e),this.k4=new Array(e),this.k5=new Array(e),this.k6=new Array(e),this.history=new u.HistoryElement(this.order*e)}return t.prototype.step=function(t,e){var r=this.n,n=this.y,o=this.yNext,i=this.k1,a=this.k2,s=this.k3,u=this.k4,l=this.k5,c=this.k6,f=this.history.data,h=0;for(h=0;h<r;++h)o[h]=n[h]+.2*e*i[h];for(this.rhs(t+.2*e,o,a),h=0;h<r;++h)o[h]=n[h]+e*(.075*i[h]+.225*a[h]);for(this.rhs(t+.3*e,o,s),h=0;h<r;++h)o[h]=n[h]+e*(.9777777777777777*i[h]+-3.7333333333333334*a[h]+3.5555555555555554*s[h]);for(this.rhs(t+.8*e,o,u),h=0;h<r;++h)o[h]=n[h]+e*(2.9525986892242035*i[h]+-11.595793324188385*a[h]+9.822892851699436*s[h]+-.2908093278463649*u[h]);for(this.rhs(t+.8888888888888888*e,o,l),h=0;h<r;++h)this.yStiff[h]=n[h]+e*(2.8462752525252526*i[h]+-10.757575757575758*a[h]+8.906422717743473*s[h]+.2784090909090909*u[h]+-.2735313036020583*l[h]);var p=t+e;for(this.rhs(p,this.yStiff,c),h=0;h<r;++h)o[h]=n[h]+e*(.09114583333333333*i[h]+.44923629829290207*s[h]+.6510416666666666*u[h]+-.322376179245283*l[h]+.13095238095238096*c[h]);this.rhs(p,o,a);var d=4*r;for(h=0;h<r;++h)f[d++]=e*(-1.1270175653862835*i[h]+2.675424484351598*s[h]+-5.685526961588504*u[h]+3.5219323679207912*l[h]+-1.7672812570757455*c[h]+2.382468931778144*a[h]);for(h=0;h<r;++h)u[h]=e*(.0012326388888888888*i[h]+-.0042527702905061394*s[h]+.03697916666666667*u[h]+-.05086379716981132*l[h]+.0419047619047619*c[h]+-.025*a[h]);this.nEval+=6},t.prototype.stepComplete=function(t,e){this.saveHistory(t,e),a.copyArray(this.k1,this.k2),a.copyArray(this.y,this.yNext)},t.prototype.error=function(t,e){var r=0,n=0;for(n=0;n<this.n;++n){var o=t+e*Math.max(Math.abs(this.y[n]),Math.abs(this.yNext[n]));r+=a.square(this.k4[n]/o)}return Math.sqrt(r/this.n)},t.prototype.interpolate=function(t,e){for(var r=e.data,n=(t-e.t)/e.h,o=1-n,i=this.n,a=new Array(i),s=0;s<i;++s)a[s]=r[s]+n*(r[i+s]+o*(r[2*i+s]+n*(r[3*i+s]+o*r[4*i+s])));return a},t.prototype.isStiff=function(t){for(var e=0,r=0,n=0;n<this.n;++n)e+=a.square(this.k2[n]-this.k6[n]),r+=a.square(this.yNext[n]-this.yStiff[n]);return r>0&&Math.abs(t)*Math.sqrt(e/r)>3.25},t.prototype.reset=function(t,e){for(var r=0;r<this.n;++r)this.y[r]=e[r];this.rhs(t,e,this.k1),this.nEval=1},t.prototype.saveHistory=function(t,e){for(var r=this.history,n=this.n,o=0;o<n;++o){var i=this.yNext[o]-this.y[o],a=e*this.k1[o]-i;r.data[o]=this.y[o],r.data[n+o]=i,r.data[2*n+o]=a,r.data[3*n+o]=-e*this.k2[o]+i-a}r.t=t,r.h=e},t}();e.Dopri5=l},606:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.HistoryElement=void 0;var r=function(){function t(t){this.t=0,this.h=0,this.data=new Array(t)}return t.prototype.clone=function(){var e=new t(this.data.length);return e.t=this.t,e.h=this.h,e.data=this.data.slice(0),e},t}();e.HistoryElement=r},245:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.integrate=e.integrateDDE=e.DDE=e.integrateDopri=e.Dopri=void 0;var n=r(223);Object.defineProperty(e,"Dopri",{enumerable:!0,get:function(){return n.Dopri}}),Object.defineProperty(e,"integrateDopri",{enumerable:!0,get:function(){return n.integrateDopri}});var o=r(742);Object.defineProperty(e,"DDE",{enumerable:!0,get:function(){return o.DDE}}),Object.defineProperty(e,"integrateDDE",{enumerable:!0,get:function(){return o.integrateDDE}});var i=r(961);Object.defineProperty(e,"integrate",{enumerable:!0,get:function(){return i.integrate}})},961:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.integrate=void 0;var n=r(742),o=r(223);e.integrate=function(t,e,r,i,a,s){if(void 0===a&&(a={}),void 0===s&&(s=null),function(t){return 3===t.length}(t)){if(!function(t){return null===t||2===t.length}(s))throw new Error("Can't used delayed output with non-delayed rhs");return(0,o.integrateDopri)(t,e,r,i,a,s)}if(!function(t){return null===t||3===t.length}(s))throw new Error("Can't used non-delayed output with delayed rhs");return(0,n.integrateDDE)(t,e,r,i,a,s)}},216:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.interpolator=void 0,e.interpolator=function(t,e,r){return function(n){return function(t,e,r,n){for(var o=[],i=e,a=0,s=0,u=t;s<u.length;s++){for(var l=u[s];i[a].t+i[a].h<l;)a++;var c=r.interpolate(l,i[a]);null!==n&&(c=c.concat(n(l,c))),o.push(c)}return o}(n,t,e,r)}}},430:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.HistoryElement=void 0;var n=r(606);Object.defineProperty(e,"HistoryElement",{enumerable:!0,get:function(){return n.HistoryElement}})},660:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.search=e.seqLen=e.approxEqualArray=e.approxEqual=e.zeros=e.copyArray=e.constrain=e.square=void 0;var r=Math.pow(2,-26);e.square=function(t){return t*t},e.constrain=function(t,e,r){return Math.max(Math.min(t,r),e)},e.copyArray=function(t,e){for(var r=t.length,n=0;n<r;n++)t[n]=e[n]},e.zeros=function(t){for(var e=Array(t),r=0;r<t;++r)e[r]=0;return e},e.approxEqual=function(t,e,n){void 0===n&&(n=r);var o=Math.abs(t-e),i=Math.abs(t);return i>n&&(o/=i),o<n},e.approxEqualArray=function(t,e,n){if(void 0===n&&(n=r),e.length!==t.length)throw Error("Incompatible arrays");for(var o=0,i=0,a=0,s=0;s<t.length;++s)t[s]!==e[s]&&(o+=Math.abs(t[s]),i+=Math.abs(t[s]-e[s]),a++);return 0===a||(i/=a,(o/=a)>n&&(i/=o),i<n)},e.seqLen=function(t,e,r){for(var n=[],o=0;o<r;++o){var i=o/(r-1);n.push((1-i)*t+i*e)}return n},e.search=function(t,e){var r=0,n=t.length-1;if(0===t.length||e(t[r]))return-1;if(!e(t[n]))return n;for(;n-r>1;){var o=Math.floor((r+n)/2);e(t[o])?n=o:r=o}return r}},863:t=>{t.exports=function(t,e,r,n,o){var i,a;for(i=1;i<t;i++){if(0===r[i-1])return console.log("tridiagonalSolve: failed due to lack of diagonal dominance"),!1;a=e[i]/r[i-1],r[i]-=a*n[i-1],o[i]-=a*o[i-1]}if(0===r[t-1])return console.log("tridiagonalSolve: failed due to singular matrix"),!1;for(o[t-1]/=r[t-1],i=t-2;i>=0;i--){if(0===r[i])return console.log("tridiagonalSolve: failed due to singular matrix"),!1;o[i]=(o[i]-n[i]*o[i+1])/r[i]}return!0}},653:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.base=void 0;const n=r(458),o=r(166),i=r(650),a=r(536);e.base={delay:n.delay,interpolate:o.interpolate,maths:i.maths,user:a.user}},401:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.updatePars=e.batchParsDisplace=e.batchParsRange=e.batchRun=e.Batch=void 0;const n=r(822),o=r(882),i=r(915);class a{constructor(t,e,r,n,o){this.pars=e,this.tStart=r,this.tEnd=n,this.solutions=e.values.map((a=>{const s=u(e.base,e.name,a);return(0,i.wodinRun)(t,s,r,n,o)}))}valueAtTime(t){const e=this.solutions.map((e=>e({mode:n.TimeMode.Given,times:[t]})));return{names:e[0].names,x:this.pars.values,y:e[0].names.map(((t,r)=>e.map((t=>t.y[r][0]))))}}extreme(t){return this.findExtremes()[t]}findExtremes(){if(void 0===this._extremes){const t=501,e={mode:n.TimeMode.Grid,nPoints:t,tEnd:this.tEnd,tStart:this.tStart},r=this.solutions.map((t=>t(e))),i=r[0].names,a=this.pars.values,s=(0,o.loop)(i.length,(t=>r.map((e=>function(t,e){const r=e.x,n=e.y[t],i=(0,o.whichMin)(n),a=(0,o.whichMax)(n),s=r[i],u=r[a],l=n[i];return{tMax:u,tMin:s,yMax:n[a],yMin:l}}(t,e)))));this._extremes={tMax:c("tMax",i,a,s),tMin:c("tMin",i,a,s),yMax:c("yMax",i,a,s),yMin:c("yMin",i,a,s)}}return this._extremes}}function s(t,e,r,n,i,a){const s=l(t,e);if(i>s)throw Error(`Expected lower bound to be no greater than ${s}`);if(a<s)throw Error(`Expected upper bound to be no less than ${s}`);if(i>=a)throw Error("Expected upper bound to be greater than lower bound");if(r<2)throw Error("Must include at least 2 traces in the batch");if(n&&i<=0)throw Error("Lower bound must be greater than 0 for logarithmic scale");return{base:t,name:e,values:n?(0,o.gridLog)(i,a,r):(0,o.grid)(i,a,r)}}function u(t,e,r){const n=Object.assign({},t);return n[e]=r,n}function l(t,e){const r=t[e];if(void 0===r)throw Error(`Expected a value for '${e}'`);if("number"!=typeof r)throw Error(`Expected a number for '${e}'`);return r}function c(t,e,r,n){return{names:e,x:r,y:(0,o.loop)(e.length,(e=>n[e].map((e=>e[t]))))}}e.Batch=a,e.batchRun=function(t,e,r,n,o){return new a(t,e,r,n,o)},e.batchParsRange=s,e.batchParsDisplace=function(t,e,r,n,o){const i=l(t,e),a=o/100;return s(t,e,r,n,i*(1-a),i*(1+a))},e.updatePars=u},458:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.delay=void 0,e.delay=function(t,e,r,n){const o=t(e);for(let t=0;t<r.length;++t)n[t]=o[r[t]]}},410:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.sumOfSquares=e.updatePars=e.fitTarget=void 0;const n=r(653),o=r(134),i=r(822);function a(t,e){const r=Object.assign({},t.base);for(let n=0;n<t.vary.length;++n)r[t.vary[n]]=e[n];return r}function s(t,e){let r=0;for(let n=0;n<t.length;++n){const o=t[n];isNaN(o)||(r+=Math.pow(o-e[n],2))}return r}e.fitTarget=function(t,e,r,u,l){if(e.time[0]<0)throw Error("Expected the first time to be at least 0");const c=e.time[e.time.length-1];return f=>{const h=a(r,f),p=new t(n.base,h,"error"),d=(0,o.runModel)(p,null,0,c,l).solution,v=p.names(),_=v.indexOf(u),y=d(e.time).map((t=>t[_]));return{data:{endTime:c,names:v,pars:h,solution:(0,i.interpolatedSolution)(d,v,0,c),startTime:0},value:s(e.value,y)}}},e.updatePars=a,e.sumOfSquares=s},166:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.interpolate=e.interpolateTimes=e.interpolateCheckT=e.interpolateCheckY=e.interpolateAlloc=void 0;const n=r(673);function o(t,e,r){const o=function(t,e){if(t.length===e)return[t];const r=[],n=t.length/e;for(let o=0;o<n;++o)r.push(t.slice(e*o,e*(o+1)));return r}(r,e.length);if("constant"===t)return new n.InterpolatorConstant(e,o);if("linear"===t)return new n.InterpolatorLinear(e,o);if("spline"===t)return new n.InterpolatorSpline(e,o);throw Error(`Invalid interpolation method '${t}'`)}function i(t,e,r,n){if(1===e.length){if(t[0]!==e[0])throw Error(`Expected ${r} to have length ${t[0]} (for ${n})`)}else for(let o=0;o<e.length;++o)if(t[o]!==e[o])throw Error(`Expected dimension ${o+1} of ${r} to have size ${t[o]} (for ${n})`)}function a(t,e){const r=Math.max(...t);return{max:Math.min(...e),min:r}}e.interpolateAlloc=o,e.interpolateCheckY=i,e.interpolateCheckT=function(t,e,r,n){if(void 0===r)return 1/0;if(t<r.min)throw Error(`Integration times do not span interpolation range; min: ${r.min}`);if(e>r.max)throw Error(`Integration times do not span interpolation range; max: ${r.max}`);return void 0===n&&(n=r.max),n},e.interpolateTimes=a,e.interpolate={alloc:o,checkY:i,times:a}},650:(t,e)=>{function r(t){return.5===o(t,1)?2*Math.round(t/2):Math.round(t)}function n(t,e){if(void 0===e||0===e)return r(t);{const n=Math.pow(10,e);return r(t*n)/n}}function o(t,e){let r=t%e;return r*e<0&&(r+=e),r}function i(t,e){return Math.floor(t/e)}function a(t,e,r){let n=0;for(let o=e;o<r;++o)n+=t[o];return n}function s(t,e,r,n,o,i){let a=0;for(let s=n;s<o;++s){const n=s*i;for(let o=e;o<r;++o)a+=t[o+n]}return a}function u(t,e,r,n,o,i,a,s,u){let l=0;for(let c=i;c<a;++c){const i=c*u;for(let a=n;a<o;++a){const n=a*s+i;for(let o=e;o<r;++o)l+=t[o+n]}}return l}Object.defineProperty(e,"__esModule",{value:!0}),e.maths=e.odinSum3=e.odinSum2=e.odinSum1=e.intdivr=e.modr=e.round2=void 0,e.round2=n,e.modr=o,e.intdivr=i,e.odinSum1=a,e.odinSum2=s,e.odinSum3=u,e.maths={intdivr:i,modr:o,odinSum1:a,odinSum2:s,odinSum3:u,round2:n}},134:function(t,e,r){var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var o=Object.getOwnPropertyDescriptor(e,r);o&&!("get"in o?!e.__esModule:o.writable||o.configurable)||(o={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,o)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),o=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),i=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return o(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.runModel=e.isODEModel=e.isDDEModel=void 0;const a=i(r(245)),s=r(166);function u(t){return 4===t.rhs.length}e.isDDEModel=u,e.isODEModel=function(t){return 3===t.rhs.length},e.runModel=function(t,e,r,n,o){const i=t.getMetadata().interpolateTimes;return o.tcrit=(0,s.interpolateCheckT)(r,n,i,o.tcrit),u(t)?function(t,e,r,n,o){let i=null;"function"==typeof t.output&&(i=(e,r,n)=>t.output(e,r,n)),null===e&&(e=t.initial(r)),t.getInternal().initial_t=r;const s=new a.DDE((function(e,r,n,o){t.rhs(e,r,n,o)}),e.length,o,i);return s.initialise(r,e),{solution:s.run(n),statistics:s.statistics()}}(t,e,r,n,o):function(t,e,r,n,o){let i=null;"function"==typeof t.output&&(i=(e,r)=>t.output(e,r)),null===e&&(e=t.initial(r));const s=new a.Dopri((function(e,r,n){t.rhs(e,r,n)}),e.length,o,i);return s.initialise(r,e),{solution:s.run(n),statistics:s.statistics()}}(t,e,r,n,o)}},586:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.PkgWrapper=void 0;const n=r(653),o=r(134);e.PkgWrapper=class{constructor(t,e,r){this.model=new t(n.base,e,r)}initial(t){return this.model.initial(t)}rhs(t,e){const r=new Array(e.length).fill(0);let n=null;if(!(0,o.isODEModel)(this.model))throw Error("Can't use rhs() with delay models");return this.model,this.model.rhs(t,e,r),this.model.output&&(n=this.model.output(t,e)),{output:n,state:r}}getMetadata(){return this.model.getMetadata()}getInternal(){return this.model.getInternal()}setUser(t,e){this.model.setUser(t,e)}run(t,e,r){const n=t[0],i=t[t.length-1],a=(0,o.runModel)(this.model,e,n,i,r);return{names:this.model.names(),statistics:a.statistics,y:a.solution(t)}}}},822:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.interpolatedSolution=e.TimeMode=void 0;const n=r(882);var o;!function(t){t.Grid="grid",t.Given="given"}(o=e.TimeMode||(e.TimeMode={})),e.interpolatedSolution=function(t,e,r,i){return a=>{const s=function(t,e,r){switch(t.mode){case o.Grid:return(0,n.grid)(Math.max(e,t.tStart),Math.min(r,t.tEnd),t.nPoints);case o.Given:return t.times}}(a,r,i),u=t(s),l=u[0].map(((t,e)=>u.map((t=>t[e]))));return{names:e,x:s,y:l}}}},536:(t,e)=>{function r(t,e,r){if("ignore"===r)return;const n=[];for(const r of Object.keys(t))e.includes(r)||n.push(r);if(n.length>0){const t="Unknown user parameters: "+n.join(", ");if("message"===r)console.log(t);else{if("warning"!==r)throw"stop"===r?Error(t):Error(t+" (and invalid value for unusedUserAction)");console.warn(t)}}}function n(t,e,r,n,o,i,a){const s=t[e];if(void 0===s){if(void 0!==r[e])return;if(null===n)throw Error(`Expected a value for '${e}'`);r[e]=n}else{if("number"!=typeof s)throw Error(`Expected a number for '${e}'`);l(s,o,i,a,e),r[e]=s}}function o(t,e,r,n,o,i,l){let c=t[e];if(void 0===c){if(void 0!==r[e])return;throw Error(`Expected a value for '${e}'`)}{const t=n.length-1;c=a(c),s(t,c,e),function(t,e,r){const n=t.length-1;for(let o=0;o<n;++o){const i=t[o+1];if(e.dim[o]!==i)throw 1===n?Error(`Expected length ${i} value for '${r}'`):Error(`Incorrect size of dimension ${o+1} of '${r}' (expected ${i})`)}}(n,c,e),u(c,o,i,l,e),r[e]=c.data.slice()}}function i(t,e,r,n,o,i,l){let c=t[e];if(void 0===c){if(void 0!==r[e])return;throw Error("Expected a value for '"+e+"'")}{const t=n.length-1;c=a(c),s(t,c,e),u(c,o,i,l,e),n[0]=c.data.length;for(let e=0;e<t;++e)n[e+1]=c.dim[e];r[e]=c.data.slice()}}function a(t,e){return Array.isArray(t)?t={data:t,dim:[t.length]}:"number"==typeof t&&(t={data:[t],dim:[1]}),t}function s(t,e,r){if(e.dim.length!==t)throw 1===t?Error(`Expected a numeric vector for '${r}'`):2===t?Error(`Expected a numeric matrix for '${r}'`):Error(`Expected a numeric array of rank ${t} for '${r}'`)}function u(t,e,r,n,o){for(const i of t.data){if(null===i)throw Error(`'${o}' must not contain any NA values`);l(i,e,r,n,o)}}function l(t,e,r,n,o){if("number"!=typeof t)throw Error(`Expected a number for '${o}'`);if(t<e)throw Error(`Expected '${o}' to be at least ${e}`);if(t>r)throw Error(`Expected '${o}' to be at most ${r}`);if(n&&!Number.isInteger(t))throw Error(`Expected '${o}' to be integer-like`)}Object.defineProperty(e,"__esModule",{value:!0}),e.user=e.setUserArrayVariable=e.setUserArrayFixed=e.setUserScalar=e.checkUser=void 0,e.checkUser=r,e.setUserScalar=n,e.setUserArrayFixed=o,e.setUserArrayVariable=i,e.user={checkUser:r,setUserArrayFixed:o,setUserArrayVariable:i,setUserScalar:n}},882:(t,e)=>{function r(t,e,r){const n=(e-t)/(r-1),o=[];for(let e=0;e<r-1;++e)o.push(t+e*n);return o.push(e),o}Object.defineProperty(e,"__esModule",{value:!0}),e.loop=e.whichMax=e.whichMin=e.gridLog=e.grid=void 0,e.grid=r,e.gridLog=function(t,e,n){return r(Math.log(t),Math.log(e),n).map(Math.exp)},e.whichMin=function(t){let e=-1,r=1/0;for(let n=0;n<t.length;++n)t[n]<r&&(e=n,r=t[n]);return e},e.whichMax=function(t){let e=-1,r=-1/0;for(let n=0;n<t.length;++n)t[n]>r&&(e=n,r=t[n]);return e},e.loop=function(t,e){const r=[];for(let n=0;n<t;++n)r.push(e(n));return r}},556:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.versions=void 0,e.versions=function(){return{dfoptim:"0.0.5",dopri:"0.0.12",odinjs:"0.0.15"}}},915:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.wodinFitValue=e.wodinFit=e.wodinRun=void 0;const n=r(153),o=r(653),i=r(410),a=r(134),s=r(822);e.wodinRun=function(t,e,r,n,i){const u=new t(o.base,e,"error"),l=(0,a.runModel)(u,null,r,n,i).solution,c=u.names();return(0,s.interpolatedSolution)(l,c,r,n)},e.wodinFit=function(t,e,r,o,a,s){const u=(0,i.fitTarget)(t,e,r,o,a),l=r.vary.map((t=>r.base[t]));return new n.Simplex(u,l,s)},e.wodinFitValue=function(t,e,r){const{names:n,y:o}=t({mode:s.TimeMode.Given,times:e.time}),a=n.indexOf(r);return(0,i.sumOfSquares)(e.value,o[a])}}},e={};function r(n){var o=e[n];if(void 0!==o)return o.exports;var i=e[n]={exports:{}};return t[n].call(i.exports,i,i.exports,r),i.exports}var n={};(()=>{var t=n;Object.defineProperty(t,"__esModule",{value:!0}),t.versions=t.batchRun=t.batchParsRange=t.batchParsDisplace=t.Batch=t.base=t.PkgWrapper=t.wodinRun=t.wodinFitValue=t.wodinFit=void 0;var e=r(915);Object.defineProperty(t,"wodinFit",{enumerable:!0,get:function(){return e.wodinFit}}),Object.defineProperty(t,"wodinFitValue",{enumerable:!0,get:function(){return e.wodinFitValue}}),Object.defineProperty(t,"wodinRun",{enumerable:!0,get:function(){return e.wodinRun}});var o=r(586);Object.defineProperty(t,"PkgWrapper",{enumerable:!0,get:function(){return o.PkgWrapper}});var i=r(653);Object.defineProperty(t,"base",{enumerable:!0,get:function(){return i.base}});var a=r(401);Object.defineProperty(t,"Batch",{enumerable:!0,get:function(){return a.Batch}}),Object.defineProperty(t,"batchParsDisplace",{enumerable:!0,get:function(){return a.batchParsDisplace}}),Object.defineProperty(t,"batchParsRange",{enumerable:!0,get:function(){return a.batchParsRange}}),Object.defineProperty(t,"batchRun",{enumerable:!0,get:function(){return a.batchRun}});var s=r(556);Object.defineProperty(t,"versions",{enumerable:!0,get:function(){return s.versions}})})(),odinjs=n})();
+var odinjs;(()=>{"use strict";var t={6:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.interpolateSearch=e.InterpolatorBase=void 0;var r=function(){function t(t,e){(function(t){return Array.isArray(t[0])})(e)||(e=[e]),this._i=0,this._x=t,this._y=e,this.nX=this._x.length,this.nY=this._y.length;for(var r=0;r<this.nY;++r)if(this._y[r].length!==this.nX)throw Error("Invalid length for 'y', expected ".concat(this.nX))}return t.prototype.evalAll=function(t){for(var e=Array(this.nY),r=0;r<this.nY;++r)e[r]=this.eval(t,r);return e},t.prototype.search=function(t,e){var r=n(t,this._x,this._i);if(r<0||!e&&r===this.nX)throw Error("Interpolation failed as 'x' is out of range");return this._i=r,r},t}();function n(t,e,r){var n=r,o=r,i=1,a=e.length;if(e[n]<=t){if(n>=a-1)return a;for(o=n+i;e[o]<t;)if(n=o,(o+=i*=2)>=a){if(e[o=a-1]<t)return a;break}}else{if(0===n)return-1;for(n-=i;e[n]>t;){if(o=n,n<(i*=2)){if(e[n=0]>t)return-1;break}n-=i}}for(;o-n>1;){var s=Math.floor((o+n)/2);e[s]<t?n=s:o=s}return n}e.InterpolatorBase=r,e.interpolateSearch=n},731:function(t,e,r){var n,o=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorConstant=void 0;var i=function(t){function e(){return null!==t&&t.apply(this,arguments)||this}return o(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=Math.min(this.search(t,!0),this.nX-1);return r!==this.nX-1&&this._x[r+1]===t&&++r,this._y[e][r]},e}(r(6).InterpolatorBase);e.InterpolatorConstant=i},673:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorSpline=e.InterpolatorLinear=e.InterpolatorConstant=e.InterpolatorBase=void 0;var n=r(6);Object.defineProperty(e,"InterpolatorBase",{enumerable:!0,get:function(){return n.InterpolatorBase}});var o=r(731);Object.defineProperty(e,"InterpolatorConstant",{enumerable:!0,get:function(){return o.InterpolatorConstant}});var i=r(477);Object.defineProperty(e,"InterpolatorLinear",{enumerable:!0,get:function(){return i.InterpolatorLinear}});var a=r(869);Object.defineProperty(e,"InterpolatorSpline",{enumerable:!0,get:function(){return a.InterpolatorSpline}})},477:function(t,e,r){var n,o=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.InterpolatorLinear=void 0;var i=function(t){function e(){return null!==t&&t.apply(this,arguments)||this}return o(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=this.search(t,!1),n=this._x[r],o=(t-n)/(this._x[r+1]-n),i=this._y[e],a=i[r];return a+(i[r+1]-a)*o},e}(r(6).InterpolatorBase);e.InterpolatorLinear=i},869:function(t,e,r){var n,o=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)}),i=this&&this.__importDefault||function(t){return t&&t.__esModule?t:{default:t}};Object.defineProperty(e,"__esModule",{value:!0}),e.solve=e.InterpolatorSpline=void 0;var a=i(r(863)),s=function(t){function e(e,r){var n=t.call(this,e,r)||this,o=function(t){var e=t.length,r=Array(e),n=Array(e),o=Array(e);r[0]=0,n[0]=2/(t[1]-t[0]),o[0]=1/(t[1]-t[0]);for(var i=1;i<e-1;++i){var a=t[i-1],s=t[i],u=t[i+1];r[i]=1/(s-a),n[i]=2*(1/(s-a)+1/(u-s)),o[i]=1/(u-s)}return r[e-1]=1/(t[e-1]-t[e-2]),n[e-1]=2/(t[e-1]-t[e-2]),o[e-1]=0,[r,n,o]}(n._x),i=function(t,e){for(var r=t.length,n=e.length,o=r-1,i=[],a=0;a<n;++a){var s=Array(r),u=e[a];s[0]=3*(u[1]-u[0])/Math.pow(t[1]-t[0],2);for(var l=1;l<o;++l){var c=t[l-1],f=t[l],h=t[l+1],p=u[l-1],d=u[l],v=u[l+1];s[l]=3*((d-p)/Math.pow(f-c,2)+(v-d)/Math.pow(h-f,2))}s[o]=3*(u[o]-u[o-1])/Math.pow(t[o]-t[o-1],2),i.push(s)}return i}(n._x,n._y);return n._k=function(t,e){for(var r=t[0],n=t[1],o=t[2],i=r.length,a=0,s=e;a<s.length;a++){var l=s[a];u(i,r,n.slice(),o,l)}return e}(o,i),n}return o(e,t),e.prototype.eval=function(t,e){void 0===e&&(e=0);var r=this.search(t,!1);return Array(this.nY),function(t,e,r,n,o){var i=(e-r[t])/(r[t+1]-r[t]),a=o[t]*(r[t+1]-r[t])-(n[t+1]-n[t]),s=-o[t+1]*(r[t+1]-r[t])+(n[t+1]-n[t]);return(1-i)*n[t]+i*n[t+1]+i*(1-i)*(a*(1-i)+s*i)}(r,t,this._x,this._y[e],this._k[e])},e}(r(6).InterpolatorBase);function u(t,e,r,n,o){if(!(0,a.default)(t,e,r,n,o))throw Error("solve failed: singular matrix?")}e.InterpolatorSpline=s,e.solve=u},394:(t,e)=>{function r(t,e){return void 0===t?e:t}Object.defineProperty(e,"__esModule",{value:!0}),e.simplexControl=void 0,e.simplexControl=function(t){void 0===t&&(t={});var e={deltaNonZero:r(t.deltaNonZero,.05),deltaZero:r(t.deltaZero,.001),errorOnFailure:r(t.errorOnFailure,!1),tolerance:r(t.tolerance,1e-5)};if(e.tolerance<=0)throw new Error("Invalid control parameter: 'tolerance' must be strictly positive");return e}},153:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0});var n=r(169);Object.defineProperty(e,"Simplex",{enumerable:!0,get:function(){return n.Simplex}}),Object.defineProperty(e,"fitSimplex",{enumerable:!0,get:function(){return n.fitSimplex}})},169:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.Simplex=e.fitSimplex=void 0;var n=r(394),o=r(3);function i(t,e,r){for(var n=e.length,o=new Array(n),i=1+t,a=-t,s=0;s<o.length;++s)o[s]=i*e[s]+a*r[s];return o}e.fitSimplex=function(t,e,r,n){return void 0===r&&(r={}),void 0===n&&(n=200),new a(t,e,r).run(n*e.length)};var a=function(){function t(t,e,r){void 0===r&&(r={}),this.rho=1,this.chi=2,this.psi=-.5,this.sigma=.5,this._iterations=0,this._id=0,this._converged=!1,this._target=t,this._control=n.simplexControl(r),this._n=e.length,this._simplex=[],this._simplex.push(this._point(e.slice())),r.errorOnFailure||(this._target=function(t){return function(e){try{return t(e)}catch(t){return 1/0}}}(t));for(var o=0;o<this._n;++o){var i=e.slice();i[o]?i[o]*=1+this._control.deltaNonZero:i[o]=this._control.deltaZero,this._simplex.push(this._point(i))}this._sort()}return t.prototype.step=function(){if(this._iterations++,this._isConverged())return this._converged=!0,!0;var t=this._n,e=this._simplex[0],r=this._simplex[t],n=this._centroid(),o=this._reflect(r,n);if(o.value<e.value){var i=this._expand(r,n);this._update(i.value<o.value?i:o)}else if(o.value>=this._simplex[t-1].value){var a=o.value>r.value?this._contractInside(r,n):this._contractOutside(r,n);if(a.value<r.value)this._update(a);else{for(var s=1;s<=t;++s)this._simplex[s]=this._shrink(this._simplex[s],e.location);this._sort()}}else this._update(o);return!1},t.prototype.run=function(t){if(!this._converged)for(var e=0;e<t&&!this.step();++e);return this.result()},t.prototype.result=function(){var t=this._simplex[0];return{converged:this._converged,data:t.data,evaluations:this._id,iterations:this._iterations,location:t.location,value:t.value}},t.prototype.simplex=function(){return this._simplex.map((function(t){return{location:t.location,value:t.value}}))},t.prototype._point=function(t){var e=o.checkResult(this._target(t));return{data:e.data,id:this._id++,location:t,value:e.value}},t.prototype._sort=function(){this._simplex.sort((function(t,e){return t.value-e.value}))},t.prototype._centroid=function(){for(var t=this._n,e=new Array(t),r=0;r<t;++r){e[r]=0;for(var n=0;n<t;++n)e[r]+=this._simplex[n].location[r];e[r]/=t}return e},t.prototype._update=function(t){this._simplex[this._n]=t,this._sort()},t.prototype._reflect=function(t,e){return this._point(i(this.rho,e,t.location))},t.prototype._expand=function(t,e){return this._point(i(this.chi,e,t.location))},t.prototype._contractInside=function(t,e){return this._point(i(this.psi,e,t.location))},t.prototype._contractOutside=function(t,e){return this._point(i(-this.psi*this.rho,e,t.location))},t.prototype._shrink=function(t,e){return this._point(i(-this.sigma,e,t.location))},t.prototype._isConverged=function(){for(var t=this._control.tolerance,e=this._simplex[0],r=this._simplex[this._n],n=r.value-e.value<t||1-e.value/r.value<t,o=0,i=0;i<this._n;++i)o=Math.max(o,Math.abs(e.location[i]-r.location[i]));return o<t&&n},t}();e.Simplex=a},3:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.checkResult=void 0,e.checkResult=function(t){return"number"==typeof t?{data:null,value:t}:t}},886:(t,e)=>{function r(t,e){return new Error("Invalid control parameter: '".concat(t,"' ").concat(e))}function n(t,e){return void 0===t?e:t}Object.defineProperty(e,"__esModule",{value:!0}),e.dopriControl=void 0,e.dopriControl=function(t){void 0===t&&(t={});var e=1e-6,o=1e4,i=1e-6,a=1/0,s=1e-8,u=!1,l=0,c=1/0,f={atol:n(t.atol,e),maxSteps:n(t.maxSteps,o),rtol:n(t.rtol,i),stepSizeMax:n(t.stepSizeMax,a),stepSizeMin:n(t.stepSizeMin,s),stepSizeMinAllow:n(t.stepSizeMinAllow,u),stiffCheck:n(t.stiffCheck,l),tcrit:n(t.tcrit,c)};if(f.maxSteps<1)throw r("maxSteps","must be at least 1");if(f.atol<=0)throw r("atol","must be strictly positive");if(f.rtol<=0)throw r("rtol","must be strictly positive");return f}},742:function(t,e,r){var n,o=this&&this.__extends||(n=function(t,e){return n=Object.setPrototypeOf||{__proto__:[]}instanceof Array&&function(t,e){t.__proto__=e}||function(t,e){for(var r in e)Object.prototype.hasOwnProperty.call(e,r)&&(t[r]=e[r])},n(t,e)},function(t,e){if("function"!=typeof e&&null!==e)throw new TypeError("Class extends value "+String(e)+" is not a constructor or null");function r(){this.constructor=t}n(t,e),t.prototype=null===e?Object.create(e):(r.prototype=e.prototype,new r)});Object.defineProperty(e,"__esModule",{value:!0}),e.DDE=e.integrateDDE=void 0;var i=r(223),a=r(660);e.integrateDDE=function(t,e,r,n,o,i){void 0===o&&(o={}),void 0===i&&(i=null);var a=new s(t,e.length,o,i);return a.initialise(r,e),a.run(n)};var s=function(t){function e(e,r,n,o){void 0===n&&(n={}),void 0===o&&(o=null);var i=this,a=function(t){return i._interpolate(t)},s=null===o?null:function(t,e){return o(t,e,a)};return(i=t.call(this,(function(t,r,n){return e(t,r,n,a)}),r,n,s)||this)._y0=new Array(r),i}return o(e,t),e.prototype.initialise=function(e,r){return this._y0=r,t.prototype.initialise.call(this,e,r),this},e.prototype._interpolate=function(t){var e=this._findHistory(t);return e<0?this._y0.slice(0):this._stepper.interpolate(t,this._history[e])},e.prototype._findHistory=function(t){return(0,a.search)(this._history,(function(e){return e.t>t}))},e}(i.Dopri);e.DDE=s},223:function(t,e,r){var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var o=Object.getOwnPropertyDescriptor(e,r);o&&!("get"in o?!e.__esModule:o.writable||o.configurable)||(o={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,o)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),o=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),i=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return o(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri=e.integrateDopri=void 0;var a=r(886),s=i(r(922)),u=r(216),l=i(r(660)),c=Math.pow(2,-52);function f(t,e){return new Error("Integration failure: ".concat(t," at ").concat(e))}e.integrateDopri=function(t,e,r,n,o,i){void 0===o&&(o={}),void 0===i&&(i=null);var a=new h(t,e.length,o,i);return a.initialise(r,e),a.run(n)};var h=function(){function t(t,e,r,n){void 0===r&&(r={}),void 0===n&&(n=null),this._history=[],this._t=0,this._h=0,this._nSteps=0,this._nStepsAccepted=0,this._nStepsRejected=0,this._stiffNStiff=0,this._stiffNNonstiff=0,this._lastError=0,this._stepper=new s.Dopri5(t,e),this._control=(0,a.dopriControl)(r),this._output=n}return t.prototype.initialise=function(t,e){var r=this._stepper.n;if(e.length!==r)throw Error("Invalid size 'y' - expected a length ".concat(r," array"));return this._stepper.reset(t,e),this._reset(),this._h=function(t,e,r,n,o,i){var a=new Array(t.n),s=new Array(t.n),u=new Array(t.n);t.rhs(e,r,a),t.nEval++;var c=0,f=0,h=0;for(h=0;h<t.n;++h){var p=n+o*Math.abs(r[h]);c+=l.square(a[h]/p),f+=l.square(r[h]/p)}var d=c<=1e-10||f<=1e-10?1e-6:.01*Math.sqrt(f/c);for(d=Math.min(d,i),h=0;h<t.n;++h)u[h]=r[h]+d*a[h];t.rhs(e+d,u,s),t.nEval++;var v=0;for(h=0;h<t.n;++h)p=n+o*Math.abs(r[h]),v+=l.square((s[h]-a[h])/p);v=Math.sqrt(v)/d;var _=Math.max(Math.abs(v),Math.sqrt(c)),y=_<=1e-15?Math.max(1e-6,.001*Math.abs(d)):Math.pow(.01/_,1/t.order);return Math.min(Math.min(100*Math.abs(d),y),i)}(this._stepper,t,e,this._control.atol,this._control.rtol,this._control.stepSizeMax),this._t=t,this._history=[],this},t.prototype.run=function(t){for(;this._t<t;)this._step(),this._history.push(this._stepper.history.clone());return(0,u.interpolator)(this._history.slice(0),this._stepper,this._output)},t.prototype.statistics=function(){return{lastError:this._lastError,nEval:this._stepper.nEval,nSteps:this._nSteps,nStepsAccepted:this._nStepsAccepted,nStepsRejected:this._nStepsRejected,stiffNNonstiff:this._stiffNNonstiff,stiffNStiff:this._stiffNStiff}},t.prototype._reset=function(){this._nSteps=0,this._nStepsAccepted=0,this._nStepsRejected=0,this._stiffNStiff=0,this._stiffNNonstiff=0,this._lastError=0},t.prototype._step=function(){for(var t=this._t,e=this._h,r=!1,n=!1,o=Math.max(this._lastError,1e-4),i=this._stepper.stepControl,a=this._control;!r;){if(this._nSteps>a.maxSteps)throw f("too many steps",t);if(e<a.stepSizeMin){if(!a.stepSizeMinAllow)throw f("step too small",t);e=a.stepSizeMin}if(e<=Math.abs(t)*c)throw f("step size vanished",t);t+e>a.tcrit&&(e=a.tcrit-t),this._stepper.step(t,e),this._nSteps++;var s=this._stepper.error(a.atol,a.rtol),u=Math.pow(s,i.constant),h=1/i.factorMin,p=1/i.factorMax;if(s<=1){if(r=!0,this._nStepsAccepted++,this._isStiff(e))throw f("problem became stiff",t);this._stepper.stepComplete(t,e);var d=u/Math.pow(o,i.beta),v=e/(d=l.constrain(d/i.factorSafe,p,h));this._t+=e,this._h=n?Math.min(v,e):Math.min(v,a.stepSizeMax),this._lastError=s}else n=!0,this._nStepsAccepted>=1&&this._nStepsRejected++,e/=Math.min(h,u/i.factorSafe)}return this._t},t.prototype._isStiff=function(t){if(this._stiffNStiff>0||this._nStepsAccepted%this._control.stiffCheck==0)if(this._stepper.isStiff(t)){if(this._stiffNNonstiff=0,this._stiffNStiff++>=15)return!0}else this._stiffNStiff>0&&this._stiffNNonstiff++>=6&&(this._stiffNStiff=0,this._stiffNNonstiff=0);return!1},t}();e.Dopri=h},167:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri5StepControl=void 0;e.Dopri5StepControl=function(){this.factorSafe=.9,this.factorMin=.2,this.factorMax=10,this.beta=.04,this.constant=.2-.75*this.beta}},922:function(t,e,r){var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var o=Object.getOwnPropertyDescriptor(e,r);o&&!("get"in o?!e.__esModule:o.writable||o.configurable)||(o={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,o)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),o=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),i=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return o(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.Dopri5=void 0;var a=i(r(660)),s=i(r(167)),u=r(430),l=function(){function t(t,e){this.order=5,this.stepControl=new s.Dopri5StepControl,this.nEval=0,this.rhs=t,this.n=e,this.y=new Array(e),this.yNext=new Array(e),this.yStiff=new Array(e),this.k1=new Array(e),this.k2=new Array(e),this.k3=new Array(e),this.k4=new Array(e),this.k5=new Array(e),this.k6=new Array(e),this.history=new u.HistoryElement(this.order*e)}return t.prototype.step=function(t,e){var r=this.n,n=this.y,o=this.yNext,i=this.k1,a=this.k2,s=this.k3,u=this.k4,l=this.k5,c=this.k6,f=this.history.data,h=0;for(h=0;h<r;++h)o[h]=n[h]+.2*e*i[h];for(this.rhs(t+.2*e,o,a),h=0;h<r;++h)o[h]=n[h]+e*(.075*i[h]+.225*a[h]);for(this.rhs(t+.3*e,o,s),h=0;h<r;++h)o[h]=n[h]+e*(.9777777777777777*i[h]+-3.7333333333333334*a[h]+3.5555555555555554*s[h]);for(this.rhs(t+.8*e,o,u),h=0;h<r;++h)o[h]=n[h]+e*(2.9525986892242035*i[h]+-11.595793324188385*a[h]+9.822892851699436*s[h]+-.2908093278463649*u[h]);for(this.rhs(t+.8888888888888888*e,o,l),h=0;h<r;++h)this.yStiff[h]=n[h]+e*(2.8462752525252526*i[h]+-10.757575757575758*a[h]+8.906422717743473*s[h]+.2784090909090909*u[h]+-.2735313036020583*l[h]);var p=t+e;for(this.rhs(p,this.yStiff,c),h=0;h<r;++h)o[h]=n[h]+e*(.09114583333333333*i[h]+.44923629829290207*s[h]+.6510416666666666*u[h]+-.322376179245283*l[h]+.13095238095238096*c[h]);this.rhs(p,o,a);var d=4*r;for(h=0;h<r;++h)f[d++]=e*(-1.1270175653862835*i[h]+2.675424484351598*s[h]+-5.685526961588504*u[h]+3.5219323679207912*l[h]+-1.7672812570757455*c[h]+2.382468931778144*a[h]);for(h=0;h<r;++h)u[h]=e*(.0012326388888888888*i[h]+-.0042527702905061394*s[h]+.03697916666666667*u[h]+-.05086379716981132*l[h]+.0419047619047619*c[h]+-.025*a[h]);this.nEval+=6},t.prototype.stepComplete=function(t,e){this.saveHistory(t,e),a.copyArray(this.k1,this.k2),a.copyArray(this.y,this.yNext)},t.prototype.error=function(t,e){var r=0,n=0;for(n=0;n<this.n;++n){var o=t+e*Math.max(Math.abs(this.y[n]),Math.abs(this.yNext[n]));r+=a.square(this.k4[n]/o)}return Math.sqrt(r/this.n)},t.prototype.interpolate=function(t,e){for(var r=e.data,n=(t-e.t)/e.h,o=1-n,i=this.n,a=new Array(i),s=0;s<i;++s)a[s]=r[s]+n*(r[i+s]+o*(r[2*i+s]+n*(r[3*i+s]+o*r[4*i+s])));return a},t.prototype.isStiff=function(t){for(var e=0,r=0,n=0;n<this.n;++n)e+=a.square(this.k2[n]-this.k6[n]),r+=a.square(this.yNext[n]-this.yStiff[n]);return r>0&&Math.abs(t)*Math.sqrt(e/r)>3.25},t.prototype.reset=function(t,e){for(var r=0;r<this.n;++r)this.y[r]=e[r];this.rhs(t,e,this.k1),this.nEval=1},t.prototype.saveHistory=function(t,e){for(var r=this.history,n=this.n,o=0;o<n;++o){var i=this.yNext[o]-this.y[o],a=e*this.k1[o]-i;r.data[o]=this.y[o],r.data[n+o]=i,r.data[2*n+o]=a,r.data[3*n+o]=-e*this.k2[o]+i-a}r.t=t,r.h=e},t}();e.Dopri5=l},606:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.HistoryElement=void 0;var r=function(){function t(t){this.t=0,this.h=0,this.data=new Array(t)}return t.prototype.clone=function(){var e=new t(this.data.length);return e.t=this.t,e.h=this.h,e.data=this.data.slice(0),e},t}();e.HistoryElement=r},245:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.integrate=e.integrateDDE=e.DDE=e.integrateDopri=e.Dopri=void 0;var n=r(223);Object.defineProperty(e,"Dopri",{enumerable:!0,get:function(){return n.Dopri}}),Object.defineProperty(e,"integrateDopri",{enumerable:!0,get:function(){return n.integrateDopri}});var o=r(742);Object.defineProperty(e,"DDE",{enumerable:!0,get:function(){return o.DDE}}),Object.defineProperty(e,"integrateDDE",{enumerable:!0,get:function(){return o.integrateDDE}});var i=r(961);Object.defineProperty(e,"integrate",{enumerable:!0,get:function(){return i.integrate}})},961:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.integrate=void 0;var n=r(742),o=r(223);e.integrate=function(t,e,r,i,a,s){if(void 0===a&&(a={}),void 0===s&&(s=null),function(t){return 3===t.length}(t)){if(!function(t){return null===t||2===t.length}(s))throw new Error("Can't used delayed output with non-delayed rhs");return(0,o.integrateDopri)(t,e,r,i,a,s)}if(!function(t){return null===t||3===t.length}(s))throw new Error("Can't used non-delayed output with delayed rhs");return(0,n.integrateDDE)(t,e,r,i,a,s)}},216:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.interpolator=void 0,e.interpolator=function(t,e,r){return function(n){return function(t,e,r,n){for(var o=[],i=e,a=0,s=0,u=t;s<u.length;s++){for(var l=u[s];i[a].t+i[a].h<l;)a++;var c=r.interpolate(l,i[a]);null!==n&&(c=c.concat(n(l,c))),o.push(c)}return o}(n,t,e,r)}}},430:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.HistoryElement=void 0;var n=r(606);Object.defineProperty(e,"HistoryElement",{enumerable:!0,get:function(){return n.HistoryElement}})},660:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.search=e.seqLen=e.approxEqualArray=e.approxEqual=e.zeros=e.copyArray=e.constrain=e.square=void 0;var r=Math.pow(2,-26);e.square=function(t){return t*t},e.constrain=function(t,e,r){return Math.max(Math.min(t,r),e)},e.copyArray=function(t,e){for(var r=t.length,n=0;n<r;n++)t[n]=e[n]},e.zeros=function(t){for(var e=Array(t),r=0;r<t;++r)e[r]=0;return e},e.approxEqual=function(t,e,n){void 0===n&&(n=r);var o=Math.abs(t-e),i=Math.abs(t);return i>n&&(o/=i),o<n},e.approxEqualArray=function(t,e,n){if(void 0===n&&(n=r),e.length!==t.length)throw Error("Incompatible arrays");for(var o=0,i=0,a=0,s=0;s<t.length;++s)t[s]!==e[s]&&(o+=Math.abs(t[s]),i+=Math.abs(t[s]-e[s]),a++);return 0===a||(i/=a,(o/=a)>n&&(i/=o),i<n)},e.seqLen=function(t,e,r){for(var n=[],o=0;o<r;++o){var i=o/(r-1);n.push((1-i)*t+i*e)}return n},e.search=function(t,e){var r=0,n=t.length-1;if(0===t.length||e(t[r]))return-1;if(!e(t[n]))return n;for(;n-r>1;){var o=Math.floor((r+n)/2);e(t[o])?n=o:r=o}return r}},863:t=>{t.exports=function(t,e,r,n,o){var i,a;for(i=1;i<t;i++){if(0===r[i-1])return console.log("tridiagonalSolve: failed due to lack of diagonal dominance"),!1;a=e[i]/r[i-1],r[i]-=a*n[i-1],o[i]-=a*o[i-1]}if(0===r[t-1])return console.log("tridiagonalSolve: failed due to singular matrix"),!1;for(o[t-1]/=r[t-1],i=t-2;i>=0;i--){if(0===r[i])return console.log("tridiagonalSolve: failed due to singular matrix"),!1;o[i]=(o[i]-n[i]*o[i+1])/r[i]}return!0}},653:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.base=void 0;const n=r(458),o=r(166),i=r(650),a=r(536);e.base={delay:n.delay,interpolate:o.interpolate,maths:i.maths,user:a.user}},401:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.updatePars=e.batchParsDisplace=e.batchParsRange=e.batchRun=e.Batch=void 0;const n=r(822),o=r(882),i=r(915);class a{constructor(t,e,r,n,o){this.pars=e,this.tStart=r,this.tEnd=n;const a=[],s=[],l=[];if(e.values.forEach((c=>{const f=u(e.base,e.name,c);try{a.push((0,i.wodinRun)(t,f,r,n,o)),l.push(c)}catch(t){s.push({value:c,error:t.message})}})),0===a.length)throw Error(`All solutions failed; first error: ${s[0].error}`);this.pars=Object.assign(Object.assign({},e),{values:l}),this.solutions=a,this.errors=s}valueAtTime(t){const e=this.solutions.map((e=>e({mode:n.TimeMode.Given,times:[t]}))),r=this.pars.values;return{values:e[0].values.map(((t,r)=>(t=>({name:e[0].values[t].name,y:e.map((e=>e.values[t].y[0]))}))(r))),x:r}}extreme(t){return this.findExtremes()[t]}findExtremes(){if(void 0===this._extremes){const t=501,e={mode:n.TimeMode.Grid,nPoints:t,tEnd:this.tEnd,tStart:this.tStart},r=this.solutions.map((t=>t(e))),i=r[0].x,a=r[0].values.map((t=>t.name)),s=(0,o.loop)(a.length,(t=>r.map((e=>function(t,e){const r=(0,o.whichMin)(e),n=(0,o.whichMax)(e),i=t[r],a=t[n],s=e[r];return{tMax:a,tMin:i,yMax:e[n],yMin:s}}(i,e.values[t].y))))),u=this.pars.values;this._extremes={tMax:c("tMax",a,u,s),tMin:c("tMin",a,u,s),yMax:c("yMax",a,u,s),yMin:c("yMin",a,u,s)}}return this._extremes}}function s(t,e,r,n,i,a){const s=l(t,e);if(i>s)throw Error(`Expected lower bound to be no greater than ${s}`);if(a<s)throw Error(`Expected upper bound to be no less than ${s}`);if(i>=a)throw Error("Expected upper bound to be greater than lower bound");if(r<2)throw Error("Must include at least 2 traces in the batch");if(n&&i<=0)throw Error("Lower bound must be greater than 0 for logarithmic scale");return{base:t,name:e,values:n?(0,o.gridLog)(i,a,r):(0,o.grid)(i,a,r)}}function u(t,e,r){const n=Object.assign({},t);return n[e]=r,n}function l(t,e){const r=t[e];if(void 0===r)throw Error(`Expected a value for '${e}'`);if("number"!=typeof r)throw Error(`Expected a number for '${e}'`);return r}function c(t,e,r,n){return{x:r,values:(0,o.loop)(e.length,(r=>({name:e[r],y:n[r].map((e=>e[t]))})))}}e.Batch=a,e.batchRun=function(t,e,r,n,o){return new a(t,e,r,n,o)},e.batchParsRange=s,e.batchParsDisplace=function(t,e,r,n,o){const i=l(t,e),a=o/100;return s(t,e,r,n,i*(1-a),i*(1+a))},e.updatePars=u},458:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.delay=void 0,e.delay=function(t,e,r,n){const o=t(e);for(let t=0;t<r.length;++t)n[t]=o[r[t]]}},410:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.sumOfSquares=e.updatePars=e.fitTarget=void 0;const n=r(653),o=r(134),i=r(822);function a(t,e){const r=Object.assign({},t.base);for(let n=0;n<t.vary.length;++n)r[t.vary[n]]=e[n];return r}function s(t,e){let r=0;for(let n=0;n<t.length;++n){const o=t[n];isNaN(o)||(r+=Math.pow(o-e[n],2))}return r}e.fitTarget=function(t,e,r,u,l){if(e.time[0]<0)throw Error("Expected the first time to be at least 0");const c=e.time[e.time.length-1];return f=>{const h=a(r,f),p=new t(n.base,h,"error"),d=(0,o.runModel)(p,null,0,c,l).solution,v=p.names(),_=v.indexOf(u),y=d(e.time).map((t=>t[_]));return{data:{endTime:c,names:v,pars:h,solution:(0,i.interpolatedSolution)(d,v,0,c),startTime:0},value:s(e.value,y)}}},e.updatePars=a,e.sumOfSquares=s},166:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.interpolate=e.interpolateTimes=e.interpolateCheckT=e.interpolateCheckY=e.interpolateAlloc=void 0;const n=r(673);function o(t,e,r){const o=function(t,e){if(t.length===e)return[t];const r=[],n=t.length/e;for(let o=0;o<n;++o)r.push(t.slice(e*o,e*(o+1)));return r}(r,e.length);if("constant"===t)return new n.InterpolatorConstant(e,o);if("linear"===t)return new n.InterpolatorLinear(e,o);if("spline"===t)return new n.InterpolatorSpline(e,o);throw Error(`Invalid interpolation method '${t}'`)}function i(t,e,r,n){if(1===e.length){if(t[0]!==e[0])throw Error(`Expected ${r} to have length ${t[0]} (for ${n})`)}else for(let o=0;o<e.length;++o)if(t[o]!==e[o])throw Error(`Expected dimension ${o+1} of ${r} to have size ${t[o]} (for ${n})`)}function a(t,e){const r=Math.max(...t);return{max:Math.min(...e),min:r}}e.interpolateAlloc=o,e.interpolateCheckY=i,e.interpolateCheckT=function(t,e,r,n){if(void 0===r)return 1/0;if(t<r.min)throw Error(`Integration times do not span interpolation range; min: ${r.min}`);if(e>r.max)throw Error(`Integration times do not span interpolation range; max: ${r.max}`);return void 0===n&&(n=r.max),n},e.interpolateTimes=a,e.interpolate={alloc:o,checkY:i,times:a}},650:(t,e)=>{function r(t){return.5===o(t,1)?2*Math.round(t/2):Math.round(t)}function n(t,e){if(void 0===e||0===e)return r(t);{const n=Math.pow(10,e);return r(t*n)/n}}function o(t,e){let r=t%e;return r*e<0&&(r+=e),r}function i(t,e){return Math.floor(t/e)}function a(t,e,r){let n=0;for(let o=e;o<r;++o)n+=t[o];return n}function s(t,e,r,n,o,i){let a=0;for(let s=n;s<o;++s){const n=s*i;for(let o=e;o<r;++o)a+=t[o+n]}return a}function u(t,e,r,n,o,i,a,s,u){let l=0;for(let c=i;c<a;++c){const i=c*u;for(let a=n;a<o;++a){const n=a*s+i;for(let o=e;o<r;++o)l+=t[o+n]}}return l}Object.defineProperty(e,"__esModule",{value:!0}),e.maths=e.odinSum3=e.odinSum2=e.odinSum1=e.intdivr=e.modr=e.round2=void 0,e.round2=n,e.modr=o,e.intdivr=i,e.odinSum1=a,e.odinSum2=s,e.odinSum3=u,e.maths={intdivr:i,modr:o,odinSum1:a,odinSum2:s,odinSum3:u,round2:n}},134:function(t,e,r){var n=this&&this.__createBinding||(Object.create?function(t,e,r,n){void 0===n&&(n=r);var o=Object.getOwnPropertyDescriptor(e,r);o&&!("get"in o?!e.__esModule:o.writable||o.configurable)||(o={enumerable:!0,get:function(){return e[r]}}),Object.defineProperty(t,n,o)}:function(t,e,r,n){void 0===n&&(n=r),t[n]=e[r]}),o=this&&this.__setModuleDefault||(Object.create?function(t,e){Object.defineProperty(t,"default",{enumerable:!0,value:e})}:function(t,e){t.default=e}),i=this&&this.__importStar||function(t){if(t&&t.__esModule)return t;var e={};if(null!=t)for(var r in t)"default"!==r&&Object.prototype.hasOwnProperty.call(t,r)&&n(e,t,r);return o(e,t),e};Object.defineProperty(e,"__esModule",{value:!0}),e.runModel=e.isODEModel=e.isDDEModel=void 0;const a=i(r(245)),s=r(166);function u(t){return 4===t.rhs.length}e.isDDEModel=u,e.isODEModel=function(t){return 3===t.rhs.length},e.runModel=function(t,e,r,n,o){const i=t.getMetadata().interpolateTimes;return o.tcrit=(0,s.interpolateCheckT)(r,n,i,o.tcrit),u(t)?function(t,e,r,n,o){let i=null;"function"==typeof t.output&&(i=(e,r,n)=>t.output(e,r,n)),null===e&&(e=t.initial(r)),t.getInternal().initial_t=r;const s=new a.DDE((function(e,r,n,o){t.rhs(e,r,n,o)}),e.length,o,i);return s.initialise(r,e),{solution:s.run(n),statistics:s.statistics()}}(t,e,r,n,o):function(t,e,r,n,o){let i=null;"function"==typeof t.output&&(i=(e,r)=>t.output(e,r)),null===e&&(e=t.initial(r));const s=new a.Dopri((function(e,r,n){t.rhs(e,r,n)}),e.length,o,i);return s.initialise(r,e),{solution:s.run(n),statistics:s.statistics()}}(t,e,r,n,o)}},586:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.PkgWrapper=void 0;const n=r(653),o=r(134);e.PkgWrapper=class{constructor(t,e,r){this.model=new t(n.base,e,r)}initial(t){return this.model.initial(t)}rhs(t,e){const r=new Array(e.length).fill(0);let n=null;if(!(0,o.isODEModel)(this.model))throw Error("Can't use rhs() with delay models");return this.model,this.model.rhs(t,e,r),this.model.output&&(n=this.model.output(t,e)),{output:n,state:r}}getMetadata(){return this.model.getMetadata()}getInternal(){return this.model.getInternal()}setUser(t,e){this.model.setUser(t,e)}run(t,e,r){const n=t[0],i=t[t.length-1],a=(0,o.runModel)(this.model,e,n,i,r);return{names:this.model.names(),statistics:a.statistics,y:a.solution(t)}}}},822:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.interpolatedSolution=e.TimeMode=void 0;const n=r(882);var o;!function(t){t.Grid="grid",t.Given="given"}(o=e.TimeMode||(e.TimeMode={})),e.interpolatedSolution=function(t,e,r,i){return a=>{const s=function(t,e,r){switch(t.mode){case o.Grid:return(0,n.grid)(Math.max(e,t.tStart),Math.min(r,t.tEnd),t.nPoints);case o.Given:return t.times}}(a,r,i),u=t(s);return{values:e.map(((t,r)=>({name:e[r],y:u.map((t=>t[r]))}))),x:s}}}},536:(t,e)=>{function r(t,e,r){if("ignore"===r)return;const n=[];for(const r of Object.keys(t))e.includes(r)||n.push(r);if(n.length>0){const t="Unknown user parameters: "+n.join(", ");if("message"===r)console.log(t);else{if("warning"!==r)throw"stop"===r?Error(t):Error(t+" (and invalid value for unusedUserAction)");console.warn(t)}}}function n(t,e,r,n,o,i,a){const s=t[e];if(void 0===s){if(void 0!==r[e])return;if(null===n)throw Error(`Expected a value for '${e}'`);r[e]=n}else{if("number"!=typeof s)throw Error(`Expected a number for '${e}'`);l(s,o,i,a,e),r[e]=s}}function o(t,e,r,n,o,i,l){let c=t[e];if(void 0===c){if(void 0!==r[e])return;throw Error(`Expected a value for '${e}'`)}{const t=n.length-1;c=a(c),s(t,c,e),function(t,e,r){const n=t.length-1;for(let o=0;o<n;++o){const i=t[o+1];if(e.dim[o]!==i)throw 1===n?Error(`Expected length ${i} value for '${r}'`):Error(`Incorrect size of dimension ${o+1} of '${r}' (expected ${i})`)}}(n,c,e),u(c,o,i,l,e),r[e]=c.data.slice()}}function i(t,e,r,n,o,i,l){let c=t[e];if(void 0===c){if(void 0!==r[e])return;throw Error("Expected a value for '"+e+"'")}{const t=n.length-1;c=a(c),s(t,c,e),u(c,o,i,l,e),n[0]=c.data.length;for(let e=0;e<t;++e)n[e+1]=c.dim[e];r[e]=c.data.slice()}}function a(t,e){return Array.isArray(t)?t={data:t,dim:[t.length]}:"number"==typeof t&&(t={data:[t],dim:[1]}),t}function s(t,e,r){if(e.dim.length!==t)throw 1===t?Error(`Expected a numeric vector for '${r}'`):2===t?Error(`Expected a numeric matrix for '${r}'`):Error(`Expected a numeric array of rank ${t} for '${r}'`)}function u(t,e,r,n,o){for(const i of t.data){if(null===i)throw Error(`'${o}' must not contain any NA values`);l(i,e,r,n,o)}}function l(t,e,r,n,o){if("number"!=typeof t)throw Error(`Expected a number for '${o}'`);if(t<e)throw Error(`Expected '${o}' to be at least ${e}`);if(t>r)throw Error(`Expected '${o}' to be at most ${r}`);if(n&&!Number.isInteger(t))throw Error(`Expected '${o}' to be integer-like`)}Object.defineProperty(e,"__esModule",{value:!0}),e.user=e.setUserArrayVariable=e.setUserArrayFixed=e.setUserScalar=e.checkUser=void 0,e.checkUser=r,e.setUserScalar=n,e.setUserArrayFixed=o,e.setUserArrayVariable=i,e.user={checkUser:r,setUserArrayFixed:o,setUserArrayVariable:i,setUserScalar:n}},882:(t,e)=>{function r(t,e,r){const n=(e-t)/(r-1),o=[];for(let e=0;e<r-1;++e)o.push(t+e*n);return o.push(e),o}Object.defineProperty(e,"__esModule",{value:!0}),e.loop=e.whichMax=e.whichMin=e.gridLog=e.grid=void 0,e.grid=r,e.gridLog=function(t,e,n){return r(Math.log(t),Math.log(e),n).map(Math.exp)},e.whichMin=function(t){let e=-1,r=1/0;for(let n=0;n<t.length;++n)t[n]<r&&(e=n,r=t[n]);return e},e.whichMax=function(t){let e=-1,r=-1/0;for(let n=0;n<t.length;++n)t[n]>r&&(e=n,r=t[n]);return e},e.loop=function(t,e){const r=[];for(let n=0;n<t;++n)r.push(e(n));return r}},556:(t,e)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.versions=void 0,e.versions=function(){return{dfoptim:"0.0.5",dopri:"0.0.12",odinjs:"0.0.17"}}},915:(t,e,r)=>{Object.defineProperty(e,"__esModule",{value:!0}),e.wodinFitValue=e.wodinFit=e.wodinRun=void 0;const n=r(153),o=r(653),i=r(410),a=r(134),s=r(822);e.wodinRun=function(t,e,r,n,i){const u=new t(o.base,e,"error"),l=(0,a.runModel)(u,null,r,n,i).solution,c=u.names();return(0,s.interpolatedSolution)(l,c,r,n)},e.wodinFit=function(t,e,r,o,a,s){const u=(0,i.fitTarget)(t,e,r,o,a),l=r.vary.map((t=>r.base[t]));return new n.Simplex(u,l,s)},e.wodinFitValue=function(t,e,r){const n=t({mode:s.TimeMode.Given,times:e.time}),o=n.values.findIndex((t=>t.name===r));return(0,i.sumOfSquares)(e.value,n.values[o].y)}}},e={};function r(n){var o=e[n];if(void 0!==o)return o.exports;var i=e[n]={exports:{}};return t[n].call(i.exports,i,i.exports,r),i.exports}var n={};(()=>{var t=n;Object.defineProperty(t,"__esModule",{value:!0}),t.versions=t.batchRun=t.batchParsRange=t.batchParsDisplace=t.Batch=t.base=t.PkgWrapper=t.wodinRun=t.wodinFitValue=t.wodinFit=void 0;var e=r(915);Object.defineProperty(t,"wodinFit",{enumerable:!0,get:function(){return e.wodinFit}}),Object.defineProperty(t,"wodinFitValue",{enumerable:!0,get:function(){return e.wodinFitValue}}),Object.defineProperty(t,"wodinRun",{enumerable:!0,get:function(){return e.wodinRun}});var o=r(586);Object.defineProperty(t,"PkgWrapper",{enumerable:!0,get:function(){return o.PkgWrapper}});var i=r(653);Object.defineProperty(t,"base",{enumerable:!0,get:function(){return i.base}});var a=r(401);Object.defineProperty(t,"Batch",{enumerable:!0,get:function(){return a.Batch}}),Object.defineProperty(t,"batchParsDisplace",{enumerable:!0,get:function(){return a.batchParsDisplace}}),Object.defineProperty(t,"batchParsRange",{enumerable:!0,get:function(){return a.batchParsRange}}),Object.defineProperty(t,"batchRun",{enumerable:!0,get:function(){return a.batchRun}});var s=r(556);Object.defineProperty(t,"versions",{enumerable:!0,get:function(){return s.versions}})})(),odinjs=n})();
diff --git a/inst/js/test-discrete.js b/inst/js/test-discrete.js
index 7250ebae..05e8f702 100644
--- a/inst/js/test-discrete.js
+++ b/inst/js/test-discrete.js
@@ -1,3 +1,4 @@
 function call_odin_bundle(Model, pars, tStart, tEnd, dt, nParticles) {
-    return dust.wodinRunDiscrete(Model, pars, tStart, tEnd, dt, nParticles)
+    const solution = dust.wodinRunDiscrete(Model, pars, tStart, tEnd, dt, nParticles);
+    return solution({mode: "grid", tStart, tEnd, nPoints: Infinity});
 }
diff --git a/tests/testthat/helper-js.R b/tests/testthat/helper-js.R
index 8dc74d50..5515ab39 100644
--- a/tests/testthat/helper-js.R
+++ b/tests/testthat/helper-js.R
@@ -14,9 +14,8 @@ call_odin_bundle_continuous <- function(bundle, user, t0, t1, tn,
   } else {
     control_js <- V8::JS(jsonlite::toJSON(control, auto_unbox = TRUE))
   }
-  res <- ct$call("call_odin_bundle", odin_js, user_js, t0, t1, tn, control_js)
-  res$y <- t(res$y)
-  res
+
+  ct$call("call_odin_bundle", odin_js, user_js, t0, t1, tn, control_js)
 }
 
 
@@ -74,3 +73,7 @@ skip_if_no_js <- function() {
   ## behaves poorly.
   skip_on_cran()
 }
+
+list_to_matrix <- function(x) {
+  matrix(unlist(x), ncol = length(x))
+}
diff --git a/tests/testthat/test-js-bundle.R b/tests/testthat/test-js-bundle.R
index 5c58e487..dab0d02b 100644
--- a/tests/testthat/test-js-bundle.R
+++ b/tests/testthat/test-js-bundle.R
@@ -22,9 +22,9 @@ test_that("bundle works", {
   sol <- with(list(K = 100, r = 0.5, y0 = 1),
               K / (1 + (K / y0 - 1) * exp(-r * t)))
 
-  expect_equal(res$names, "N")
   expect_equal(res$x, t)
-  expect_equal(drop(res$y), sol, tolerance = 1e-6)
+  expect_equal(res$values$name, "N")
+  expect_equal(res$values$y[[1]], sol, tolerance = 1e-6)
 })
 
 
@@ -52,9 +52,9 @@ test_that("include interpolate", {
 
   tt <- seq(t0, t1, length.out = tn)
 
-  expect_equal(res$names, c("y", "p"))
+  expect_equal(res$values$name, c("y", "p"))
   zz <- ifelse(tt < 1, 0, ifelse(tt > 2, 1, tt - 1))
-  expect_equal(res$y[, 1], zz, tolerance = 2e-5)
+  expect_equal(res$values$y[[1]], zz, tolerance = 2e-5)
 })
 
 
@@ -78,12 +78,12 @@ test_that("include sum", {
 
   res <- call_odin_bundle_continuous(bundle, NULL, t0, t1, tn)
   expect_equal(
-    res$names,
+    res$values$name,
     c("y[1]", "y[2]", "y[3]", "ytot", "y2[1]", "y2[2]", "y2[3]"))
 
-  y <- res$y[, 1:3]
-  ytot <- res$y[, 4]
-  y2 <- res$y[, 5:7]
+  y <- list_to_matrix(res$values$y[1:3])
+  ytot <- res$values$y[[4]]
+  y2 <- list_to_matrix(res$values$y[5:7])
 
   expect_equal(ytot, rowSums(y))
   expect_equal(y2, y * 2)
@@ -127,8 +127,10 @@ test_that("include fancy sum", {
   tt <- seq(t0, t1, length.out = tn)
   cmp <- odin::odin_(code, target = "r")$new(user = user)$run(tt)
 
-  expect_equal(res$names, c("y[1]", "y[2]", "y[3]", "y[4]"))
-  expect_equal(res$y, unname(cmp[, -1]), tolerance = 1e-5)
+  expect_equal(res$values$name, c("y[1]", "y[2]", "y[3]", "y[4]"))
+
+  expect_equal(list_to_matrix(res$values$y),
+               unname(cmp[, -1]), tolerance = 1e-5)
 })