universal-physics-tensor 9.0.0 → 11.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.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"geodesic-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/geodesic-integrator.ts"],"names":[],"mappings":"AA2BA;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC;;;;;;;;;;;;OAYG;IACH,aAAa,EAAE,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EAAE,GAAG,CAAC,EAAE,YAAY,KAAK,YAAY,CAAC;IAE9E,iDAAiD;IACjD,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE9C;;;;;;;;OAQG;IACH,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE9C,yDAAyD;IACzD,QAAQ,EAAE,MAAM,CAAC;IAEjB,uDAAuD;IACvD,MAAM,EAAE,MAAM,CAAC;IAEf,8DAA8D;IAC9D,KAAK,EAAE,MAAM,CAAC;IAEd;;;;;;OAMG;IACH,eAAe,CAAC,EAAE,MAAM,CAAC;CAC1B;AAED;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE3D,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE3D;;;OAGG;IACH,QAAQ,CAAC,UAAU,EAAE,aAAa,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;CAC/E;AAuFD;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,wBAAwB,GAC/B,wBAAwB,CAuD1B"}
+{"version":3,"file":"geodesic-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/geodesic-integrator.ts"],"names":[],"mappings":"AA4BA;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC;;;;;;;;;;;;OAYG;IACH,aAAa,EAAE,CAAC,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,EAAE,GAAG,CAAC,EAAE,YAAY,KAAK,YAAY,CAAC;IAE9E,iDAAiD;IACjD,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE9C;;;;;;;;OAQG;IACH,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE9C,yDAAyD;IACzD,QAAQ,EAAE,MAAM,CAAC;IAEjB,uDAAuD;IACvD,MAAM,EAAE,MAAM,CAAC;IAEf,8DAA8D;IAC9D,KAAK,EAAE,MAAM,CAAC;IAEd;;;;;;OAMG;IACH,eAAe,CAAC,EAAE,MAAM,CAAC;CAC1B;AAED;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE3D,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IAE3D;;;OAGG;IACH,QAAQ,CAAC,UAAU,EAAE,aAAa,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;CAC/E;AAuFD;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,wBAAwB,GAC/B,wBAAwB,CA8D1B"}

package/dist/numerical/geodesic-integrator.js

@@ -14,16 +14,17 @@
  * The Christoffel symbol Γ^μ_{νρ} is supplied as a plain-JS closure that
  * maps a 4-coordinate array to a [4][4][4] tensor.  No TensorEngine
  * dependency — the integrator is self-contained. The fixed steps are
- * MathTS `solveODESystem` with `dt` set, which is that library's classical
- * RK4. Adaptive `solveODE` is not this step.
+ * `rk4Trajectory` (MathTS `solveODESystem` with `dt` set, that library's
+ * classical RK4), which takes exactly `steps` steps. Adaptive `solveODE` is
+ * not this step.
  *
  * Validated against the cycloid parametric form for radial Schwarzschild
  * infall to ±1e-6 (Task 14 test, η = 0.5).
  *
  * @module numerical/geodesic-integrator
  */
-import { solveODESystem } from '@danielsimonjr/mathts-functions';
 import { NumericalBackendError } from './errors.js';
+import { rk4Trajectory } from './rk4-trajectory.js';
 /**
  * Evaluate the Christoffel acceleration −Γ^μ_{νρ} v^ν v^ρ at position x
  * with 4-velocity v.  Returns (dx/dτ, dv/dτ) = (v, accel).
@@ -149,28 +150,35 @@ export function integrateGeodesic(inputs) {
     if (!Number.isInteger(steps) || steps <= 0) {
         throw new NumericalBackendError(`integrateGeodesic: steps must be a positive integer, got ${steps}`);
     }
-    const h = (tauEnd - tauStart) / steps;
+    // A NaN or infinite end made every step NaN and the result the start state; an
+    // end equal to the start is a zero step, which is not an integration either.
+    if (!Number.isFinite(tauStart)) {
+        throw new NumericalBackendError(`integrateGeodesic: tauStart must be finite, got ${tauStart}`);
+    }
+    if (!Number.isFinite(tauEnd)) {
+        throw new NumericalBackendError(`integrateGeodesic: tauEnd must be finite, got ${tauEnd}`);
+    }
+    if (tauEnd === tauStart) {
+        throw new NumericalBackendError(`integrateGeodesic: tauEnd must differ from tauStart, both are ${tauStart}`);
+    }
     const y0 = [x0[0], x0[1], x0[2], x0[3], v0[0], v0[1], v0[2], v0[3]];
     // One scratch buffer reused across all Christoffel evaluations (4 per step).
     // Each geodesicRHS consumes its Γ before the next call, so reuse is safe and
     // avoids ~4·steps Float64Array(64) allocations per integration.
     const scratch = new Float64Array(64);
-    const sol = solveODESystem((_tau, y) => {
+    const states = rk4Trajectory((_tau, y) => {
         const x = [y[0], y[1], y[2], y[3]];
         const v = [y[4], y[5], y[6], y[7]];
         const s = geodesicRHS(christoffelFn, x, v, scratch);
         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]];
-    }, y0, [tauStart, tauEnd], { dt: h });
+    }, y0, tauStart, tauEnd, steps);
     const asVec4 = (row) => [row[0], row[1], row[2], row[3]];
     const sampleEvery = Math.max(1, Math.floor(steps / 100));
-    const trajectory = [asVec4(sol.y[0] ?? y0)];
-    for (let i = 0; i < steps; i++) {
-        if ((i + 1) % sampleEvery === 0) {
-            const row = sol.y[i + 1] ?? sol.y[sol.y.length - 1] ?? y0;
-            trajectory.push(asVec4(row));
-        }
-    }
-    const last = sol.y[sol.y.length - 1] ?? y0;
+    // `rk4Trajectory` returns exactly steps + 1 states, so state i is the one after i steps.
+    const trajectory = [asVec4(states[0] ?? y0)];
+    for (let i = sampleEvery; i <= steps; i += sampleEvery)
+        trajectory.push(asVec4(states[i] ?? y0));
+    const last = states[steps] ?? y0;
     return {
         xFinal: asVec4(last),
         vFinal: [last[4], last[5], last[6], last[7]],

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;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,aAAa,EAAE,MAAM,qBAAqB,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;IACD,+EAA+E;IAC/E,6EAA6E;IAC7E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,EAAE,CAAC;QAC/B,MAAM,IAAI,qBAAqB,CAAC,mDAAmD,QAAQ,EAAE,CAAC,CAAC;IACjG,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC7B,MAAM,IAAI,qBAAqB,CAAC,iDAAiD,MAAM,EAAE,CAAC,CAAC;IAC7F,CAAC;IACD,IAAI,MAAM,KAAK,QAAQ,EAAE,CAAC;QACxB,MAAM,IAAI,qBAAqB,CAAC,iEAAiE,QAAQ,EAAE,CAAC,CAAC;IAC/G,CAAC;IAED,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,MAAM,GAAG,aAAa,CAC1B,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,QAAQ,EACR,MAAM,EACN,KAAK,CACN,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,yFAAyF;IACzF,MAAM,UAAU,GAAW,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACrD,KAAK,IAAI,CAAC,GAAG,WAAW,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,IAAI,WAAW;QAAE,UAAU,CAAC,IAAI,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACjG,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC;IAEjC,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 * `rk4Trajectory` (MathTS `solveODESystem` with `dt` set, that library's\n * classical RK4), which takes exactly `steps` steps. Adaptive `solveODE` is\n * 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 { NumericalBackendError } from './errors.js';\nimport { rk4Trajectory } from './rk4-trajectory.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  // A NaN or infinite end made every step NaN and the result the start state; an\n  // end equal to the start is a zero step, which is not an integration either.\n  if (!Number.isFinite(tauStart)) {\n    throw new NumericalBackendError(`integrateGeodesic: tauStart must be finite, got ${tauStart}`);\n  }\n  if (!Number.isFinite(tauEnd)) {\n    throw new NumericalBackendError(`integrateGeodesic: tauEnd must be finite, got ${tauEnd}`);\n  }\n  if (tauEnd === tauStart) {\n    throw new NumericalBackendError(`integrateGeodesic: tauEnd must differ from tauStart, both are ${tauStart}`);\n  }\n\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 states = rk4Trajectory(\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,\n    tauEnd,\n    steps,\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  // `rk4Trajectory` returns exactly steps + 1 states, so state i is the one after i steps.\n  const trajectory: Vec4[] = [asVec4(states[0] ?? y0)];\n  for (let i = sampleEvery; i <= steps; i += sampleEvery) trajectory.push(asVec4(states[i] ?? y0));\n  const last = states[steps] ?? 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.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"geometrized.d.ts","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAEH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAGpD;;;;;GAKG;AACH,qBAAa,8BAA+B,SAAQ,QAAQ;IAC1D,YAAY,OAAO,EAAE,MAAM,EAI1B;CACF;AAED;;;;;;;;GAQG;AACH,wBAAgB,iBAAiB,CAAC,GAAG,EAAE,SAAS,GAAG,MAAM,CASxD;AAED,4CAA4C;AAC5C,wBAAgB,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAErE;AAED,4CAA4C;AAC5C,wBAAgB,eAAe,CAAC,SAAS,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAEzE"}
+{"version":3,"file":"geometrized.d.ts","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAEH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAGpD;;;;;GAKG;AACH,qBAAa,8BAA+B,SAAQ,QAAQ;IAC1D,YAAY,OAAO,EAAE,MAAM,EAG1B;CACF;AAED;;;;;;;;GAQG;AACH,wBAAgB,iBAAiB,CAAC,GAAG,EAAE,SAAS,GAAG,MAAM,CASxD;AAED,4CAA4C;AAC5C,wBAAgB,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAErE;AAED,4CAA4C;AAC5C,wBAAgB,eAAe,CAAC,SAAS,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAEzE"}

package/dist/numerical/geometrized.js

@@ -35,7 +35,6 @@ export class NonGeometrizableDimensionError extends UPTError {
     constructor(message) {
         super(message);
         this.name = 'NonGeometrizableDimensionError';
-        Object.setPrototypeOf(this, NonGeometrizableDimensionError.prototype);
     }
 }
 /**

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;IAC/C,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  }\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++;
                 }