diff-grok 1.1.0 → 1.2.0

package/src/solver-tools/lsoda/stoda.js.map

@@ -1 +1 @@
-{"version":3,"file":"stoda.js","sourceRoot":"","sources":["stoda.ts"],"names":[],"mappings":"AAAA,OAAO,EAAe,KAAK,EAAE,IAAI,EAAC,MAAM,UAAU,CAAC;AACnD,OAAO,EAAC,MAAM,EAAC,MAAM,QAAQ,CAAC;AAC9B,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAC9B,OAAO,EAAC,MAAM,EAAC,MAAM,UAAU,CAAC;AAChC,OAAO,EAAC,UAAU,EAAkB,MAAM,cAAc,CAAC;AACzD,OAAO,EAAC,WAAW,EAAoB,MAAM,eAAe,CAAC;AAC7D,OAAO,EAAC,YAAY,EAAqB,MAAM,gBAAgB,CAAC;AAEhE;;;;;GAKG;AACH,MAAM,UAAU,KAAK,CAAC,GAAiB,EAAE,CAAe,EAAE,MAAc;IACtE,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,IAAI,GAAG,GAAG,CAAC,GAAI,CAAC,IAAI,CAAC;IAC3B,MAAM,MAAM,GAAG,GAAG,CAAC,GAAI,CAAC,MAAM,CAAC;IAC/B,MAAM,MAAM,GAAG,GAAG,CAAC,GAAI,CAAC,MAAM,CAAC;IAC/B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IAEpB,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,MAAM,IAAI,GAAG,CAAC,CAAC,EAAE,CAAC;IAClB,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IACV,IAAI,IAAI,GAAG,EAAE,CAAC;IAEd,IAAI,MAAM,GAAG,MAAM,CAAC;IACpB,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;QAAE,MAAM,GAAG,MAAM,CAAC;IAElC,+CAA+C;IAC/C,SAAS,QAAQ;QACf,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAC3B,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QACjB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;IACf,CAAC;IAED,6DAA6D;IAC7D,SAAS,UAAU;QACjB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE;YAChC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5B,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAC9B,CAAC;IAED,yCAAyC;IAEzC,IAAI,MAAM,KAAK,CAAC,EAAE,CAAC;QACjB,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACT,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;QACZ,CAAC,CAAC,IAAI,GAAG,MAAM,CAAC;QAChB,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC;QACV,CAAC,CAAC,KAAK,GAAG,GAAG,CAAC;QACd,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;QACb,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;QACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;QACjB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,CAAC,CAAC,MAAM,GAAG,EAAE,CAAC;QACd,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;QACb,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC;QACb,CAAC,CAAC,MAAM,GAAG,EAAE,CAAC;QACd,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACd,UAAU,EAAE,CAAC;IACf,CAAC;IAED,IAAI,MAAM,KAAK,CAAC,CAAC,EAAE,CAAC;QAClB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;QACjB,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC;YAAE,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;QAC/B,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,EAAE,CAAC;YACvB,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;YACnB,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;YACnB,UAAU,EAAE,CAAC;QACf,CAAC;QACD,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC;YACnB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACxB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACb,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAClB,CAAC;IACH,CAAC;IAED,IAAI,MAAM,KAAK,CAAC,CAAC,EAAE,CAAC;QAClB,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC;YACnB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACxB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACb,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAClB,CAAC;IACH,CAAC;IAED,gCAAgC;IAEhC,IAAI,GAAG,GAAG,GAAG,CAAC;IACd,MAAM,EAAE,GAAoB,EAAC,GAAG,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAC,CAAC;IACvD,MAAM,QAAQ,GAAsB,EAAC,EAAE,EAAE,CAAC,EAAC,CAAC;IAC5C,MAAM,QAAQ,GAAuB,EAAC,EAAE,EAAE,CAAC,EAAC,CAAC;IAC7C,IAAI,EAAU,CAAC;IAEf,KAAK,EACL,OAAO,IAAI,EAAE,CAAC;QACZ,8BAA8B;QAC9B,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;QAEX,KAAK,EACL,OAAO,IAAI,EAAE,CAAC;YACZ,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,KAAK;gBAC7B,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;YACnB,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,IAAI,GAAG,IAAI;gBACxB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;YAEnB,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;YAEZ,oDAAoD;YACpD,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;oBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnC,CAAC;YACH,CAAC;YAED,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;YAE1C,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;YACX,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC;YACf,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;YACT,MAAM,OAAO,GAAG,UAAU,CAAC,GAAG,EAAE,CAAC,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;YACpD,IAAI,GAAG,EAAE,CAAC,IAAI,CAAC;YAEf,IAAI,OAAO,KAAK,CAAC;gBAAE,MAAM,KAAK,CAAC;YAE/B,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;gBAClB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1C,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;gBAChB,SAAS,KAAK,CAAC;YACjB,CAAC;YAED,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;gBAClB,KAAK,GAAG,CAAC,CAAC,CAAC;gBACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;gBACb,MAAM,GAAG,CAAC,CAAC;gBACX,OAAO,KAAK,CAAC;YACf,CAAC;QACH,CAAC;QAED,mBAAmB;QACnB,IAAI,EAAE,CAAC,CAAC,KAAK,CAAC;YACZ,GAAG,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAClC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC;YACV,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAEtD,IAAI,GAAG,IAAI,EAAE,EAAE,CAAC;YACd,kBAAkB;YAClB,KAAK,GAAG,CAAC,CAAC;YACV,CAAC,CAAC,GAAG,EAAE,CAAC;YACR,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACX,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC;YACb,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,IAAI,CAAC;YAEjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACnC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;oBAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAChC,CAAC;YAED,CAAC,CAAC,MAAM,EAAE,CAAC;YACX,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;gBACjB,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,YAAY,CAAC,GAAG,EAAE,GAAG,EAAE,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,CAAC;gBAC9D,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,EAAE,CAAC;oBACvB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;oBAChB,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC;oBACb,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;YACH,CAAC;YAED,6CAA6C;YAC7C,CAAC,CAAC,KAAK,EAAE,CAAC;YACV,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;gBAClB,IAAI,IAAI,GAAG,EAAE,CAAC;gBACd,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,KAAK,MAAM,GAAG,CAAC,EAAE,CAAC;oBAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC9C,MAAM,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC1D,MAAM,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;oBAC7B,IAAI,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC;gBACtD,CAAC;gBAED,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,MAAM,SAAS,GAAG,WAAW,CAAC,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;gBAEvE,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;oBAChB,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC;oBACb,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,UAAU,EAAE,CAAC;oBACb,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;oBAChB,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC;oBACb,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;YACH,CAAC;YAED,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,KAAK,MAAM,GAAG,CAAC,EAAE,CAAC;gBAC3C,QAAQ,EAAE,CAAC;gBACX,MAAM,KAAK,CAAC;YACd,CAAC;YAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;gBAC3B,CAAC,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAClC,QAAQ,EAAE,CAAC;YACX,MAAM,KAAK,CAAC;QACd,CAAC;aAAM,CAAC;YACN,oBAAoB;YACpB,KAAK,EAAE,CAAC;YACR,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC;YAEZ,mBAAmB;YACnB,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;oBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnC,CAAC;YACH,CAAC;YAED,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC;YACZ,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,OAAO,EAAE,CAAC;gBACpC,KAAK,GAAG,CAAC,CAAC,CAAC;gBACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;gBACb,MAAM,GAAG,CAAC,CAAC;gBACX,MAAM,KAAK,CAAC;YACd,CAAC;YAED,IAAI,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC;gBACf,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,MAAM,SAAS,GAAG,WAAW,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;gBACrE,IAAI,SAAS,KAAK,CAAC,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACvC,IAAI,SAAS,KAAK,CAAC;wBACjB,QAAQ,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;oBAC3C,QAAQ,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC1D,MAAM,CAAC,GAAG,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC;gBAC3B,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,UAAU,EAAE,CAAC;oBACb,QAAQ,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC1D,MAAM,CAAC,GAAG,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC;gBAC3B,CAAC;gBACD,SAAS,KAAK,CAAC;YACjB,CAAC;iBAAM,CAAC;gBACN,qBAAqB;gBACrB,IAAI,KAAK,KAAK,CAAC,EAAE,EAAE,CAAC;oBAClB,KAAK,GAAG,CAAC,CAAC,CAAC;oBACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;oBACb,MAAM,GAAG,CAAC,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;qBAAM,CAAC;oBACN,EAAE,GAAG,GAAG,CAAC;oBACT,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;oBACxC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;oBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;oBAC5D,CAAC,CAAC,GAAG,EAAE,CAAC;oBACR,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;oBAC/B,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;oBACjB,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;oBACZ,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC;wBAAE,SAAS,KAAK,CAAC;oBAC/B,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;oBACT,UAAU,EAAE,CAAC;oBACb,SAAS,KAAK,CAAC;gBACjB,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,KAAK,CAAC;AACf,CAAC"}
+{"version":3,"file":"stoda.js","sourceRoot":"","sources":["stoda.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AAEH,OAAO,EAAe,KAAK,EAAE,IAAI,EAAC,MAAM,UAAU,CAAC;AACnD,OAAO,EAAC,MAAM,EAAC,MAAM,QAAQ,CAAC;AAC9B,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAC9B,OAAO,EAAC,MAAM,EAAC,MAAM,UAAU,CAAC;AAChC,OAAO,EAAC,UAAU,EAAkB,MAAM,cAAc,CAAC;AACzD,OAAO,EAAC,WAAW,EAAoB,MAAM,eAAe,CAAC;AAC7D,OAAO,EAAC,YAAY,EAAqB,MAAM,gBAAgB,CAAC;AAEhE;;;;;GAKG;AACH,MAAM,UAAU,KAAK,CAAC,GAAiB,EAAE,CAAe,EAAE,MAAc;IACtE,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,IAAI,GAAG,GAAG,CAAC,GAAI,CAAC,IAAI,CAAC;IAC3B,MAAM,MAAM,GAAG,GAAG,CAAC,GAAI,CAAC,MAAM,CAAC;IAC/B,MAAM,MAAM,GAAG,GAAG,CAAC,GAAI,CAAC,MAAM,CAAC;IAC/B,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IAEpB,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,MAAM,IAAI,GAAG,CAAC,CAAC,EAAE,CAAC;IAClB,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IACV,IAAI,IAAI,GAAG,EAAE,CAAC;IAEd,IAAI,MAAM,GAAG,MAAM,CAAC;IACpB,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;QAAE,MAAM,GAAG,MAAM,CAAC;IAElC,+CAA+C;IAC/C,SAAS,QAAQ;QACf,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAC3B,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;QACjB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;IACf,CAAC;IAED,6DAA6D;IAC7D,SAAS,UAAU;QACjB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE;YAChC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAC5B,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IAC9B,CAAC;IAED,yCAAyC;IAEzC,IAAI,MAAM,KAAK,CAAC,EAAE,CAAC;QACjB,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACT,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;QACZ,CAAC,CAAC,IAAI,GAAG,MAAM,CAAC;QAChB,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC;QACV,CAAC,CAAC,KAAK,GAAG,GAAG,CAAC;QACd,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;QACb,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;QACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;QACjB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,CAAC,CAAC,MAAM,GAAG,EAAE,CAAC;QACd,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;QACb,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC;QACb,CAAC,CAAC,MAAM,GAAG,EAAE,CAAC;QACd,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACd,UAAU,EAAE,CAAC;IACf,CAAC;IAED,IAAI,MAAM,KAAK,CAAC,CAAC,EAAE,CAAC;QAClB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;QACjB,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC;YAAE,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;QAC/B,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,EAAE,CAAC;YACvB,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC;YACnB,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;YACnB,UAAU,EAAE,CAAC;QACf,CAAC;QACD,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC;YACnB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACxB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACb,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAClB,CAAC;IACH,CAAC;IAED,IAAI,MAAM,KAAK,CAAC,CAAC,EAAE,CAAC;QAClB,IAAI,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC;YACnB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACxB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC;YACb,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAClB,CAAC;IACH,CAAC;IAED,gCAAgC;IAEhC,IAAI,GAAG,GAAG,GAAG,CAAC;IACd,MAAM,EAAE,GAAoB,EAAC,GAAG,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAC,CAAC;IACvD,MAAM,QAAQ,GAAsB,EAAC,EAAE,EAAE,CAAC,EAAC,CAAC;IAC5C,MAAM,QAAQ,GAAuB,EAAC,EAAE,EAAE,CAAC,EAAC,CAAC;IAC7C,IAAI,EAAU,CAAC;IAEf,KAAK,EACL,OAAO,IAAI,EAAE,CAAC;QACZ,8BAA8B;QAC9B,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;QAEX,KAAK,EACL,OAAO,IAAI,EAAE,CAAC;YACZ,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,KAAK;gBAC7B,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;YACnB,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,IAAI,GAAG,IAAI;gBACxB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;YAEnB,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;YAEZ,oDAAoD;YACpD,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;oBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnC,CAAC;YACH,CAAC;YAED,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;YAE1C,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC;YACX,EAAE,CAAC,IAAI,GAAG,IAAI,CAAC;YACf,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;YACT,MAAM,OAAO,GAAG,UAAU,CAAC,GAAG,EAAE,CAAC,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;YACpD,IAAI,GAAG,EAAE,CAAC,IAAI,CAAC;YAEf,IAAI,OAAO,KAAK,CAAC;gBAAE,MAAM,KAAK,CAAC;YAE/B,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;gBAClB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1C,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;gBAChB,SAAS,KAAK,CAAC;YACjB,CAAC;YAED,IAAI,OAAO,KAAK,CAAC,EAAE,CAAC;gBAClB,KAAK,GAAG,CAAC,CAAC,CAAC;gBACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;gBACb,MAAM,GAAG,CAAC,CAAC;gBACX,OAAO,KAAK,CAAC;YACf,CAAC;QACH,CAAC;QAED,mBAAmB;QACnB,IAAI,EAAE,CAAC,CAAC,KAAK,CAAC;YACZ,GAAG,GAAG,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAClC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC;YACV,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QAEtD,IAAI,GAAG,IAAI,EAAE,EAAE,CAAC;YACd,kBAAkB;YAClB,KAAK,GAAG,CAAC,CAAC;YACV,CAAC,CAAC,GAAG,EAAE,CAAC;YACR,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACX,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC;YACb,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,IAAI,CAAC;YAEjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACnC,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;oBAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAChC,CAAC;YAED,CAAC,CAAC,MAAM,EAAE,CAAC;YACX,IAAI,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;gBACjB,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,YAAY,CAAC,GAAG,EAAE,GAAG,EAAE,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,QAAQ,CAAC,CAAC;gBAC9D,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,EAAE,CAAC;oBACvB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;oBAChB,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC;oBACb,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;YACH,CAAC;YAED,6CAA6C;YAC7C,CAAC,CAAC,KAAK,EAAE,CAAC;YACV,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;gBAClB,IAAI,IAAI,GAAG,EAAE,CAAC;gBACd,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,KAAK,MAAM,GAAG,CAAC,EAAE,CAAC;oBAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC9C,MAAM,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC1D,MAAM,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;oBAC7B,IAAI,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,SAAS,CAAC,CAAC;gBACtD,CAAC;gBAED,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,MAAM,SAAS,GAAG,WAAW,CAAC,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;gBAEvE,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;oBAChB,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC;oBACb,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,UAAU,EAAE,CAAC;oBACb,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjD,MAAM,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;oBAChB,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC;oBACb,QAAQ,EAAE,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;YACH,CAAC;YAED,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC,KAAK,MAAM,GAAG,CAAC,EAAE,CAAC;gBAC3C,QAAQ,EAAE,CAAC;gBACX,MAAM,KAAK,CAAC;YACd,CAAC;YAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;gBAC3B,CAAC,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YAClC,QAAQ,EAAE,CAAC;YACX,MAAM,KAAK,CAAC;QACd,CAAC;aAAM,CAAC;YACN,oBAAoB;YACpB,KAAK,EAAE,CAAC;YACR,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC;YAEZ,mBAAmB;YACnB,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC/B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;oBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnC,CAAC;YACH,CAAC;YAED,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC;YACZ,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,OAAO,EAAE,CAAC;gBACpC,KAAK,GAAG,CAAC,CAAC,CAAC;gBACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;gBACb,MAAM,GAAG,CAAC,CAAC;gBACX,MAAM,KAAK,CAAC;YACd,CAAC;YAED,IAAI,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC;gBACf,QAAQ,CAAC,EAAE,GAAG,CAAC,CAAC;gBAChB,MAAM,SAAS,GAAG,WAAW,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,QAAQ,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;gBACrE,IAAI,SAAS,KAAK,CAAC,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACvC,IAAI,SAAS,KAAK,CAAC;wBACjB,QAAQ,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;oBAC3C,QAAQ,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC1D,MAAM,CAAC,GAAG,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC;gBAC3B,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,UAAU,EAAE,CAAC;oBACb,QAAQ,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAE,EAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC1D,MAAM,CAAC,GAAG,EAAE,QAAQ,CAAC,EAAE,CAAC,CAAC;gBAC3B,CAAC;gBACD,SAAS,KAAK,CAAC;YACjB,CAAC;iBAAM,CAAC;gBACN,qBAAqB;gBACrB,IAAI,KAAK,KAAK,CAAC,EAAE,EAAE,CAAC;oBAClB,KAAK,GAAG,CAAC,CAAC,CAAC;oBACX,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;oBACb,MAAM,GAAG,CAAC,CAAC;oBACX,MAAM,KAAK,CAAC;gBACd,CAAC;qBAAM,CAAC;oBACN,EAAE,GAAG,GAAG,CAAC;oBACT,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;oBACxC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;oBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;oBAC5D,CAAC,CAAC,GAAG,EAAE,CAAC;oBACR,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;wBAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;oBAC/B,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC;oBACjB,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;oBACZ,IAAI,CAAC,CAAC,EAAE,KAAK,CAAC;wBAAE,SAAS,KAAK,CAAC;oBAC/B,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;oBACT,UAAU,EAAE,CAAC;oBACb,SAAS,KAAK,CAAC;gBACjB,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,OAAO,KAAK,CAAC;AACf,CAAC"}

package/src/solver-tools/lsoda/stoda.ts

@@ -1,3 +1,15 @@
+/**
+ * TypeScript port of liblsoda (https://github.com/sdwfrost/liblsoda)
+ * Original LSODA algorithm by Linda Petzold & Alan Hindmarsh (1983)
+ *
+ * Copyright (c) 2011 Yu Feng, McWilliam Cosmology Center, Carnegie Mellon University
+ * Copyright (c) 2016 Simon Frost
+ * Copyright (c) 2026 Datagrok
+ *
+ * Licensed under the MIT License. See THIRD_PARTY_LICENSES for the full
+ * license text and provenance chain of the original code.
+ */
+
 import {LsodaContext, CCMAX, MSBP} from './common';
 import {vmnorm} from './blas';
 import {cfode} from './cfode';

package/src/solver-tools/lsoda-method.js

@@ -77,23 +77,37 @@ export function lsoda(odes, callback) {
     // 3. Integrate through checkpoints, collecting Nordsieck snapshots for dense output.
     //    Checkpoints are spaced coarsely (up to 1000) so the solver handles adaptive
     //    stepping internally; mxstep per call stays within bounds.
-    if (callback)
-        callback.onIterationStart();
     let y = [...odes.initial];
     let t = t0;
-    // Warmup: a tiny initial step bootstraps the BDF method for extremely stiff problems
-    const warmupTout = t0 + Math.min(step, 1.0) * 1e-5;
-    if (warmupTout > t0 && warmupTout < t1) {
+    // Warmup: a tiny initial step bootstraps the BDF method for extremely stiff problems.
+    // Try 1e-k for k from 4 to 10, picking the largest step that succeeds.
+    const base = Math.min(step, 1.0);
+    let warmupOk = false;
+    for (let k = 5; k <= 10; k++) {
+        const warmupTout = t0 + base * Math.pow(10, -k);
+        if (warmupTout <= t0 || warmupTout >= t1)
+            continue;
         const wr = solver.solve(y, t, warmupTout);
-        y = wr.y;
-        t = wr.t;
-        if (solver.state <= 0)
-            throw new Error(ERROR_MSG.LSODA_FAILS);
+        if (solver.state > 0) {
+            y = wr.y;
+            t = wr.t;
+            warmupOk = true;
+            break;
+        }
+        // Reset solver to retry with a smaller warmup step
+        solver.reset();
+        y = [...odes.initial];
+        t = t0;
     }
+    if (!warmupOk && t === t0)
+        throw new Error(ERROR_MSG.LSODA_FAILS);
     const gridPoints = Math.trunc((t1 - t0) / step) + 1;
     const numCheckpoints = Math.min(gridPoints, 1000);
     const cpStep = (t1 - t0) / numCheckpoints;
     for (let i = 1; i <= numCheckpoints; i++) {
+        // check whether to go on computations
+        if (callback)
+            callback.onIterationStart();
         const tout = (i < numCheckpoints) ? t0 + i * cpStep : t1;
         const result = solver.solve(y, t, tout);
         y = result.y;

package/src/solver-tools/lsoda-method.js.map

@@ -1 +1 @@
-{"version":3,"file":"lsoda-method.js","sourceRoot":"","sources":["lsoda-method.ts"],"names":[],"mappings":"AAAA;;;;;;;;;EASE;AAEF,OAAO,EAAO,SAAS,EAAC,MAAM,eAAe,CAAC;AAE9C,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAE9B,2FAA2F;AAC3F,SAAS,QAAQ,CAAC,CAAe;IAC/B,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE;QACpB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QACd,OAAO,CAAC,CAAC;IACX,CAAC,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EA0CE;AACF,MAAM,UAAU,KAAK,CAAC,IAAU,EAAE,QAAmB;IACnD,gCAAgC;IAChC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;IAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC;IAC3B,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC;IAC3B,MAAM,SAAS,GAAG,IAAI,CAAC,SAAS,CAAC;IACjC,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,CAAC,MAAM,CAAC;IAChC,MAAM,MAAM,GAAG,KAAK,CAAC;IAErB,mDAAmD;IACnD,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IACnD,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IAEnD,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE;QACjD,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM;KAC1C,CAAC,CAAC;IAEH,qFAAqF;IACrF,iFAAiF;IACjF,+DAA+D;IAC/D,IAAI,QAAQ;QACV,QAAQ,CAAC,gBAAgB,EAAE,CAAC;IAE9B,IAAI,CAAC,GAAsB,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC;IAC7C,IAAI,CAAC,GAAG,EAAE,CAAC;IAEX,qFAAqF;IACrF,MAAM,UAAU,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC;IACnD,IAAI,UAAU,GAAG,EAAE,IAAI,UAAU,GAAG,EAAE,EAAE,CAAC;QACvC,MAAM,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QAC1C,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACT,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACT,IAAI,MAAM,CAAC,KAAK,IAAI,CAAC;YACnB,MAAM,IAAI,KAAK,CAAC,SAAS,CAAC,WAAW,CAAC,CAAC;IAC3C,CAAC;IAED,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;IACpD,MAAM,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC;IAClD,MAAM,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,cAAc,CAAC;IAE1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,cAAc,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;QACzD,MAAM,MAAM,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QACxC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QACb,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QAEb,IAAI,MAAM,CAAC,KAAK,IAAI,CAAC;YACnB,MAAM,IAAI,KAAK,CAAC,SAAS,CAAC,WAAW,CAAC,CAAC;IAC3C,CAAC;IAED,gEAAgE;IAChE,MAAM,KAAK,GAAG,MAAM,CAAC,cAAc,EAAE,CAAC;IACtC,MAAM,IAAI,GAAG,KAAK,CAAC,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;IAE7C,IAAI,QAAQ;QACV,QAAQ,CAAC,uBAAuB,EAAE,CAAC;IAErC,gDAAgD;IAChD,MAAM,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC;IAC/B,MAAM,QAAQ,GAAG,KAAK,CAAe,GAAG,GAAG,CAAC,CAAC,CAAC;IAC9C,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;IACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAE9B,8DAA8D;IAC9D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAEvC,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC,QAAQ"}
+{"version":3,"file":"lsoda-method.js","sourceRoot":"","sources":["lsoda-method.ts"],"names":[],"mappings":"AAAA;;;;;;;;;EASE;AAEF,OAAO,EAAO,SAAS,EAAC,MAAM,eAAe,CAAC;AAE9C,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAE9B,2FAA2F;AAC3F,SAAS,QAAQ,CAAC,CAAe;IAC/B,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,EAAE;QACpB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QACd,OAAO,CAAC,CAAC;IACX,CAAC,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EA0CE;AACF,MAAM,UAAU,KAAK,CAAC,IAAU,EAAE,QAAmB;IACnD,gCAAgC;IAChC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;IAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC;IAC3B,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC;IAC3B,MAAM,SAAS,GAAG,IAAI,CAAC,SAAS,CAAC;IACjC,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,CAAC,MAAM,CAAC;IAChC,MAAM,MAAM,GAAG,KAAK,CAAC;IAErB,mDAAmD;IACnD,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IACnD,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IAEnD,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE;QACjD,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM;KAC1C,CAAC,CAAC;IAEH,qFAAqF;IACrF,iFAAiF;IACjF,+DAA+D;IAE/D,IAAI,CAAC,GAAsB,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC;IAC7C,IAAI,CAAC,GAAG,EAAE,CAAC;IAEX,sFAAsF;IACtF,uEAAuE;IACvE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;IACjC,IAAI,QAAQ,GAAG,KAAK,CAAC;IACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,UAAU,GAAG,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QAChD,IAAI,UAAU,IAAI,EAAE,IAAI,UAAU,IAAI,EAAE;YACtC,SAAS;QAEX,MAAM,EAAE,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,UAAU,CAAC,CAAC;QAC1C,IAAI,MAAM,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC;YACrB,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YACT,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;YACT,QAAQ,GAAG,IAAI,CAAC;YAChB,MAAM;QACR,CAAC;QAED,mDAAmD;QACnD,MAAM,CAAC,KAAK,EAAE,CAAC;QACf,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,CAAC;QACtB,CAAC,GAAG,EAAE,CAAC;IACT,CAAC;IAED,IAAI,CAAC,QAAQ,IAAI,CAAC,KAAK,EAAE;QACvB,MAAM,IAAI,KAAK,CAAC,SAAS,CAAC,WAAW,CAAC,CAAC;IAEzC,MAAM,UAAU,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;IACpD,MAAM,cAAc,GAAG,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC;IAClD,MAAM,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,cAAc,CAAC;IAE1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,cAAc,EAAE,CAAC,EAAE,EAAE,CAAC;QACzC,sCAAsC;QACtC,IAAI,QAAQ;YACV,QAAQ,CAAC,gBAAgB,EAAE,CAAC;QAE9B,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC;QACzD,MAAM,MAAM,GAAG,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QACxC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QACb,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC;QAEb,IAAI,MAAM,CAAC,KAAK,IAAI,CAAC;YACnB,MAAM,IAAI,KAAK,CAAC,SAAS,CAAC,WAAW,CAAC,CAAC;IAC3C,CAAC;IAED,gEAAgE;IAChE,MAAM,KAAK,GAAG,MAAM,CAAC,cAAc,EAAE,CAAC;IACtC,MAAM,IAAI,GAAG,KAAK,CAAC,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,CAAC;IAE7C,IAAI,QAAQ;QACV,QAAQ,CAAC,uBAAuB,EAAE,CAAC;IAErC,gDAAgD;IAChD,MAAM,QAAQ,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC;IAC/B,MAAM,QAAQ,GAAG,KAAK,CAAe,GAAG,GAAG,CAAC,CAAC,CAAC;IAC9C,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;IACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAE9B,8DAA8D;IAC9D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAEvC,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC,QAAQ"}

package/src/solver-tools/lsoda-method.ts

@@ -84,27 +84,45 @@ export function lsoda(odes: ODEs, callback?: Callback): Float64Array[] {
   // 3. Integrate through checkpoints, collecting Nordsieck snapshots for dense output.
   //    Checkpoints are spaced coarsely (up to 1000) so the solver handles adaptive
   //    stepping internally; mxstep per call stays within bounds.
-  if (callback)
-    callback.onIterationStart();
 
   let y: ArrayLike<number> = [...odes.initial];
   let t = t0;
 
-  // Warmup: a tiny initial step bootstraps the BDF method for extremely stiff problems
-  const warmupTout = t0 + Math.min(step, 1.0) * 1e-5;
-  if (warmupTout > t0 && warmupTout < t1) {
+  // Warmup: a tiny initial step bootstraps the BDF method for extremely stiff problems.
+  // Try 1e-k for k from 4 to 10, picking the largest step that succeeds.
+  const base = Math.min(step, 1.0);
+  let warmupOk = false;
+  for (let k = 5; k <= 10; k++) {
+    const warmupTout = t0 + base * Math.pow(10, -k);
+    if (warmupTout <= t0 || warmupTout >= t1)
+      continue;
+
     const wr = solver.solve(y, t, warmupTout);
-    y = wr.y;
-    t = wr.t;
-    if (solver.state <= 0)
-      throw new Error(ERROR_MSG.LSODA_FAILS);
+    if (solver.state > 0) {
+      y = wr.y;
+      t = wr.t;
+      warmupOk = true;
+      break;
+    }
+
+    // Reset solver to retry with a smaller warmup step
+    solver.reset();
+    y = [...odes.initial];
+    t = t0;
   }
 
+  if (!warmupOk && t === t0)
+    throw new Error(ERROR_MSG.LSODA_FAILS);
+
   const gridPoints = Math.trunc((t1 - t0) / step) + 1;
   const numCheckpoints = Math.min(gridPoints, 1000);
   const cpStep = (t1 - t0) / numCheckpoints;
 
   for (let i = 1; i <= numCheckpoints; i++) {
+    // check whether to go on computations
+    if (callback)
+      callback.onIterationStart();
+
     const tout = (i < numCheckpoints) ? t0 + i * cpStep : t1;
     const result = solver.solve(y, t, tout);
     y = result.y;

package/src/solver-tools/solver-defs.js

@@ -49,6 +49,7 @@ export var ERROR_MSG;
     ERROR_MSG["RKDP_FAILS"] = "The Dormand-Prince 5(4) method fails";
     ERROR_MSG["RK3_FAILS"] = "The Bogacki-Shampine 3(2) method fails";
     ERROR_MSG["LSODA_FAILS"] = "The LSODA method fails";
+    ERROR_MSG["CVODE_FAILS"] = "The CVODE method fails";
 })(ERROR_MSG || (ERROR_MSG = {}));
 ;
 /** Callback action

package/src/solver-tools/solver-defs.js.map

@@ -1 +1 @@
-{"version":3,"file":"solver-defs.js","sourceRoot":"","sources":["solver-defs.ts"],"names":[],"mappings":"AAAA,qBAAqB;AAmGrB,uBAAuB;AACvB,MAAM,CAAC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAEnD,uBAAuB;AACvB,MAAM,CAAC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAE7D,uBAAuB;AACvB,MAAM,CAAC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAE7D,gDAAgD;AAChD,MAAM,CAAC,MAAM,MAAM,GAAG,GAAG,CAAC;AAC1B,MAAM,CAAC,MAAM,MAAM,GAAG,CAAC,IAAI,CAAC;AAC5B,MAAM,CAAC,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC;AAC3B,MAAM,CAAC,MAAM,WAAW,GAAG,IAAI,CAAC;AAChC,MAAM,CAAC,MAAM,SAAS,GAAG,GAAG,CAAC;AAC7B,MAAM,CAAC,MAAM,SAAS,GAAG,OAAO,CAAC;AAEjC,WAAW;AACX,MAAM,CAAC,MAAM,IAAI,GAAG,KAAK,CAAC;AAC1B,MAAM,CAAC,MAAM,GAAG,GAAG,OAAO,CAAC;AAE3B,4CAA4C;AAC5C,MAAM,UAAU,WAAW,CAAC,CAAS,EAAE,CAAe,EAAE,CAAO,EAAE,GAAW,EAC1E,KAAmB,EAAE,KAAmB,EAAE,MAAoB;IAC9D,MAAM,IAAI,GAAG,CAAC,CAAC,MAAM,CAAC;IACtB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACf,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IAErB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;AAC5C,CAAC;AAED,wBAAwB;AACxB,MAAM,UAAU,QAAQ,CAAC,CAAS,EAAE,CAAe,EAAE,CAAO,EAAE,GAAW,EACvE,KAAmB,EAAE,KAAmB,EAAE,MAAoB;IAC9D,MAAM,IAAI,GAAG,CAAC,CAAC,MAAM,CAAC;IACtB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IAEf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,EAAE,CAAC,EAAE,CAAC;QAC9B,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;QACZ,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QAEf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAErD,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;IACd,CAAC;AACH,CAAC;AAED;;EAEE;AACF,MAAM,CAAN,IAAY,SAUX;AAVD,WAAY,SAAS;IACnB,sEAAyD,CAAA;IACzD,uDAA0C,CAAA;IAC1C,yDAA4C,CAAA;IAC5C,qEAAwD,CAAA;IACxD,qEAAwD,CAAA;IACxD,qEAAwD,CAAA;IACxD,gEAAmD,CAAA;IACnD,iEAAoD,CAAA;IACpD,mDAAsC,CAAA;AACxC,CAAC,EAVW,SAAS,KAAT,SAAS,QAUpB;AAAA,CAAC;AAEF;;GAEG;AACH,MAAM,OAAO,cAAe,SAAQ,KAAK;IACvC,YAAY,GAAW;QACrB,KAAK,CAAC,GAAG,CAAC,CAAC;IACb,CAAC;CACF;AAED;;GAEG;AACH,MAAM,CAAN,IAAY,eAGX;AAHD,WAAY,eAAe;IACzB,mDAAgC,CAAA;IAChC,oCAAiB,CAAA;AACnB,CAAC,EAHW,eAAe,KAAf,eAAe,QAG1B"}
+{"version":3,"file":"solver-defs.js","sourceRoot":"","sources":["solver-defs.ts"],"names":[],"mappings":"AAAA,qBAAqB;AAmGrB,uBAAuB;AACvB,MAAM,CAAC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAEnD,uBAAuB;AACvB,MAAM,CAAC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAE7D,uBAAuB;AACvB,MAAM,CAAC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAE7D,gDAAgD;AAChD,MAAM,CAAC,MAAM,MAAM,GAAG,GAAG,CAAC;AAC1B,MAAM,CAAC,MAAM,MAAM,GAAG,CAAC,IAAI,CAAC;AAC5B,MAAM,CAAC,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC;AAC3B,MAAM,CAAC,MAAM,WAAW,GAAG,IAAI,CAAC;AAChC,MAAM,CAAC,MAAM,SAAS,GAAG,GAAG,CAAC;AAC7B,MAAM,CAAC,MAAM,SAAS,GAAG,OAAO,CAAC;AAEjC,WAAW;AACX,MAAM,CAAC,MAAM,IAAI,GAAG,KAAK,CAAC;AAC1B,MAAM,CAAC,MAAM,GAAG,GAAG,OAAO,CAAC;AAE3B,4CAA4C;AAC5C,MAAM,UAAU,WAAW,CAAC,CAAS,EAAE,CAAe,EAAE,CAAO,EAAE,GAAW,EAC1E,KAAmB,EAAE,KAAmB,EAAE,MAAoB;IAC9D,MAAM,IAAI,GAAG,CAAC,CAAC,MAAM,CAAC;IACtB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACf,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IAErB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,EAAE,CAAC;QAC3B,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;AAC5C,CAAC;AAED,wBAAwB;AACxB,MAAM,UAAU,QAAQ,CAAC,CAAS,EAAE,CAAe,EAAE,CAAO,EAAE,GAAW,EACvE,KAAmB,EAAE,KAAmB,EAAE,MAAoB;IAC9D,MAAM,IAAI,GAAG,CAAC,CAAC,MAAM,CAAC;IACtB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IAEf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,EAAE,CAAC,EAAE,CAAC;QAC9B,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;QACZ,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;QAEf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAErD,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;IACd,CAAC;AACH,CAAC;AAED;;EAEE;AACF,MAAM,CAAN,IAAY,SAWX;AAXD,WAAY,SAAS;IACnB,sEAAyD,CAAA;IACzD,uDAA0C,CAAA;IAC1C,yDAA4C,CAAA;IAC5C,qEAAwD,CAAA;IACxD,qEAAwD,CAAA;IACxD,qEAAwD,CAAA;IACxD,gEAAmD,CAAA;IACnD,iEAAoD,CAAA;IACpD,mDAAsC,CAAA;IACtC,mDAAsC,CAAA;AACxC,CAAC,EAXW,SAAS,KAAT,SAAS,QAWpB;AAAA,CAAC;AAEF;;GAEG;AACH,MAAM,OAAO,cAAe,SAAQ,KAAK;IACvC,YAAY,GAAW;QACrB,KAAK,CAAC,GAAG,CAAC,CAAC;IACb,CAAC;CACF;AAED;;GAEG;AACH,MAAM,CAAN,IAAY,eAGX;AAHD,WAAY,eAAe;IACzB,mDAAgC,CAAA;IAChC,oCAAiB,CAAA;AACnB,CAAC,EAHW,eAAe,KAAf,eAAe,QAG1B"}

package/src/solver-tools/solver-defs.ts

@@ -159,6 +159,7 @@ export enum ERROR_MSG {
   RKDP_FAILS = 'The Dormand-Prince 5(4) method fails',
   RK3_FAILS = 'The Bogacki-Shampine 3(2) method fails',
   LSODA_FAILS = 'The LSODA method fails',
+  CVODE_FAILS = 'The CVODE method fails',
 };
 
 /** Callback action

package/src/tests/cvode-lib-correctness.test.js

@@ -0,0 +1,404 @@
+import { Cvode, cvodeCreate, cvodeInit, cvode, cvodeSStolerances, cvodeSVtolerances, cvodeSetLinearSolver, CV_BDF, CV_NORMAL, CV_SUCCESS, CV_ROOT_RETURN, CV_TOO_MUCH_WORK, } from '../solver-tools/cvode';
+import { corrProbs } from '../examples/corr-probs';
+// ─── Robertson chemical kinetics (original test) ────────────────────────────
+function robertsonRhs(t, y, ydot) {
+    ydot[0] = 1.0e4 * y[1] * y[2] - 0.04 * y[0];
+    ydot[2] = 3.0e7 * y[1] * y[1];
+    ydot[1] = -(ydot[0] + ydot[2]);
+}
+/** Low-level RHS for Robertson (returns int). */
+function robertsonRhsLowLevel(t, y, ydot, _userData) {
+    robertsonRhs(t, y, ydot);
+    return 0;
+}
+const robertsonExpected = [
+    { t: 4.0000e-01, y: [9.851712e-01, 3.386380e-05, 1.479493e-02] },
+    { t: 4.0000e+00, y: [9.055333e-01, 2.240655e-05, 9.444430e-02] },
+    { t: 4.0000e+01, y: [7.158403e-01, 9.186334e-06, 2.841505e-01] },
+    { t: 4.0000e+02, y: [4.505250e-01, 3.222964e-06, 5.494717e-01] },
+    { t: 4.0000e+03, y: [1.831976e-01, 8.941773e-07, 8.168015e-01] },
+    { t: 4.0000e+04, y: [3.898729e-02, 1.621940e-07, 9.610125e-01] },
+    { t: 4.0000e+05, y: [4.936362e-03, 1.984221e-08, 9.950636e-01] },
+    { t: 4.0000e+06, y: [5.161833e-04, 2.065787e-09, 9.994838e-01] },
+    { t: 4.0000e+07, y: [5.179804e-05, 2.072027e-10, 9.999482e-01] },
+    { t: 4.0000e+08, y: [5.283675e-06, 2.113481e-11, 9.999947e-01] },
+    { t: 4.0000e+09, y: [4.658667e-07, 1.863468e-12, 9.999995e-01] },
+    { t: 4.0000e+10, y: [1.431100e-08, 5.724404e-14, 1.000000e+00] },
+];
+// Robertson solver tolerances for CVODE.
+const ROBERTSON_SOLVER_RTOL = 1e-8;
+const ROBERTSON_SOLVER_ATOL = new Float64Array([1e-10, 1e-14, 1e-10]);
+// Number of Robertson output points to check with strict relative tolerance.
+// At very late times (t>=4e8) the solution components become extremely small
+// (y1~1e-5, y2~1e-10) and accumulated integration error over many orders of
+// magnitude of time makes pure relative error checks impractical. We check
+// the first 9 points (up to t=4e7) with strict tolerance and use qualitative
+// checks for the remaining points.
+const ROBERTSON_STRICT_POINTS = 9;
+const ROBERTSON_CHECK_TOL = 5e-3;
+// ─── Helpers ─────────────────────────────────────────────────────────────────
+function relError(computed, reference) {
+    if (reference === 0)
+        return Math.abs(computed);
+    return Math.abs((computed - reference) / reference);
+}
+/**
+ * Solves an ODE with the Cvode class and checks against exact solution at the final time.
+ * Uses combined absolute + relative tolerance: |computed - exact| < checkTol * |exact| + checkAbsTol.
+ */
+function solveAndCheck(ode, exactSolution, lmm, checkTol = 1e-3, checkAbsTol = 1e-6) {
+    const neq = ode.initial.length;
+    const tol = ode.tolerance;
+    const solver = new Cvode(ode.func, neq, ode.arg.start, Float64Array.from(ode.initial), {
+        lmm,
+        rtol: tol,
+        atol: tol,
+    });
+    let t = ode.arg.start;
+    const step = ode.arg.step;
+    const finish = ode.arg.finish;
+    let result;
+    let tout = t + step;
+    while (tout <= finish + step / 2) {
+        result = solver.solve(tout);
+        expect(result.flag).toBe(CV_SUCCESS);
+        t = result.t;
+        tout += step;
+    }
+    // Check final solution against exact
+    const exact = exactSolution(t);
+    const y = result.y;
+    for (let j = 0; j < neq; j++) {
+        const absErr = Math.abs(y[j] - exact[j]);
+        const allowed = checkTol * Math.abs(exact[j]) + checkAbsTol;
+        expect(absErr).toBeLessThan(allowed);
+    }
+}
+/**
+ * Dense output helper: solve an ODE to evenly spaced output points, then use getDky
+ * to verify the interpolation matches the exact solution.
+ */
+function denseExactCheck(ode, exactSolution, lmm, checkTol = 1e-3, checkAbsTol = 1e-6) {
+    const neq = ode.initial.length;
+    const tol = ode.tolerance;
+    const solver = new Cvode(ode.func, neq, ode.arg.start, Float64Array.from(ode.initial), {
+        lmm,
+        rtol: tol,
+        atol: tol,
+    });
+    const start = ode.arg.start;
+    const finish = ode.arg.finish;
+    const nOutputPts = 100;
+    let nChecked = 0;
+    for (let i = 1; i <= nOutputPts; i++) {
+        const tout = start + (finish - start) * i / nOutputPts;
+        const result = solver.solve(tout);
+        expect(result.flag).toBe(CV_SUCCESS);
+        // getDky at the output time should match the solve result
+        const dky = solver.getDky(result.t, 0);
+        const exact = exactSolution(result.t);
+        for (let j = 0; j < neq; j++) {
+            // Dense output should match the direct solve result
+            expect(Math.abs(dky[j] - result.y[j])).toBeLessThan(1e-12);
+            // Dense output should match the exact solution
+            const absErr = Math.abs(dky[j] - exact[j]);
+            const allowed = checkTol * Math.abs(exact[j]) + checkAbsTol;
+            expect(absErr).toBeLessThan(allowed);
+        }
+        nChecked++;
+    }
+    // Ensure we actually checked points
+    expect(nChecked).toBe(nOutputPts);
+}
+// ─── Tests ───────────────────────────────────────────────────────────────────
+const nonStiffProbs = corrProbs.filter((p) => p.odes.name.startsWith('non-stiff'));
+const stiffProbs = corrProbs.filter((p) => p.odes.name.startsWith('stiff'));
+describe('CVODE LIB', () => {
+    it('should solve Robertson chemical kinetics with Cvode class', () => {
+        const neq = 3;
+        const solver = new Cvode(robertsonRhs, neq, 0.0, Float64Array.from([1.0, 0.0, 0.0]), {
+            lmm: 'bdf',
+            rtol: ROBERTSON_SOLVER_RTOL,
+            atol: ROBERTSON_SOLVER_ATOL,
+        });
+        let tout = 0.4;
+        for (let iout = 0; iout < 12; iout++) {
+            const result = solver.solve(tout);
+            expect(result.flag).toBe(CV_SUCCESS);
+            const exp = robertsonExpected[iout];
+            if (iout < ROBERTSON_STRICT_POINTS) {
+                for (let j = 0; j < 3; j++) {
+                    const err = relError(result.y[j], exp.y[j]);
+                    expect(err).toBeLessThanOrEqual(ROBERTSON_CHECK_TOL);
+                }
+            }
+            else {
+                const sum = result.y[0] + result.y[1] + result.y[2];
+                expect(Math.abs(sum - 1.0)).toBeLessThan(1e-4);
+                expect(result.y[2]).toBeGreaterThan(0.999);
+            }
+            tout *= 10.0;
+        }
+    });
+    it('should solve Robertson chemical kinetics with the low-level API', () => {
+        const neq = 3;
+        const y0 = Float64Array.from([1.0, 0.0, 0.0]);
+        const mem = cvodeCreate(CV_BDF);
+        cvodeInit(mem, robertsonRhsLowLevel, 0.0, y0);
+        cvodeSVtolerances(mem, ROBERTSON_SOLVER_RTOL, ROBERTSON_SOLVER_ATOL);
+        cvodeSetLinearSolver(mem);
+        const yout = new Float64Array(neq);
+        let tout = 0.4;
+        for (let iout = 0; iout < 12; iout++) {
+            const { flag, t } = cvode(mem, tout, yout, CV_NORMAL);
+            expect(flag).toBe(CV_SUCCESS);
+            const exp = robertsonExpected[iout];
+            if (iout < ROBERTSON_STRICT_POINTS) {
+                for (let j = 0; j < 3; j++) {
+                    const err = relError(yout[j], exp.y[j]);
+                    expect(err).toBeLessThanOrEqual(ROBERTSON_CHECK_TOL);
+                }
+            }
+            else {
+                const sum = yout[0] + yout[1] + yout[2];
+                expect(Math.abs(sum - 1.0)).toBeLessThan(1e-4);
+                expect(yout[2]).toBeGreaterThan(0.999);
+            }
+            tout *= 10.0;
+        }
+    });
+    describe('non-stiff problems (Adams)', () => {
+        nonStiffProbs.forEach(({ odes, exact }) => {
+            it(`should solve ${odes.name} problem`, () => {
+                solveAndCheck(odes, exact, 'adams');
+            });
+        });
+    });
+    describe('stiff problems (BDF)', () => {
+        stiffProbs.forEach(({ odes, exact }) => {
+            it(`should solve ${odes.name} problem`, () => {
+                solveAndCheck(odes, exact, 'bdf');
+            });
+        });
+    });
+    // ─── Rootfinding ─────────────────────────────────────────────────────────
+    describe('rootfinding', () => {
+        it('should detect roots in Robertson system', () => {
+            const neq = 3;
+            // Root functions: g1 = y1 - 1e-4, g2 = y3 - 0.01
+            const rootFn = (t, y, gout) => {
+                gout[0] = y[0] - 1e-4;
+                gout[1] = y[2] - 0.01;
+            };
+            const solver = new Cvode(robertsonRhs, neq, 0.0, Float64Array.from([1.0, 0.0, 0.0]), {
+                lmm: 'bdf',
+                rtol: ROBERTSON_SOLVER_RTOL,
+                atol: ROBERTSON_SOLVER_ATOL,
+                rootFn,
+                nRootFns: 2,
+            });
+            const rootTimes = [];
+            const rootIndices = [];
+            // Advance through the problem, collecting root returns
+            let tout = 0.4;
+            const tEnd = 4e10;
+            while (tout <= tEnd) {
+                const result = solver.solve(tout);
+                if (result.flag === CV_ROOT_RETURN) {
+                    rootTimes.push(result.t);
+                    rootIndices.push(result.rootsFound);
+                    // Continue from the root time; do not advance tout yet
+                    continue;
+                }
+                expect(result.flag).toBe(CV_SUCCESS);
+                tout *= 10.0;
+            }
+            // Expect at least 2 roots: y3 reaches 0.01 first, then y1 drops to 1e-4
+            expect(rootTimes.length).toBeGreaterThanOrEqual(2);
+            // First root: y3 = 0.01, near t ~ 2.6391e-1
+            const tRoot1 = rootTimes[0];
+            expect(relError(tRoot1, 2.6391e-1)).toBeLessThan(2e-2);
+            // rootsFound should indicate root function index 1 (g2)
+            expect(rootIndices[0][1]).not.toBe(0);
+            // Second root: y1 = 1e-4, near t ~ 2.0790e+7
+            const tRoot2 = rootTimes[1];
+            expect(relError(tRoot2, 2.0790e7)).toBeLessThan(2e-2);
+            // rootsFound should indicate root function index 0 (g1)
+            expect(rootIndices[1][0]).not.toBe(0);
+        });
+    });
+    // ─── Dense output tests ─────────────────────────────────────────────────
+    describe('dense output (getDky)', () => {
+        it('consistency: getDky at output points should match solve results (Robertson)', () => {
+            const neq = 3;
+            const solver = new Cvode(robertsonRhs, neq, 0.0, Float64Array.from([1.0, 0.0, 0.0]), {
+                lmm: 'bdf',
+                rtol: ROBERTSON_SOLVER_RTOL,
+                atol: ROBERTSON_SOLVER_ATOL,
+            });
+            let tout = 0.4;
+            for (let iout = 0; iout < 12; iout++) {
+                const result = solver.solve(tout);
+                expect(result.flag).toBe(CV_SUCCESS);
+                // getDky at the output time should match the solve result
+                const dky = solver.getDky(result.t, 0);
+                for (let j = 0; j < neq; j++)
+                    expect(Math.abs(dky[j] - result.y[j])).toBeLessThan(1e-12);
+                tout *= 10.0;
+            }
+        });
+        it('monotonicity: getDky at successive output points should produce finite values', () => {
+            const neq = 3;
+            const solver = new Cvode(robertsonRhs, neq, 0.0, Float64Array.from([1.0, 0.0, 0.0]), {
+                lmm: 'bdf',
+                rtol: ROBERTSON_SOLVER_RTOL,
+                atol: ROBERTSON_SOLVER_ATOL,
+            });
+            // Solve to 1000 uniformly spaced points across [0.01, 1.0]
+            const n = 1000;
+            const tMin = 0.01;
+            const tMax = 1.0;
+            let prevT = 0.0;
+            for (let i = 0; i < n; i++) {
+                const tout = tMin + (tMax - tMin) * i / (n - 1);
+                const result = solver.solve(tout);
+                expect(result.flag).toBe(CV_SUCCESS);
+                // Time should advance monotonically
+                expect(result.t).toBeGreaterThan(prevT);
+                prevT = result.t;
+                // getDky at the output point should return finite values
+                const dky = solver.getDky(result.t, 0);
+                for (let j = 0; j < neq; j++)
+                    expect(Number.isFinite(dky[j])).toBe(true);
+            }
+        });
+        it('boundary: getDky at endpoints of the current step interval', () => {
+            const neq = 3;
+            const solver = new Cvode(robertsonRhs, neq, 0.0, Float64Array.from([1.0, 0.0, 0.0]), {
+                lmm: 'bdf',
+                rtol: ROBERTSON_SOLVER_RTOL,
+                atol: ROBERTSON_SOLVER_ATOL,
+            });
+            // Advance to a point to establish state
+            const result = solver.solve(0.4);
+            expect(result.flag).toBe(CV_SUCCESS);
+            // getDky at the last output time should succeed
+            const dkyEnd = solver.getDky(result.t, 0);
+            for (let j = 0; j < neq; j++)
+                expect(Number.isFinite(dkyEnd[j])).toBe(true);
+            // Query derivative (k=1) - first derivative at output time
+            const dkyDeriv = solver.getDky(result.t, 1);
+            for (let j = 0; j < neq; j++)
+                expect(Number.isFinite(dkyDeriv[j])).toBe(true);
+        });
+        corrProbs.forEach(({ odes, exact }) => {
+            const lmm = odes.name.startsWith('stiff') ? 'bdf' : 'adams';
+            it(`dense output on ${odes.name} problem should match exact solution`, () => {
+                denseExactCheck(odes, exact, lmm);
+            });
+        });
+    });
+    // ─── Error handling ─────────────────────────────────────────────────────
+    describe('error handling', () => {
+        it('should return error for negative tolerances (low-level API)', () => {
+            const mem = cvodeCreate(CV_BDF);
+            cvodeInit(mem, robertsonRhsLowLevel, 0.0, Float64Array.from([1.0, 0.0, 0.0]));
+            const flag = cvodeSStolerances(mem, -1e-4, 1e-6);
+            expect(flag).toBeLessThan(0);
+        });
+        it('should return CV_TOO_MUCH_WORK with maxSteps: 1', () => {
+            const neq = 3;
+            const solver = new Cvode(robertsonRhs, neq, 0.0, Float64Array.from([1.0, 0.0, 0.0]), {
+                lmm: 'bdf',
+                rtol: ROBERTSON_SOLVER_RTOL,
+                atol: ROBERTSON_SOLVER_ATOL,
+                maxSteps: 1,
+            });
+            const result = solver.solve(4e10);
+            expect(result.flag).toBe(CV_TOO_MUCH_WORK);
+        });
+        it('should return error flag when RHS function fails', () => {
+            const failingRhs = (_t, _y, ydot) => {
+                ydot[0] = NaN;
+            };
+            const solver = new Cvode(failingRhs, 1, 0.0, Float64Array.from([1.0]), {
+                lmm: 'bdf',
+                rtol: 1e-4,
+                atol: 1e-6,
+            });
+            const result = solver.solve(1.0);
+            expect(result.flag).toBeLessThan(0);
+        });
+    });
+    // ─── Cvode class features ──────────────────────────────────────────────
+    describe('Cvode class features', () => {
+        it('getStats should return valid integrator statistics', () => {
+            const neq = 3;
+            const solver = new Cvode(robertsonRhs, neq, 0.0, Float64Array.from([1.0, 0.0, 0.0]), {
+                lmm: 'bdf',
+                rtol: ROBERTSON_SOLVER_RTOL,
+                atol: ROBERTSON_SOLVER_ATOL,
+            });
+            // Solve to first output point
+            solver.solve(0.4);
+            const stats = solver.getStats();
+            expect(stats.nSteps).toBeGreaterThan(0);
+            expect(stats.nRhsEvals).toBeGreaterThan(0);
+            expect(stats.nLinSolvSetups).toBeGreaterThanOrEqual(0);
+            expect(stats.nErrTestFails).toBeGreaterThanOrEqual(0);
+            expect(stats.nNonlinSolvIters).toBeGreaterThan(0);
+            expect(stats.nNonlinSolvConvFails).toBeGreaterThanOrEqual(0);
+            expect(stats.nJacEvals).toBeGreaterThanOrEqual(0);
+            expect(stats.lastOrder).toBeGreaterThanOrEqual(1);
+            expect(stats.lastStep).toBeGreaterThan(0);
+            expect(stats.currentTime).toBeGreaterThan(0);
+        });
+        it('reInit should allow solving with new initial conditions', () => {
+            const neq = 3;
+            const solver = new Cvode(robertsonRhs, neq, 0.0, Float64Array.from([1.0, 0.0, 0.0]), {
+                lmm: 'bdf',
+                rtol: ROBERTSON_SOLVER_RTOL,
+                atol: ROBERTSON_SOLVER_ATOL,
+            });
+            // Solve to first output point
+            const result1 = solver.solve(0.4);
+            expect(result1.flag).toBe(CV_SUCCESS);
+            const y1 = new Float64Array(result1.y);
+            // Re-initialize with the same initial conditions and solve again
+            solver.reInit(0.0, Float64Array.from([1.0, 0.0, 0.0]));
+            const result2 = solver.solve(0.4);
+            expect(result2.flag).toBe(CV_SUCCESS);
+            // Results should match closely
+            for (let j = 0; j < neq; j++) {
+                const err = relError(result2.y[j], y1[j]);
+                expect(err).toBeLessThanOrEqual(1e-6);
+            }
+        });
+        it('reInit with different initial conditions should produce different results', () => {
+            const neq = 1;
+            const simpleRhs = (_t, y, ydot) => {
+                ydot[0] = -y[0]; // y' = -y, solution: y(t) = y0 * exp(-t)
+            };
+            const solver = new Cvode(simpleRhs, neq, 0.0, Float64Array.from([1.0]), {
+                lmm: 'adams',
+                rtol: 1e-10,
+                atol: 1e-12,
+            });
+            // Solve with y0 = 1
+            const result1 = solver.solve(1.0);
+            expect(result1.flag).toBe(CV_SUCCESS);
+            const y1 = result1.y[0]; // should be ~ exp(-1) ~ 0.3679
+            // Re-initialize with y0 = 2
+            solver.reInit(0.0, Float64Array.from([2.0]));
+            const result2 = solver.solve(1.0);
+            expect(result2.flag).toBe(CV_SUCCESS);
+            const y2 = result2.y[0]; // should be ~ 2*exp(-1) ~ 0.7358
+            // y2 should be approximately 2 * y1
+            expect(relError(y2, 2 * y1)).toBeLessThan(1e-6);
+            // Check against exact solutions
+            expect(relError(y1, Math.exp(-1))).toBeLessThan(1e-6);
+            expect(relError(y2, 2 * Math.exp(-1))).toBeLessThan(1e-6);
+        });
+    });
+});
+//# sourceMappingURL=cvode-lib-correctness.test.js.map