universal-physics-tensor 9.0.0 → 10.0.0

package/dist/numerical/formula-mathts.js.map

@@ -1 +1 @@
-{"version":3,"file":"formula-mathts.js","sourceRoot":"","sources":["../../src/numerical/formula-mathts.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;GAkBG;AAEH,OAAO,EAAE,WAAW,EAAE,KAAK,IAAI,WAAW,EAAE,MAAM,iCAAiC,CAAC;AAEpF,OAAO,EAAE,sBAAsB,EAAE,mBAAmB,EAAE,gBAAgB,EAAE,YAAY,EAAE,sBAAsB,EAAE,MAAM,uBAAuB,CAAC;AAE5I,0DAA0D;AAC1D,MAAM,cAAc,GAAG,EAAE,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,QAAQ,EAAW,CAAC;AAkBpE;;;GAGG;AACH,SAAS,yBAAyB,CAChC,GAA0B;IAE1B,MAAM,YAAY,GAAG,IAAI,GAAG,EAAmB,CAAC;IAChD,MAAM,SAAS,GAAG,CAAC,IAAY,EAAW,EAAE;QAC1C,MAAM,MAAM,GAAG,YAAY,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACtC,IAAI,MAAM,KAAK,SAAS;YAAE,OAAO,MAAM,CAAC;QACxC,IAAI,OAAgB,CAAC;QACrB,IAAI,CAAC;YACH,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;YAC7B,OAAO,GAAG,IAAI,CAAC,CAAC,+CAA+C;QACjE,CAAC;QAAC,MAAM,CAAC;YACP,OAAO,GAAG,KAAK,CAAC,CAAC,2BAA2B;QAC9C,CAAC;QACD,YAAY,CAAC,GAAG,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;QAChC,OAAO,OAAO,CAAC;IACjB,CAAC,CAAC;IAEF,4EAA4E;IAC5E,4EAA4E;IAC5E,2EAA2E;IAC3E,uEAAuE;IACvE,sEAAsE;IACtE,6EAA6E;IAC7E,+EAA+E;IAC/E,yFAAyF;IACzF,uEAAuE;IACvE,mEAAmE;IACnE,MAAM,kBAAkB,GAAG,IAAI,GAAG,EAAmB,CAAC;IACtD,MAAM,eAAe,GAAG,CAAC,IAAY,EAAW,EAAE;QAChD,MAAM,MAAM,GAAG,kBAAkB,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC5C,IAAI,MAAM,KAAK,SAAS;YAAE,OAAO,MAAM,CAAC;QACxC,IAAI,QAAQ,GAAG,KAAK,CAAC;QACrB,IAAI,CAAC;YACH,QAAQ,GAAG,OAAO,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,EAAE,CAAC,KAAK,QAAQ,CAAC;QAC9D,CAAC;QAAC,MAAM,CAAC;YACP,QAAQ,GAAG,KAAK,CAAC;QACnB,CAAC;QACD,kBAAkB,CAAC,GAAG,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QACvC,OAAO,QAAQ,CAAC;IAClB,CAAC,CAAC;IAEF,6EAA6E;IAC7E,sEAAsE;IACtE,4DAA4D;IAC5D,MAAM,KAAK,GAAkD,EAAE,CAAC;IAChE,KAAK,MAAM,IAAI,IAAI,sBAAsB,EAAE,CAAC;QAC1C,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC;YAAE,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,EAAE,EAAE,CAAC,mBAAmB,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IACnF,CAAC;IAED,OAAO;QACL,KAAK,CAAC,IAAY;YAChB,IAAI,CAAC,IAAI,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE;gBAAE,MAAM,IAAI,YAAY,CAAC,eAAe,CAAC,CAAC;YACnE,IAAI,IAAc,CAAC;YACnB,IAAI,CAAC;gBACH,IAAI,GAAG,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;YACzB,CAAC;YAAC,OAAO,GAAG,EAAE,CAAC;gBACb,MAAM,IAAI,YAAY,CACpB,gBAAgB,GAAG,YAAY,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,CACnE,CAAC;YACJ,CAAC;YAED,wEAAwE;YACxE,wEAAwE;YACxE,wDAAwD;YACxD,MAAM,OAAO,GAAG,IAAI,GAAG,CACrB,IAAI;iBACD,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,cAAc,KAAK,IAAI,CAAC;iBACxC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC;iBACtB,MAAM,CAAC,CAAC,CAAC,EAAe,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,CAAC,CACrD,CAAC;YACF,6FAA6F;YAC7F,uFAAuF;YACvF,MAAM,WAAW,GAAG,IAAI;iBACrB,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,KAAK,IAAI,CAAC;iBACtC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC;iBAClB,MAAM,CAAC,CAAC,CAAC,EAAe,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,CAAC,CAAC;YACrD,IAAI,WAAW,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,OAAO,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC;gBAC1D,MAAM,IAAI,gBAAgB,EAAE,CAAC;YAC/B,CAAC;YACD,MAAM,aAAa,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/E,MAAM,SAAS,GAAG;gBAChB,GAAG,IAAI,GAAG,CACR,IAAI;qBACD,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,KAAK,IAAI,CAAC;qBACtC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC;qBAClB,MAAM,CAAC,CAAC,CAAC,EAAe,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,CAAC;qBACjD,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,CACzD;aACF,CAAC,IAAI,EAAE,CAAC;YAET,sEAAsE;YACtE,8EAA8E;YAC9E,MAAM,QAAQ,GAAG,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;YACnD,OAAO;gBACL,MAAM,EAAE,IAAI;gBACZ,SAAS;gBACT,QAAQ,CAAC,KAA6B;oBACpC,IAAI,aAAa,KAAK,SAAS;wBAAE,MAAM,IAAI,YAAY,CAAC,sBAAsB,CAAC,aAAa,CAAC,CAAC,CAAC;oBAC/F,IAAI,MAAe,CAAC;oBACpB,IAAI,CAAC;wBACH,MAAM,GAAG,QAAQ,CAAC,QAAQ,CAAC,EAAE,GAAG,KAAK,EAAE,GAAG,KAAK,EAAE,CAAC,CAAC;oBACrD,CAAC;oBAAC,OAAO,GAAG,EAAE,CAAC;wBACb,gEAAgE;wBAChE,2DAA2D;wBAC3D,sFAAsF;wBACtF,0FAA0F;wBAC1F,mGAAmG;wBACnG,+EAA+E;wBAC/E,MAAM,OAAO,GAAG,GAAG,YAAY,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;wBACjE,MAAM,IAAI,YAAY,CAAC,OAAO,CAAC,OAAO,CAAC,8BAA8B,EAAE,qBAAqB,CAAC,CAAC,CAAC;oBACjG,CAAC;oBACD,IAAI,OAAO,MAAM,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;wBAC3D,oEAAoE;wBACpE,2EAA2E;wBAC3E,MAAM,IAAI,YAAY,CAAC,oDAAoD,iBAAiB,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;oBAC3G,CAAC;oBACD,OAAO,MAAM,CAAC;gBAChB,CAAC;aACF,CAAC;QACJ,CAAC;KACF,CAAC;AACJ,CAAC;AAED,sEAAsE;AACtE,SAAS,iBAAiB,CAAC,MAAe;IACxC,IAAI,OAAO,MAAM,KAAK,QAAQ;QAAE,OAAO,MAAM,CAAC,MAAM,CAAC,CAAC;IACtD,IAAI,MAAM,KAAK,IAAI,IAAI,OAAO,MAAM,KAAK,QAAQ,EAAE,CAAC;QAClD,MAAM,CAAC,GAAG,MAAiE,CAAC;QAC5E,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,EAAE,CAAC;YACzD,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,QAAQ,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC;YAClF,OAAO,qBAAqB,IAAI,GAAG,CAAC;QACtC,CAAC;IACH,CAAC;IACD,OAAO,OAAO,MAAM,CAAC;AACvB,CAAC;AAED,MAAM,YAAY,GAA0B;IAC1C,KAAK,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,WAAW,CAAC,IAAI,CAAwB;CAC1D,CAAC;AAEF,2CAA2C;AAC3C,MAAM,CAAC,MAAM,mBAAmB,GAAkB,yBAAyB,CAAC,YAAY,CAAC,CAAC;AAE1F,wCAAwC;AACxC,MAAM,UAAU,YAAY,CAAC,IAAY;IACvC,OAAO,mBAAmB,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;AACzC,CAAC"}
+{"version":3,"file":"formula-mathts.js","sourceRoot":"","sources":["../../src/numerical/formula-mathts.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;GAkBG;AAEH,OAAO,EAAE,WAAW,EAAE,KAAK,IAAI,WAAW,EAAE,MAAM,iCAAiC,CAAC;AAEpF,OAAO,EAAE,sBAAsB,EAAE,mBAAmB,EAAE,gBAAgB,EAAE,YAAY,EAAE,sBAAsB,EAAE,iBAAiB,EAAE,MAAM,uBAAuB,CAAC;AAE/J,0DAA0D;AAC1D,MAAM,cAAc,GAAG,EAAE,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,QAAQ,EAAW,CAAC;AAkBpE;;;GAGG;AACH,SAAS,yBAAyB,CAChC,GAA0B;IAE1B,MAAM,YAAY,GAAG,IAAI,GAAG,EAAmB,CAAC;IAChD,MAAM,SAAS,GAAG,CAAC,IAAY,EAAW,EAAE;QAC1C,MAAM,MAAM,GAAG,YAAY,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACtC,IAAI,MAAM,KAAK,SAAS;YAAE,OAAO,MAAM,CAAC;QACxC,IAAI,OAAgB,CAAC;QACrB,IAAI,CAAC;YACH,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;YAC7B,OAAO,GAAG,IAAI,CAAC,CAAC,+CAA+C;QACjE,CAAC;QAAC,MAAM,CAAC;YACP,OAAO,GAAG,KAAK,CAAC,CAAC,2BAA2B;QAC9C,CAAC;QACD,YAAY,CAAC,GAAG,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;QAChC,OAAO,OAAO,CAAC;IACjB,CAAC,CAAC;IAEF,4EAA4E;IAC5E,4EAA4E;IAC5E,2EAA2E;IAC3E,uEAAuE;IACvE,sEAAsE;IACtE,6EAA6E;IAC7E,+EAA+E;IAC/E,yFAAyF;IACzF,uEAAuE;IACvE,mEAAmE;IACnE,MAAM,kBAAkB,GAAG,IAAI,GAAG,EAAmB,CAAC;IACtD,MAAM,eAAe,GAAG,CAAC,IAAY,EAAW,EAAE;QAChD,MAAM,MAAM,GAAG,kBAAkB,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC5C,IAAI,MAAM,KAAK,SAAS;YAAE,OAAO,MAAM,CAAC;QACxC,IAAI,QAAQ,GAAG,KAAK,CAAC;QACrB,IAAI,CAAC;YACH,QAAQ,GAAG,OAAO,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,EAAE,CAAC,KAAK,QAAQ,CAAC;QAC9D,CAAC;QAAC,MAAM,CAAC;YACP,QAAQ,GAAG,KAAK,CAAC;QACnB,CAAC;QACD,kBAAkB,CAAC,GAAG,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;QACvC,OAAO,QAAQ,CAAC;IAClB,CAAC,CAAC;IAEF,6EAA6E;IAC7E,sEAAsE;IACtE,4DAA4D;IAC5D,MAAM,KAAK,GAAkD,EAAE,CAAC;IAChE,KAAK,MAAM,IAAI,IAAI,sBAAsB,EAAE,CAAC;QAC1C,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC;YAAE,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,EAAE,EAAE,CAAC,mBAAmB,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IACnF,CAAC;IAED,OAAO;QACL,KAAK,CAAC,IAAY;YAChB,IAAI,CAAC,IAAI,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE;gBAAE,MAAM,IAAI,YAAY,CAAC,eAAe,CAAC,CAAC;YACnE,IAAI,IAAc,CAAC;YACnB,IAAI,CAAC;gBACH,IAAI,GAAG,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;YACzB,CAAC;YAAC,OAAO,GAAG,EAAE,CAAC;gBACb,MAAM,IAAI,YAAY,CACpB,gBAAgB,GAAG,YAAY,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,CACnE,CAAC;YACJ,CAAC;YAED,wEAAwE;YACxE,wEAAwE;YACxE,wDAAwD;YACxD,MAAM,OAAO,GAAG,IAAI,GAAG,CACrB,IAAI;iBACD,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,cAAc,KAAK,IAAI,CAAC;iBACxC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC;iBACtB,MAAM,CAAC,CAAC,CAAC,EAAe,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,CAAC,CACrD,CAAC;YACF,6FAA6F;YAC7F,uFAAuF;YACvF,MAAM,WAAW,GAAG,IAAI;iBACrB,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,KAAK,IAAI,CAAC;iBACtC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC;iBAClB,MAAM,CAAC,CAAC,CAAC,EAAe,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,CAAC,CAAC;YACrD,IAAI,WAAW,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,OAAO,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC;gBAC1D,MAAM,IAAI,gBAAgB,EAAE,CAAC;YAC/B,CAAC;YACD,MAAM,aAAa,GAAG,CAAC,GAAG,OAAO,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/E,MAAM,SAAS,GAAG;gBAChB,GAAG,IAAI,GAAG,CACR,IAAI;qBACD,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,KAAK,IAAI,CAAC;qBACtC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC;qBAClB,MAAM,CAAC,CAAC,CAAC,EAAe,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,CAAC;qBACjD,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,CACzD;aACF,CAAC,IAAI,EAAE,CAAC;YAET,sEAAsE;YACtE,8EAA8E;YAC9E,MAAM,QAAQ,GAAG,WAAW,CAAC,IAAI,EAAE,cAAc,CAAC,CAAC;YACnD,OAAO;gBACL,MAAM,EAAE,IAAI;gBACZ,SAAS;gBACT,QAAQ,CAAC,KAA6B;oBACpC,IAAI,aAAa,KAAK,SAAS;wBAAE,MAAM,IAAI,YAAY,CAAC,sBAAsB,CAAC,aAAa,CAAC,CAAC,CAAC;oBAC/F,IAAI,MAAe,CAAC;oBACpB,IAAI,CAAC;wBACH,MAAM,GAAG,QAAQ,CAAC,QAAQ,CAAC,EAAE,GAAG,KAAK,EAAE,GAAG,KAAK,EAAE,CAAC,CAAC;oBACrD,CAAC;oBAAC,OAAO,GAAG,EAAE,CAAC;wBACb,gEAAgE;wBAChE,2DAA2D;wBAC3D,sFAAsF;wBACtF,0FAA0F;wBAC1F,mGAAmG;wBACnG,+EAA+E;wBAC/E,MAAM,OAAO,GAAG,GAAG,YAAY,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;wBACjE,MAAM,IAAI,YAAY,CAAC,OAAO,CAAC,OAAO,CAAC,8BAA8B,EAAE,qBAAqB,CAAC,CAAC,CAAC;oBACjG,CAAC;oBACD,IAAI,OAAO,MAAM,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;wBAC3D,oEAAoE;wBACpE,2EAA2E;wBAC3E,MAAM,IAAI,iBAAiB,CAAC,oDAAoD,iBAAiB,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;oBAChH,CAAC;oBACD,OAAO,MAAM,CAAC;gBAChB,CAAC;aACF,CAAC;QACJ,CAAC;KACF,CAAC;AACJ,CAAC;AAED,sEAAsE;AACtE,SAAS,iBAAiB,CAAC,MAAe;IACxC,IAAI,OAAO,MAAM,KAAK,QAAQ;QAAE,OAAO,MAAM,CAAC,MAAM,CAAC,CAAC;IACtD,IAAI,MAAM,KAAK,IAAI,IAAI,OAAO,MAAM,KAAK,QAAQ,EAAE,CAAC;QAClD,MAAM,CAAC,GAAG,MAAiE,CAAC;QAC5E,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,IAAI,OAAO,CAAC,CAAC,EAAE,KAAK,QAAQ,EAAE,CAAC;YACzD,MAAM,IAAI,GAAG,OAAO,CAAC,CAAC,QAAQ,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,GAAG,CAAC;YAClF,OAAO,qBAAqB,IAAI,GAAG,CAAC;QACtC,CAAC;IACH,CAAC;IACD,OAAO,OAAO,MAAM,CAAC;AACvB,CAAC;AAED,MAAM,YAAY,GAA0B;IAC1C,0EAA0E;IAC1E,6EAA6E;IAC7E,gEAAgE;IAChE,KAAK,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,WAAW,CAAC,IAAI,CAAwB;CAC1D,CAAC;AAEF,2CAA2C;AAC3C,MAAM,CAAC,MAAM,mBAAmB,GAAkB,yBAAyB,CAAC,YAAY,CAAC,CAAC;AAE1F,wCAAwC;AACxC,MAAM,UAAU,YAAY,CAAC,IAAY;IACvC,OAAO,mBAAmB,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;AACzC,CAAC","sourcesContent":["/**\n * MathTS-backed scalar-formula parser.\n *\n * Backed by `@danielsimonjr/mathts-functions` (`parse(expr) → Node`,\n * `node.evaluate(scope)`). The package is a required dependency.\n *\n * SCALAR-ONLY guard: the CLI/inference contract returns a `number`. If a\n * formula evaluates to a non-number (matrix, complex, unit, function), this\n * throws a {@link FormulaError} rather than leaking MathTS types through the\n * seam.\n *\n * Bare `e` is the elementary charge. MathTS `{ physics: true, charge: 'scalar' }`\n * binds that symbol to the SI magnitude `1.602176634e-19`, so `1-e^2` is\n * ordinary arithmetic. `{ physics: true }` without `charge: 'scalar'` is a\n * Unit, and subtract does not accept it. Euler's number is `exp(x)`. The name\n * `euler` is refused. `E` stays unbound.\n *\n * @module numerical/formula-mathts\n */\n\nimport { compileExpr, parse as parseMathTs } from '@danielsimonjr/mathts-functions';\nimport type { CompiledFormula, FormulaParser } from './formula-contract.js';\nimport { BUILTIN_FUNCTION_NAMES, callBuiltinFunction, EulerNumberError, FormulaError, unknownFunctionMessage, FormulaValueError } from './formula-contract.js';\n\n/** Bare `e` is the SI magnitude, not the coulomb Unit. */\nconst PHYSICS_SCALAR = { physics: true, charge: 'scalar' } as const;\n\n/** Minimal structural shape of a MathTS AST node (the bits we use). */\ninterface MathNode {\n  evaluate(scope: Record<string, number>): unknown;\n  filter(predicate: (node: MathNode) => boolean): MathNode[];\n  toString(): string;\n  readonly isSymbolNode?: boolean;\n  readonly isFunctionNode?: boolean;\n  readonly name?: string;\n  readonly fn?: { readonly name?: string };\n}\n\n/** Minimal structural shape of the @danielsimonjr/mathts-functions module. */\ninterface MathtsFunctionsModule {\n  parse(expr: string): MathNode;\n}\n\n/**\n * Build a {@link FormulaParser} bound to an already-loaded mathts-functions\n * module. Pure (no I/O) — the dynamic import lives in the registry.\n */\nfunction createMathtsFormulaParser(\n  mod: MathtsFunctionsModule,\n): FormulaParser {\n  const builtinCache = new Map<string, boolean>();\n  const isBuiltin = (name: string): boolean => {\n    const cached = builtinCache.get(name);\n    if (cached !== undefined) return cached;\n    let builtin: boolean;\n    try {\n      mod.parse(name).evaluate({});\n      builtin = true; // pi / e / sin / … resolve with an empty scope\n    } catch {\n      builtin = false; // a free variable does not\n    }\n    builtinCache.set(name, builtin);\n    return builtin;\n  };\n\n  // A MathTS *function* (`gamma`, `distance`, `zeta`, `sin`, …) also resolves\n  // with an empty scope: it returns the function. That is not a value. A bare\n  // symbol with that name is a quantity (the adiabatic index, a length), and\n  // treating it as a built-in dropped it from the free-variable list, so\n  // `sqrt(gamma*pressure/density)` died as \"undeclared symbol 'gamma'\".\n  // Only a name that evaluates to a number is a MathTS constant (`pi`, `tau`).\n  // Bare `e` is one of those under MathTS's default (Euler's number). Evaluation\n  // below reads it as the elementary charge through `{ physics: true, charge: 'scalar' }`.\n  // An explicit scope `e` still wins. The name `euler` is refused below.\n  // A call such as `gamma(5)` is still a callee, decided separately.\n  const valueConstantCache = new Map<string, boolean>();\n  const isValueConstant = (name: string): boolean => {\n    const cached = valueConstantCache.get(name);\n    if (cached !== undefined) return cached;\n    let constant = false;\n    try {\n      constant = typeof mod.parse(name).evaluate({}) === 'number';\n    } catch {\n      constant = false;\n    }\n    valueConstantCache.set(name, constant);\n    return constant;\n  };\n\n  // Functions the built-in parser documents but MathTS does not define (`ln`),\n  // supplied through the scope so a formula means the same under either\n  // parser. Anything outside the documented list still fails.\n  const shims: Record<string, (...args: number[]) => number> = {};\n  for (const name of BUILTIN_FUNCTION_NAMES) {\n    if (!isBuiltin(name)) shims[name] = (...args) => callBuiltinFunction(name, args);\n  }\n\n  return {\n    parse(expr: string): CompiledFormula {\n      if (!expr || !expr.trim()) throw new FormulaError('empty formula');\n      let node: MathNode;\n      try {\n        node = mod.parse(expr);\n      } catch (err) {\n        throw new FormulaError(\n          `parse error: ${err instanceof Error ? err.message : String(err)}`,\n        );\n      }\n\n      // Free variables = symbol names − function callees − numeric constants.\n      // A function name is excluded only where it is called (`gamma(5)`), not\n      // where it is a quantity (`gamma` the adiabatic index).\n      const callees = new Set(\n        node\n          .filter((n) => n.isFunctionNode === true)\n          .map((n) => n.fn?.name)\n          .filter((n): n is string => typeof n === 'string'),\n      );\n      // A callee MathTS does not resolve and no shim supplies: evaluation would fail with MathTS's\n      // own \"Undefined function\", so it fails with the shared diagnostic instead (audit I4).\n      const symbolNames = node\n        .filter((n) => n.isSymbolNode === true)\n        .map((n) => n.name)\n        .filter((n): n is string => typeof n === 'string');\n      if (symbolNames.includes('euler') || callees.has('euler')) {\n        throw new EulerNumberError();\n      }\n      const unknownCallee = [...callees].find((n) => !(n in shims) && !isBuiltin(n));\n      const variables = [\n        ...new Set(\n          node\n            .filter((n) => n.isSymbolNode === true)\n            .map((n) => n.name)\n            .filter((n): n is string => typeof n === 'string')\n            .filter((n) => !callees.has(n) && !isValueConstant(n)),\n        ),\n      ].sort();\n\n      // `compileExpr` applies `{ physics: true, charge: 'scalar' }` on each\n      // evaluation, so bare `e` is the SI magnitude and a scope entry for `e` wins.\n      const compiled = compileExpr(expr, PHYSICS_SCALAR);\n      return {\n        source: expr,\n        variables,\n        evaluate(scope: Record<string, number>): number {\n          if (unknownCallee !== undefined) throw new FormulaError(unknownFunctionMessage(unknownCallee));\n          let result: unknown;\n          try {\n            result = compiled.evaluate({ ...shims, ...scope });\n          } catch (err) {\n            // `compileExpr` quotes the name (`Undefined symbol \"process\"`).\n            // The interpreter path reports `Undefined symbol process`.\n            // The one place src/ reads MathTS message text: MathTS throws a plain `Error` with no\n            // symbol field, and a pre-check against the scope cannot stand in because MathTS resolves\n            // some names itself (`hbar` reads as a unit). It only normalizes the text; nothing branches on it.\n            // Tracked in ACTIVE.md as a MathTS follow-up (a typed undefined-symbol error).\n            const message = err instanceof Error ? err.message : String(err);\n            throw new FormulaError(message.replace(/^Undefined symbol \"([^\"]+)\"$/, 'Undefined symbol $1'));\n          }\n          if (typeof result !== 'number' || !Number.isFinite(result)) {\n            // typeof Infinity and NaN is \"number\", which reads as a type error.\n            // A Complex is an object; \"got object\" hides that sqrt(-1) left the reals.\n            throw new FormulaValueError(`formula did not evaluate to a finite number (got ${describeNonFinite(result)})`);\n          }\n          return result;\n        },\n      };\n    },\n  };\n}\n\n/** What a non-finite or non-number result was, for the error text. */\nfunction describeNonFinite(result: unknown): string {\n  if (typeof result === 'number') return String(result);\n  if (result !== null && typeof result === 'object') {\n    const c = result as { re?: unknown; im?: unknown; toString?: () => string };\n    if (typeof c.re === 'number' && typeof c.im === 'number') {\n      const text = typeof c.toString === 'function' ? c.toString() : `${c.re}+${c.im}i`;\n      return `a complex number (${text})`;\n    }\n  }\n  return typeof result;\n}\n\nconst mathTsModule: MathtsFunctionsModule = {\n  // MathTS's Node type is not imported; `MathNode` is the structural subset\n  // this parser reads (`evaluate`, `filter`, `isSymbolNode`, `isFunctionNode`,\n  // `name`, `fn`), every member of which exists on a MathTS node.\n  parse: (expr) => parseMathTs(expr) as unknown as MathNode,\n};\n\n/** The MathTS formula parser. @internal */\nexport const mathtsFormulaParser: FormulaParser = createMathtsFormulaParser(mathTsModule);\n\n/** Parse a scalar formula. @internal */\nexport function parseFormula(expr: string): CompiledFormula {\n  return mathtsFormulaParser.parse(expr);\n}\n"]}

package/dist/numerical/formula-registry.d.ts

@@ -10,10 +10,20 @@
 import type { FormulaParser } from './formula-contract.js';
 import type { Dimension } from '../dimensional/types.js';
 import type { FormulaDimensionChecker, ParsedPhysics } from './formula-dimension.js';
-type FormulaParserKind = 'mathts' | 'builtin';
-/** Resolve the active formula parser. It is the MathTS parser. @internal */
+/**
+ * The one formula parser kind. The `builtin` recursive-descent parser was
+ * removed with Path B, so the type has no second member; the CLI records
+ * this value in an experiment record's environment.
+ */
+type FormulaParserKind = 'mathts';
+/**
+ * Resolve the active formula parser. It is the MathTS parser. The function
+ * stays async because the CLI awaits it on a code path that once loaded an
+ * optional peer; the value is constant.
+ * @internal
+ */
 export declare function getFormulaParser(): Promise<FormulaParser>;
-/** Which parser is active. Always `mathts`. @internal */
+/** Which parser is active: `mathts`, the one value of {@link FormulaParserKind}. @internal */
 export declare function getFormulaParserKind(): Promise<FormulaParserKind>;
 /**
  * The dimensional checker for user formulas. It transpiles the MathTS AST.

package/dist/numerical/formula-registry.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"formula-registry.d.ts","sourceRoot":"","sources":["../../src/numerical/formula-registry.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AAE3D,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,KAAK,EAAE,uBAAuB,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AAGrF,KAAK,iBAAiB,GAAG,QAAQ,GAAG,SAAS,CAAC;AAE9C,4EAA4E;AAC5E,wBAAsB,gBAAgB,IAAI,OAAO,CAAC,aAAa,CAAC,CAE/D;AAED,yDAAyD;AACzD,wBAAsB,oBAAoB,IAAI,OAAO,CAAC,iBAAiB,CAAC,CAEvE;AAED;;;GAGG;AACH,wBAAsB,0BAA0B,IAAI,OAAO,CAAC,uBAAuB,CAAC,CAEnF;AAED;;;;;;GAMG;AACH,wBAAsB,YAAY,CAChC,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC,GACxC,OAAO,CAAC,aAAa,CAAC,CAExB"}
+{"version":3,"file":"formula-registry.d.ts","sourceRoot":"","sources":["../../src/numerical/formula-registry.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AAE3D,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,KAAK,EAAE,uBAAuB,EAAE,aAAa,EAAE,MAAM,wBAAwB,CAAC;AAGrF;;;;GAIG;AACH,KAAK,iBAAiB,GAAG,QAAQ,CAAC;AAElC;;;;;GAKG;AACH,wBAAsB,gBAAgB,IAAI,OAAO,CAAC,aAAa,CAAC,CAE/D;AAED,8FAA8F;AAC9F,wBAAsB,oBAAoB,IAAI,OAAO,CAAC,iBAAiB,CAAC,CAEvE;AAED;;;GAGG;AACH,wBAAsB,0BAA0B,IAAI,OAAO,CAAC,uBAAuB,CAAC,CAEnF;AAED;;;;;;GAMG;AACH,wBAAsB,YAAY,CAChC,IAAI,EAAE,MAAM,EACZ,IAAI,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC,GACxC,OAAO,CAAC,aAAa,CAAC,CAExB"}

package/dist/numerical/formula-registry.js

@@ -9,11 +9,16 @@
  */
 import { mathtsFormulaParser } from './formula-mathts.js';
 import { builtinFormulaDimensionChecker } from './formula-dimension.js';
-/** Resolve the active formula parser. It is the MathTS parser. @internal */
+/**
+ * Resolve the active formula parser. It is the MathTS parser. The function
+ * stays async because the CLI awaits it on a code path that once loaded an
+ * optional peer; the value is constant.
+ * @internal
+ */
 export async function getFormulaParser() {
     return mathtsFormulaParser;
 }
-/** Which parser is active. Always `mathts`. @internal */
+/** Which parser is active: `mathts`, the one value of {@link FormulaParserKind}. @internal */
 export async function getFormulaParserKind() {
     return 'mathts';
 }

package/dist/numerical/formula-registry.js.map

@@ -1 +1 @@
-{"version":3,"file":"formula-registry.js","sourceRoot":"","sources":["../../src/numerical/formula-registry.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAGH,OAAO,EAAE,mBAAmB,EAAE,MAAM,qBAAqB,CAAC;AAG1D,OAAO,EAAE,8BAA8B,EAAE,MAAM,wBAAwB,CAAC;AAIxE,4EAA4E;AAC5E,MAAM,CAAC,KAAK,UAAU,gBAAgB;IACpC,OAAO,mBAAmB,CAAC;AAC7B,CAAC;AAED,yDAAyD;AACzD,MAAM,CAAC,KAAK,UAAU,oBAAoB;IACxC,OAAO,QAAQ,CAAC;AAClB,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,0BAA0B;IAC9C,OAAO,8BAA8B,EAAE,CAAC;AAC1C,CAAC;AAED;;;;;;GAMG;AACH,MAAM,CAAC,KAAK,UAAU,YAAY,CAChC,IAAY,EACZ,IAAyC;IAEzC,OAAO,8BAA8B,EAAE,CAAC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;AAC5D,CAAC"}
+{"version":3,"file":"formula-registry.js","sourceRoot":"","sources":["../../src/numerical/formula-registry.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAGH,OAAO,EAAE,mBAAmB,EAAE,MAAM,qBAAqB,CAAC;AAG1D,OAAO,EAAE,8BAA8B,EAAE,MAAM,wBAAwB,CAAC;AASxE;;;;;GAKG;AACH,MAAM,CAAC,KAAK,UAAU,gBAAgB;IACpC,OAAO,mBAAmB,CAAC;AAC7B,CAAC;AAED,8FAA8F;AAC9F,MAAM,CAAC,KAAK,UAAU,oBAAoB;IACxC,OAAO,QAAQ,CAAC;AAClB,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,0BAA0B;IAC9C,OAAO,8BAA8B,EAAE,CAAC;AAC1C,CAAC;AAED;;;;;;GAMG;AACH,MAAM,CAAC,KAAK,UAAU,YAAY,CAChC,IAAY,EACZ,IAAyC;IAEzC,OAAO,8BAA8B,EAAE,CAAC,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;AAC5D,CAAC","sourcesContent":["/**\n * Formula-parser registry.\n *\n * `getFormulaParser()` returns the MathTS-backed parser. The MathTS packages\n * are required dependencies, so there is no absent-peer branch and no Path B\n * parser.\n *\n * @module numerical/formula-registry\n */\n\nimport type { FormulaParser } from './formula-contract.js';\nimport { mathtsFormulaParser } from './formula-mathts.js';\nimport type { Dimension } from '../dimensional/types.js';\nimport type { FormulaDimensionChecker, ParsedPhysics } from './formula-dimension.js';\nimport { builtinFormulaDimensionChecker } from './formula-dimension.js';\n\n/**\n * The one formula parser kind. The `builtin` recursive-descent parser was\n * removed with Path B, so the type has no second member; the CLI records\n * this value in an experiment record's environment.\n */\ntype FormulaParserKind = 'mathts';\n\n/**\n * Resolve the active formula parser. It is the MathTS parser. The function\n * stays async because the CLI awaits it on a code path that once loaded an\n * optional peer; the value is constant.\n * @internal\n */\nexport async function getFormulaParser(): Promise<FormulaParser> {\n  return mathtsFormulaParser;\n}\n\n/** Which parser is active: `mathts`, the one value of {@link FormulaParserKind}. @internal */\nexport async function getFormulaParserKind(): Promise<FormulaParserKind> {\n  return 'mathts';\n}\n\n/**\n * The dimensional checker for user formulas. It transpiles the MathTS AST.\n * @internal\n */\nexport async function getFormulaDimensionChecker(): Promise<FormulaDimensionChecker> {\n  return builtinFormulaDimensionChecker();\n}\n\n/**\n * Parse physics text to a dimensional `ExprNode` + its inferred dimension.\n * The single public string→`ExprNode` entry point.\n *\n * @throws {FormulaDimensionError} on a parse error or non-homogeneity.\n * @public\n */\nexport async function parsePhysics(\n  text: string,\n  dims: Readonly<Record<string, Dimension>>,\n): Promise<ParsedPhysics> {\n  return builtinFormulaDimensionChecker().parse(text, dims);\n}\n"]}

package/dist/numerical/geodesic-integrator.js.map

@@ -1 +1 @@
-{"version":3,"file":"geodesic-integrator.js","sourceRoot":"","sources":["../../src/numerical/geodesic-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,OAAO,EAAE,cAAc,EAAE,MAAM,iCAAiC,CAAC;AACjE,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AA+EpD;;;;;;GAMG;AACH,SAAS,WAAW,CAClB,aAA6E,EAC7E,CAAO,EACP,CAAO;AACP,yEAAyE;AACzE,4EAA4E;AAC5E,8DAA8D;AAC9D,OAAsB;IAEtB,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;IAEvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;QAC9B,6EAA6E;QAC7E,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;QACrB,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,OAAO,GAAG,CAAC,CAAC;QAChB,IAAI,CAAC,GAAG,CAAC,CAAC;QAEV,6EAA6E;QAC7E,sFAAsF;QAEtF,SAAS;QACT,IAAI,MAAM,GAAG,IAAI,CAAC;QAClB,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;aACpB,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;aACzB,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;;YACzB,GAAG,GAAG,CAAC,GAAG,CAAC;IAClB,CAAC;IAED,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC;AAC5D,CAAC;AAED,8EAA8E;AAC9E,aAAa;AACb,8EAA8E;AAE9E;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAC/B,MAAgC;IAEhC,MAAM,EAAE,aAAa,EAAE,EAAE,EAAE,EAAE,EAAE,QAAQ,EAAE,MAAM,EAAE,KAAK,EAAE,eAAe,EAAE,GAAG,MAAM,CAAC;IACnF,mFAAmF;IAEnF,2EAA2E;IAC3E,yEAAyE;IACzE,IAAI,OAAO,eAAe,KAAK,QAAQ,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACnE,MAAM,IAAI,qBAAqB,CAC7B,wCAAwC,EAAE,CAAC,CAAC,CAAC,sBAAsB,eAAe,IAAI;YACtF,oDAAoD;YACpD,2EAA2E,CAC5E,CAAC;IACJ,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QAC3C,MAAM,IAAI,qBAAqB,CAC7B,4DAA4D,KAAK,EAAE,CACpE,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,KAAK,CAAC;IACtC,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAEpE,6EAA6E;IAC7E,6EAA6E;IAC7E,gEAAgE;IAChE,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,CAAC;IACrC,MAAM,GAAG,GAAG,cAAc,CACxB,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;QACV,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,WAAW,CAAC,aAAa,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;QACpD,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAClF,CAAC,EACD,EAAE,EACF,CAAC,QAAQ,EAAE,MAAM,CAAC,EAClB,EAAE,EAAE,EAAE,CAAC,EAAE,CACV,CAAC;IAEF,MAAM,MAAM,GAAG,CAAC,GAAsB,EAAQ,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAClF,MAAM,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC;IACzD,MAAM,UAAU,GAAW,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACpD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,WAAW,KAAK,CAAC,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;YAC1D,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/B,CAAC;IACH,CAAC;IACD,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;IAE3C,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,IAAI,CAAC;QACpB,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAC5C,UAAU;KACX,CAAC;AACJ,CAAC"}
+{"version":3,"file":"geodesic-integrator.js","sourceRoot":"","sources":["../../src/numerical/geodesic-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,OAAO,EAAE,cAAc,EAAE,MAAM,iCAAiC,CAAC;AACjE,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AA+EpD;;;;;;GAMG;AACH,SAAS,WAAW,CAClB,aAA6E,EAC7E,CAAO,EACP,CAAO;AACP,yEAAyE;AACzE,4EAA4E;AAC5E,8DAA8D;AAC9D,OAAsB;IAEtB,MAAM,CAAC,GAAG,aAAa,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;IACpC,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACjD,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;IAEvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;QAC9B,6EAA6E;QAC7E,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;QACrB,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,OAAO,GAAG,CAAC,CAAC;QAChB,IAAI,CAAC,GAAG,CAAC,CAAC;QAEV,6EAA6E;QAC7E,sFAAsF;QAEtF,SAAS;QACT,IAAI,MAAM,GAAG,IAAI,CAAC;QAClB,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,SAAS;QACT,MAAM,IAAI,CAAC,CAAC;QACZ,OAAO,GAAG,CAAC,CAAC;QACZ,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAC9C,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,CAAC,GAAG,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;QAAC,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,IAAI,CAAC,GAAG,EAAE,CAAC;QAChD,IAAI,OAAO,KAAK,CAAC;YAAE,GAAG,IAAI,OAAO,GAAG,EAAE,CAAC;QAEvC,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;aACpB,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;aACzB,IAAI,EAAE,KAAK,CAAC;YAAE,GAAG,GAAG,CAAC,GAAG,CAAC;;YACzB,GAAG,GAAG,CAAC,GAAG,CAAC;IAClB,CAAC;IAED,OAAO,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC;AAC5D,CAAC;AAED,8EAA8E;AAC9E,aAAa;AACb,8EAA8E;AAE9E;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAC/B,MAAgC;IAEhC,MAAM,EAAE,aAAa,EAAE,EAAE,EAAE,EAAE,EAAE,QAAQ,EAAE,MAAM,EAAE,KAAK,EAAE,eAAe,EAAE,GAAG,MAAM,CAAC;IACnF,mFAAmF;IAEnF,2EAA2E;IAC3E,yEAAyE;IACzE,IAAI,OAAO,eAAe,KAAK,QAAQ,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACnE,MAAM,IAAI,qBAAqB,CAC7B,wCAAwC,EAAE,CAAC,CAAC,CAAC,sBAAsB,eAAe,IAAI;YACtF,oDAAoD;YACpD,2EAA2E,CAC5E,CAAC;IACJ,CAAC;IAED,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QAC3C,MAAM,IAAI,qBAAqB,CAC7B,4DAA4D,KAAK,EAAE,CACpE,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,GAAG,CAAC,MAAM,GAAG,QAAQ,CAAC,GAAG,KAAK,CAAC;IACtC,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAEpE,6EAA6E;IAC7E,6EAA6E;IAC7E,gEAAgE;IAChE,MAAM,OAAO,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,CAAC;IACrC,MAAM,GAAG,GAAG,cAAc,CACxB,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE;QACV,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAS,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACzC,MAAM,CAAC,GAAG,WAAW,CAAC,aAAa,EAAE,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;QACpD,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAClF,CAAC,EACD,EAAE,EACF,CAAC,QAAQ,EAAE,MAAM,CAAC,EAClB,EAAE,EAAE,EAAE,CAAC,EAAE,CACV,CAAC;IAEF,MAAM,MAAM,GAAG,CAAC,GAAsB,EAAQ,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAClF,MAAM,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,CAAC;IACzD,MAAM,UAAU,GAAW,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACpD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,WAAW,KAAK,CAAC,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;YAC1D,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/B,CAAC;IACH,CAAC;IACD,MAAM,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;IAE3C,OAAO;QACL,MAAM,EAAE,MAAM,CAAC,IAAI,CAAC;QACpB,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAC5C,UAAU;KACX,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Fixed-step RK4 integrator for the geodesic equation in an arbitrary\n * Lorentzian manifold (v0.4.0, Task 14 [U]).\n *\n * Solves the second-order ODE system\n *\n *   d²x^μ/dτ² + Γ^μ_{νρ}(x) (dx^ν/dτ)(dx^ρ/dτ) = 0\n *\n * by casting it as a first-order system on the (x, v) phase space:\n *\n *   dx^μ/dτ  = v^μ\n *   dv^μ/dτ  = −Γ^μ_{νρ}(x) v^ν v^ρ    (Einstein summation over ν, ρ)\n *\n * The Christoffel symbol Γ^μ_{νρ} is supplied as a plain-JS closure that\n * maps a 4-coordinate array to a [4][4][4] tensor.  No TensorEngine\n * dependency — the integrator is self-contained. The fixed steps are\n * MathTS `solveODESystem` with `dt` set, which is that library's classical\n * RK4. Adaptive `solveODE` is not this step.\n *\n * Validated against the cycloid parametric form for radial Schwarzschild\n * infall to ±1e-6 (Task 14 test, η = 0.5).\n *\n * @module numerical/geodesic-integrator\n */\nimport { solveODESystem } from '@danielsimonjr/mathts-functions';\nimport { NumericalBackendError } from './errors.js';\n\n/**\n * Input bundle for `integrateGeodesic`.\n * @public\n */\nexport interface GeodesicIntegratorInputs {\n  /**\n   * Christoffel-symbol closure.  Maps a 4-coordinate array x^μ to a\n   * Float64Array(64) containing Γ^μ_{νρ}(x) in λ-major layout:\n   * index (λ, μ, ν) → 16·λ + 4·μ + ν.\n   *\n   * BR-2 (v0.6.0 Task 2.9): migrated from nested number[][][] to flat\n   * Float64Array(64). Consumers call `arr[16*mu + 4*nu + rho]`.\n   *\n   * The optional `out` buffer lets the integrator pass a reused scratch\n   * Float64Array(64) so the hot loop avoids per-call allocation. A function\n   * that ignores `out` (e.g. allocates fresh each call) still satisfies this\n   * type — the optimization is opt-in.\n   */\n  christoffelFn: (x: ReadonlyArray<number>, out?: Float64Array) => Float64Array;\n\n  /** Initial 4-position x^μ(τ₀) = [t, r, θ, φ]. */\n  x0: readonly [number, number, number, number];\n\n  /**\n   * Initial 4-velocity dx^μ/dτ at τ₀.\n   *\n   * I2 fix: for a Schwarzschild radial timelike geodesic starting at\n   * coordinate radius r₀, use\n   *   v0 = [1/√(1 − r_s/r₀), 0, 0, 0]\n   * so that g_{μν} u^μ u^ν = −c² holds.  At r₀ = 100·r_s this gives\n   * dt/dτ ≈ 1.005, not 1.\n   */\n  v0: readonly [number, number, number, number];\n\n  /** Proper-time parameter at the start of integration. */\n  tauStart: number;\n\n  /** Proper-time parameter at the end of integration. */\n  tauEnd: number;\n\n  /** Number of fixed RK4 steps.  Must be a positive integer. */\n  steps: number;\n\n  /**\n   * E11 fix: explicit domain-validity option in place of the\n   * `(christoffelFn as any).r_s` monkey-patch.  If supplied, the\n   * integrator throws `NumericalBackendError` when `x0[1] < domainMinRadius`\n   * (e.g. for Schwarzschild: caller passes `3 * r_s` to enforce the\n   * v0.4.0 domain restriction r ≥ 3·r_Schwarz).\n   */\n  domainMinRadius?: number;\n}\n\n/**\n * Result of `integrateGeodesic`.\n * @public\n */\nexport interface GeodesicIntegratorResult {\n  /** Final 4-position x^μ(τEnd). */\n  readonly xFinal: readonly [number, number, number, number];\n\n  /** Final 4-velocity v^μ(τEnd). */\n  readonly vFinal: readonly [number, number, number, number];\n\n  /**\n   * Sampled trajectory positions (approximately 100 evenly-spaced snapshots\n   * plus the initial point).  Useful for visualisation and diagnostics.\n   */\n  readonly trajectory: ReadonlyArray<readonly [number, number, number, number]>;\n}\n\n// ---------------------------------------------------------------------------\n// Internal helpers\n// ---------------------------------------------------------------------------\n\ntype Vec4 = [number, number, number, number];\n\n/**\n * Evaluate the Christoffel acceleration −Γ^μ_{νρ} v^ν v^ρ at position x\n * with 4-velocity v.  Returns (dx/dτ, dv/dτ) = (v, accel).\n *\n * BR-2 (v0.6.0 Task 2.9): christoffelFn now returns Float64Array(64).\n * Index access: G[16*mu + 4*nu + rho] (λ-major layout).\n */\nfunction geodesicRHS(\n  christoffelFn: (x: ReadonlyArray<number>, out?: Float64Array) => Float64Array,\n  x: Vec4,\n  v: Vec4,\n  // Reused scratch buffer for the Christoffel evaluation. Safe: G is fully\n  // consumed into `dv` below before geodesicRHS returns, so the next call may\n  // overwrite it. `dx`/`dv` are fresh arrays and never alias G.\n  scratch?: Float64Array,\n): { dx: Vec4; dv: Vec4 } {\n  const G = christoffelFn(x, scratch);\n  const v0 = v[0], v1 = v[1], v2 = v[2], v3 = v[3];\n  let dv0 = 0, dv1 = 0, dv2 = 0, dv3 = 0;\n\n  for (let mu = 0; mu < 4; mu++) {\n    // Bolt: Hoist the invariant 16 * mu multiplication outside the nu, rho loops\n    const mu16 = mu * 16;\n    let acc = 0;\n    let accTerm = 0;\n    let g = 0;\n\n    // Bolt: Fully unroll the inner `nu` and `rho` loops for a fixed 4D spacetime\n    // to eliminate loop overhead and index arithmetic, caching `v` into scalar variables.\n\n    // nu = 0\n    let offset = mu16;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v0;\n\n    // nu = 1\n    offset += 4;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v1;\n\n    // nu = 2\n    offset += 4;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v2;\n\n    // nu = 3\n    offset += 4;\n    accTerm = 0;\n    g = G[offset]; if (g !== 0) accTerm += g * v0;\n    g = G[offset+1]; if (g !== 0) accTerm += g * v1;\n    g = G[offset+2]; if (g !== 0) accTerm += g * v2;\n    g = G[offset+3]; if (g !== 0) accTerm += g * v3;\n    if (accTerm !== 0) acc += accTerm * v3;\n\n    if (mu === 0) dv0 = -acc;\n    else if (mu === 1) dv1 = -acc;\n    else if (mu === 2) dv2 = -acc;\n    else dv3 = -acc;\n  }\n\n  return { dx: [v0, v1, v2, v3], dv: [dv0, dv1, dv2, dv3] };\n}\n\n// ---------------------------------------------------------------------------\n// Public API\n// ---------------------------------------------------------------------------\n\n/**\n * Numerically integrate the geodesic equation for a timelike (or null)\n * worldline in an arbitrary curved spacetime using a fixed-step classical\n * fourth-order Runge-Kutta scheme.\n *\n * @public\n */\nexport function integrateGeodesic(\n  inputs: GeodesicIntegratorInputs,\n): GeodesicIntegratorResult {\n  const { christoffelFn, x0, v0, tauStart, tauEnd, steps, domainMinRadius } = inputs;\n  // BR-2: christoffelFn now returns Float64Array(64); geodesicRHS reads flat layout.\n\n  // E11 fix: explicit option in place of monkey-patched .r_s on the closure.\n  // Callers (e.g. the Schwarzschild test fixture) pass `3 * r_s` directly.\n  if (typeof domainMinRadius === 'number' && x0[1] < domainMinRadius) {\n    throw new NumericalBackendError(\n      `integrateGeodesic: initial radius r0=${x0[1]} < domainMinRadius=${domainMinRadius}. ` +\n      `Position is outside the valid integration domain. ` +\n      `For Schwarzschild geodesics this restriction is r ≥ 3·r_Schwarz (v0.4.0).`,\n    );\n  }\n\n  if (!Number.isInteger(steps) || steps <= 0) {\n    throw new NumericalBackendError(\n      `integrateGeodesic: steps must be a positive integer, got ${steps}`,\n    );\n  }\n\n  const h = (tauEnd - tauStart) / steps;\n  const y0 = [x0[0], x0[1], x0[2], x0[3], v0[0], v0[1], v0[2], v0[3]];\n\n  // One scratch buffer reused across all Christoffel evaluations (4 per step).\n  // Each geodesicRHS consumes its Γ before the next call, so reuse is safe and\n  // avoids ~4·steps Float64Array(64) allocations per integration.\n  const scratch = new Float64Array(64);\n  const sol = solveODESystem(\n    (_tau, y) => {\n      const x: Vec4 = [y[0], y[1], y[2], y[3]];\n      const v: Vec4 = [y[4], y[5], y[6], y[7]];\n      const s = geodesicRHS(christoffelFn, x, v, scratch);\n      return [s.dx[0], s.dx[1], s.dx[2], s.dx[3], s.dv[0], s.dv[1], s.dv[2], s.dv[3]];\n    },\n    y0,\n    [tauStart, tauEnd],\n    { dt: h },\n  );\n\n  const asVec4 = (row: readonly number[]): Vec4 => [row[0], row[1], row[2], row[3]];\n  const sampleEvery = Math.max(1, Math.floor(steps / 100));\n  const trajectory: Vec4[] = [asVec4(sol.y[0] ?? y0)];\n  for (let i = 0; i < steps; i++) {\n    if ((i + 1) % sampleEvery === 0) {\n      const row = sol.y[i + 1] ?? sol.y[sol.y.length - 1] ?? y0;\n      trajectory.push(asVec4(row));\n    }\n  }\n  const last = sol.y[sol.y.length - 1] ?? y0;\n\n  return {\n    xFinal: asVec4(last),\n    vFinal: [last[4], last[5], last[6], last[7]],\n    trajectory,\n  };\n}\n"]}

package/dist/numerical/geometrized.js.map

@@ -1 +1 @@
-{"version":3,"file":"geometrized.js","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAGH,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AACpD,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,sBAAsB,CAAC;AAElD;;;;;GAKG;AACH,MAAM,OAAO,8BAA+B,SAAQ,QAAQ;IAC1D,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,gCAAgC,CAAC;QAC7C,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,8BAA8B,CAAC,SAAS,CAAC,CAAC;IACxE,CAAC;CACF;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,iBAAiB,CAAC,GAAc;IAC9C,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;QACjE,MAAM,IAAI,8BAA8B,CACtC,qEAAqE;YACnE,wCAAwC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,KAAK,IAAI;YACjE,KAAK,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,CAAC,iDAAiD,CAC1E,CAAC;IACJ,CAAC;IACD,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;AACnE,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,aAAa,CAAC,OAAe,EAAE,GAAc;IAC3D,OAAO,OAAO,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,eAAe,CAAC,SAAiB,EAAE,GAAc;IAC/D,OAAO,SAAS,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC5C,CAAC"}
+{"version":3,"file":"geometrized.js","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAGH,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AACpD,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,sBAAsB,CAAC;AAElD;;;;;GAKG;AACH,MAAM,OAAO,8BAA+B,SAAQ,QAAQ;IAC1D,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,gCAAgC,CAAC;QAC7C,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,8BAA8B,CAAC,SAAS,CAAC,CAAC;IACxE,CAAC;CACF;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,iBAAiB,CAAC,GAAc;IAC9C,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;QACjE,MAAM,IAAI,8BAA8B,CACtC,qEAAqE;YACnE,wCAAwC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,KAAK,IAAI;YACjE,KAAK,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,CAAC,iDAAiD,CAC1E,CAAC;IACJ,CAAC;IACD,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;AACnE,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,aAAa,CAAC,OAAe,EAAE,GAAc;IAC3D,OAAO,OAAO,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,eAAe,CAAC,SAAiB,EAAE,GAAc;IAC/D,OAAO,SAAS,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC5C,CAAC","sourcesContent":["/**\n * Geometrized-units boundary adapters (G-9 increment 1).\n *\n * Converts a scalar between SI and geometrized (G = c = 1) units, driven\n * MECHANICALLY by the `Dimension` exponent vector — the dimension functor\n * earning its keep. For an SI quantity of dimension L^l M^m T^t the single\n * conversion factor is `G^m · c^(t − 2m)` (each kg → G/c² metres, each second\n * → c metres); the geometrized quantity is a pure length L^(l+m+t).\n *\n * This is the self-contained FOUNDATION of the units-normalization layer\n * (design: docs/planning/v0.10.0-Units-Normalization-Design-Note.md, Adam-\n * vetted r2; plan: docs/planning/v0.13-G9-Adapters-Plan.md, Eve-vetted). Only c\n * and G are in scope: a nonzero electromagnetic / thermal / molar / luminous\n * exponent (I/Θ/N/J) throws.\n *\n * PUBLIC boundary API (G-9 increment 2): a consumer running the GR pipeline in\n * geometrized (G = c = 1) units converts scalars at the boundary with these\n * adapters. The consumer-wide DEFAULT-pipeline migration (increment 3) was\n * DECLINED after adversarial review — it is measured to give no precision win\n * (and far worse far-field; see docs/planning/v0.14-G9-Increment3-Disposition.md).\n * The SI pipeline stays the default; this geometrized layer is the optional\n * boundary API.\n *\n * @module numerical/geometrized\n */\n\nimport type { Dimension } from '../dimensional/types.js';\nimport { UPTError } from '../dimensional/errors.js';\nimport { C_SI, G_SI } from '../core/constants.js';\n\n/**\n * A dimension carries a non-c/G base (I/Θ/N/J ≠ 0) and cannot be geometrized\n * by c and G alone (k_B / ε₀ extensions are out of scope for this increment).\n *\n * @public\n */\nexport class NonGeometrizableDimensionError extends UPTError {\n  constructor(message: string) {\n    super(message);\n    this.name = 'NonGeometrizableDimensionError';\n    Object.setPrototypeOf(this, NonGeometrizableDimensionError.prototype);\n  }\n}\n\n/**\n * The single SI→geometrized conversion factor for `dim`: `G^M · c^(T − 2M)`\n * (only the mass and time exponents enter). Shared by `toGeometrized` /\n * `fromGeometrized` so the round-trip is a clean inverse. A nonzero I/Θ/N/J\n * exponent throws `NonGeometrizableDimensionError` BEFORE any factor is\n * computed.\n *\n * @public\n */\nexport function geometrizedFactor(dim: Dimension): number {\n  if (dim.I !== 0 || dim.Theta !== 0 || dim.N !== 0 || dim.J !== 0) {\n    throw new NonGeometrizableDimensionError(\n      `geometrizedFactor: only c and G are in scope; this dimension has a ` +\n        `nonzero non-mechanizable exponent (I=${dim.I}, Θ=${dim.Theta}, ` +\n        `N=${dim.N}, J=${dim.J}). Geometrization handles length / mass / time.`,\n    );\n  }\n  return Math.pow(G_SI, dim.M) * Math.pow(C_SI, dim.T - 2 * dim.M);\n}\n\n/** SI → geometrized (G = c = 1). @public */\nexport function toGeometrized(valueSI: number, dim: Dimension): number {\n  return valueSI * geometrizedFactor(dim);\n}\n\n/** Geometrized (G = c = 1) → SI. @public */\nexport function fromGeometrized(valueGeom: number, dim: Dimension): number {\n  return valueGeom / geometrizedFactor(dim);\n}\n"]}

package/dist/numerical/gl4-integrator.d.ts

@@ -26,8 +26,9 @@ export interface GL4State {
 }
 /**
  * Per-step snapshot recorded by the integrator. `v` (the contravariant
- * 4-velocity v^μ = g^{μν} p_ν) is optional — emitted when the caller asks
- * for it, since it requires an extra metric-inverse contraction.
+ * 4-velocity v^μ = g^{μν} p_ν) is declared optional and is never emitted:
+ * no option asks for it, and nothing in `src/` reads it. It stays so a caller's
+ * own snapshot type can carry one.
  *
  * @public
  */
@@ -90,34 +91,27 @@ export interface GL4Options {
     }) => void;
 }
 /**
- * Result of `solveGL4Stage` — the two converged stage values plus the
- * iteration count actually consumed. Consumed by the
- * `integrateGeodesicGL4` step driver.
- *
- * v0.6.1: dropped export — internal-only result shape (was already
- * @internal-tagged but had no external consumer).
+ * The state after one GL4 step, with the Picard iteration count that step
+ * consumed. MathTS `gaussLegendre4` does not expose its two internal stage
+ * values, so this is the advanced state and nothing else.
+ * @internal
  */
-interface StageSolveResult {
-    readonly stageX: readonly [readonly number[], readonly number[]];
-    readonly stageP: readonly [readonly number[], readonly number[]];
-    readonly stageDx: readonly [readonly number[], readonly number[]];
-    readonly stageDp: readonly [readonly number[], readonly number[]];
+export interface GL4StepResult {
+    readonly x: readonly number[];
+    readonly p: readonly number[];
     readonly iterations: number;
 }
 /**
- * One GL4 step, reported in the stage-result shape the stage tests read.
- *
- * The step itself is MathTS `gaussLegendre4`. In flat space `p` is constant,
- * so both stage momenta equal that advanced `p`. Throws `GL4ConvergenceError`
- * with a message matching `/Picard iteration did not converge/` when the
- * Picard cap is exhausted.
- *
+ * One GL4 step of the geodesic Hamiltonian, via MathTS `gaussLegendre4`.
+ * Throws `GL4ConvergenceError` with a message matching
+ * `/Picard iteration did not converge/` when the Picard cap is exhausted;
+ * any other error from the metric closures propagates unchanged.
  * @internal
  */
-export declare function solveGL4Stage(state: GL4State, h: number, gInverseFn: (x: readonly number[]) => Float64Array, dgInverseFn: (x: readonly number[]) => Float64Array, opts: {
+export declare function gl4Step(state: GL4State, h: number, gInverseFn: (x: readonly number[]) => Float64Array, dgInverseFn: (x: readonly number[]) => Float64Array, opts: {
     picardTol: number;
     picardMaxIter: number;
-}): StageSolveResult;
+}): GL4StepResult;
 /**
  * GL4 symplectic integrator on the canonical (x, p) geodesic Hamiltonian.
  *
@@ -163,10 +157,11 @@ export declare function solveGL4Stage(state: GL4State, h: number, gInverseFn: (x
  * @param initialState — canonical (x, p) at τ = 0. Units follow the
  *   `gInverseFn` convention (see above).
  * @param options — see {@link GL4Options}:
- *   - `steps` — integer step count (dimensionless).
- *   - `tauMax` — affine-parameter extent (seconds in canonical SI).
- *   - `gInverseFn(x)[μ][ν]` — inverse metric g^{μν}(x).
- *   - `dgInverseFn(x)[λ][μ][ν]` — ∂_λ g^{μν}(x).
+ *   - `steps` — positive integer step count (dimensionless).
+ *   - `tauMax` — positive finite affine-parameter extent (seconds in
+ *     canonical SI). The integrator runs forward from τ = 0 only.
+ *   - `gInverseFn(x)[μ*dim + ν]` — inverse metric g^{μν}(x), row-major flat.
+ *   - `dgInverseFn(x)[λ*dim² + μ*dim + ν]` — ∂_λ g^{μν}(x), flat.
  *   - `picardTol` — convergence tolerance (dimensionless, default 1e-12).
  *   - `picardMaxIter` — fixed-point iteration cap (dimensionless, default 50).
  *   - `hMin` — step-halving floor (same units as `tauMax / steps`).
@@ -174,8 +169,14 @@ export declare function solveGL4Stage(state: GL4State, h: number, gInverseFn: (x
  * @returns `steps + 1` snapshots: index `0` is the initial state, index `n`
  *   is the state after `n` steps (τ = n · h). Each snapshot carries `tau`,
  *   `x`, `p` (and optional `v` = g^{μν} p_ν) in the units chosen above.
- * @throws NumericalBackendError if `initialState.x[1] < domainMinRadius`.
+ * @throws NumericalBackendError if `steps` is not a positive integer, if
+ *   `tauMax` is not a positive finite number, or if
+ *   `initialState.x[1] < domainMinRadius`. The two option checks run before
+ *   any metric evaluation.
  * @throws GL4ConvergenceError if Picard fails even after step-halving to h_min.
+ *   Any other error thrown by `gInverseFn` or `dgInverseFn` (a domain
+ *   crossing the caller signals by throwing, or a programming error) is
+ *   rethrown unchanged; it is not halved away.
  *
  * @public
  */

package/dist/numerical/gl4-integrator.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"gl4-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AA4BA;;;;;;GAMG;AACH,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,iCAAiC,CAAC;AAEtE;;;;;;;;;;;GAWG;AACH,MAAM,WAAW,QAAQ;IACvB,+BAA+B;IAC/B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,yCAAyC;IACzC,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAC/B;AAED;;;;;;GAMG;AACH,MAAM,WAAW,WAAW;IAC1B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAChC;AAED;;;;GAIG;AACH,MAAM,WAAW,UAAU;IACzB;qCACiC;IACjC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,+CAA+C;IAC/C,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB;gFAC4E;IAC5E,QAAQ,CAAC,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC5D;;;;;OAKG;IACH,QAAQ,CAAC,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC7D,+DAA+D;IAC/D,QAAQ,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC;IAC5B,mEAAmE;IACnE,QAAQ,CAAC,aAAa,CAAC,EAAE,MAAM,CAAC;IAChC,+FAA+F;IAC/F,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;IACvB;;qCAEiC;IACjC,QAAQ,CAAC,eAAe,CAAC,EAAE,MAAM,CAAC;IAClC;;;;;;;;;;;;;;OAcG;IACH,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC,KAAK,EAAE;QAAE,IAAI,EAAE,MAAM,CAAC;QAAC,UAAU,EAAE,MAAM,CAAC;QAAC,QAAQ,EAAE,MAAM,CAAA;KAAE,KAAK,IAAI,CAAC;CAC3F;AAED;;;;;;;GAOG;AACH,UAAU,gBAAgB;IACxB,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IACjE,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IACjE,QAAQ,CAAC,OAAO,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IAClE,QAAQ,CAAC,OAAO,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IAClE,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AA2FD;;;;;;;;;GASG;AACH,wBAAgB,aAAa,CAC3B,KAAK,EAAE,QAAQ,EACf,CAAC,EAAE,MAAM,EACT,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EAClD,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EACnD,IAAI,EAAE;IAAE,SAAS,EAAE,MAAM,CAAC;IAAC,aAAa,EAAE,MAAM,CAAA;CAAE,GACjD,gBAAgB,CAYlB;AAMD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4DG;AACH,wBAAgB,oBAAoB,CAClC,YAAY,EAAE,QAAQ,EACtB,OAAO,EAAE,UAAU,GAClB,SAAS,WAAW,EAAE,CAmGxB"}
+{"version":3,"file":"gl4-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AA4BA;;;;;;GAMG;AACH,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,iCAAiC,CAAC;AAEtE;;;;;;;;;;;GAWG;AACH,MAAM,WAAW,QAAQ;IACvB,+BAA+B;IAC/B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,yCAAyC;IACzC,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAC/B;AAED;;;;;;;GAOG;AACH,MAAM,WAAW,WAAW;IAC1B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAChC;AAED;;;;GAIG;AACH,MAAM,WAAW,UAAU;IACzB;qCACiC;IACjC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,+CAA+C;IAC/C,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB;gFAC4E;IAC5E,QAAQ,CAAC,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC5D;;;;;OAKG;IACH,QAAQ,CAAC,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC7D,+DAA+D;IAC/D,QAAQ,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC;IAC5B,mEAAmE;IACnE,QAAQ,CAAC,aAAa,CAAC,EAAE,MAAM,CAAC;IAChC,+FAA+F;IAC/F,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;IACvB;;qCAEiC;IACjC,QAAQ,CAAC,eAAe,CAAC,EAAE,MAAM,CAAC;IAClC;;;;;;;;;;;;;;OAcG;IACH,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC,KAAK,EAAE;QAAE,IAAI,EAAE,MAAM,CAAC;QAAC,UAAU,EAAE,MAAM,CAAC;QAAC,QAAQ,EAAE,MAAM,CAAA;KAAE,KAAK,IAAI,CAAC;CAC3F;AAgDD;;;;;GAKG;AACH,MAAM,WAAW,aAAa;IAC5B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AAED;;;;;;GAMG;AACH,wBAAgB,OAAO,CACrB,KAAK,EAAE,QAAQ,EACf,CAAC,EAAE,MAAM,EACT,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EAClD,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EACnD,IAAI,EAAE;IAAE,SAAS,EAAE,MAAM,CAAC;IAAC,aAAa,EAAE,MAAM,CAAA;CAAE,GACjD,aAAa,CA6Bf;AAMD;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAmEG;AACH,wBAAgB,oBAAoB,CAClC,YAAY,EAAE,QAAQ,EACtB,OAAO,EAAE,UAAU,GAClB,SAAS,WAAW,EAAE,CA0GxB"}

package/dist/numerical/gl4-integrator.js

@@ -17,9 +17,9 @@
  * p_μ = g_μν dx^ν/dτ — not (x, v). This is what makes the flow symplectic
  * on T*M.
  *
- * This module holds the types, the Butcher constants, the implicit Picard
- * stage solver (`solveGL4Stage`, internal), and the integrator entry point
- * `integrateGeodesicGL4`.
+ * This module holds the types, the Butcher constants, the one-step driver
+ * (`gl4Step`, internal; the implicit Picard stages are MathTS
+ * `gaussLegendre4`'s), and the integrator entry point `integrateGeodesicGL4`.
  *
  * @module numerical/gl4-integrator
  */
@@ -41,14 +41,15 @@ export { GL4_A, GL4_B, GL4_C } from '@danielsimonjr/mathts-functions';
  *
  * `dg[μ*dim² + ν*dim + ρ] = ∂_μ g^{νρ}`. A coefficient that is exactly 0 is
  * skipped, so `0 * NaN` does not poison a component the term does not enter.
+ * `xBuf` is a reused coordinate buffer of length `dim`; it is the array the
+ * metric closures receive, so it is a plain `number[]` and needs no cast.
  */
 function geodesicDeriv(y, gInverseFn, dgInverseFn, xBuf) {
     const dim = xBuf.length;
     for (let i = 0; i < dim; i++)
         xBuf[i] = y[i];
-    const coords = xBuf;
-    const gInv = gInverseFn(coords);
-    const dgInv = dgInverseFn(coords);
+    const gInv = gInverseFn(xBuf);
+    const dgInv = dgInverseFn(xBuf);
     const out = new Array(dim * 2);
     for (let mu = 0; mu < dim; mu++) {
         let dx = 0;
@@ -76,15 +77,21 @@ function geodesicDeriv(y, gInverseFn, dgInverseFn, xBuf) {
     }
     return out;
 }
-/** One GL4 step of the geodesic Hamiltonian, via MathTS `gaussLegendre4`. */
-function advanceGl4(state, h, gInverseFn, dgInverseFn, opts) {
+/**
+ * One GL4 step of the geodesic Hamiltonian, via MathTS `gaussLegendre4`.
+ * Throws `GL4ConvergenceError` with a message matching
+ * `/Picard iteration did not converge/` when the Picard cap is exhausted;
+ * any other error from the metric closures propagates unchanged.
+ * @internal
+ */
+export function gl4Step(state, h, gInverseFn, dgInverseFn, opts) {
     const dim = state.x.length;
     const y0 = new Array(dim * 2);
     for (let i = 0; i < dim; i++) {
         y0[i] = state.x[i];
         y0[dim + i] = state.p[i];
     }
-    const xBuf = new Float64Array(dim);
+    const xBuf = new Array(dim).fill(0);
     let solved;
     try {
         solved = gaussLegendre4((_t, y) => geodesicDeriv(y, gInverseFn, dgInverseFn, xBuf), y0, [0, h], { steps: 1, picardTol: opts.picardTol, picardMaxIter: opts.picardMaxIter });
@@ -102,29 +109,6 @@ function advanceGl4(state, h, gInverseFn, dgInverseFn, opts) {
     }
     return { x: y1.slice(0, dim), p: y1.slice(dim), iterations };
 }
-/**
- * One GL4 step, reported in the stage-result shape the stage tests read.
- *
- * The step itself is MathTS `gaussLegendre4`. In flat space `p` is constant,
- * so both stage momenta equal that advanced `p`. Throws `GL4ConvergenceError`
- * with a message matching `/Picard iteration did not converge/` when the
- * Picard cap is exhausted.
- *
- * @internal
- */
-export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
-    const step = advanceGl4(state, h, gInverseFn, dgInverseFn, opts);
-    const dim = step.p.length;
-    const zeros = new Array(dim).fill(0);
-    // In flat space p is constant, so both stage momenta equal the advanced p.
-    return {
-        stageX: [step.x.slice(), step.x.slice()],
-        stageP: [step.p.slice(), step.p.slice()],
-        stageDx: [zeros, zeros.slice()],
-        stageDp: [zeros.slice(), zeros.slice()],
-        iterations: step.iterations,
-    };
-}
 // ---------------------------------------------------------------------------
 // Integrator entry-point (Task 3, Phase 1a-iii)
 // ---------------------------------------------------------------------------
@@ -173,10 +157,11 @@ export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
  * @param initialState — canonical (x, p) at τ = 0. Units follow the
  *   `gInverseFn` convention (see above).
  * @param options — see {@link GL4Options}:
- *   - `steps` — integer step count (dimensionless).
- *   - `tauMax` — affine-parameter extent (seconds in canonical SI).
- *   - `gInverseFn(x)[μ][ν]` — inverse metric g^{μν}(x).
- *   - `dgInverseFn(x)[λ][μ][ν]` — ∂_λ g^{μν}(x).
+ *   - `steps` — positive integer step count (dimensionless).
+ *   - `tauMax` — positive finite affine-parameter extent (seconds in
+ *     canonical SI). The integrator runs forward from τ = 0 only.
+ *   - `gInverseFn(x)[μ*dim + ν]` — inverse metric g^{μν}(x), row-major flat.
+ *   - `dgInverseFn(x)[λ*dim² + μ*dim + ν]` — ∂_λ g^{μν}(x), flat.
  *   - `picardTol` — convergence tolerance (dimensionless, default 1e-12).
  *   - `picardMaxIter` — fixed-point iteration cap (dimensionless, default 50).
  *   - `hMin` — step-halving floor (same units as `tauMax / steps`).
@@ -184,21 +169,31 @@ export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
  * @returns `steps + 1` snapshots: index `0` is the initial state, index `n`
  *   is the state after `n` steps (τ = n · h). Each snapshot carries `tau`,
  *   `x`, `p` (and optional `v` = g^{μν} p_ν) in the units chosen above.
- * @throws NumericalBackendError if `initialState.x[1] < domainMinRadius`.
+ * @throws NumericalBackendError if `steps` is not a positive integer, if
+ *   `tauMax` is not a positive finite number, or if
+ *   `initialState.x[1] < domainMinRadius`. The two option checks run before
+ *   any metric evaluation.
  * @throws GL4ConvergenceError if Picard fails even after step-halving to h_min.
+ *   Any other error thrown by `gInverseFn` or `dgInverseFn` (a domain
+ *   crossing the caller signals by throwing, or a programming error) is
+ *   rethrown unchanged; it is not halved away.
  *
  * @public
  */
 export function integrateGeodesicGL4(initialState, options) {
     const { steps, tauMax, gInverseFn, dgInverseFn, picardTol = 1e-12, picardMaxIter = 50, hMin, domainMinRadius, onStep, } = options;
+    if (!Number.isInteger(steps) || steps <= 0) {
+        throw new NumericalBackendError(`GL4 integrator: steps must be a positive integer, got ${steps}`);
+    }
+    if (!Number.isFinite(tauMax) || tauMax <= 0) {
+        throw new NumericalBackendError(`GL4 integrator: tauMax must be a positive finite number, got ${tauMax}`);
+    }
     if (domainMinRadius !== undefined && initialState.x[1] < domainMinRadius) {
         throw new NumericalBackendError(`GL4 integrator: initial r=${initialState.x[1]} < domainMinRadius=${domainMinRadius} (domain violation)`);
     }
     const h = tauMax / steps;
     const hFloor = hMin ?? h * 1e-9;
-    const snapshots = [
-        { tau: 0, x: initialState.x.slice(), p: initialState.p.slice() },
-    ];
+    const snapshots = [{ tau: 0, x: initialState.x.slice(), p: initialState.p.slice() }];
     let x = initialState.x.slice();
     let p = initialState.p.slice();
     const stateDim = x.length;
@@ -225,14 +220,18 @@ export function integrateGeodesicGL4(initialState, options) {
             let subH = trialH;
             while (subH >= hFloor) {
                 try {
-                    advanced = advanceGl4({ x, p }, subH, gInverseFn, dgInverseFn, {
+                    advanced = gl4Step({ x, p }, subH, gInverseFn, dgInverseFn, {
                         picardTol,
                         picardMaxIter,
                     });
                     stepSucceeded = true;
                     break;
                 }
-                catch {
+                catch (err) {
+                    // Only a Picard failure is a reason to halve. A throw from the metric
+                    // closures (a domain crossing the caller signals, or a bug) is that error.
+                    if (!(err instanceof GL4ConvergenceError))
+                        throw err;
                     subH /= 2;
                     halvings++;
                 }

package/dist/numerical/gl4-integrator.js.map

@@ -1 +1 @@
-{"version":3,"file":"gl4-integrator.js","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,OAAO,EAAE,cAAc,EAAE,mBAAmB,IAAI,yBAAyB,EAAE,MAAM,iCAAiC,CAAC;AACnH,OAAO,EAAE,mBAAmB,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEzE;;;;;;GAMG;AACH,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,iCAAiC,CAAC;AAoGtE;;;;;;;;GAQG;AACH,SAAS,aAAa,CACpB,CAAoB,EACpB,UAAkD,EAClD,WAAmD,EACnD,IAAkB;IAElB,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC;IACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;QAAE,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;IAC9C,MAAM,MAAM,GAAG,IAAoC,CAAC;IACpD,MAAM,IAAI,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IAChC,MAAM,KAAK,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC;IAClC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,GAAG,CAAC,CAAC,CAAC;IACvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;QAChC,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,MAAM,IAAI,GAAG,EAAE,GAAG,GAAG,CAAC;QACtB,MAAM,MAAM,GAAG,IAAI,GAAG,GAAG,CAAC;QAC1B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAE,CAAC;YACzB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,GAAG,EAAE,CAAE,CAAC;YAC3B,wFAAwF;YACxF,IAAI,CAAC,KAAK,CAAC;gBAAE,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC;YAC3B,IAAI,IAAI,GAAG,CAAC,CAAC;YACb,MAAM,MAAM,GAAG,MAAM,GAAG,EAAE,GAAG,GAAG,CAAC;YACjC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;gBACnC,MAAM,EAAE,GAAG,KAAK,CAAC,MAAM,GAAG,GAAG,CAAE,CAAC;gBAChC,IAAI,EAAE,KAAK,CAAC;oBAAE,IAAI,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAE,CAAC;YAC3C,CAAC;YACD,IAAI,IAAI,KAAK,CAAC;gBAAE,EAAE,IAAI,IAAI,GAAG,GAAG,CAAC;QACnC,CAAC;QACD,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;QACb,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;IAC3C,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAQD,6EAA6E;AAC7E,SAAS,UAAU,CACjB,KAAe,EACf,CAAS,EACT,UAAkD,EAClD,WAAmD,EACnD,IAAkD;IAElD,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;IAC3B,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,GAAG,GAAG,CAAC,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;QACpB,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;IAC5B,CAAC;IACD,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACnC,IAAI,MAAyC,CAAC;IAC9C,IAAI,CAAC;QACH,MAAM,GAAG,cAAc,CACrB,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,aAAa,CAAC,CAAC,EAAE,UAAU,EAAE,WAAW,EAAE,IAAI,CAAC,EAC1D,EAAE,EACF,CAAC,CAAC,EAAE,CAAC,CAAC,EACN,EAAE,KAAK,EAAE,CAAC,EAAE,SAAS,EAAE,IAAI,CAAC,SAAS,EAAE,aAAa,EAAE,IAAI,CAAC,aAAa,EAAE,CAC3E,CAAC;IACJ,CAAC;IAAC,OAAO,GAAG,EAAE,CAAC;QACb,4FAA4F;QAC5F,IAAI,GAAG,YAAY,yBAAyB;YAAE,MAAM,IAAI,mBAAmB,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QACzF,MAAM,GAAG,CAAC;IACZ,CAAC;IACD,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACzC,MAAM,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;IACxC,IAAI,EAAE,KAAK,SAAS,IAAI,UAAU,KAAK,SAAS,EAAE,CAAC;QACjD,MAAM,IAAI,mBAAmB,CAC3B,wCAAwC,IAAI,CAAC,aAAa,yCAAyC,IAAI,CAAC,SAAS,GAAG,CACrH,CAAC;IACJ,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,UAAU,EAAE,CAAC;AAC/D,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,aAAa,CAC3B,KAAe,EACf,CAAS,EACT,UAAkD,EAClD,WAAmD,EACnD,IAAkD;IAElD,MAAM,IAAI,GAAG,UAAU,CAAC,KAAK,EAAE,CAAC,EAAE,UAAU,EAAE,WAAW,EAAE,IAAI,CAAC,CAAC;IACjE,MAAM,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC;IAC1B,MAAM,KAAK,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC7C,2EAA2E;IAC3E,OAAO;QACL,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;QACxC,MAAM,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,EAAE,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;QACxC,OAAO,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC,KAAK,EAAE,CAAC;QAC/B,OAAO,EAAE,CAAC,KAAK,CAAC,KAAK,EAAE,EAAE,KAAK,CAAC,KAAK,EAAE,CAAC;QACvC,UAAU,EAAE,IAAI,CAAC,UAAU;KAC5B,CAAC;AACJ,CAAC;AAED,8EAA8E;AAC9E,gDAAgD;AAChD,8EAA8E;AAE9E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4DG;AACH,MAAM,UAAU,oBAAoB,CAClC,YAAsB,EACtB,OAAmB;IAEnB,MAAM,EACJ,KAAK,EACL,MAAM,EACN,UAAU,EACV,WAAW,EACX,SAAS,GAAG,KAAK,EACjB,aAAa,GAAG,EAAE,EAClB,IAAI,EACJ,eAAe,EACf,MAAM,GACP,GAAG,OAAO,CAAC;IAEZ,IAAI,eAAe,KAAK,SAAS,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACzE,MAAM,IAAI,qBAAqB,CAC7B,6BAA6B,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,sBAAsB,eAAe,qBAAqB,CACzG,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC;IACzB,MAAM,MAAM,GAAG,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC;IAChC,MAAM,SAAS,GAAkB;QAC/B,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,EAAE;KACzF,CAAC;IACF,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,CAAC;IAC3C,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,CAAC;IAE3C,MAAM,QAAQ,GAAG,CAAC,CAAC,MAAM,CAAC;IAC1B,IAAI,IAAI,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC;IACvC,IAAI,IAAI,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC;IAEvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,8DAA8D;QAC9D,uEAAuE;QACvE,yEAAyE;QACzE,uEAAuE;QACvE,wEAAwE;QACxE,sEAAsE;QACtE,4DAA4D;QAC5D,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,cAAc,GAAG,CAAC,CAAC;QACvB,+EAA+E;QAC/E,MAAM,SAAS,GAAG,CAAC,GAAG,KAAK,CAAC;QAC5B,OAAO,SAAS,GAAG,SAAS,EAAE,CAAC;YAC7B,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;YAC1C,IAAI,QAAgC,CAAC;YACrC,IAAI,aAAa,GAAG,KAAK,CAAC;YAC1B,IAAI,IAAI,GAAG,MAAM,CAAC;YAClB,OAAO,IAAI,IAAI,MAAM,EAAE,CAAC;gBACtB,IAAI,CAAC;oBACH,QAAQ,GAAG,UAAU,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,IAAI,EAAE,UAAU,EAAE,WAAW,EAAE;wBAC7D,SAAS;wBACT,aAAa;qBACd,CAAC,CAAC;oBACH,aAAa,GAAG,IAAI,CAAC;oBACrB,MAAM;gBACR,CAAC;gBAAC,MAAM,CAAC;oBACP,IAAI,IAAI,CAAC,CAAC;oBACV,QAAQ,EAAE,CAAC;gBACb,CAAC;YACH,CAAC;YACD,IAAI,CAAC,aAAa,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;gBAC7C,MAAM,IAAI,mBAAmB,CAC3B,mEAAmE,MAAM,UAAU,CAAC,8CAA8C,CACnI,CAAC;YACJ,CAAC;YACD,IAAI,QAAQ,KAAK,CAAC,EAAE,CAAC;gBACnB,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;YACrD,CAAC;iBAAM,CAAC;gBACN,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC;oBACrC,IAAI,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAE,CAAC;oBAC3B,IAAI,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAE,CAAC;gBAC7B,CAAC;YACH,CAAC;YAED,MAAM,IAAI,GAAG,CAAC,CAAC;YACf,CAAC,GAAG,IAAI,CAAC;YACT,IAAI,GAAG,IAAI,CAAC;YAEZ,MAAM,IAAI,GAAG,CAAC,CAAC;YACf,CAAC,GAAG,IAAI,CAAC;YACT,IAAI,GAAG,IAAI,CAAC;YAEZ,SAAS,IAAI,IAAI,CAAC;YAClB,KAAK,GAAG,IAAI,CAAC,CAAC,2DAA2D;YACzE,cAAc,GAAG,QAAQ,CAAC,UAAU,CAAC;QACvC,CAAC;QACD,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACjE,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;YACzB,MAAM,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,cAAc,EAAE,QAAQ,EAAE,CAAC,CAAC;QAC5D,CAAC;IACH,CAAC;IAED,OAAO,SAAS,CAAC;AACnB,CAAC"}
+{"version":3,"file":"gl4-integrator.js","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,OAAO,EAAE,cAAc,EAAE,mBAAmB,IAAI,yBAAyB,EAAE,MAAM,iCAAiC,CAAC;AACnH,OAAO,EAAE,mBAAmB,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEzE;;;;;;GAMG;AACH,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,iCAAiC,CAAC;AAqFtE;;;;;;;;;;GAUG;AACH,SAAS,aAAa,CACpB,CAAoB,EACpB,UAAkD,EAClD,WAAmD,EACnD,IAAc;IAEd,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC;IACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;QAAE,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAE,CAAC;IAC9C,MAAM,IAAI,GAAG,UAAU,CAAC,IAAI,CAAC,CAAC;IAC9B,MAAM,KAAK,GAAG,WAAW,CAAC,IAAI,CAAC,CAAC;IAChC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,GAAG,CAAC,CAAC,CAAC;IACvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;QAChC,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,MAAM,IAAI,GAAG,EAAE,GAAG,GAAG,CAAC;QACtB,MAAM,MAAM,GAAG,IAAI,GAAG,GAAG,CAAC;QAC1B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAE,CAAC;YACzB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,GAAG,EAAE,CAAE,CAAC;YAC3B,wFAAwF;YACxF,IAAI,CAAC,KAAK,CAAC;gBAAE,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC;YAC3B,IAAI,IAAI,GAAG,CAAC,CAAC;YACb,MAAM,MAAM,GAAG,MAAM,GAAG,EAAE,GAAG,GAAG,CAAC;YACjC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;gBACnC,MAAM,EAAE,GAAG,KAAK,CAAC,MAAM,GAAG,GAAG,CAAE,CAAC;gBAChC,IAAI,EAAE,KAAK,CAAC;oBAAE,IAAI,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAE,CAAC;YAC3C,CAAC;YACD,IAAI,IAAI,KAAK,CAAC;gBAAE,EAAE,IAAI,IAAI,GAAG,GAAG,CAAC;QACnC,CAAC;QACD,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC;QACb,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;IAC3C,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAcD;;;;;;GAMG;AACH,MAAM,UAAU,OAAO,CACrB,KAAe,EACf,CAAS,EACT,UAAkD,EAClD,WAAmD,EACnD,IAAkD;IAElD,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;IAC3B,MAAM,EAAE,GAAG,IAAI,KAAK,CAAS,GAAG,GAAG,CAAC,CAAC,CAAC;IACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;QACpB,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;IAC5B,CAAC;IACD,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC5C,IAAI,MAAyC,CAAC;IAC9C,IAAI,CAAC;QACH,MAAM,GAAG,cAAc,CACrB,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,aAAa,CAAC,CAAC,EAAE,UAAU,EAAE,WAAW,EAAE,IAAI,CAAC,EAC1D,EAAE,EACF,CAAC,CAAC,EAAE,CAAC,CAAC,EACN,EAAE,KAAK,EAAE,CAAC,EAAE,SAAS,EAAE,IAAI,CAAC,SAAS,EAAE,aAAa,EAAE,IAAI,CAAC,aAAa,EAAE,CAC3E,CAAC;IACJ,CAAC;IAAC,OAAO,GAAG,EAAE,CAAC;QACb,4FAA4F;QAC5F,IAAI,GAAG,YAAY,yBAAyB;YAAE,MAAM,IAAI,mBAAmB,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;QACzF,MAAM,GAAG,CAAC;IACZ,CAAC;IACD,MAAM,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACzC,MAAM,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;IACxC,IAAI,EAAE,KAAK,SAAS,IAAI,UAAU,KAAK,SAAS,EAAE,CAAC;QACjD,MAAM,IAAI,mBAAmB,CAC3B,wCAAwC,IAAI,CAAC,aAAa,yCAAyC,IAAI,CAAC,SAAS,GAAG,CACrH,CAAC;IACJ,CAAC;IACD,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,UAAU,EAAE,CAAC;AAC/D,CAAC;AAED,8EAA8E;AAC9E,gDAAgD;AAChD,8EAA8E;AAE9E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAmEG;AACH,MAAM,UAAU,oBAAoB,CAClC,YAAsB,EACtB,OAAmB;IAEnB,MAAM,EACJ,KAAK,EACL,MAAM,EACN,UAAU,EACV,WAAW,EACX,SAAS,GAAG,KAAK,EACjB,aAAa,GAAG,EAAE,EAClB,IAAI,EACJ,eAAe,EACf,MAAM,GACP,GAAG,OAAO,CAAC;IAEZ,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QAC3C,MAAM,IAAI,qBAAqB,CAAC,yDAAyD,KAAK,EAAE,CAAC,CAAC;IACpG,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,MAAM,IAAI,CAAC,EAAE,CAAC;QAC5C,MAAM,IAAI,qBAAqB,CAAC,gEAAgE,MAAM,EAAE,CAAC,CAAC;IAC5G,CAAC;IACD,IAAI,eAAe,KAAK,SAAS,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACzE,MAAM,IAAI,qBAAqB,CAC7B,6BAA6B,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,sBAAsB,eAAe,qBAAqB,CACzG,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC;IACzB,MAAM,MAAM,GAAG,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC;IAChC,MAAM,SAAS,GAAkB,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;IACpG,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;IAC/B,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC;IAE/B,MAAM,QAAQ,GAAG,CAAC,CAAC,MAAM,CAAC;IAC1B,IAAI,IAAI,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC;IACvC,IAAI,IAAI,GAAG,IAAI,KAAK,CAAS,QAAQ,CAAC,CAAC;IAEvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,8DAA8D;QAC9D,uEAAuE;QACvE,yEAAyE;QACzE,uEAAuE;QACvE,wEAAwE;QACxE,sEAAsE;QACtE,4DAA4D;QAC5D,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,cAAc,GAAG,CAAC,CAAC;QACvB,+EAA+E;QAC/E,MAAM,SAAS,GAAG,CAAC,GAAG,KAAK,CAAC;QAC5B,OAAO,SAAS,GAAG,SAAS,EAAE,CAAC;YAC7B,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;YAC1C,IAAI,QAAmC,CAAC;YACxC,IAAI,aAAa,GAAG,KAAK,CAAC;YAC1B,IAAI,IAAI,GAAG,MAAM,CAAC;YAClB,OAAO,IAAI,IAAI,MAAM,EAAE,CAAC;gBACtB,IAAI,CAAC;oBACH,QAAQ,GAAG,OAAO,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,IAAI,EAAE,UAAU,EAAE,WAAW,EAAE;wBAC1D,SAAS;wBACT,aAAa;qBACd,CAAC,CAAC;oBACH,aAAa,GAAG,IAAI,CAAC;oBACrB,MAAM;gBACR,CAAC;gBAAC,OAAO,GAAG,EAAE,CAAC;oBACb,sEAAsE;oBACtE,2EAA2E;oBAC3E,IAAI,CAAC,CAAC,GAAG,YAAY,mBAAmB,CAAC;wBAAE,MAAM,GAAG,CAAC;oBACrD,IAAI,IAAI,CAAC,CAAC;oBACV,QAAQ,EAAE,CAAC;gBACb,CAAC;YACH,CAAC;YACD,IAAI,CAAC,aAAa,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;gBAC7C,MAAM,IAAI,mBAAmB,CAC3B,mEAAmE,MAAM,UAAU,CAAC,8CAA8C,CACnI,CAAC;YACJ,CAAC;YACD,IAAI,QAAQ,KAAK,CAAC,EAAE,CAAC;gBACnB,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBACnD,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;YACrD,CAAC;iBAAM,CAAC;gBACN,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,QAAQ,EAAE,EAAE,EAAE,EAAE,CAAC;oBACrC,IAAI,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAE,CAAC;oBAC3B,IAAI,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAE,CAAC;gBAC7B,CAAC;YACH,CAAC;YAED,MAAM,IAAI,GAAG,CAAC,CAAC;YACf,CAAC,GAAG,IAAI,CAAC;YACT,IAAI,GAAG,IAAI,CAAC;YAEZ,MAAM,IAAI,GAAG,CAAC,CAAC;YACf,CAAC,GAAG,IAAI,CAAC;YACT,IAAI,GAAG,IAAI,CAAC;YAEZ,SAAS,IAAI,IAAI,CAAC;YAClB,KAAK,GAAG,IAAI,CAAC,CAAC,2DAA2D;YACzE,cAAc,GAAG,QAAQ,CAAC,UAAU,CAAC;QACvC,CAAC;QACD,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACjE,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;YACzB,MAAM,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,cAAc,EAAE,QAAQ,EAAE,CAAC,CAAC;QAC5D,CAAC;IACH,CAAC;IAED,OAAO,SAAS,CAAC;AACnB,CAAC","sourcesContent":["/**\n * Gauss-Legendre 4th-order (GL4) symplectic integrator — types + Butcher\n * tableau scaffold (v0.5.0 Task 1, Phase 1a-i).\n *\n * GL4 is a 2-stage implicit Runge-Kutta method built on the roots of the\n * shifted Legendre polynomial P₂, with order p = 4 and stage order s = 2.\n * It is symplectic for non-separable Hamiltonians — the property that\n * justifies its selection over Ruth-4 for the geodesic Hamiltonian\n *\n *   H(x, p) = ½ g^{μν}(x) p_μ p_ν,\n *\n * which is non-separable because g^{μν} depends on x. (v0.5.0 Decision #2,\n * post-adversarial reconciliation; Sanz-Serna 1988; Hairer/Lubich/Wanner\n * \"Geometric Numerical Integration\" §II.1.)\n *\n * The state is canonical (x, p) (Decision #3) — covariant momentum\n * p_μ = g_μν dx^ν/dτ — not (x, v). This is what makes the flow symplectic\n * on T*M.\n *\n * This module holds the types, the Butcher constants, the one-step driver\n * (`gl4Step`, internal; the implicit Picard stages are MathTS\n * `gaussLegendre4`'s), and the integrator entry point `integrateGeodesicGL4`.\n *\n * @module numerical/gl4-integrator\n */\nimport { gaussLegendre4, GL4ConvergenceError as MathTsGL4ConvergenceError } from '@danielsimonjr/mathts-functions';\nimport { GL4ConvergenceError, NumericalBackendError } from './errors.js';\n\n/**\n * Gauss–Legendre order-4 Butcher tableau from MathTS. `GL4_A` is the 2×2\n * stage matrix, `GL4_B` the weights (both ½), and `GL4_C` the two nodes.\n * `integrateGeodesicGL4` does not read these names; the tableau tests do.\n *\n * @internal\n */\nexport { GL4_A, GL4_B, GL4_C } from '@danielsimonjr/mathts-functions';\n\n/**\n * Canonical (x, p) phase-space state for the geodesic flow on T*M.\n *\n *   x^μ        — coordinate 4-vector\n *   p_μ        — covariant momentum, p_μ = g_μν dx^ν/dτ\n *\n * Decision #3 (v0.5.0): the canonical state is (x, p), not (x, v). The\n * symplectic 2-form ω = dp_μ ∧ dx^μ is preserved by the GL4 flow only on\n * this representation.\n *\n * @public\n */\nexport interface GL4State {\n  /** Coordinate 4-vector x^μ. */\n  readonly x: readonly number[];\n  /** Covariant momentum p_μ = g_μν v^ν. */\n  readonly p: readonly number[];\n}\n\n/**\n * Per-step snapshot recorded by the integrator. `v` (the contravariant\n * 4-velocity v^μ = g^{μν} p_ν) is declared optional and is never emitted:\n * no option asks for it, and nothing in `src/` reads it. It stays so a caller's\n * own snapshot type can carry one.\n *\n * @public\n */\nexport interface GL4Snapshot {\n  readonly tau: number;\n  readonly x: readonly number[];\n  readonly p: readonly number[];\n  readonly v?: readonly number[];\n}\n\n/**\n * Options for `integrateGeodesicGL4`.\n *\n * @public\n */\nexport interface GL4Options {\n  /** Number of integration steps (uniform-step baseline; may be subdivided\n   *  by adaptive step-halving). */\n  readonly steps: number;\n  /** Final proper time τ_max (initial τ = 0). */\n  readonly tauMax: number;\n  /** Inverse-metric closure (v0.9.0 O-1 flat layout):\n   *  `gInverseFn(x)[μ*dim + ν] = g^{μν}(x)`, row-major Float64Array(dim²). */\n  readonly gInverseFn: (x: readonly number[]) => Float64Array;\n  /**\n   * Partial derivatives of the inverse metric.\n   * Index order (v0.9.0 O-1 flat layout):\n   * `dgInverseFn(x)[lambda*dim² + mu*dim + nu] = ∂_lambda g^{mu nu}` at coords x.\n   * (I2: axis semantics pinned here to prevent silent transposition bugs.)\n   */\n  readonly dgInverseFn: (x: readonly number[]) => Float64Array;\n  /** Picard fixed-point tolerance (default chosen in Task 2). */\n  readonly picardTol?: number;\n  /** Picard fixed-point iteration cap (default chosen in Task 2). */\n  readonly picardMaxIter?: number;\n  /** Adaptive step-halving floor (I4). If step-halving reaches h_min, throws with diagnostic. */\n  readonly hMin?: number;\n  /** Minimum radial coordinate (or domain analog) — abort integration if\n   *  the trajectory crosses inside this radius (e.g., the Schwarzschild\n   *  event horizon at r = r_s). */\n  readonly domainMinRadius?: number;\n  /**\n   * v0.5.1 PD-4: opt-in per-step diagnostics callback. Fires once per\n   * successful step with the Picard iteration count consumed and whether\n   * adaptive step-halving had to subdivide (an \"exhaustion\" event from the\n   * caller's perspective: the original h failed Picard and was halved).\n   *\n   * Used by the gated `GL4_LONG=1` Mercury 100-orbit Picard-convergence\n   * test to measure the failure fraction across millions of steps without\n   * polluting the integrator's return shape for normal callers.\n   *\n   * - `iterations`: Picard iteration count actually consumed at the\n   *   successful step size (always 1..picardMaxIter).\n   * - `halvings`: number of step-halvings the step required before\n   *   succeeding (0 = first try; ≥1 means original h hit picardMaxIter).\n   */\n  readonly onStep?: (event: { step: number; iterations: number; halvings: number }) => void;\n}\n\n/**\n * Geodesic Hamiltonian derivative, packed as `y = [x..., p...]`.\n *\n *   dx^μ = g^{μν} p_ν\n *   dp_μ = −½ (∂_μ g^{νρ}) p_ν p_ρ\n *\n * `dg[μ*dim² + ν*dim + ρ] = ∂_μ g^{νρ}`. A coefficient that is exactly 0 is\n * skipped, so `0 * NaN` does not poison a component the term does not enter.\n * `xBuf` is a reused coordinate buffer of length `dim`; it is the array the\n * metric closures receive, so it is a plain `number[]` and needs no cast.\n */\nfunction geodesicDeriv(\n  y: readonly number[],\n  gInverseFn: (x: readonly number[]) => Float64Array,\n  dgInverseFn: (x: readonly number[]) => Float64Array,\n  xBuf: number[],\n): number[] {\n  const dim = xBuf.length;\n  for (let i = 0; i < dim; i++) xBuf[i] = y[i]!;\n  const gInv = gInverseFn(xBuf);\n  const dgInv = dgInverseFn(xBuf);\n  const out = new Array<number>(dim * 2);\n  for (let mu = 0; mu < dim; mu++) {\n    let dx = 0;\n    let dp = 0;\n    const gRow = mu * dim;\n    const dgBase = gRow * dim;\n    for (let nu = 0; nu < dim; nu++) {\n      const pNu = y[dim + nu]!;\n      const g = gInv[gRow + nu]!;\n      // Skip an exact zero so `0 * NaN` does not poison a component this term does not enter.\n      if (g !== 0) dx += g * pNu;\n      let pDot = 0;\n      const offset = dgBase + nu * dim;\n      for (let rho = 0; rho < dim; rho++) {\n        const dg = dgInv[offset + rho]!;\n        if (dg !== 0) pDot += dg * y[dim + rho]!;\n      }\n      if (pDot !== 0) dp += pDot * pNu;\n    }\n    out[mu] = dx;\n    out[dim + mu] = dp === 0 ? 0 : -0.5 * dp;\n  }\n  return out;\n}\n\n/**\n * The state after one GL4 step, with the Picard iteration count that step\n * consumed. MathTS `gaussLegendre4` does not expose its two internal stage\n * values, so this is the advanced state and nothing else.\n * @internal\n */\nexport interface GL4StepResult {\n  readonly x: readonly number[];\n  readonly p: readonly number[];\n  readonly iterations: number;\n}\n\n/**\n * One GL4 step of the geodesic Hamiltonian, via MathTS `gaussLegendre4`.\n * Throws `GL4ConvergenceError` with a message matching\n * `/Picard iteration did not converge/` when the Picard cap is exhausted;\n * any other error from the metric closures propagates unchanged.\n * @internal\n */\nexport function gl4Step(\n  state: GL4State,\n  h: number,\n  gInverseFn: (x: readonly number[]) => Float64Array,\n  dgInverseFn: (x: readonly number[]) => Float64Array,\n  opts: { picardTol: number; picardMaxIter: number },\n): GL4StepResult {\n  const dim = state.x.length;\n  const y0 = new Array<number>(dim * 2);\n  for (let i = 0; i < dim; i++) {\n    y0[i] = state.x[i]!;\n    y0[dim + i] = state.p[i]!;\n  }\n  const xBuf = new Array<number>(dim).fill(0);\n  let solved: ReturnType<typeof gaussLegendre4>;\n  try {\n    solved = gaussLegendre4(\n      (_t, y) => geodesicDeriv(y, gInverseFn, dgInverseFn, xBuf),\n      y0,\n      [0, h],\n      { steps: 1, picardTol: opts.picardTol, picardMaxIter: opts.picardMaxIter },\n    );\n  } catch (err) {\n    // MathTS throws its own typed GL4ConvergenceError for a Picard step that does not converge.\n    if (err instanceof MathTsGL4ConvergenceError) throw new GL4ConvergenceError(err.message);\n    throw err;\n  }\n  const y1 = solved.y[solved.y.length - 1];\n  const iterations = solved.iterations[0];\n  if (y1 === undefined || iterations === undefined) {\n    throw new GL4ConvergenceError(\n      `Picard iteration did not converge in ${opts.picardMaxIter} iterations (maxDelta above picardTol=${opts.picardTol})`,\n    );\n  }\n  return { x: y1.slice(0, dim), p: y1.slice(dim), iterations };\n}\n\n// ---------------------------------------------------------------------------\n// Integrator entry-point (Task 3, Phase 1a-iii)\n// ---------------------------------------------------------------------------\n\n/**\n * GL4 symplectic integrator on the canonical (x, p) geodesic Hamiltonian.\n *\n *   H(x, p) = ½ g^{μν}(x) p_μ p_ν\n *\n * Each accepted step is MathTS `gaussLegendre4` on the geodesic Hamiltonian.\n * Adaptive step-halving stays here. `onStep` reads that step's `iterations`.\n * with **adaptive step-halving on Picard non-convergence** (Adam+Eve I4,\n * replaces the single-retry R8): if Picard fails at step size h, retry at\n * h/2, h/4, … down to `hMin` (default `h · 1e-9`); throw\n * `GL4ConvergenceError` with a diagnostic message only when h_min is also\n * exhausted.\n *\n * Symplecticity (preservation of ω = dp_μ ∧ dx^μ) is a property of the\n * Butcher tableau, not of the inner Picard solver — see Sanz-Serna 1988,\n * Hairer/Lubich/Wanner §II.1. Hamiltonian drift over long integrations is\n * bounded; for non-resonant systems it remains O(h^p) over exponentially\n * long times (`p = 4` for GL4).\n *\n * **Domain guard.** If `domainMinRadius` is provided and `initialState.x[1]`\n * (radial coordinate) is below the bound, throws `NumericalBackendError`\n * synchronously with a `/domain/i`-matching message. The mid-trajectory\n * domain crossing is not checked here — callers needing that supply a\n * `gInverseFn` that throws on out-of-domain input.\n *\n * **Units.** The integrator is metric-agnostic — units follow the units of\n * the supplied `gInverseFn` and `initialState`. For UPT's canonical SI\n * Schwarzschild applications (the covariant-eikonal method used by\n * catalog relation be-37, and the perihelion method used by catalog\n * relation be-52):\n *   - `initialState.x` — `(t, r, θ, φ)` in **(s, m, rad, rad)** (SI).\n *   - `initialState.p` — covariant 4-momentum `p_μ = g_μν v^ν` in\n *     **(J·s, kg·m, kg·m², kg·m²)** under the affine normalization\n *     `p_t = −c²` used by `evaluateCovariantEikonalNumerical`.\n *   - `tauMax` — affine-parameter (proper-time for timelike, coordinate-\n *     time-like for the null normalization) extent in **seconds** under\n *     the covariant-eikonal method's convention (catalog relation be-37);\n *     **dimensionless** if the caller chose\n *     geometric units. The integrator does not enforce a choice.\n *   - `domainMinRadius` — radial coordinate lower bound in the same length\n *     units as `initialState.x[1]` (typically **metres** for SI).\n *\n * @param initialState — canonical (x, p) at τ = 0. Units follow the\n *   `gInverseFn` convention (see above).\n * @param options — see {@link GL4Options}:\n *   - `steps` — positive integer step count (dimensionless).\n *   - `tauMax` — positive finite affine-parameter extent (seconds in\n *     canonical SI). The integrator runs forward from τ = 0 only.\n *   - `gInverseFn(x)[μ*dim + ν]` — inverse metric g^{μν}(x), row-major flat.\n *   - `dgInverseFn(x)[λ*dim² + μ*dim + ν]` — ∂_λ g^{μν}(x), flat.\n *   - `picardTol` — convergence tolerance (dimensionless, default 1e-12).\n *   - `picardMaxIter` — fixed-point iteration cap (dimensionless, default 50).\n *   - `hMin` — step-halving floor (same units as `tauMax / steps`).\n *   - `domainMinRadius` — radial cutoff (same units as `initialState.x[1]`).\n * @returns `steps + 1` snapshots: index `0` is the initial state, index `n`\n *   is the state after `n` steps (τ = n · h). Each snapshot carries `tau`,\n *   `x`, `p` (and optional `v` = g^{μν} p_ν) in the units chosen above.\n * @throws NumericalBackendError if `steps` is not a positive integer, if\n *   `tauMax` is not a positive finite number, or if\n *   `initialState.x[1] < domainMinRadius`. The two option checks run before\n *   any metric evaluation.\n * @throws GL4ConvergenceError if Picard fails even after step-halving to h_min.\n *   Any other error thrown by `gInverseFn` or `dgInverseFn` (a domain\n *   crossing the caller signals by throwing, or a programming error) is\n *   rethrown unchanged; it is not halved away.\n *\n * @public\n */\nexport function integrateGeodesicGL4(\n  initialState: GL4State,\n  options: GL4Options,\n): readonly GL4Snapshot[] {\n  const {\n    steps,\n    tauMax,\n    gInverseFn,\n    dgInverseFn,\n    picardTol = 1e-12,\n    picardMaxIter = 50,\n    hMin,\n    domainMinRadius,\n    onStep,\n  } = options;\n\n  if (!Number.isInteger(steps) || steps <= 0) {\n    throw new NumericalBackendError(`GL4 integrator: steps must be a positive integer, got ${steps}`);\n  }\n  if (!Number.isFinite(tauMax) || tauMax <= 0) {\n    throw new NumericalBackendError(`GL4 integrator: tauMax must be a positive finite number, got ${tauMax}`);\n  }\n  if (domainMinRadius !== undefined && initialState.x[1] < domainMinRadius) {\n    throw new NumericalBackendError(\n      `GL4 integrator: initial r=${initialState.x[1]} < domainMinRadius=${domainMinRadius} (domain violation)`,\n    );\n  }\n\n  const h = tauMax / steps;\n  const hFloor = hMin ?? h * 1e-9;\n  const snapshots: GL4Snapshot[] = [{ tau: 0, x: initialState.x.slice(), p: initialState.p.slice() }];\n  let x = initialState.x.slice();\n  let p = initialState.p.slice();\n\n  const stateDim = x.length;\n  let newX = new Array<number>(stateDim);\n  let newP = new Array<number>(stateDim);\n\n  for (let n = 0; n < steps; n++) {\n    // I4: adaptive step-halving loop (not single-retry) on Picard\n    // non-convergence. When the full step fails, we advance by the SMALLER\n    // converging step and sub-step until the macro-step `h` is covered — the\n    // stages are solved at `stepH`, so the state/τ update MUST use `stepH`\n    // too (advancing by full `h` after a halved solve destroys accuracy and\n    // symplecticity). τ is reported on the fixed `(n+1)·h` grid; only the\n    // residual-to-cover accumulator tracks intra-step progress.\n    let remaining = h;\n    let stepH = h;\n    let halvings = 0;\n    let lastIterations = 0;\n    // Guard against FP residue: treat anything below a tiny fraction of h as done.\n    const remainEps = h * 1e-12;\n    while (remaining > remainEps) {\n      const trialH = Math.min(stepH, remaining);\n      let advanced: GL4StepResult | undefined;\n      let stepSucceeded = false;\n      let subH = trialH;\n      while (subH >= hFloor) {\n        try {\n          advanced = gl4Step({ x, p }, subH, gInverseFn, dgInverseFn, {\n            picardTol,\n            picardMaxIter,\n          });\n          stepSucceeded = true;\n          break;\n        } catch (err) {\n          // Only a Picard failure is a reason to halve. A throw from the metric\n          // closures (a domain crossing the caller signals, or a bug) is that error.\n          if (!(err instanceof GL4ConvergenceError)) throw err;\n          subH /= 2;\n          halvings++;\n        }\n      }\n      if (!stepSucceeded || advanced === undefined) {\n        throw new GL4ConvergenceError(\n          `GL4 integrator: Picard iteration did not converge even at h_min=${hFloor} (step ${n}). Diagnose step-size or metric singularity.`,\n        );\n      }\n      if (stateDim === 4) {\n        newX[0] = advanced.x[0]!; newP[0] = advanced.p[0]!;\n        newX[1] = advanced.x[1]!; newP[1] = advanced.p[1]!;\n        newX[2] = advanced.x[2]!; newP[2] = advanced.p[2]!;\n        newX[3] = advanced.x[3]!; newP[3] = advanced.p[3]!;\n      } else {\n        for (let mu = 0; mu < stateDim; mu++) {\n          newX[mu] = advanced.x[mu]!;\n          newP[mu] = advanced.p[mu]!;\n        }\n      }\n\n      const tmpX = x;\n      x = newX;\n      newX = tmpX;\n\n      const tmpP = p;\n      p = newP;\n      newP = tmpP;\n\n      remaining -= subH;\n      stepH = subH; // keep the converging size for the rest of this macro-step\n      lastIterations = advanced.iterations;\n    }\n    snapshots.push({ tau: (n + 1) * h, x: x.slice(), p: p.slice() });\n    if (onStep !== undefined) {\n      onStep({ step: n, iterations: lastIterations, halvings });\n    }\n  }\n\n  return snapshots;\n}\n"]}

package/dist/numerical/grid-field.js.map

@@ -1 +1 @@
-{"version":3,"file":"grid-field.js","sourceRoot":"","sources":["../../src/numerical/grid-field.ts"],"names":[],"mappings":""}
+{"version":3,"file":"grid-field.js","sourceRoot":"","sources":["../../src/numerical/grid-field.ts"],"names":[],"mappings":"","sourcesContent":["/**\n * GridField — a sampled field on a regular grid, for the 'grid' numericalForm\n * finite-difference path. See v0.3.5-Design.md §5.\n *\n * The interface now lives in the leaf `./types.js` module (it only depends on\n * `NestedArray`), which removes the former `types ↔ grid-field` type cycle. This\n * module re-exports it so the established `numerical/grid-field` import path —\n * used by `pderiv.ts`, `numerical/index.ts`, and the public `src/index.ts` — keeps\n * working unchanged.\n *\n * @module numerical/grid-field\n */\nexport type { GridField } from './types.js';\n"]}

package/dist/numerical/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/numerical/index.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAGvD,OAAO,EAAE,SAAS,EAAE,MAAM,eAAe,CAAC;AAC1C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,qBAAqB,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAI/E,cAAc;AACd,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,oBAAoB,CAAC;AAK/E,cAAc;AACd,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACxE,cAAc;AACd,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD;;0EAE0E;AAC1E,OAAO,EAAE,0BAA0B,EAAE,MAAM,0BAA0B,CAAC;AACtE;;;;GAIG;AACH,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAC5D,mEAAmE;AACnE,OAAO,EACL,iCAAiC,GAClC,MAAM,wBAAwB,CAAC;AAOhC,kDAAkD;AAClD,OAAO,EAAE,oBAAoB,EAAE,MAAM,qBAAqB,CAAC;AAI3D,6DAA6D;AAC7D,OAAO,EAAE,cAAc,EAAE,MAAM,wBAAwB,CAAC;AAoCxD,wEAAwE;AACxE,KAAK,UAAU,OAAO,CAAC,IAAc,EAAE,OAAyB;IAC9D,MAAM,EAAE,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC;IAC1B,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACX,MAAM,IAAI,qBAAqB,CAC7B,0DAA0D;cACxD,GAAG,EAAE,CAAC,UAAU,CAAC,MAAM,iBAAiB;cACxC,EAAE,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,QAAQ,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CACnE,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GAAG,OAAO,EAAE,MAAM,IAAI,MAAM,eAAe,EAAE,CAAC;IAC1D,MAAM,QAAQ,GAAG,EAAE,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC;IACvE,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC,iBAA8B,EAAE,WAAW,EAAE,EAAE,CAAC,WAAW,EAAE,QAAQ,EAAE,CAAC;AACnG,CAAC;AAED;;mDAEmD;AACnD,KAAK,UAAU,4BAA4B,CACzC,IAAc,EACd,MAAuB,EACvB,MAAoB;IAEpB,MAAM,IAAI,GAAG,iBAAiB,CAAC,IAAI,CAAC,CAAC;IACrC,IAAI,CAAC,IAAI;QAAE,OAAO,EAAE,CAAC;IACrB,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM,qBAAqB,CAC7C,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,MAAM,EAAE,SAAS,EAAE,EAAE,MAAM,EAAE,CACxD,CAAC;IACF,OAAO,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;AAClC,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,iBAAiB,CACrC,IAAc,EACd,MAAuB,EACvB,OAAyB;IAEzB,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,GAAG,MAAM,OAAO,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAC5E,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,qBAAqB,GAAG,MAAM,4BAA4B,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACvF,OAAO;QACL,KAAK,EAAE,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;QAC9B,GAAG;QACH,WAAW;QACX,QAAQ,EAAE,CAAC,GAAG,QAAQ,EAAE,GAAG,qBAAqB,CAAC;KAClD,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,oBAAoB,CACxC,IAAc,EACd,MAAuB,EACvB,OAAyB;IAEzB,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,GAAG,MAAM,OAAO,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAC5E,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,qBAAqB,GAAG,MAAM,4BAA4B,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACvF,OAAO;QACL,KAAK,EAAE,MAAM;QACb,GAAG;QACH,WAAW;QACX,QAAQ,EAAE,CAAC,GAAG,QAAQ,EAAE,GAAG,qBAAqB,CAAC;QACjD,OAAO,EAAE,GAAG,EAAE,GAAG,MAAM,CAAC,OAAO,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;KAC7C,CAAC;AACJ,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/numerical/index.ts"],"names":[],"mappings":"AAQA,OAAO,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAGvD,OAAO,EAAE,SAAS,EAAE,MAAM,eAAe,CAAC;AAC1C,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,qBAAqB,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAI/E,cAAc;AACd,OAAO,EAAE,kBAAkB,EAAE,qBAAqB,EAAE,MAAM,oBAAoB,CAAC;AAK/E,cAAc;AACd,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACxE,cAAc;AACd,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD;;0EAE0E;AAC1E,OAAO,EAAE,0BAA0B,EAAE,MAAM,0BAA0B,CAAC;AACtE;;;;GAIG;AACH,OAAO,EAAE,qBAAqB,EAAE,MAAM,qBAAqB,CAAC;AAC5D,mEAAmE;AACnE,OAAO,EACL,iCAAiC,GAClC,MAAM,wBAAwB,CAAC;AAOhC,kDAAkD;AAClD,OAAO,EAAE,oBAAoB,EAAE,MAAM,qBAAqB,CAAC;AAI3D,6DAA6D;AAC7D,OAAO,EAAE,cAAc,EAAE,MAAM,wBAAwB,CAAC;AAoCxD,wEAAwE;AACxE,KAAK,UAAU,OAAO,CAAC,IAAc,EAAE,OAAyB;IAC9D,MAAM,EAAE,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC;IAC1B,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACX,MAAM,IAAI,qBAAqB,CAC7B,0DAA0D;cACxD,GAAG,EAAE,CAAC,UAAU,CAAC,MAAM,iBAAiB;cACxC,EAAE,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC,QAAQ,KAAK,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CACnE,CAAC;IACJ,CAAC;IACD,MAAM,MAAM,GAAG,OAAO,EAAE,MAAM,IAAI,MAAM,eAAe,EAAE,CAAC;IAC1D,MAAM,QAAQ,GAAG,EAAE,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC;IACvE,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC,iBAA8B,EAAE,WAAW,EAAE,EAAE,CAAC,WAAW,EAAE,QAAQ,EAAE,CAAC;AACnG,CAAC;AAED;;mDAEmD;AACnD,KAAK,UAAU,4BAA4B,CACzC,IAAc,EACd,MAAuB,EACvB,MAAoB;IAEpB,MAAM,IAAI,GAAG,iBAAiB,CAAC,IAAI,CAAC,CAAC;IACrC,IAAI,CAAC,IAAI;QAAE,OAAO,EAAE,CAAC;IACrB,MAAM,EAAE,OAAO,EAAE,GAAG,MAAM,qBAAqB,CAC7C,IAAI,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,MAAM,EAAE,SAAS,EAAE,EAAE,MAAM,EAAE,CACxD,CAAC;IACF,OAAO,OAAO,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;AAClC,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,iBAAiB,CACrC,IAAc,EACd,MAAuB,EACvB,OAAyB;IAEzB,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,GAAG,MAAM,OAAO,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAC5E,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,qBAAqB,GAAG,MAAM,4BAA4B,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACvF,OAAO;QACL,KAAK,EAAE,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC;QAC9B,GAAG;QACH,WAAW;QACX,QAAQ,EAAE,CAAC,GAAG,QAAQ,EAAE,GAAG,qBAAqB,CAAC;KAClD,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,MAAM,CAAC,KAAK,UAAU,oBAAoB,CACxC,IAAc,EACd,MAAuB,EACvB,OAAyB;IAEzB,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,WAAW,EAAE,QAAQ,EAAE,GAAG,MAAM,OAAO,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAC5E,MAAM,MAAM,GAAG,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAC/C,MAAM,qBAAqB,GAAG,MAAM,4BAA4B,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACvF,OAAO;QACL,KAAK,EAAE,MAAM;QACb,GAAG;QACH,WAAW;QACX,QAAQ,EAAE,CAAC,GAAG,QAAQ,EAAE,GAAG,qBAAqB,CAAC;QACjD,OAAO,EAAE,GAAG,EAAE,GAAG,MAAM,CAAC,OAAO,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;KAC7C,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Public surface of the UPT numerical-contraction backend (v0.3.5).\n * See docs/planning/v0.3.5-Design.md §5.\n *\n * @module numerical\n */\nimport type { Dimension } from '../dimensional/types.js';\nimport type { ExprNode, Violation } from '../dimensional/validator.js';\nimport { validate } from '../dimensional/validator.js';\nimport type { EngineTensor, TensorEngine } from './tensor-engine.js';\nimport type { NumericalInputs, NestedArray } from './types.js';\nimport { lowerNode } from './lowering.js';\nimport { getActiveEngine } from './engine-registry.js';\nimport { NumericalBackendError } from './errors.js';\nimport { evaluateMetricInverse, scanForMetricPair } from './metric-inverse.js';\n\n/** @public */\nexport type { TensorEngine, EngineTensor, EinsumSpec, ForwardGradResult, ReverseGradResult } from './tensor-engine.js';\n/** @public */\nexport { hasAutogradSupport, EngineCapabilityError } from './tensor-engine.js';\n/** @public */\nexport type { NumericalInputs, NestedArray } from './types.js';\n/** @public — part of the `NumericalInputs.grids` public contract. */\nexport type { GridField } from './grid-field.js';\n/** @public */\nexport { getActiveEngine, setActiveEngine } from './engine-registry.js';\n/** @public */\nexport { NumericalBackendError } from './errors.js';\n/** @public — re-exported from dimensional/errors to keep numerical/index.ts\n *  as the single public API surface, without creating a dimensional→numerical\n *  import cycle. See v0.4.0-Implementation-Plan Task 13 for rationale. */\nexport { DuplicateCoordinateWarning } from '../dimensional/errors.js';\n/**\n * Evaluate a metric inverse from a detected pair.\n *\n * @public\n */\nexport { evaluateMetricInverse } from './metric-inverse.js';\n/** Integrate the covariant eikonal on a numerical grid. @public */\nexport {\n  evaluateCovariantEikonalNumerical,\n} from './covariant-eikonal.js';\n/** Inputs and the result of the covariant eikonal integrator. */\nexport type {\n  CovariantEikonalInputs,\n  CovariantEikonalResult,\n} from './covariant-eikonal.js';\n\n/** @public — v0.5.0 GL4 symplectic integrator. */\nexport { integrateGeodesicGL4 } from './gl4-integrator.js';\n/** @public — v0.5.0 GL4 types. */\nexport type { GL4State, GL4Snapshot, GL4Options } from './gl4-integrator.js';\n\n/** @public — v0.5.0 perihelion finder (Task 4, Phase 1b). */\nexport { findPerihelion } from './perihelion-finder.js';\n/** @public — v0.5.0 perihelion-finder types. */\nexport type { PerihelionResult, FindPerihelionOptions } from './perihelion-finder.js';\n\n/**\n * Plain-JS result of `evaluateNumerical`.\n * @public\n */\nexport interface NumericalResult {\n  readonly value: NestedArray;\n  readonly dim: Dimension;\n  readonly freeIndices: ReadonlyMap<string, { upper: number; lower: number }>;\n  readonly warnings: ReadonlyArray<Violation>;\n}\n\n/**\n * Result of `evaluateNumericalRaw` — carries a live `EngineTensor` for\n * chaining workloads; the caller must `dispose()` it.\n * @public\n */\nexport interface NumericalRawResult {\n  readonly value: EngineTensor;\n  readonly dim: Dimension;\n  readonly freeIndices: ReadonlyMap<string, { upper: number; lower: number }>;\n  readonly warnings: ReadonlyArray<Violation>;\n  dispose(): void;\n}\n\n/**\n * Per-call options for the `evaluateNumerical*` entry points.\n * @public\n */\nexport interface EvaluateOptions {\n  readonly engine?: TensorEngine;\n}\n\n/** Validate the AST and return (engine, validation result) or throw. */\nasync function prepare(node: ExprNode, options?: EvaluateOptions) {\n  const vr = validate(node);\n  if (!vr.ok) {\n    throw new NumericalBackendError(\n      `evaluateNumerical: the AST is not dimensionally valid — `\n      + `${vr.violations.length} violation(s): `\n      + vr.violations.map((v) => `[${v.location}] ${v.note}`).join('; '),\n    );\n  }\n  const engine = options?.engine ?? await getActiveEngine();\n  const warnings = vr.violations.filter((v) => v.severity === 'warning');\n  return { engine, dim: vr.inferredDimension as Dimension, freeIndices: vr.freeIndices, warnings };\n}\n\n/** If the AST contains an identifiable lower/upper metric pair, run the\n *  numerical inverse-metric check and return any warning (deduplicated —\n *  scanForMetricPair returns at most one pair). */\nasync function collectInverseMetricWarnings(\n  node: ExprNode,\n  inputs: NumericalInputs,\n  engine: TensorEngine,\n): Promise<Violation[]> {\n  const pair = scanForMetricPair(node);\n  if (!pair) return [];\n  const { warning } = await evaluateMetricInverse(\n    pair.gUpper, pair.gLower, inputs, undefined, { engine },\n  );\n  return warning ? [warning] : [];\n}\n\n/**\n * Evaluate a validated AST to plain JS.\n * @public\n */\nexport async function evaluateNumerical(\n  node: ExprNode,\n  inputs: NumericalInputs,\n  options?: EvaluateOptions,\n): Promise<NumericalResult> {\n  const { engine, dim, freeIndices, warnings } = await prepare(node, options);\n  const tensor = lowerNode(node, inputs, engine);\n  const inverseMetricWarnings = await collectInverseMetricWarnings(node, inputs, engine);\n  return {\n    value: engine.toNested(tensor),\n    dim,\n    freeIndices,\n    warnings: [...warnings, ...inverseMetricWarnings],\n  };\n}\n\n/**\n * Evaluate to a live EngineTensor for chaining workloads.\n * @public\n */\nexport async function evaluateNumericalRaw(\n  node: ExprNode,\n  inputs: NumericalInputs,\n  options?: EvaluateOptions,\n): Promise<NumericalRawResult> {\n  const { engine, dim, freeIndices, warnings } = await prepare(node, options);\n  const tensor = lowerNode(node, inputs, engine);\n  const inverseMetricWarnings = await collectInverseMetricWarnings(node, inputs, engine);\n  return {\n    value: tensor,\n    dim,\n    freeIndices,\n    warnings: [...warnings, ...inverseMetricWarnings],\n    dispose: () => { engine.dispose?.(tensor); },\n  };\n}\n"]}

package/dist/numerical/input-validation.d.ts

@@ -1,12 +1,9 @@
 /**
  * Runtime input validation for numeric evaluators — `validateFiniteInputs`.
  *
- * A generic, dependency-free finite-and-in-range check shared by the numerical
- * evaluators (`klein-gordon.ts`) and the 42 bridge-equation evaluators
- * (`bridges/equations/be-*.ts`). It lived in `bridges/equations/_be-helpers.ts`,
- * which forced the lower numerical layer to import upward from `bridges`; moving
- * it here (a leaf module with no imports) removes that dependency while
- * `_be-helpers.ts` re-exports it so every bridge evaluator is unchanged.
+ * A generic, dependency-free finite-and-in-range check for a numeric
+ * evaluator's input record (`klein-gordon.ts` reads it). It is a leaf module
+ * with no imports, so any layer may read it without an upward edge.
  *
  * @module numerical/input-validation
  */

package/dist/numerical/input-validation.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"input-validation.d.ts","sourceRoot":"","sources":["../../src/numerical/input-validation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH;;;;;;;;;;GAUG;AACH,MAAM,WAAW,SAAS;IACxB,yDAAyD;IACzD,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,mDAAmD;IACnD,QAAQ,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC;IACtB,kDAAkD;IAClD,QAAQ,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC;IACtB,oDAAoD;IACpD,QAAQ,CAAC,UAAU,CAAC,EAAE,OAAO,CAAC;IAC9B,oDAAoD;IACpD,QAAQ,CAAC,UAAU,CAAC,EAAE,OAAO,CAAC;IAC9B,wEAAwE;IACxE,QAAQ,CAAC,SAAS,CAAC,EAAE,OAAO,CAAC;IAC7B;;;;;OAKG;IACH,QAAQ,CAAC,WAAW,CAAC,EAAE,MAAM,CAAC;CAC/B;AAgCD;;;;;;;;;GASG;AACH,wBAAgB,oBAAoB,CAAC,CAAC,SAAS,MAAM,EACnD,KAAK,EAAE,CAAC,EACR,UAAU,EAAE,aAAa,CAAC,SAAS,CAAC,EACpC,aAAa,EAAE,MAAM,GACpB,IAAI,CAgCN"}
+{"version":3,"file":"input-validation.d.ts","sourceRoot":"","sources":["../../src/numerical/input-validation.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH;;;;;;;;;;GAUG;AACH,MAAM,WAAW,SAAS;IACxB,yDAAyD;IACzD,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,mDAAmD;IACnD,QAAQ,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC;IACtB,kDAAkD;IAClD,QAAQ,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC;IACtB,oDAAoD;IACpD,QAAQ,CAAC,UAAU,CAAC,EAAE,OAAO,CAAC;IAC9B,oDAAoD;IACpD,QAAQ,CAAC,UAAU,CAAC,EAAE,OAAO,CAAC;IAC9B,wEAAwE;IACxE,QAAQ,CAAC,SAAS,CAAC,EAAE,OAAO,CAAC;IAC7B;;;;;OAKG;IACH,QAAQ,CAAC,WAAW,CAAC,EAAE,MAAM,CAAC;CAC/B;AAgCD;;;;;;;;;GASG;AACH,wBAAgB,oBAAoB,CAAC,CAAC,SAAS,MAAM,EACnD,KAAK,EAAE,CAAC,EACR,UAAU,EAAE,aAAa,CAAC,SAAS,CAAC,EACpC,aAAa,EAAE,MAAM,GACpB,IAAI,CA+BN"}