@johnhenry/math-sample 0.0.0-stage → 0.1.0

package/dist/mesh.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"mesh.js","sourceRoot":"","sources":["../src/mesh.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,WAAW,EAAE,WAAW,EAAiC,KAAK,EAAiB,QAAQ,EAAE,MAAM,aAAa,CAAC;AAuCtH,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;AAEhH,SAAS,KAAK,CACZ,MAAmE,EACnE,IAAY,EACZ,IAAY,EACZ,OAAoB,EACpB,MAA6B,EAC7B,IAAa;IAEb,MAAM,OAAO,GAAG,OAAO,CAAC,OAAO,IAAI,EAAE,CAAC;IACtC,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI,IAAI,EAAE,CAAC;IAChC,WAAW,CAAC,SAAS,EAAE,OAAO,CAAC,CAAC;IAChC,WAAW,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC;IAC1B,WAAW,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC;IAC9B,WAAW,CAAC,UAAU,EAAE,IAAI,CAAC,CAAC;IAC9B,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAC1C,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,OAAO,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAClE,MAAM,EAAE,GAAG,IAAI,UAAU,CAAC,OAAO,GAAG,IAAI,CAAC,CAAC;IAC1C,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,MAAoC,CAAC;IACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9B,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;QAC9B,IAAI,MAAM;YAAE,MAAM;QAClB,MAAM,CAAC,GAAG,CAAC,KAAK,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;QAC1F,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,MAAM,CAAC,GAAG,CAAC,KAAK,OAAO,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC,CAAC;YAChG,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACvB,MAAM,CAAC,GAAG,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC;YAC1B,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,MAAM,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;gBACnF,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClB,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACtB,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACtB,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACZ,CAAC;QACH,CAAC;QACD,KAAK,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC;QACrB,IAAI,EAAE,CAAC;IACT,CAAC;IACD,MAAM,QAAQ,GAAG,IAAI,CAAC;IACtB,MAAM,WAAW,GAAG,OAAO,GAAG,QAAQ,CAAC;IACvC,MAAM,SAAS,GAAG,IAAI,YAAY,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC;IACpD,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC;IAClD,MAAM,KAAK,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC;IACvC,MAAM,GAAG,GAA6B,CAAC,QAAQ,EAAE,QAAQ,EAAE,QAAQ,CAAC,CAAC;IACrE,MAAM,GAAG,GAA6B,CAAC,CAAC,QAAQ,EAAE,CAAC,QAAQ,EAAE,CAAC,QAAQ,CAAC,CAAC;IACxE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAW,CAAC;YACrC,SAAS,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;YAC3B,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;gBACb,IAAI,GAAG,GAAI,GAAG,CAAC,CAAC,CAAY;oBAAE,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;gBAC3C,IAAI,GAAG,GAAI,GAAG,CAAC,CAAC,CAAY;oBAAE,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YAC7C,CAAC;QACH,CAAC;IACH,CAAC;IACD,wEAAwE;IACxE,MAAM,EAAE,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC,CAAC;IAC3D,MAAM,IAAI,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,EAAE,CACpC,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,QAAQ,IAAI,KAAK,CAAC,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC;IAClF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC;gBAAE,SAAS;YAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAY,GAAI,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAY,CAAC,CAAC;YACjG,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAE,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAY,GAAI,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAY,CAAC,CAAC;YACjG,IAAI,EAAE,GAAI,EAAE,CAAC,CAAC,CAAY,GAAI,EAAE,CAAC,CAAC,CAAY,GAAI,EAAE,CAAC,CAAC,CAAY,GAAI,EAAE,CAAC,CAAC,CAAY,CAAC;YACvF,IAAI,EAAE,GAAI,EAAE,CAAC,CAAC,CAAY,GAAI,EAAE,CAAC,CAAC,CAAY,GAAI,EAAE,CAAC,CAAC,CAAY,GAAI,EAAE,CAAC,CAAC,CAAY,CAAC;YACvF,IAAI,EAAE,GAAI,EAAE,CAAC,CAAC,CAAY,GAAI,EAAE,CAAC,CAAC,CAAY,GAAI,EAAE,CAAC,CAAC,CAAY,GAAI,EAAE,CAAC,CAAC,CAAY,CAAC;YACvF,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YACnC,IAAI,GAAG,GAAG,CAAC,IAAI,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,EAAE,CAAC;gBACpC,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;gBAChC,EAAE,IAAI,CAAC,CAAC;gBACR,EAAE,IAAI,CAAC,CAAC;gBACR,EAAE,IAAI,CAAC,CAAC;gBACR,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACnB,OAAO,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;gBAChB,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;gBACpB,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;YACtB,CAAC;QACH,CAAC;IACH,CAAC;IACD,4GAA4G;IAC5G,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACnB,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC;gBAAE,SAAS;YAC9F,KAAK,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,IAAI;gBACrB,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;gBACP,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;gBACP,CAAC,CAAC,EAAE,CAAC,CAAC;gBACN,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;aACC,EAAE,CAAC;gBACX,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC;oBAAE,SAAS;gBACpC,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;gBAC7B,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;oBACrE,OAAO,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAW,CAAC;oBAClC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAW,CAAC;oBAC1C,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,GAAG,CAAC,CAAW,CAAC;oBAC1C,MAAM;gBACR,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IACD,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;YACpB,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC;QACtB,CAAC;IACH,CAAC;IACD,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;QACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACrC,MAAM,CAAC,GAAG,CAAC,GAAG,OAAO,GAAG,CAAC,CAAC;YAC1B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAChB,MAAM,CAAC,GAAG,CAAC,GAAG,OAAO,CAAC;YACtB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAChB,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAClE,mDAAmD;YACnD,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE;gBAAE,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACtC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE;gBAAE,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACtC,qGAAqG;YACrG,IAAI,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC;gBAC3C,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE;oBAAE,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;qBACjC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE;oBAAE,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAC7C,CAAC;QACH,CAAC;IACH,CAAC;IACD,OAAO;QACL,SAAS;QACT,OAAO;QACP,OAAO,EAAE,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC;QAC9B,KAAK;QACL,WAAW;QACX,OAAO;QACP,IAAI,EAAE,QAAQ;QACd,MAAM,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE;QACpB,QAAQ,EAAE,QAAQ,CAAC,QAAQ,EAAE,IAAI,EAAE,MAAM,CAAC;KAC3C,CAAC;AACJ,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,aAAa,CAAC,CAAM,EAAE,OAAe,EAAE,OAAe,EAAE,UAAuB,EAAE;IAC/F,OAAO,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;AAC5F,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,uBAAuB,CACrC,CAAsB,EACtB,OAAe,EACf,OAAe,EACf,UAAuB,EAAE;IAEzB,OAAO,KAAK,CAAC,CAAC,EAAE,OAAO,EAAE,OAAO,EAAE,EAAE,GAAG,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,EAAE,SAAS,EAAE,KAAK,CAAC,CAAC;AACzF,CAAC;AAkBD,+FAA+F;AAC/F,MAAM,UAAU,gBAAgB,CAAC,CAAW,EAAE,OAAe,EAAE,UAAwB,EAAE;IACvF,MAAM,CAAC,GAAG,OAAO,CAAC,OAAO,IAAI,GAAG,CAAC;IACjC,WAAW,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;IAC1B,WAAW,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAC1C,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,MAAM,EAAE,GAAa,EAAE,CAAC;IACxB,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,MAAoC,CAAC;IACzC,MAAM,KAAK,GAAG,GAAG,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,IAAI,CAAC,GAAG,KAAK,KAAK,CAAC,EAAE,CAAC;YACpB,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC7C,IAAI,MAAM;gBAAE,MAAM;QACpB,CAAC;QACD,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,GAAG,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAChG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QACX,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACjE,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;QACf,IAAI,EAAE,CAAC;IACT,CAAC;IACD,2EAA2E;IAC3E,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,IAAI,IAAwB,CAAC,CAAC,wCAAwC;IACtE,IAAI,UAAU,GAAG,KAAK,CAAC;IACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAC9B,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAW,CAAC;QAChC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC;YACpB,UAAU,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;YAC7B,IAAI,GAAG,SAAS,CAAC;YACjB,SAAS;QACX,CAAC;QACD,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAW,CAAC;QACpC,MAAM,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAW,CAAC;QACpC,IAAI,GAAG,GAAG,UAAU,CAAC;QACrB,IAAI,IAAI,KAAK,SAAS,IAAI,OAAO,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;YACxD,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAClB,CAAC,GAAI,IAAI,CAAC,CAAC,GAAG,IAAI,CAAY,EAC9B,CAAC,GAAI,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAY,EAClC,CAAC,GAAI,IAAI,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAY,CACnC,CAAC;YACF,IAAI,CAAC,GAAG,OAAO,CAAC,OAAO;gBAAE,GAAG,GAAG,IAAI,CAAC;QACtC,CAAC;QACD,IAAI,GAAG,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YAC3B,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC;YAC9C,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;QACxB,CAAC;QACD,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACnB,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAW,CAAC,CAAC;QAC3B,IAAI,GAAG,CAAC,CAAC;QACT,UAAU,GAAG,KAAK,CAAC;IACrB,CAAC;IACD,OAAO;QACL,SAAS,EAAE,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC;QAClC,CAAC,EAAE,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC;QAC1B,KAAK,EAAE,IAAI,CAAC,MAAM;QAClB,QAAQ,EAAE,QAAQ,CAAC,IAAI,EAAE,CAAC,EAAE,MAAM,CAAC;KACpC,CAAC;AACJ,CAAC;AAED,gFAAgF;AAChF,MAAM,UAAU,UAAU,CAAC,CAAS;IAClC,MAAM,IAAI,GAA4B,EAAE,CAAC;IACzC,IAAI,KAAK,GAAG,CAAC,CAAC,CAAC;IACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QACjC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAW,CAAC,EAAE,CAAC;YAC/C,IAAI,KAAK,IAAI,CAAC;gBAAE,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;YACtC,KAAK,GAAG,CAAC,CAAC,CAAC;QACb,CAAC;aAAM,IAAI,KAAK,GAAG,CAAC;YAAE,KAAK,GAAG,CAAC,CAAC;IAClC,CAAC;IACD,IAAI,KAAK,IAAI,CAAC;QAAE,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC;IAC5C,OAAO,IAAI,CAAC;AACd,CAAC;AAqBD,+GAA+G;AAC/G,MAAM,UAAU,kBAAkB,CAChC,CAAY,EACZ,OAAe,EACf,OAAe,EACf,OAAe,EACf,UAA+B,EAAE;IAEjC,MAAM,CAAC,GAAG,OAAO,CAAC,IAAI,IAAI,CAAC,CAAC;IAC5B,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,OAAO,CAAC,KAAK,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3D,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;IACxB,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;IACxB,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAC;IACxB,WAAW,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IACjC,WAAW,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IACjC,WAAW,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC;IAC1C,MAAM,CAAC,GAAa,EAAE,CAAC;IACvB,MAAM,EAAE,GAAa,EAAE,CAAC;IACxB,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,IAAI,OAAO,GAAG,CAAC,CAAC;IAChB,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,MAAoC,CAAC;IACzC,MAAM,IAAI,GAAG,CAAC,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAClH,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;QAC/B,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QAC9B,IAAI,MAAM;YAAE,MAAM;QAClB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;YAC/B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC/B,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBAChC,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBAChC,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;gBAChC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;gBACrB,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;oBAC3B,OAAO,EAAE,CAAC;oBACV,SAAS;gBACX,CAAC;gBACD,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;gBAChB,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1B,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC1C,CAAC;QACH,CAAC;QACD,KAAK,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACrB,KAAK,EAAE,CAAC;IACV,CAAC;IACD,IAAI,YAAY,GAAG,CAAC,CAAC;IACrB,KAAK,MAAM,CAAC,IAAI,IAAI;QAAE,IAAI,CAAC,GAAG,YAAY;YAAE,YAAY,GAAG,CAAC,CAAC;IAC7D,OAAO;QACL,OAAO,EAAE,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC;QAC7B,OAAO,EAAE,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC;QAC9B,UAAU,EAAE,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC;QACnC,YAAY;QACZ,KAAK,EAAE,IAAI,CAAC,MAAM;QAClB,OAAO;QACP,QAAQ,EAAE,QAAQ,CAAC,KAAK,EAAE,EAAE,EAAE,MAAM,CAAC;KACtC,CAAC;AACJ,CAAC"}

package/dist/ode.d.ts

@@ -0,0 +1,147 @@
+import { type Budget, type Cancellable, type Progress } from "./budget.ts";
+import { type Domain, type Samples } from "./sample.ts";
+/** `dy/dx = f(x, y)` for a scalar equation. */
+export type Slope = (x: number, y: number) => number;
+/** `dy/dt = f(t, y)` for a system; returns a new array (or the same-length output) of derivatives. */
+export type System = (t: number, y: readonly number[]) => readonly number[];
+/** Autonomous planar system `[dx/dt, dy/dt] = F(x, y)`. */
+export type Field2 = (x: number, y: number) => readonly [number, number];
+export interface SlopeField {
+    columns: number;
+    rows: number;
+    /** Grid points (cell centres when `centered`, else grid nodes), `columns * rows` each, row-major from the smallest y. */
+    x: Float64Array;
+    y: Float64Array;
+    slope: Float64Array;
+    /**
+     * Unit direction `(dirX, dirY)` of each segment in data space (`(1, m) / |(1, m)|`; a vertical slope `+-Infinity`
+     * is `(0, +-1)`). `NaN` where the slope is `NaN` (undefined there).
+     */
+    dirX: Float64Array;
+    dirY: Float64Array;
+    progress: Progress;
+}
+export interface SlopeFieldOptions extends Budget, Cancellable {
+    columns?: number;
+    rows?: number;
+}
+/** Evaluate the slope field of `dy/dx = f(x, y)` on a `columns` x `rows` grid including both edges. */
+export declare function sampleSlopeField(f: Slope, xDomain: Domain, yDomain: Domain, options?: SlopeFieldOptions): SlopeField;
+export type StopReason = "end" | "blowup" | "maxSteps" | "stopped" | "stepTooSmall";
+export interface Trajectory {
+    /** Times of the accepted steps, including the start. */
+    t: Float64Array;
+    /** `steps * dim` states, row-major. */
+    y: Float64Array;
+    dim: number;
+    steps: number;
+    reason: StopReason;
+}
+export interface OdeOptions extends Cancellable {
+    /** Stop (reason `"stopped"`) when this returns true for an accepted state. */
+    stop?: (t: number, y: readonly number[]) => boolean;
+    /** Abort with `"blowup"` when any component is non-finite or its magnitude exceeds this. Default 1e12. */
+    blowup?: number;
+    /** Hard cap on accepted steps. Default 100000. */
+    maxSteps?: number;
+}
+export interface Rk4Options extends OdeOptions {
+    /** Number of equal steps from `t0` to `tEnd`. Default 200. */
+    steps?: number;
+}
+export interface AdaptiveOdeOptions extends OdeOptions {
+    /** Absolute tolerance per step. Default 1e-8. */
+    atol?: number;
+    /** Relative tolerance per step. Default 1e-6. */
+    rtol?: number;
+    /** First step size; default `|tEnd - t0| / 100`. */
+    h0?: number;
+    /** Smallest step; below it the integration stops with `"stepTooSmall"`. Default `|tEnd - t0| * 1e-12`. */
+    hMin?: number;
+    /** Largest step. Default `|tEnd - t0|`. */
+    hMax?: number;
+}
+/**
+ * Fixed-step classical Runge-Kutta (order 4) for a system `y' = f(t, y)` from `t0` to `tEnd` (which may be
+ * smaller than `t0`: integrates backward). Stops early, with a reason, on blow-up or `options.stop`.
+ */
+export declare function integrateRk4(f: System, t0: number, y0: readonly number[], tEnd: number, options?: Rk4Options): Trajectory;
+/**
+ * Adaptive Dormand-Prince (RK45) with error control, forward or backward in `t`. Stops with `"blowup"` near a
+ * finite-time singularity, or `"stepTooSmall"` when the step would fall below `hMin`.
+ */
+export declare function integrateAdaptive(f: System, t0: number, y0: readonly number[], tEnd: number, options?: AdaptiveOdeOptions): Trajectory;
+export interface SolutionCurveOptions extends AdaptiveOdeOptions {
+    /** `"adaptive"` (default) or `"rk4"` with `steps` per direction. */
+    method?: "adaptive" | "rk4";
+    steps?: number;
+    /** Curve points with |y| beyond this end the curve (a blow-up asymptote). Default 1e6. */
+    yLimit?: number;
+}
+/**
+ * The solution curve of `dy/dx = f(x, y)` through `(x0, y0)`, extended both ways across `xDomain`, as
+ * `Samples` (x ascending). It ends where the solution blows up, so a vertical asymptote leaves a clean
+ * end, not a spike. `reasonBefore` / `reasonAfter` say why each side stopped.
+ */
+export declare function sampleSolutionCurve(f: Slope, x0: number, y0: number, xDomain: Domain, options?: SolutionCurveOptions): Samples & {
+    reasonBefore: StopReason;
+    reasonAfter: StopReason;
+};
+export type FixedPointType = "stable node" | "unstable node" | "saddle" | "stable spiral" | "unstable spiral" | "center" | "stable star or degenerate node" | "unstable star or degenerate node" | "non-hyperbolic";
+export interface FixedPoint {
+    x: number;
+    y: number;
+    /** `[[dFx/dx, dFx/dy], [dFy/dx, dFy/dy]]` by central differences. */
+    jacobian: [[number, number], [number, number]];
+    trace: number;
+    determinant: number;
+    /** Eigenvalues as `[re, im]` pairs (complex conjugates for a spiral or center). */
+    eigenvalues: [[number, number], [number, number]];
+    type: FixedPointType;
+}
+/** Numerical Jacobian of a planar field at a point, by central differences. */
+export declare function jacobian2(F: Field2, x: number, y: number): [[number, number], [number, number]];
+/** Classify an equilibrium from its 2x2 Jacobian (linearisation; `tolerance` is the zero threshold, default 1e-7). */
+export declare function classifyFixedPoint(J: [[number, number], [number, number]], tolerance?: number): Pick<FixedPoint, "trace" | "determinant" | "eigenvalues" | "type">;
+export interface PhasePortraitOptions extends Budget, Cancellable {
+    /** Field arrows per axis (`[columns, rows]` or one number). Default 15. */
+    grid?: number | [number, number];
+    /** Trajectory start points. Default: a 5 x 5 interior grid. */
+    seeds?: Array<[number, number]>;
+    /** Integrate each seed forward to `+tMax` and backward to `-tMax`. Default: 10. */
+    tMax?: number;
+    /** Search for fixed points with Newton from this many starts per axis. Default 8. 0 disables. */
+    fixedPointStarts?: number;
+    /** Trajectories are cut when they leave the window grown by this factor on each side. Default 0.5. */
+    margin?: number;
+}
+export interface PhasePortrait {
+    field: {
+        columns: number;
+        rows: number;
+        x: Float64Array;
+        y: Float64Array;
+        /** The vector at each node; `NaN` where `F` is non-finite. */
+        u: Float64Array;
+        v: Float64Array;
+        magnitude: Float64Array;
+    };
+    /** One entry per seed: backward half reversed then forward half, joined at the seed. */
+    trajectories: Array<{
+        seed: [number, number];
+        x: Float64Array;
+        y: Float64Array;
+        reasonBefore: StopReason;
+        reasonAfter: StopReason;
+    }>;
+    fixedPoints: FixedPoint[];
+    progress: Progress;
+}
+/**
+ * Phase portrait of an autonomous planar system: the vector field on a grid, trajectories from seed points
+ * (forward and backward, cut at the window's margin or at blow-up), and fixed points classified from the Jacobian.
+ * With a budget, `progress.rows` counts field rows and the trajectories and fixed points are only computed
+ * when the field finished.
+ */
+export declare function samplePhasePortrait(F: Field2, xDomain: Domain, yDomain: Domain, options?: PhasePortraitOptions): PhasePortrait;
+//# sourceMappingURL=ode.d.ts.map

package/dist/ode.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ode.d.ts","sourceRoot":"","sources":["../src/ode.ts"],"names":[],"mappings":"AAAA,OAAO,EAA4B,KAAK,MAAM,EAAE,KAAK,WAAW,EAAS,KAAK,QAAQ,EAAY,MAAM,aAAa,CAAC;AACtH,OAAO,EAAE,KAAK,MAAM,EAAwB,KAAK,OAAO,EAAE,MAAM,aAAa,CAAC;AAE9E,+CAA+C;AAC/C,MAAM,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,MAAM,CAAC;AACrD,sGAAsG;AACtG,MAAM,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,SAAS,MAAM,EAAE,CAAC;AAC5E,2DAA2D;AAC3D,MAAM,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,KAAK,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;AAEzE,MAAM,WAAW,UAAU;IACzB,OAAO,EAAE,MAAM,CAAC;IAChB,IAAI,EAAE,MAAM,CAAC;IACb,yHAAyH;IACzH,CAAC,EAAE,YAAY,CAAC;IAChB,CAAC,EAAE,YAAY,CAAC;IAChB,KAAK,EAAE,YAAY,CAAC;IACpB;;;OAGG;IACH,IAAI,EAAE,YAAY,CAAC;IACnB,IAAI,EAAE,YAAY,CAAC;IACnB,QAAQ,EAAE,QAAQ,CAAC;CACpB;AAED,MAAM,WAAW,iBAAkB,SAAQ,MAAM,EAAE,WAAW;IAC5D,OAAO,CAAC,EAAE,MAAM,CAAC;IACjB,IAAI,CAAC,EAAE,MAAM,CAAC;CACf;AAED,uGAAuG;AACvG,wBAAgB,gBAAgB,CAC9B,CAAC,EAAE,KAAK,EACR,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,OAAO,GAAE,iBAAsB,GAC9B,UAAU,CAmDZ;AAED,MAAM,MAAM,UAAU,GAAG,KAAK,GAAG,QAAQ,GAAG,UAAU,GAAG,SAAS,GAAG,cAAc,CAAC;AAEpF,MAAM,WAAW,UAAU;IACzB,wDAAwD;IACxD,CAAC,EAAE,YAAY,CAAC;IAChB,uCAAuC;IACvC,CAAC,EAAE,YAAY,CAAC;IAChB,GAAG,EAAE,MAAM,CAAC;IACZ,KAAK,EAAE,MAAM,CAAC;IACd,MAAM,EAAE,UAAU,CAAC;CACpB;AAED,MAAM,WAAW,UAAW,SAAQ,WAAW;IAC7C,8EAA8E;IAC9E,IAAI,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,OAAO,CAAC;IACpD,0GAA0G;IAC1G,MAAM,CAAC,EAAE,MAAM,CAAC;IAChB,kDAAkD;IAClD,QAAQ,CAAC,EAAE,MAAM,CAAC;CACnB;AAED,MAAM,WAAW,UAAW,SAAQ,UAAU;IAC5C,8DAA8D;IAC9D,KAAK,CAAC,EAAE,MAAM,CAAC;CAChB;AAED,MAAM,WAAW,kBAAmB,SAAQ,UAAU;IACpD,iDAAiD;IACjD,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,iDAAiD;IACjD,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,oDAAoD;IACpD,EAAE,CAAC,EAAE,MAAM,CAAC;IACZ,0GAA0G;IAC1G,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,2CAA2C;IAC3C,IAAI,CAAC,EAAE,MAAM,CAAC;CACf;AAuCD;;;GAGG;AACH,wBAAgB,YAAY,CAC1B,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,SAAS,MAAM,EAAE,EACrB,IAAI,EAAE,MAAM,EACZ,OAAO,GAAE,UAAe,GACvB,UAAU,CAkCZ;AAgBD;;;GAGG;AACH,wBAAgB,iBAAiB,CAC/B,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,SAAS,MAAM,EAAE,EACrB,IAAI,EAAE,MAAM,EACZ,OAAO,GAAE,kBAAuB,GAC/B,UAAU,CAkFZ;AAED,MAAM,WAAW,oBAAqB,SAAQ,kBAAkB;IAC9D,oEAAoE;IACpE,MAAM,CAAC,EAAE,UAAU,GAAG,KAAK,CAAC;IAC5B,KAAK,CAAC,EAAE,MAAM,CAAC;IACf,0FAA0F;IAC1F,MAAM,CAAC,EAAE,MAAM,CAAC;CACjB;AAED;;;;GAIG;AACH,wBAAgB,mBAAmB,CACjC,CAAC,EAAE,KAAK,EACR,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,OAAO,EAAE,MAAM,EACf,OAAO,GAAE,oBAAyB,GACjC,OAAO,GAAG;IAAE,YAAY,EAAE,UAAU,CAAC;IAAC,WAAW,EAAE,UAAU,CAAA;CAAE,CA4CjE;AAED,MAAM,MAAM,cAAc,GACtB,aAAa,GACb,eAAe,GACf,QAAQ,GACR,eAAe,GACf,iBAAiB,GACjB,QAAQ,GACR,gCAAgC,GAChC,kCAAkC,GAClC,gBAAgB,CAAC;AAErB,MAAM,WAAW,UAAU;IACzB,CAAC,EAAE,MAAM,CAAC;IACV,CAAC,EAAE,MAAM,CAAC;IACV,qEAAqE;IACrE,QAAQ,EAAE,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;IAC/C,KAAK,EAAE,MAAM,CAAC;IACd,WAAW,EAAE,MAAM,CAAC;IACpB,mFAAmF;IACnF,WAAW,EAAE,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;IAClD,IAAI,EAAE,cAAc,CAAC;CACtB;AAED,+EAA+E;AAC/E,wBAAgB,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAU/F;AAED,sHAAsH;AACtH,wBAAgB,kBAAkB,CAChC,CAAC,EAAE,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,EACvC,SAAS,SAAO,GACf,IAAI,CAAC,UAAU,EAAE,OAAO,GAAG,aAAa,GAAG,aAAa,GAAG,MAAM,CAAC,CA4BpE;AAED,MAAM,WAAW,oBAAqB,SAAQ,MAAM,EAAE,WAAW;IAC/D,2EAA2E;IAC3E,IAAI,CAAC,EAAE,MAAM,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IACjC,+DAA+D;IAC/D,KAAK,CAAC,EAAE,KAAK,CAAC,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;IAChC,mFAAmF;IACnF,IAAI,CAAC,EAAE,MAAM,CAAC;IACd,iGAAiG;IACjG,gBAAgB,CAAC,EAAE,MAAM,CAAC;IAC1B,sGAAsG;IACtG,MAAM,CAAC,EAAE,MAAM,CAAC;CACjB;AAED,MAAM,WAAW,aAAa;IAC5B,KAAK,EAAE;QACL,OAAO,EAAE,MAAM,CAAC;QAChB,IAAI,EAAE,MAAM,CAAC;QACb,CAAC,EAAE,YAAY,CAAC;QAChB,CAAC,EAAE,YAAY,CAAC;QAChB,8DAA8D;QAC9D,CAAC,EAAE,YAAY,CAAC;QAChB,CAAC,EAAE,YAAY,CAAC;QAChB,SAAS,EAAE,YAAY,CAAC;KACzB,CAAC;IACF,wFAAwF;IACxF,YAAY,EAAE,KAAK,CAAC;QAClB,IAAI,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;QACvB,CAAC,EAAE,YAAY,CAAC;QAChB,CAAC,EAAE,YAAY,CAAC;QAChB,YAAY,EAAE,UAAU,CAAC;QACzB,WAAW,EAAE,UAAU,CAAC;KACzB,CAAC,CAAC;IACH,WAAW,EAAE,UAAU,EAAE,CAAC;IAC1B,QAAQ,EAAE,QAAQ,CAAC;CACpB;AAED;;;;;GAKG;AACH,wBAAgB,mBAAmB,CACjC,CAAC,EAAE,MAAM,EACT,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,OAAO,GAAE,oBAAyB,GACjC,aAAa,CAkIf"}

package/dist/ode.js

@@ -0,0 +1,466 @@
+import { assertCount, assertRange, Meter, progress } from "./budget.js";
+import { SampleCancelledError } from "./sample.js";
+/** Evaluate the slope field of `dy/dx = f(x, y)` on a `columns` x `rows` grid including both edges. */
+export function sampleSlopeField(f, xDomain, yDomain, options = {}) {
+    const columns = options.columns ?? 20;
+    const rows = options.rows ?? 20;
+    assertCount("columns", columns);
+    assertCount("rows", rows);
+    assertRange("x domain", xDomain);
+    assertRange("y domain", yDomain);
+    const meter = new Meter(options, options);
+    const n = columns * rows;
+    const x = new Float64Array(n);
+    const y = new Float64Array(n);
+    const slope = new Float64Array(n).fill(Number.NaN);
+    const dirX = new Float64Array(n).fill(Number.NaN);
+    const dirY = new Float64Array(n).fill(Number.NaN);
+    let done = 0;
+    let reason;
+    for (let j = 0; j < rows; j++) {
+        reason = meter.check(columns);
+        if (reason)
+            break;
+        const yy = j === rows - 1 ? yDomain.max : yDomain.min + ((yDomain.max - yDomain.min) * j) / (rows - 1);
+        for (let i = 0; i < columns; i++) {
+            const xx = i === columns - 1 ? xDomain.max : xDomain.min + ((xDomain.max - xDomain.min) * i) / (columns - 1);
+            const k = j * columns + i;
+            x[k] = xx;
+            y[k] = yy;
+            const m = f(xx, yy);
+            slope[k] = m;
+            if (Number.isNaN(m))
+                continue;
+            if (!Number.isFinite(m)) {
+                dirX[k] = 0;
+                dirY[k] = m > 0 ? 1 : -1;
+            }
+            else {
+                const len = Math.hypot(1, m);
+                dirX[k] = 1 / len;
+                dirY[k] = m / len;
+            }
+        }
+        meter.spend(columns);
+        done++;
+    }
+    const m = columns * done;
+    return {
+        columns,
+        rows: done,
+        x: x.slice(0, m),
+        y: y.slice(0, m),
+        slope: slope.slice(0, m),
+        dirX: dirX.slice(0, m),
+        dirY: dirY.slice(0, m),
+        progress: progress(done, rows, reason),
+    };
+}
+function bad(y, limit) {
+    for (const v of y)
+        if (!Number.isFinite(v) || Math.abs(v) > limit)
+            return true;
+    return false;
+}
+function collect(ts, ys, dim, reason) {
+    return { t: Float64Array.from(ts), y: Float64Array.from(ys), dim, steps: ts.length, reason };
+}
+function check(opts) {
+    if (opts.signal?.aborted || opts.isCancelled?.())
+        throw new SampleCancelledError();
+}
+/** One classical RK4 step. */
+function rk4Step(f, t, y, h) {
+    const n = y.length;
+    const k1 = f(t, y);
+    const k2 = f(t + h / 2, y.map((v, i) => v + (h / 2) * k1[i]));
+    const k3 = f(t + h / 2, y.map((v, i) => v + (h / 2) * k2[i]));
+    const k4 = f(t + h, y.map((v, i) => v + h * k3[i]));
+    const out = new Array(n);
+    for (let i = 0; i < n; i++)
+        out[i] =
+            y[i] +
+                (h / 6) * (k1[i] + 2 * k2[i] + 2 * k3[i] + k4[i]);
+    return out;
+}
+/**
+ * Fixed-step classical Runge-Kutta (order 4) for a system `y' = f(t, y)` from `t0` to `tEnd` (which may be
+ * smaller than `t0`: integrates backward). Stops early, with a reason, on blow-up or `options.stop`.
+ */
+export function integrateRk4(f, t0, y0, tEnd, options = {}) {
+    const steps = options.steps ?? 200;
+    assertCount("steps", steps, 1);
+    if (!Number.isFinite(t0) || !Number.isFinite(tEnd))
+        throw new RangeError("t0 and tEnd must be finite");
+    const limit = options.blowup ?? 1e12;
+    const maxSteps = Math.max(options.maxSteps ?? 100000, steps + 1); // never truncate a run the caller asked for
+    const dim = y0.length;
+    const h = (tEnd - t0) / steps;
+    const ts = [t0];
+    const ys = [...y0];
+    let y = y0;
+    let reason = "end";
+    for (let k = 0; k < steps; k++) {
+        if (ts.length >= maxSteps) {
+            reason = "maxSteps";
+            break;
+        }
+        if (k % 64 === 0)
+            check(options);
+        const t = t0 + k * h;
+        const next = rk4Step(f, t, y, h);
+        if (bad(next, limit)) {
+            reason = "blowup";
+            break;
+        }
+        y = next;
+        const tn = k === steps - 1 ? tEnd : t0 + (k + 1) * h;
+        ts.push(tn);
+        ys.push(...y);
+        if (options.stop?.(tn, y)) {
+            reason = "stopped";
+            break;
+        }
+    }
+    return collect(ts, ys, dim, reason);
+}
+// Dormand-Prince 5(4)
+const C = [0, 1 / 5, 3 / 10, 4 / 5, 8 / 9, 1, 1];
+const A = [
+    [],
+    [1 / 5],
+    [3 / 40, 9 / 40],
+    [44 / 45, -56 / 15, 32 / 9],
+    [19372 / 6561, -25360 / 2187, 64448 / 6561, -212 / 729],
+    [9017 / 3168, -355 / 33, 46732 / 5247, 49 / 176, -5103 / 18656],
+    [35 / 384, 0, 500 / 1113, 125 / 192, -2187 / 6784, 11 / 84],
+];
+const B5 = [35 / 384, 0, 500 / 1113, 125 / 192, -2187 / 6784, 11 / 84, 0];
+const B4 = [5179 / 57600, 0, 7571 / 16695, 393 / 640, -92097 / 339200, 187 / 2100, 1 / 40];
+/**
+ * Adaptive Dormand-Prince (RK45) with error control, forward or backward in `t`. Stops with `"blowup"` near a
+ * finite-time singularity, or `"stepTooSmall"` when the step would fall below `hMin`.
+ */
+export function integrateAdaptive(f, t0, y0, tEnd, options = {}) {
+    if (!Number.isFinite(t0) || !Number.isFinite(tEnd))
+        throw new RangeError("t0 and tEnd must be finite");
+    const span = Math.abs(tEnd - t0);
+    const dir = tEnd >= t0 ? 1 : -1;
+    const dim = y0.length;
+    const atol = options.atol ?? 1e-8;
+    const rtol = options.rtol ?? 1e-6;
+    const limit = options.blowup ?? 1e12;
+    const maxSteps = options.maxSteps ?? 100000;
+    const hMin = options.hMin ?? span * 1e-12;
+    const hMax = options.hMax ?? span;
+    let h = Math.min(options.h0 ?? span / 100, hMax);
+    const ts = [t0];
+    const ys = [...y0];
+    if (span === 0)
+        return collect(ts, ys, dim, "end");
+    let t = t0;
+    let y = y0;
+    let reason = "end";
+    let iterations = 0;
+    let badSteps = 0;
+    let sinceBad = Number.POSITIVE_INFINITY;
+    while (dir * (tEnd - t) > 0) {
+        if (ts.length >= maxSteps || iterations++ > maxSteps * 20) {
+            reason = "maxSteps";
+            break;
+        }
+        if (iterations % 64 === 0)
+            check(options);
+        if (h < hMin) {
+            reason = "stepTooSmall";
+            break;
+        }
+        const last = dir * (t + dir * h - tEnd) >= 0;
+        const hh = last ? Math.abs(tEnd - t) : h;
+        const step = dir * hh;
+        const k = [];
+        for (let s = 0; s < 7; s++) {
+            const arg = y.map((v, i) => {
+                let acc = 0;
+                for (let m = 0; m < s; m++)
+                    acc += A[s]?.[m] * k[m][i];
+                return v + step * acc;
+            });
+            k.push(f(t + C[s] * step, arg));
+        }
+        const y5 = y.map((v, i) => v + step * B5.reduce((acc, b, s) => acc + b * k[s][i], 0));
+        const y4 = y.map((v, i) => v + step * B4.reduce((acc, b, s) => acc + b * k[s][i], 0));
+        if (bad(y5, limit)) {
+            // may be a bad step rather than a blow-up: try smaller before giving up
+            // bad steps count only while they recur within 10 accepted steps (a Zeno approach to the limit);
+            // isolated ones over a long run are forgiven
+            if (sinceBad > 10)
+                badSteps = 0;
+            sinceBad = 0;
+            if (++badSteps > 12 || hh <= hMin * 2) {
+                reason = "blowup";
+                break;
+            }
+            h = hh / 4;
+            continue;
+        }
+        let err = 0;
+        for (let i = 0; i < dim; i++) {
+            const sc = atol + rtol * Math.max(Math.abs(y[i]), Math.abs(y5[i]));
+            err = Math.max(err, Math.abs(y5[i] - y4[i]) / sc);
+        }
+        if (err <= 1 || hh <= hMin) {
+            t = last ? tEnd : t + step;
+            y = y5;
+            sinceBad++;
+            ts.push(t);
+            ys.push(...y);
+            if (options.stop?.(t, y)) {
+                reason = "stopped";
+                break;
+            }
+        }
+        const factor = err === 0 ? 5 : Math.min(5, Math.max(0.2, 0.9 * err ** -0.2));
+        h = Math.min(hMax, hh * factor);
+    }
+    return collect(ts, ys, dim, reason);
+}
+/**
+ * The solution curve of `dy/dx = f(x, y)` through `(x0, y0)`, extended both ways across `xDomain`, as
+ * `Samples` (x ascending). It ends where the solution blows up, so a vertical asymptote leaves a clean
+ * end, not a spike. `reasonBefore` / `reasonAfter` say why each side stopped.
+ */
+export function sampleSolutionCurve(f, x0, y0, xDomain, options = {}) {
+    assertRange("x domain", xDomain);
+    const sys = (t, y) => [f(t, y[0])];
+    const limit = options.yLimit ?? 1e6;
+    const run = (end) => {
+        const o = { ...options, blowup: limit };
+        return options.method === "rk4"
+            ? integrateRk4(sys, x0, [y0], end, { ...o, steps: options.steps ?? 200 })
+            : integrateAdaptive(sys, x0, [y0], end, o);
+    };
+    const xs = [];
+    const ys = [];
+    let reasonBefore = "end";
+    let reasonAfter = "end";
+    // A start outside the window is first carried to the nearest edge, so the curve begins exactly at the edge.
+    if (x0 < xDomain.min || x0 > xDomain.max) {
+        const edge = x0 < xDomain.min ? xDomain.min : xDomain.max;
+        const toEdge = run(edge);
+        if (toEdge.reason !== "end") {
+            const r = toEdge.reason;
+            return { x: new Float64Array(0), y: new Float64Array(0), reasonBefore: r, reasonAfter: r };
+        }
+        x0 = edge;
+        y0 = toEdge.y[toEdge.steps - 1];
+    }
+    if (x0 > xDomain.min) {
+        const back = run(xDomain.min);
+        reasonBefore = back.reason;
+        for (let k = back.steps - 1; k >= 1; k--) {
+            xs.push(back.t[k]);
+            ys.push(back.y[k]);
+        }
+    }
+    xs.push(x0);
+    ys.push(y0);
+    if (x0 < xDomain.max) {
+        const fwd = run(xDomain.max);
+        reasonAfter = fwd.reason;
+        for (let k = 1; k < fwd.steps; k++) {
+            xs.push(fwd.t[k]);
+            ys.push(fwd.y[k]);
+        }
+    }
+    return { x: Float64Array.from(xs), y: Float64Array.from(ys), reasonBefore, reasonAfter };
+}
+/** Numerical Jacobian of a planar field at a point, by central differences. */
+export function jacobian2(F, x, y) {
+    const h = 1e-6 * Math.max(1, Math.abs(x), Math.abs(y));
+    const xp = F(x + h, y);
+    const xm = F(x - h, y);
+    const yp = F(x, y + h);
+    const ym = F(x, y - h);
+    return [
+        [(xp[0] - xm[0]) / (2 * h), (yp[0] - ym[0]) / (2 * h)],
+        [(xp[1] - xm[1]) / (2 * h), (yp[1] - ym[1]) / (2 * h)],
+    ];
+}
+/** Classify an equilibrium from its 2x2 Jacobian (linearisation; `tolerance` is the zero threshold, default 1e-7). */
+export function classifyFixedPoint(J, tolerance = 1e-7) {
+    const tr = J[0][0] + J[1][1];
+    const det = J[0][0] * J[1][1] - J[0][1] * J[1][0];
+    const disc = tr * tr - 4 * det;
+    const scale = Math.max(1, Math.abs(tr), Math.abs(det));
+    let eigenvalues;
+    if (disc >= 0) {
+        const s = Math.sqrt(disc);
+        eigenvalues = [
+            [(tr + s) / 2, 0],
+            [(tr - s) / 2, 0],
+        ];
+    }
+    else {
+        const s = Math.sqrt(-disc) / 2;
+        eigenvalues = [
+            [tr / 2, s],
+            [tr / 2, -s],
+        ];
+    }
+    let type;
+    if (Math.abs(det) <= tolerance * scale)
+        type = "non-hyperbolic";
+    else if (det < 0)
+        type = "saddle";
+    else if (Math.abs(tr) <= tolerance * scale)
+        type = "center";
+    else if (disc > tolerance * scale * scale)
+        type = tr < 0 ? "stable node" : "unstable node";
+    else if (disc >= -tolerance * scale * scale)
+        type = tr < 0 ? "stable star or degenerate node" : "unstable star or degenerate node";
+    else
+        type = tr < 0 ? "stable spiral" : "unstable spiral";
+    return { trace: tr, determinant: det, eigenvalues, type };
+}
+/**
+ * Phase portrait of an autonomous planar system: the vector field on a grid, trajectories from seed points
+ * (forward and backward, cut at the window's margin or at blow-up), and fixed points classified from the Jacobian.
+ * With a budget, `progress.rows` counts field rows and the trajectories and fixed points are only computed
+ * when the field finished.
+ */
+export function samplePhasePortrait(F, xDomain, yDomain, options = {}) {
+    assertRange("x domain", xDomain);
+    assertRange("y domain", yDomain);
+    const g = options.grid ?? 15;
+    const [columns, rows] = typeof g === "number" ? [g, g] : g;
+    assertCount("columns", columns);
+    assertCount("rows", rows);
+    const meter = new Meter(options, options);
+    const n = columns * rows;
+    const fx = new Float64Array(n);
+    const fy = new Float64Array(n);
+    const u = new Float64Array(n).fill(Number.NaN);
+    const v = new Float64Array(n).fill(Number.NaN);
+    const mag = new Float64Array(n).fill(Number.NaN);
+    let done = 0;
+    let reason;
+    for (let j = 0; j < rows; j++) {
+        reason = meter.check(columns);
+        if (reason)
+            break;
+        const y = j === rows - 1 ? yDomain.max : yDomain.min + ((yDomain.max - yDomain.min) * j) / (rows - 1);
+        for (let i = 0; i < columns; i++) {
+            const x = i === columns - 1 ? xDomain.max : xDomain.min + ((xDomain.max - xDomain.min) * i) / (columns - 1);
+            const k = j * columns + i;
+            fx[k] = x;
+            fy[k] = y;
+            const w = F(x, y);
+            if (Number.isFinite(w[0]) && Number.isFinite(w[1])) {
+                u[k] = w[0];
+                v[k] = w[1];
+                mag[k] = Math.hypot(w[0], w[1]);
+            }
+        }
+        meter.spend(columns);
+        done++;
+    }
+    const m = columns * done;
+    const out = {
+        field: {
+            columns,
+            rows: done,
+            x: fx.slice(0, m),
+            y: fy.slice(0, m),
+            u: u.slice(0, m),
+            v: v.slice(0, m),
+            magnitude: mag.slice(0, m),
+        },
+        trajectories: [],
+        fixedPoints: [],
+        progress: progress(done, rows, reason),
+    };
+    if (reason)
+        return out;
+    const w = xDomain.max - xDomain.min;
+    const hgt = yDomain.max - yDomain.min;
+    const mar = options.margin ?? 0.5;
+    const inside = (x, y) => x >= xDomain.min - mar * w &&
+        x <= xDomain.max + mar * w &&
+        y >= yDomain.min - mar * hgt &&
+        y <= yDomain.max + mar * hgt;
+    const sys = (_t, s) => {
+        const r = F(s[0], s[1]);
+        return [r[0], r[1]];
+    };
+    const tMax = options.tMax ?? 10;
+    const seeds = options.seeds ??
+        Array.from({ length: 25 }, (_, q) => [
+            xDomain.min + (w * ((q % 5) + 0.5)) / 5,
+            yDomain.min + (hgt * (Math.floor(q / 5) + 0.5)) / 5,
+        ]);
+    for (const seed of seeds) {
+        const o = {
+            ...options,
+            stop: (_t, s) => !inside(s[0], s[1]),
+            hMax: tMax / 20,
+        };
+        const fwd = integrateAdaptive(sys, 0, seed, tMax, o);
+        const bwd = integrateAdaptive(sys, 0, seed, -tMax, o);
+        const xs = [];
+        const ys = [];
+        for (let k = bwd.steps - 1; k >= 1; k--) {
+            xs.push(bwd.y[2 * k]);
+            ys.push(bwd.y[2 * k + 1]);
+        }
+        for (let k = 0; k < fwd.steps; k++) {
+            xs.push(fwd.y[2 * k]);
+            ys.push(fwd.y[2 * k + 1]);
+        }
+        out.trajectories.push({
+            seed,
+            x: Float64Array.from(xs),
+            y: Float64Array.from(ys),
+            reasonBefore: bwd.reason,
+            reasonAfter: fwd.reason,
+        });
+    }
+    const starts = options.fixedPointStarts ?? 8;
+    const eps = 1e-8 * Math.max(w, hgt);
+    for (let a = 0; a < starts && starts > 0; a++) {
+        for (let b = 0; b < starts; b++) {
+            let x = xDomain.min + (w * (a + 0.5)) / starts;
+            let y = yDomain.min + (hgt * (b + 0.5)) / starts;
+            let ok = false;
+            for (let it = 0; it < 50; it++) {
+                const r = F(x, y);
+                if (!Number.isFinite(r[0]) || !Number.isFinite(r[1]))
+                    break;
+                const J = jacobian2(F, x, y);
+                const det = J[0][0] * J[1][1] - J[0][1] * J[1][0];
+                if (!Number.isFinite(det) || Math.abs(det) < 1e-14)
+                    break;
+                const dx = (J[1][1] * r[0] - J[0][1] * r[1]) / det;
+                const dy = (-J[1][0] * r[0] + J[0][0] * r[1]) / det;
+                x -= dx;
+                y -= dy;
+                if (Math.hypot(dx, dy) < eps) {
+                    ok = true;
+                    break;
+                }
+            }
+            if (!ok || !inside(x, y) || x < xDomain.min || x > xDomain.max || y < yDomain.min || y > yDomain.max)
+                continue;
+            const r = F(x, y);
+            if (Math.hypot(r[0], r[1]) > 1e-6)
+                continue;
+            if (out.fixedPoints.some((p) => Math.hypot(p.x - x, p.y - y) < 1e-5 * Math.max(w, hgt)))
+                continue;
+            const J = jacobian2(F, x, y);
+            out.fixedPoints.push({ x, y, jacobian: J, ...classifyFixedPoint(J) });
+        }
+    }
+    out.fixedPoints.sort((p, q) => p.x - q.x || p.y - q.y);
+    return out;
+}
+//# sourceMappingURL=ode.js.map