@danielsimonjr/mathts-functions 0.3.0 → 0.4.0

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
@@ -83,7 +83,8 @@ export declare const matrixFromColumns: unknown;
83
83
  export declare const matrixFromRows: unknown;
84
84
  export declare const count: unknown;
85
85
  export declare const trace: unknown;
86
- export declare const det: number | import("complex.js").Complex | BigNumber;
86
+ declare const factoryDet: number | import("complex.js").Complex | BigNumber;
87
+ export declare const det: typeof factoryDet;
87
88
  export declare const reshape: unknown;
88
89
  export declare const factory_equal: unknown;
89
90
  export declare const isInteger: unknown;
@@ -91,7 +92,8 @@ export declare const concat: unknown;
91
92
  export declare const mapSlices: import("../types.js").TypedFunction;
92
93
  export declare const resize: unknown;
93
94
  export declare const subset: unknown;
94
- export declare const inv: number | import("complex.js").Complex | BigNumber;
95
+ declare const factoryInv: number | import("complex.js").Complex | BigNumber;
96
+ export declare const inv: typeof factoryInv;
95
97
  export declare const factory_cbrt: unknown;
96
98
  export declare const nthRoots: import("../types.js").TypedFunction;
97
99
  export declare const factory_round: unknown;
@@ -1 +1 @@
1
- {"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/factories/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAwIH,eAAO,MAAM,SAAS,SAAyE,CAAC;AAChG,eAAO,MAAM,cAAc,SAAmF,CAAC;AAC/G,eAAO,MAAM,cAAc,SAAmF,CAAC;AAe/G,eAAO,MAAM,MAAM,SAAmE,CAAC;AACvF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAC1F,eAAO,MAAM,OAAO,SAAqE,CAAC;AAC1F,eAAO,MAAM,iBAAiB,SAAyF,CAAC;AACxH,eAAO,MAAM,GAAG,SAA6D,CAAC;AAC9E,eAAO,MAAM,IAAI,SAA+D,CAAC;AACjF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAc1F,eAAO,MAAM,GAAG,SAA6D,CAAC;AAO9E,eAAO,MAAM,KAAK,qCAAiE,CAAC;AACpF,eAAO,MAAM,KAAK,qCAAiE,CAAC;AACpF,eAAO,MAAM,KAAK,qCAAiE,CAAC;AACpF,eAAO,MAAM,IAAI,qCAA+D,CAAC;AACjF,eAAO,MAAM,IAAI,qCAA+D,CAAC;AACjF,eAAO,MAAM,IAAI,qCAA+D,CAAC;AAGjF,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,KAAK,SAAiE,CAAC;AACpF,eAAO,MAAM,SAAS,qCAAyE,CAAC;AAChG,eAAO,MAAM,KAAK,SAAiE,CAAC;AACpF,eAAO,MAAM,UAAU,SAA2E,CAAC;AACnG,eAAO,MAAM,SAAS,qCAAyE,CAAC;AAChG,eAAO,MAAM,UAAU,SAA2E,CAAC;AACnG,eAAO,MAAM,OAAO,SAAqE,CAAC;AAC1F,eAAO,MAAM,OAAO,mEAAqE,CAAC;AAC1F,eAAO,MAAM,MAAM,SAAmE,CAAC;AAMvF,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,aAAa,SAAiE,CAAC;AAC5F,eAAO,MAAM,aAAa,SAAiE,CAAC;AAC5F,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,cAAc,SAAmE,CAAC;AAC/F,eAAO,MAAM,kBAAkB,SAA2E,CAAC;AAC3G,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,aAAa,qCAAiE,CAAC;AAC5F,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,aAAa,qCAAiE,CAAC;AAC5F,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,aAAa,qCAAiE,CAAC;AAC5F,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,WAAW,qCAA6D,CAAC;AACtF,eAAO,MAAM,WAAW,qCAA6D,CAAC;AACtF,eAAO,MAAM,WAAW,qCAA6D,CAAC;AACtF,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AA2EzF,eAAO,MAAM,qBAAqB,oEAAiG,CAAC;AAIpI,eAAO,MAAM,YAAY,SAA+E,CAAC;AAMzG,eAAO,MAAM,IAAI,SAA+D,CAAC;AACjF,eAAO,MAAM,IAAI,SAA+D,CAAC;AAQjF,eAAO,MAAM,eAAe,SAAqF,CAAC;AAClH,eAAO,MAAM,QAAQ,qCAAuE,CAAC;AAC7F,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,iBAAiB,SAAyE,CAAC;AACxG,eAAO,MAAM,WAAW,SAA6D,CAAC;AAwCtF,eAAO,MAAM,iBAAiB,SAAyE,CAAC;AAGxG,eAAO,MAAM,kBAAkB,SAA2E,CAAC;AAG3G,eAAO,MAAM,QAAQ,SAAuE,CAAC;AAG7F,eAAO,MAAM,KAAK,SAAiE,CAAC;AAGpF,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,kBAAkB,qCAA2F,CAAC;AAG3H,eAAO,MAAM,iBAAiB,SAAyF,CAAC;AAGxH,eAAO,MAAM,cAAc,SAAmF,CAAC;AAG/G,eAAO,MAAM,KAAK,SAAiE,CAAC;AAQpF,eAAO,MAAM,KAAK,SAAiE,CAAC;AAQpF,eAAO,MAAM,GAAG,mDAA6D,CAAC;AAK9E,eAAO,MAAM,OAAO,SAAqE,CAAC;AAuD1F,eAAO,MAAM,aAAa,SAAiE,CAAC;AAI5F,eAAO,MAAM,SAAS,SAAyE,CAAC;AAIhG,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,SAAS,qCAAyE,CAAC;AAGhG,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,GAAG,mDAA6D,CAAC;AAI9E,eAAO,MAAM,YAAY,SAA+D,CAAC;AAGzF,eAAO,MAAM,QAAQ,qCAAuE,CAAC;AAG7F,eAAO,MAAM,aAAa,SAAiE,CAAC;AAG5F,eAAO,MAAM,YAAY,qCAA+D,CAAC;AAGzF,eAAO,MAAM,WAAW,SAA6D,CAAC;AAItF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAQ1F,eAAO,MAAM,MAAM,KAAmE,CAAC;AAIvF,eAAO,MAAM,cAAc,SAAmE,CAAC;AAG/F,eAAO,MAAM,WAAW,SAA6D,CAAC;AAItF,eAAO,MAAM,QAAQ,uJAAuE,CAAC;AAG7F,eAAO,MAAM,SAAS,qKAAyE,CAAC;AAGhG,eAAO,MAAM,KAAK,qGAAiE,CAAC;AAGpF,eAAO,MAAM,KAAK,wJAAiE,CAAC;AAIpF,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,SAAS,wDAAyE,CAAC;AAGhG,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,SAAS,wDAAyE,CAAC;AAqDhG,eAAO,MAAM,eAAe,SAAqE,CAAC;AAWlG,eAAO,MAAM,cAAc,SAAmE,CAAC;AAG/F,eAAO,MAAM,gBAAgB,SAAuE,CAAC;AAGrG,eAAO,MAAM,eAAe,SAAqE,CAAC;AAGlG,eAAO,MAAM,iBAAiB,SAAyE,CAAC;AAQxG,eAAO,MAAM,SAAS,qCAAyE,CAAC;AAGhG,eAAO,MAAM,WAAW,qCAA6E,CAAC;AAGtG,eAAO,MAAM,WAAW,qCAA6D,CAAC;AAGtF,eAAO,MAAM,WAAW,qCAA6D,CAAC;AAGtF,eAAO,MAAM,aAAa,SAAiE,CAAC;AAG5F,eAAO,MAAM,WAAW,SAA6D,CAAC;AAGtF,eAAO,MAAM,eAAe,qCAAqE,CAAC;AAGlG,eAAO,MAAM,WAAW,SAA6D,CAAC;AAGtF,eAAO,MAAM,YAAY,SAA+D,CAAC;AAGzF,eAAO,MAAM,aAAa,SAAiE,CAAC;AA8B5F,eAAO,MAAM,IAAI,0FAA+D,CAAC;AAIjF,eAAO,MAAM,aAAa,SAAiE,CAAC;AAI5F,eAAO,MAAM,EAAE,SAA2D,CAAC;AAoD3E,eAAO,MAAM,MAAM,qCAAmE,CAAC;AAGvF,eAAO,MAAM,WAAW,SAA6D,CAAC;AAGtF,eAAO,MAAM,MAAM,qCAAmE,CAAC;AAWvF,eAAO,MAAM,WAAW,SAA6E,CAAC;AAItG,eAAO,MAAM,eAAe,SAAqF,CAAC;AAQlH,eAAO,MAAM,EAAE,SAA2D,CAAC;AAG3E,eAAO,MAAM,KAAK,SAAiE,CAAC;AAGpF,eAAO,MAAM,QAAQ,SAAuE,CAAC;AAI7F,eAAO,MAAM,KAAK,SAAiE,CAAC;AAIpF,eAAO,MAAM,WAAW,SAA6D,CAAC;AAGtF,eAAO,MAAM,WAAW,SAA6D,CAAC;AAItF,eAAO,MAAM,aAAa,qCAAiE,CAAC;AAW5F,eAAO,MAAM,SAAS,SAAyE,CAAC;AAOhG,eAAO,MAAM,IAAI,qCAA+D,CAAC;AAUjF,eAAO,MAAM,UAAU,SAA2E,CAAC;AAanG,eAAO,MAAM,WAAW,SAA6E,CAAC;AAgKtG,eAAO,MAAM,cAAc;;;;;;;CAAmF,CAAC;AAG/G,eAAO,MAAM,SAAS,SAAyE,CAAC;AAGhG,eAAO,MAAM,OAAO,SAAqE,CAAC;AAG1F,eAAO,MAAM,gBAAgB;;;;yCAAuF,CAAC;AAGrH,eAAO,MAAM,YAAY;;;;;;;;;;;;;;;;;;;YAn4BtB,CAAC;;;;;YAYD,CAAA;UAAgB,CAAC;;;;;YAcpB,CAAA;;;;;YAW8D,CAAC;;;;;;;;;;;;;CA81ByC,CAAC;AAGzG,eAAO,MAAM,SAAS,SAAyE,CAAC;AAGhG,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,KAAK,mDAA2E,CAAC;AA4B9F,eAAO,MAAM,gBAAgB,SAAuE,CAAC;AAGrG,eAAO,MAAM,cAAc,SAAmE,CAAC;AAO/F,eAAO,MAAM,IAAI,SAA+D,CAAC;AAMjF,eAAO,MAAM,KAAK,qCAAiE,CAAC;AAGpF,eAAO,MAAM,kBAAkB,SAA2E,CAAC;AAG3G,eAAO,MAAM,IAAI,SAA+E,CAAC;AA2CjG,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,KAAK,SAAiE,CAAC;AAGpF,eAAO,MAAM,IAAI,qCAA+D,CAAC;AAGjF,eAAO,MAAM,KAAK,2FAAiE,CAAC;AAGpF,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,MAAM,0LAAmE,CAAC;AAGvF,eAAO,MAAM,IAAI,sMAA+D,CAAC;AAGjF,eAAO,MAAM,SAAS,mNAAyE,CAAC;AAGhG,eAAO,MAAM,SAAS,SAAyE,CAAC;AAGhG,eAAO,MAAM,YAAY,SAA+D,CAAC;AAGzF,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,gBAAgB,SAAuE,CAAC;AAGrG,eAAO,MAAM,WAAW,SAA6E,CAAC;AAGtG,eAAO,MAAM,YAAY,SAA+E,CAAC;AAOzG,eAAO,MAAM,KAAK,KAAiE,CAAC;AA+BpF,eAAO,MAAM,YAAY,SAA+E,CAAC;AAGzG,eAAO,MAAM,cAAc,SAAmF,CAAC;AAG/G,eAAO,MAAM,QAAQ,KAAuE,CAAC;AAG7F,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAe1F,eAAO,MAAM,IAAI,mDAA+D,CAAC;AAGjF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAG1F,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,WAAW,SAA6D,CAAC;AAOtF,eAAO,MAAM,QAAQ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CAAuE,CAAC;AAa7F,eAAO,MAAM,UAAU;;;;;CAA2E,CAAC;AAGnG,eAAO,MAAM,YAAY,qCAA+D,CAAC;AAGzF,eAAO,MAAM,WAAW,SAA6E,CAAC;AAOtG,eAAO,MAAM,aAAa,SAAiF,CAAC;AAU5G,eAAO,MAAM,cAAc,qCAAmF,CAAC;AAG/G,eAAO,MAAM,KAAK,SAAiE,CAAC;AAUpF,eAAO,MAAM,MAAM,qCAAmE,CAAC;AAGvF,eAAO,MAAM,SAAS,SAAyE,CAAC;AAShG,eAAO,MAAM,IAAI,SAA+D,CAAC;AAkEjF,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,mBAAmB,0DAA6F,CAAC;AAC9H,eAAO,MAAM,cAAc,0DAAmF,CAAC;AAC/G,eAAO,MAAM,qBAAqB,0DAAiG,CAAC;AACpI,eAAO,MAAM,gBAAgB,0DAAuF,CAAC;AACrH,eAAO,MAAM,gBAAgB,0DAAuF,CAAC;AACrH,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,OAAO,0DAAqE,CAAC;AAC1F,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,gBAAgB,0DAAuF,CAAC;AACrH,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,kBAAkB,0DAA2F,CAAC;AAC3H,eAAO,MAAM,yBAAyB,0DAAyG,CAAC;AAChJ,eAAO,MAAM,mBAAmB,0DAA6F,CAAC;AAC9H,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,QAAQ,0DAAuE,CAAC;AAC7F,eAAO,MAAM,SAAS,0DAAyE,CAAC;AAChG,eAAO,MAAM,OAAO,0DAAqE,CAAC;AAC1F,eAAO,MAAM,aAAa,gDAAiF,CAAC;AAC5G,eAAO,MAAM,SAAS,0DAAyE,CAAC;AAChG,eAAO,MAAM,WAAW,0DAA6E,CAAC;AACtG,eAAO,MAAM,WAAW,0DAA6E,CAAC;AACtG,eAAO,MAAM,SAAS,0DAAyE,CAAC;AAChG,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,mBAAmB,0DAA6F,CAAC;AAC9H,eAAO,MAAM,QAAQ,0DAAuE,CAAC;AAC7F,eAAO,MAAM,SAAS,0DAAyE,CAAC;AAChG,eAAO,MAAM,aAAa,gDAAiF,CAAC;AAC5G,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,cAAc,0DAAmF,CAAC;AAC/G,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,gBAAgB,0DAAuF,CAAC;AACrH,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,WAAW,0DAA6E,CAAC;AACtG,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,uBAAuB,0DAAqG,CAAC;AAC1I,eAAO,MAAM,OAAO,0DAAqE,CAAC;AAC1F,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,iBAAiB,0DAAyF,CAAC;AACxH,eAAO,MAAM,aAAa,0DAAiF,CAAC;AAC5G,eAAO,MAAM,oBAAoB,0DAA+F,CAAC;AACjI,eAAO,MAAM,OAAO,0DAAqE,CAAC;AAC1F,eAAO,MAAM,mBAAmB,0DAA6F,CAAC;AAC9H,eAAO,MAAM,eAAe,gDAAqF,CAAC;AAClH,eAAO,MAAM,YAAY,gDAA+E,CAAC;AACzG,eAAO,MAAM,aAAa,0DAAiF,CAAC;AAO5G,eAAO,MAAM,OAAO,EAA2B,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,KAAK,OAAO,CAAC;AACrF,eAAO,MAAM,QAAQ,EAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,OAAO,CAAC;AACpF,eAAO,MAAM,SAAS,EAA6B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,KAAK,OAAO,CAAC;AACnF,eAAO,MAAM,MAAM,EAA0B,CAAC,IAAI,CAAC,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,MAAM,KAAK,OAAO,CAAC;AAC1F,eAAO,MAAM,MAAM,EAAsC,CAAC,IAAI,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;AACrF,eAAO,MAAM,MAAM,GAAI,IAAI,OAAO,KAAG,MAAwB,CAAC;AAC9D,eAAO,MAAM,MAAM,GAAI,IAAI,OAAO,KAAG,MAA4C,CAAC;AAClF,eAAO,MAAM,OAAO,GAAI,IAAI,OAAO,KAAG,OAAqB,CAAC;AAC5D,eAAO,MAAM,MAAM,GAAI,IAAI,OAAO,KAAG,MAAmC,CAAC;AAOzE,eAAO,MAAM,MAAM,SAAiB,CAAC;AACrC,eAAO,MAAM,UAAU,SAAqB,CAAC;AAC7C,eAAO,MAAM,UAAU,SAAqB,CAAC;AAC7C,OAAO,EAAE,SAAS,IAAI,KAAK,EAAE,OAAO,IAAI,KAAK,EAAE,CAAC"}
1
+ {"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/factories/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAwIH,eAAO,MAAM,SAAS,SAAyE,CAAC;AAChG,eAAO,MAAM,cAAc,SAAmF,CAAC;AAC/G,eAAO,MAAM,cAAc,SAAmF,CAAC;AAe/G,eAAO,MAAM,MAAM,SAAmE,CAAC;AACvF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAC1F,eAAO,MAAM,OAAO,SAAqE,CAAC;AAC1F,eAAO,MAAM,iBAAiB,SAAyF,CAAC;AACxH,eAAO,MAAM,GAAG,SAA6D,CAAC;AAC9E,eAAO,MAAM,IAAI,SAA+D,CAAC;AACjF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAc1F,eAAO,MAAM,GAAG,SAA6D,CAAC;AAO9E,eAAO,MAAM,KAAK,qCAAiE,CAAC;AACpF,eAAO,MAAM,KAAK,qCAAiE,CAAC;AACpF,eAAO,MAAM,KAAK,qCAAiE,CAAC;AACpF,eAAO,MAAM,IAAI,qCAA+D,CAAC;AACjF,eAAO,MAAM,IAAI,qCAA+D,CAAC;AACjF,eAAO,MAAM,IAAI,qCAA+D,CAAC;AAGjF,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,KAAK,SAAiE,CAAC;AACpF,eAAO,MAAM,SAAS,qCAAyE,CAAC;AAChG,eAAO,MAAM,KAAK,SAAiE,CAAC;AACpF,eAAO,MAAM,UAAU,SAA2E,CAAC;AACnG,eAAO,MAAM,SAAS,qCAAyE,CAAC;AAChG,eAAO,MAAM,UAAU,SAA2E,CAAC;AACnG,eAAO,MAAM,OAAO,SAAqE,CAAC;AAC1F,eAAO,MAAM,OAAO,mEAAqE,CAAC;AAC1F,eAAO,MAAM,MAAM,SAAmE,CAAC;AAMvF,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,aAAa,SAAiE,CAAC;AAC5F,eAAO,MAAM,aAAa,SAAiE,CAAC;AAC5F,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,cAAc,SAAmE,CAAC;AAC/F,eAAO,MAAM,kBAAkB,SAA2E,CAAC;AAC3G,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,aAAa,qCAAiE,CAAC;AAC5F,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,aAAa,qCAAiE,CAAC;AAC5F,eAAO,MAAM,YAAY,SAA+D,CAAC;AACzF,eAAO,MAAM,aAAa,qCAAiE,CAAC;AAC5F,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,WAAW,qCAA6D,CAAC;AACtF,eAAO,MAAM,WAAW,qCAA6D,CAAC;AACtF,eAAO,MAAM,WAAW,qCAA6D,CAAC;AACtF,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AACzF,eAAO,MAAM,WAAW,SAA6D,CAAC;AACtF,eAAO,MAAM,YAAY,qCAA+D,CAAC;AA2EzF,eAAO,MAAM,qBAAqB,oEAAiG,CAAC;AAIpI,eAAO,MAAM,YAAY,SAA+E,CAAC;AAMzG,eAAO,MAAM,IAAI,SAA+D,CAAC;AACjF,eAAO,MAAM,IAAI,SAA+D,CAAC;AAQjF,eAAO,MAAM,eAAe,SAAqF,CAAC;AAClH,eAAO,MAAM,QAAQ,qCAAuE,CAAC;AAC7F,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,iBAAiB,SAAyE,CAAC;AACxG,eAAO,MAAM,WAAW,SAA6D,CAAC;AAyCtF,eAAO,MAAM,iBAAiB,SAAyE,CAAC;AAGxG,eAAO,MAAM,kBAAkB,SAA2E,CAAC;AAG3G,eAAO,MAAM,QAAQ,SAAuE,CAAC;AAG7F,eAAO,MAAM,KAAK,SAAiE,CAAC;AAGpF,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,kBAAkB,qCAA2F,CAAC;AAG3H,eAAO,MAAM,iBAAiB,SAAyF,CAAC;AAGxH,eAAO,MAAM,cAAc,SAAmF,CAAC;AAG/G,eAAO,MAAM,KAAK,SAAiE,CAAC;AAQpF,eAAO,MAAM,KAAK,SAAiE,CAAC;AAQpF,QAAA,MAAM,UAAU,mDAA6D,CAAC;AAM9E,eAAO,MAAM,GAAG,EACqE,OAAO,UAAU,CAAC;AAIvG,eAAO,MAAM,OAAO,SAAqE,CAAC;AAuD1F,eAAO,MAAM,aAAa,SAAiE,CAAC;AAI5F,eAAO,MAAM,SAAS,SAAyE,CAAC;AAIhG,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,SAAS,qCAAyE,CAAC;AAGhG,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,QAAA,MAAM,UAAU,mDAA6D,CAAC;AAK9E,eAAO,MAAM,GAAG,EACqE,OAAO,UAAU,CAAC;AAGvG,eAAO,MAAM,YAAY,SAA+D,CAAC;AAGzF,eAAO,MAAM,QAAQ,qCAAuE,CAAC;AAG7F,eAAO,MAAM,aAAa,SAAiE,CAAC;AAG5F,eAAO,MAAM,YAAY,qCAA+D,CAAC;AAGzF,eAAO,MAAM,WAAW,SAA6D,CAAC;AAItF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAQ1F,eAAO,MAAM,MAAM,KAAmE,CAAC;AAIvF,eAAO,MAAM,cAAc,SAAmE,CAAC;AAG/F,eAAO,MAAM,WAAW,SAA6D,CAAC;AAItF,eAAO,MAAM,QAAQ,uJAAuE,CAAC;AAG7F,eAAO,MAAM,SAAS,qKAAyE,CAAC;AAGhG,eAAO,MAAM,KAAK,qGAAiE,CAAC;AAGpF,eAAO,MAAM,KAAK,wJAAiE,CAAC;AAIpF,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,SAAS,wDAAyE,CAAC;AAGhG,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,SAAS,wDAAyE,CAAC;AAqDhG,eAAO,MAAM,eAAe,SAAqE,CAAC;AAWlG,eAAO,MAAM,cAAc,SAAmE,CAAC;AAG/F,eAAO,MAAM,gBAAgB,SAAuE,CAAC;AAGrG,eAAO,MAAM,eAAe,SAAqE,CAAC;AAGlG,eAAO,MAAM,iBAAiB,SAAyE,CAAC;AAQxG,eAAO,MAAM,SAAS,qCAAyE,CAAC;AAGhG,eAAO,MAAM,WAAW,qCAA6E,CAAC;AAGtG,eAAO,MAAM,WAAW,qCAA6D,CAAC;AAGtF,eAAO,MAAM,WAAW,qCAA6D,CAAC;AAGtF,eAAO,MAAM,aAAa,SAAiE,CAAC;AAG5F,eAAO,MAAM,WAAW,SAA6D,CAAC;AAGtF,eAAO,MAAM,eAAe,qCAAqE,CAAC;AAGlG,eAAO,MAAM,WAAW,SAA6D,CAAC;AAGtF,eAAO,MAAM,YAAY,SAA+D,CAAC;AAGzF,eAAO,MAAM,aAAa,SAAiE,CAAC;AA8B5F,eAAO,MAAM,IAAI,0FAA+D,CAAC;AAIjF,eAAO,MAAM,aAAa,SAAiE,CAAC;AAI5F,eAAO,MAAM,EAAE,SAA2D,CAAC;AAoD3E,eAAO,MAAM,MAAM,qCAAmE,CAAC;AAGvF,eAAO,MAAM,WAAW,SAA6D,CAAC;AAGtF,eAAO,MAAM,MAAM,qCAAmE,CAAC;AAWvF,eAAO,MAAM,WAAW,SAA6E,CAAC;AAItG,eAAO,MAAM,eAAe,SAAqF,CAAC;AAQlH,eAAO,MAAM,EAAE,SAA2D,CAAC;AAG3E,eAAO,MAAM,KAAK,SAAiE,CAAC;AAGpF,eAAO,MAAM,QAAQ,SAAuE,CAAC;AAI7F,eAAO,MAAM,KAAK,SAAiE,CAAC;AAIpF,eAAO,MAAM,WAAW,SAA6D,CAAC;AAGtF,eAAO,MAAM,WAAW,SAA6D,CAAC;AAItF,eAAO,MAAM,aAAa,qCAAiE,CAAC;AAW5F,eAAO,MAAM,SAAS,SAAyE,CAAC;AAOhG,eAAO,MAAM,IAAI,qCAA+D,CAAC;AAUjF,eAAO,MAAM,UAAU,SAA2E,CAAC;AAanG,eAAO,MAAM,WAAW,SAA6E,CAAC;AAgKtG,eAAO,MAAM,cAAc;;;;;;;CAAmF,CAAC;AAG/G,eAAO,MAAM,SAAS,SAAyE,CAAC;AAGhG,eAAO,MAAM,OAAO,SAAqE,CAAC;AAG1F,eAAO,MAAM,gBAAgB;;;;yCAAuF,CAAC;AAGrH,eAAO,MAAM,YAAY;;;;;;;;;;;;;;;;;;;YA/4BtB,CAAC;;;;;YAYD,CAAA;UAAgB,CAAC;;;;;YAcpB,CAAA;;;;;YAW8D,CAAC;;;;;;;;;;;;;CA02ByC,CAAC;AAGzG,eAAO,MAAM,SAAS,SAAyE,CAAC;AAGhG,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,KAAK,mDAA2E,CAAC;AA4B9F,eAAO,MAAM,gBAAgB,SAAuE,CAAC;AAGrG,eAAO,MAAM,cAAc,SAAmE,CAAC;AAO/F,eAAO,MAAM,IAAI,SAA+D,CAAC;AAMjF,eAAO,MAAM,KAAK,qCAAiE,CAAC;AAGpF,eAAO,MAAM,kBAAkB,SAA2E,CAAC;AAG3G,eAAO,MAAM,IAAI,SAA+E,CAAC;AA2CjG,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,KAAK,SAAiE,CAAC;AAGpF,eAAO,MAAM,IAAI,qCAA+D,CAAC;AAGjF,eAAO,MAAM,KAAK,2FAAiE,CAAC;AAGpF,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,MAAM,0LAAmE,CAAC;AAGvF,eAAO,MAAM,IAAI,sMAA+D,CAAC;AAGjF,eAAO,MAAM,SAAS,mNAAyE,CAAC;AAGhG,eAAO,MAAM,SAAS,SAAyE,CAAC;AAGhG,eAAO,MAAM,YAAY,SAA+D,CAAC;AAGzF,eAAO,MAAM,MAAM,SAAmE,CAAC;AAGvF,eAAO,MAAM,gBAAgB,SAAuE,CAAC;AAGrG,eAAO,MAAM,WAAW,SAA6E,CAAC;AAGtG,eAAO,MAAM,YAAY,SAA+E,CAAC;AAOzG,eAAO,MAAM,KAAK,KAAiE,CAAC;AA+BpF,eAAO,MAAM,YAAY,SAA+E,CAAC;AAGzG,eAAO,MAAM,cAAc,SAAmF,CAAC;AAG/G,eAAO,MAAM,QAAQ,KAAuE,CAAC;AAG7F,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAe1F,eAAO,MAAM,IAAI,mDAA+D,CAAC;AAGjF,eAAO,MAAM,OAAO,SAAqE,CAAC;AAG1F,eAAO,MAAM,IAAI,SAA+D,CAAC;AAGjF,eAAO,MAAM,GAAG,SAA6D,CAAC;AAG9E,eAAO,MAAM,WAAW,SAA6D,CAAC;AAOtF,eAAO,MAAM,QAAQ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CAAuE,CAAC;AAa7F,eAAO,MAAM,UAAU;;;;;CAA2E,CAAC;AAGnG,eAAO,MAAM,YAAY,qCAA+D,CAAC;AAGzF,eAAO,MAAM,WAAW,SAA6E,CAAC;AAOtG,eAAO,MAAM,aAAa,SAAiF,CAAC;AAU5G,eAAO,MAAM,cAAc,qCAAmF,CAAC;AAG/G,eAAO,MAAM,KAAK,SAAiE,CAAC;AAUpF,eAAO,MAAM,MAAM,qCAAmE,CAAC;AAGvF,eAAO,MAAM,SAAS,SAAyE,CAAC;AAShG,eAAO,MAAM,IAAI,SAA+D,CAAC;AAkEjF,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,mBAAmB,0DAA6F,CAAC;AAC9H,eAAO,MAAM,cAAc,0DAAmF,CAAC;AAC/G,eAAO,MAAM,qBAAqB,0DAAiG,CAAC;AACpI,eAAO,MAAM,gBAAgB,0DAAuF,CAAC;AACrH,eAAO,MAAM,gBAAgB,0DAAuF,CAAC;AACrH,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,OAAO,0DAAqE,CAAC;AAC1F,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,gBAAgB,0DAAuF,CAAC;AACrH,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,kBAAkB,0DAA2F,CAAC;AAC3H,eAAO,MAAM,yBAAyB,0DAAyG,CAAC;AAChJ,eAAO,MAAM,mBAAmB,0DAA6F,CAAC;AAC9H,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,QAAQ,0DAAuE,CAAC;AAC7F,eAAO,MAAM,SAAS,0DAAyE,CAAC;AAChG,eAAO,MAAM,OAAO,0DAAqE,CAAC;AAC1F,eAAO,MAAM,aAAa,gDAAiF,CAAC;AAC5G,eAAO,MAAM,SAAS,0DAAyE,CAAC;AAChG,eAAO,MAAM,WAAW,0DAA6E,CAAC;AACtG,eAAO,MAAM,WAAW,0DAA6E,CAAC;AACtG,eAAO,MAAM,SAAS,0DAAyE,CAAC;AAChG,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,mBAAmB,0DAA6F,CAAC;AAC9H,eAAO,MAAM,QAAQ,0DAAuE,CAAC;AAC7F,eAAO,MAAM,SAAS,0DAAyE,CAAC;AAChG,eAAO,MAAM,aAAa,gDAAiF,CAAC;AAC5G,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,cAAc,0DAAmF,CAAC;AAC/G,eAAO,MAAM,eAAe,0DAAqF,CAAC;AAClH,eAAO,MAAM,gBAAgB,0DAAuF,CAAC;AACrH,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,WAAW,0DAA6E,CAAC;AACtG,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,uBAAuB,0DAAqG,CAAC;AAC1I,eAAO,MAAM,OAAO,0DAAqE,CAAC;AAC1F,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,UAAU,0DAA2E,CAAC;AACnG,eAAO,MAAM,YAAY,0DAA+E,CAAC;AACzG,eAAO,MAAM,iBAAiB,0DAAyF,CAAC;AACxH,eAAO,MAAM,aAAa,0DAAiF,CAAC;AAC5G,eAAO,MAAM,oBAAoB,0DAA+F,CAAC;AACjI,eAAO,MAAM,OAAO,0DAAqE,CAAC;AAC1F,eAAO,MAAM,mBAAmB,0DAA6F,CAAC;AAC9H,eAAO,MAAM,eAAe,gDAAqF,CAAC;AAClH,eAAO,MAAM,YAAY,gDAA+E,CAAC;AACzG,eAAO,MAAM,aAAa,0DAAiF,CAAC;AAO5G,eAAO,MAAM,OAAO,EAA2B,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,CAAC,EAAE,MAAM,KAAK,OAAO,CAAC;AACrF,eAAO,MAAM,QAAQ,EAA4B,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,OAAO,CAAC;AACpF,eAAO,MAAM,SAAS,EAA6B,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,KAAK,OAAO,CAAC;AACnF,eAAO,MAAM,MAAM,EAA0B,CAAC,IAAI,CAAC,EAAE,OAAO,EAAE,MAAM,CAAC,EAAE,MAAM,KAAK,OAAO,CAAC;AAC1F,eAAO,MAAM,MAAM,EAAsC,CAAC,IAAI,CAAC,EAAE,OAAO,KAAK,OAAO,CAAC;AACrF,eAAO,MAAM,MAAM,GAAI,IAAI,OAAO,KAAG,MAAwB,CAAC;AAC9D,eAAO,MAAM,MAAM,GAAI,IAAI,OAAO,KAAG,MAA4C,CAAC;AAClF,eAAO,MAAM,OAAO,GAAI,IAAI,OAAO,KAAG,OAAqB,CAAC;AAC5D,eAAO,MAAM,MAAM,GAAI,IAAI,OAAO,KAAG,MAAmC,CAAC;AAOzE,eAAO,MAAM,MAAM,SAAiB,CAAC;AACrC,eAAO,MAAM,UAAU,SAAqB,CAAC;AAC7C,eAAO,MAAM,UAAU,SAAqB,CAAC;AAC7C,OAAO,EAAE,SAAS,IAAI,KAAK,EAAE,OAAO,IAAI,KAAK,EAAE,CAAC"}
@@ -0,0 +1,31 @@
1
+ /**
2
+ * Forward-mode automatic differentiation over the plain functions surface.
3
+ *
4
+ * The elementary typed functions (`add`/`subtract`/`multiply`/`divide`/`pow`/
5
+ * `sin`/`cos`/`tan`/`exp`/`log`/`sqrt`/`square`/`cbrt`/`cube`/`abs`) carry `Dual`
6
+ * signatures, so a function written with them differentiates exactly when
7
+ * evaluated on a seed `Dual`:
8
+ *
9
+ * import { multiply, sin, derivativeAt } from '@danielsimonjr/mathts-functions';
10
+ * derivativeAt((x) => multiply(sin(x), x), 2) // sin(2) + 2·cos(2)
11
+ *
12
+ * This complements the symbolic `derivative` (which returns an expression) and
13
+ * the autograd package's tensor `grad`: `derivativeAt` evaluates the exact
14
+ * scalar derivative of an ordinary numeric function at a point, with no finite
15
+ * differences and no manual tape.
16
+ */
17
+ import { Dual } from '@danielsimonjr/mathts-core';
18
+ export type DualFn = (x: Dual) => Dual;
19
+ /** Exact derivative of `fn: ℝ → ℝ` at `x0` (forward-mode AD). */
20
+ export declare function derivativeAt(fn: DualFn, x0: number): number;
21
+ /** Both the value and the exact derivative of `fn` at `x0`. */
22
+ export declare function valueAndDerivativeAt(fn: DualFn, x0: number): {
23
+ value: number;
24
+ deriv: number;
25
+ };
26
+ /**
27
+ * Gradient of `fn: ℝⁿ → ℝ` at `x` by n forward passes (seed one coordinate at a
28
+ * time). `fn` receives an array of `Dual`s and must return a scalar `Dual`.
29
+ */
30
+ export declare function gradientAt(fn: (x: Dual[]) => Dual, x: readonly number[]): number[];
31
+ //# sourceMappingURL=grad-forward.d.ts.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"grad-forward.d.ts","sourceRoot":"","sources":["../src/grad-forward.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AACH,OAAO,EAAE,IAAI,EAAE,MAAM,4BAA4B,CAAC;AAElD,MAAM,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,IAAI,KAAK,IAAI,CAAC;AAEvC,iEAAiE;AACjE,wBAAgB,YAAY,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM,CAE3D;AAED,+DAA+D;AAC/D,wBAAgB,oBAAoB,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG;IAAE,KAAK,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,MAAM,CAAA;CAAE,CAG7F;AAED;;;GAGG;AACH,wBAAgB,UAAU,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,KAAK,IAAI,EAAE,CAAC,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM,EAAE,CAOlF"}
package/dist/index.d.ts CHANGED
@@ -14,4 +14,6 @@ export * from './factories/index.js';
14
14
  export { to, toBest } from './typed/unit.js';
15
15
  export { evaluate, compileExpr, parse, parser, reviver, replacer } from './factories/evaluate.js';
16
16
  export { help } from './help.js';
17
+ export { derivativeAt, valueAndDerivativeAt, gradientAt } from './grad-forward.js';
18
+ export type { DualFn } from './grad-forward.js';
17
19
  //# sourceMappingURL=index.d.ts.map
@@ -1 +1 @@
1
- {"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAGH,cAAc,kBAAkB,CAAC;AAKjC,cAAc,gBAAgB,CAAC;AAG/B,cAAc,sBAAsB,CAAC;AAOrC,OAAO,EAAE,EAAE,EAAE,MAAM,EAAE,MAAM,iBAAiB,CAAC;AAG7C,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,yBAAyB,CAAC;AAGlG,OAAO,EAAE,IAAI,EAAE,MAAM,WAAW,CAAC"}
1
+ {"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAGH,cAAc,kBAAkB,CAAC;AAKjC,cAAc,gBAAgB,CAAC;AAG/B,cAAc,sBAAsB,CAAC;AAOrC,OAAO,EAAE,EAAE,EAAE,MAAM,EAAE,MAAM,iBAAiB,CAAC;AAG7C,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,yBAAyB,CAAC;AAGlG,OAAO,EAAE,IAAI,EAAE,MAAM,WAAW,CAAC;AAGjC,OAAO,EAAE,YAAY,EAAE,oBAAoB,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AACnF,YAAY,EAAE,MAAM,EAAE,MAAM,mBAAmB,CAAC"}
package/dist/index.js CHANGED
@@ -468,7 +468,7 @@ __export(typed_exports, {
468
468
  });
469
469
 
470
470
  // src/typed/arithmetic.ts
471
- import { mathTyped, Complex, Fraction, BigNumber } from "@danielsimonjr/mathts-core";
471
+ import { mathTyped, Complex, Fraction, BigNumber, Dual } from "@danielsimonjr/mathts-core";
472
472
  import { DenseMatrix, backendManager } from "@danielsimonjr/mathts-matrix";
473
473
  import { computePool } from "@danielsimonjr/mathts-parallel";
474
474
 
@@ -510,6 +510,10 @@ var add = mathTyped("add", {
510
510
  "BigNumber, BigNumber": (a, b) => a.add(b),
511
511
  // Unit arithmetic (mathjs parity): same-dimension units add (e.g. 5 cm + 3 mm).
512
512
  "Unit, Unit": (a, b) => a.add(b),
513
+ // Dual numbers (forward-mode AD).
514
+ "Dual, Dual": (a, b) => a.add(b),
515
+ "Dual, number": (a, b) => a.add(Dual.constant(b)),
516
+ "number, Dual": (a, b) => Dual.constant(a).add(b),
513
517
  // Mixed type operations (auto-coercion via conversions)
514
518
  "number, Complex": (a, b) => Complex.fromNumber(a).add(b),
515
519
  "Complex, number": (a, b) => a.add(Complex.fromNumber(b)),
@@ -544,6 +548,10 @@ var subtract = mathTyped("subtract", {
544
548
  "BigNumber, BigNumber": (a, b) => a.subtract(b),
545
549
  // Unit subtraction (same dimension): 5 cm − 3 mm.
546
550
  "Unit, Unit": (a, b) => a.sub(b),
551
+ // Dual numbers (forward-mode AD).
552
+ "Dual, Dual": (a, b) => a.sub(b),
553
+ "Dual, number": (a, b) => a.sub(Dual.constant(b)),
554
+ "number, Dual": (a, b) => Dual.constant(a).sub(b),
547
555
  "number, Complex": (a, b) => Complex.fromNumber(a).subtract(b),
548
556
  "Complex, number": (a, b) => a.subtract(Complex.fromNumber(b)),
549
557
  "number, Fraction": (a, b) => Fraction.fromNumber(a).subtract(b),
@@ -566,6 +574,10 @@ var multiply = mathTyped("multiply", {
566
574
  "Unit, Unit": (a, b) => a.mul(b),
567
575
  "Unit, number": (a, b) => a.mul(b),
568
576
  "number, Unit": (a, b) => b.mul(a),
577
+ // Dual numbers (forward-mode AD).
578
+ "Dual, Dual": (a, b) => a.mul(b),
579
+ "Dual, number": (a, b) => a.mul(Dual.constant(b)),
580
+ "number, Dual": (a, b) => Dual.constant(a).mul(b),
569
581
  // Matrix × matrix (2-D arrays) routed through the native DenseMatrix +
570
582
  // BackendManager — i.e. WASM/GPU acceleration for large matrices (GC7).
571
583
  // Previously `multiply(2DArray, 2DArray)` threw; this is additive. Non-2-D
@@ -621,6 +633,10 @@ var divide = mathTyped("divide", {
621
633
  // Unit division / scaling (mathjs parity). Unit ÷ Unit may be dimensionless.
622
634
  "Unit, Unit": (a, b) => a.div(b),
623
635
  "Unit, number": (a, b) => a.div(b),
636
+ // Dual numbers (forward-mode AD).
637
+ "Dual, Dual": (a, b) => a.div(b),
638
+ "Dual, number": (a, b) => a.div(Dual.constant(b)),
639
+ "number, Dual": (a, b) => Dual.constant(a).div(b),
624
640
  "number, Complex": (a, b) => Complex.fromNumber(a).divide(b),
625
641
  "Complex, number": (a, b) => a.divide(Complex.fromNumber(b)),
626
642
  "number, Fraction": (a, b) => Fraction.fromNumber(a).divide(b),
@@ -660,6 +676,7 @@ var abs = mathTyped("abs", {
660
676
  Complex: (a) => a.abs(),
661
677
  Fraction: (a) => a.abs(),
662
678
  BigNumber: (a) => a.abs(),
679
+ Dual: (a) => a.abs(),
663
680
  // Parallel array abs
664
681
  Float64Array: async (a) => {
665
682
  const wasm = elementwiseUnaryDispatch("abs", a);
@@ -696,7 +713,10 @@ var pow = mathTyped("pow", {
696
713
  "Complex, Complex": (a, b) => a.pow(b),
697
714
  "Fraction, number": (a, b) => a.pow(Math.floor(b)),
698
715
  "BigNumber, number": (a, b) => a.pow(b),
699
- "BigNumber, BigNumber": (a, b) => a.pow(b.valueOf())
716
+ "BigNumber, BigNumber": (a, b) => a.pow(b.valueOf()),
717
+ // Dual numbers (forward-mode AD): constant or variable exponent.
718
+ "Dual, number": (a, b) => a.powConst(b),
719
+ "Dual, Dual": (a, b) => a.pow(b)
700
720
  });
701
721
  var sqrt = mathTyped("sqrt", {
702
722
  number: (a) => {
@@ -707,6 +727,7 @@ var sqrt = mathTyped("sqrt", {
707
727
  },
708
728
  Complex: (a) => a.sqrt(),
709
729
  BigNumber: (a) => a.sqrt(),
730
+ Dual: (a) => a.sqrt(),
710
731
  // Parallel array sqrt
711
732
  Float64Array: async (a) => {
712
733
  const result = await computePool.sqrt(a);
@@ -719,6 +740,7 @@ var square = mathTyped("square", {
719
740
  Complex: (a) => a.multiply(a),
720
741
  Fraction: (a) => a.multiply(a),
721
742
  BigNumber: (a) => a.multiply(a),
743
+ Dual: (a) => a.square(),
722
744
  // Parallel array square
723
745
  Float64Array: async (a) => {
724
746
  const result = await computePool.square(a);
@@ -731,6 +753,7 @@ var cube = mathTyped("cube", {
731
753
  Complex: (a) => a.multiply(a).multiply(a),
732
754
  Fraction: (a) => a.multiply(a).multiply(a),
733
755
  BigNumber: (a) => a.multiply(a).multiply(a),
756
+ Dual: (a) => a.powConst(3),
734
757
  // Parallel array cube
735
758
  Float64Array: async (a) => {
736
759
  if (computePool.shouldParallelize(a.length)) {
@@ -746,6 +769,7 @@ var cbrt = mathTyped("cbrt", {
746
769
  number: (a) => Math.cbrt(a),
747
770
  Complex: (a) => a.pow(1 / 3),
748
771
  BigNumber: (a) => a.cbrt(),
772
+ Dual: (a) => a.powConst(1 / 3),
749
773
  // Parallel array cbrt
750
774
  Float64Array: async (a) => {
751
775
  if (computePool.shouldParallelize(a.length)) {
@@ -766,6 +790,7 @@ var exp = mathTyped("exp", {
766
790
  number: (a) => Math.exp(a),
767
791
  Complex: (a) => a.exp(),
768
792
  BigNumber: (a) => a.exp(),
793
+ Dual: (a) => a.exp(),
769
794
  // Parallel array exp
770
795
  Float64Array: async (a) => {
771
796
  const wasm = elementwiseUnaryDispatch("exp", a);
@@ -778,6 +803,7 @@ var log = mathTyped("log", {
778
803
  number: (a) => Math.log(a),
779
804
  Complex: (a) => a.log(),
780
805
  BigNumber: (a) => a.ln(),
806
+ Dual: (a) => a.log(),
781
807
  "number, number": (a, base) => Math.log(a) / Math.log(base),
782
808
  // Parallel array log
783
809
  Float64Array: async (a) => {
@@ -1349,6 +1375,7 @@ var sin = mathTyped2("sin", {
1349
1375
  number: (a) => Math.sin(a),
1350
1376
  Complex: (a) => a.sin(),
1351
1377
  BigNumber: (a) => a.sin(),
1378
+ Dual: (a) => a.sin(),
1352
1379
  // Parallel Float64Array sin
1353
1380
  Float64Array: async (a) => {
1354
1381
  const wasm = elementwiseUnaryDispatch("sin", a);
@@ -1361,6 +1388,7 @@ var cos = mathTyped2("cos", {
1361
1388
  number: (a) => Math.cos(a),
1362
1389
  Complex: (a) => a.cos(),
1363
1390
  BigNumber: (a) => a.cos(),
1391
+ Dual: (a) => a.cos(),
1364
1392
  // Parallel Float64Array cos
1365
1393
  Float64Array: async (a) => {
1366
1394
  const wasm = elementwiseUnaryDispatch("cos", a);
@@ -1373,6 +1401,7 @@ var tan = mathTyped2("tan", {
1373
1401
  number: (a) => Math.tan(a),
1374
1402
  Complex: (a) => a.tan(),
1375
1403
  BigNumber: (a) => a.tan(),
1404
+ Dual: (a) => a.tan(),
1376
1405
  // Parallel Float64Array tan
1377
1406
  Float64Array: async (a) => {
1378
1407
  const wasm = elementwiseUnaryDispatch("tan", a);
@@ -20260,7 +20289,7 @@ var dependencies85 = ["typed", "config", "matrix", "BigNumber", "DenseMatrix", "
20260
20289
  var createIdentity = /* @__PURE__ */ factory(
20261
20290
  name85,
20262
20291
  dependencies85,
20263
- ({ typed: typed3, config, matrix: matrix2, BigNumber: BigNumber9, DenseMatrix: DenseMatrix4, SparseMatrix }) => {
20292
+ ({ typed: typed3, config, matrix: matrix2, BigNumber: BigNumber9, DenseMatrix: DenseMatrix5, SparseMatrix }) => {
20264
20293
  return typed3(name85, {
20265
20294
  "": function() {
20266
20295
  return config.matrix === "Matrix" ? matrix2([]) : [];
@@ -20323,7 +20352,7 @@ var createIdentity = /* @__PURE__ */ factory(
20323
20352
  return SparseMatrix.diagonal(size2, one, 0, defaultValue);
20324
20353
  }
20325
20354
  if (format6 === "dense") {
20326
- return DenseMatrix4.diagonal(size2, one, 0, defaultValue);
20355
+ return DenseMatrix5.diagonal(size2, one, 0, defaultValue);
20327
20356
  }
20328
20357
  throw new TypeError(`Unknown matrix type "${format6}"`);
20329
20358
  }
@@ -20431,7 +20460,7 @@ var dependencies88 = ["typed", "matrix", "DenseMatrix", "SparseMatrix"];
20431
20460
  var createDiag = /* @__PURE__ */ factory(
20432
20461
  name88,
20433
20462
  dependencies88,
20434
- ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix4, SparseMatrix }) => {
20463
+ ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix5, SparseMatrix }) => {
20435
20464
  return typed3(name88, {
20436
20465
  // FIXME: simplify this huge amount of signatures as soon as typed-function supports optional arguments
20437
20466
  Array: function(x) {
@@ -20490,7 +20519,7 @@ var createDiag = /* @__PURE__ */ factory(
20490
20519
  if (format6 && format6 !== "sparse" && format6 !== "dense") {
20491
20520
  throw new TypeError(`Unknown matrix type ${format6}"`);
20492
20521
  }
20493
- const m = format6 === "sparse" ? SparseMatrix.diagonal(ms, x, k) : DenseMatrix4.diagonal(ms, x, k);
20522
+ const m = format6 === "sparse" ? SparseMatrix.diagonal(ms, x, k) : DenseMatrix5.diagonal(ms, x, k);
20494
20523
  return format6 !== null ? m : m.valueOf();
20495
20524
  }
20496
20525
  function _getDiagonal(x, k, format6, s, kSub, kSuper) {
@@ -21052,6 +21081,79 @@ var createDet = /* @__PURE__ */ factory(
21052
21081
  }
21053
21082
  );
21054
21083
 
21084
+ // src/matrix/native-accel.ts
21085
+ import { DenseMatrix as DenseMatrix4, lu } from "@danielsimonjr/mathts-matrix";
21086
+ var NATIVE_MATRIX_THRESHOLD = 8;
21087
+ function isLargeNumericSquare(a) {
21088
+ if (!Array.isArray(a) || a.length < NATIVE_MATRIX_THRESHOLD) return false;
21089
+ const n = a.length;
21090
+ for (const row2 of a) {
21091
+ if (!Array.isArray(row2) || row2.length !== n) return false;
21092
+ for (let j = 0; j < n; j++) if (typeof row2[j] !== "number") return false;
21093
+ }
21094
+ return true;
21095
+ }
21096
+ function permutationParity(P2) {
21097
+ let sign3 = 1;
21098
+ const seen = new Array(P2.length).fill(false);
21099
+ for (let i = 0; i < P2.length; i++) {
21100
+ if (seen[i]) continue;
21101
+ let j = i;
21102
+ let len = 0;
21103
+ while (!seen[j]) {
21104
+ seen[j] = true;
21105
+ j = P2[j];
21106
+ len++;
21107
+ }
21108
+ if (len % 2 === 0) sign3 = -sign3;
21109
+ }
21110
+ return sign3;
21111
+ }
21112
+ var SINGULAR = /singular|zero pivot/i;
21113
+ function acceleratedDet(a, factoryDet2) {
21114
+ if (!isLargeNumericSquare(a)) return factoryDet2(a);
21115
+ try {
21116
+ const { U, P: P2 } = lu(DenseMatrix4.fromArray(a));
21117
+ const Ud = U.toArray();
21118
+ let d = permutationParity(P2);
21119
+ for (let i = 0; i < a.length; i++) d *= Ud[i][i];
21120
+ return d;
21121
+ } catch (e) {
21122
+ if (e instanceof Error && SINGULAR.test(e.message)) return 0;
21123
+ return factoryDet2(a);
21124
+ }
21125
+ }
21126
+ function acceleratedInv(a, factoryInv2) {
21127
+ if (!isLargeNumericSquare(a)) return factoryInv2(a);
21128
+ try {
21129
+ const { L, U, P: P2 } = lu(DenseMatrix4.fromArray(a));
21130
+ const Ld = L.toArray();
21131
+ const Ud = U.toArray();
21132
+ const n = a.length;
21133
+ const inv2 = Array.from({ length: n }, () => new Array(n).fill(0));
21134
+ for (let j = 0; j < n; j++) {
21135
+ const b = new Array(n).fill(0);
21136
+ for (let i = 0; i < n; i++) if (P2[i] === j) b[i] = 1;
21137
+ const y = new Array(n).fill(0);
21138
+ for (let i = 0; i < n; i++) {
21139
+ let s = b[i];
21140
+ for (let k = 0; k < i; k++) s -= Ld[i][k] * y[k];
21141
+ y[i] = s / Ld[i][i];
21142
+ }
21143
+ const x = new Array(n).fill(0);
21144
+ for (let i = n - 1; i >= 0; i--) {
21145
+ let s = y[i];
21146
+ for (let k = i + 1; k < n; k++) s -= Ud[i][k] * x[k];
21147
+ x[i] = s / Ud[i][i];
21148
+ }
21149
+ for (let i = 0; i < n; i++) inv2[i][j] = x[i];
21150
+ }
21151
+ return inv2;
21152
+ } catch {
21153
+ return factoryInv2(a);
21154
+ }
21155
+ }
21156
+
21055
21157
  // src/matrix/reshape.ts
21056
21158
  var name96 = "reshape";
21057
21159
  var dependencies96 = ["typed", "isInteger", "matrix"];
@@ -21958,7 +22060,7 @@ var createCsSymperm = /* @__PURE__ */ factory(
21958
22060
  );
21959
22061
 
21960
22062
  // src/algebra/solver/utils/solveValidation.ts
21961
- function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
22063
+ function createSolveValidation({ DenseMatrix: DenseMatrix5 }) {
21962
22064
  return function solveValidation(m, b, copy) {
21963
22065
  const mSize = m.size();
21964
22066
  if (mSize.length !== 2) {
@@ -21981,7 +22083,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
21981
22083
  for (let i = 0; i < rows; i++) {
21982
22084
  data[i] = [bdata[i]];
21983
22085
  }
21984
- return new DenseMatrix4({
22086
+ return new DenseMatrix5({
21985
22087
  data,
21986
22088
  size: [rows, 1],
21987
22089
  datatype: bm._datatype
@@ -21997,7 +22099,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
21997
22099
  for (let i = 0; i < rows; i++) {
21998
22100
  data[i] = [bdata[i][0]];
21999
22101
  }
22000
- return new DenseMatrix4({
22102
+ return new DenseMatrix5({
22001
22103
  data,
22002
22104
  size: [rows, 1],
22003
22105
  datatype: bm._datatype
@@ -22016,7 +22118,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
22016
22118
  const i = index[k];
22017
22119
  data[i][0] = values[k];
22018
22120
  }
22019
- return new DenseMatrix4({
22121
+ return new DenseMatrix5({
22020
22122
  data,
22021
22123
  size: [rows, 1],
22022
22124
  datatype: bm._datatype
@@ -22036,7 +22138,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
22036
22138
  for (let i = 0; i < rows; i++) {
22037
22139
  data[i] = [b[i]];
22038
22140
  }
22039
- return new DenseMatrix4({
22141
+ return new DenseMatrix5({
22040
22142
  data,
22041
22143
  size: [rows, 1]
22042
22144
  });
@@ -22048,7 +22150,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
22048
22150
  for (let i = 0; i < rows; i++) {
22049
22151
  data[i] = [b[i][0]];
22050
22152
  }
22051
- return new DenseMatrix4({
22153
+ return new DenseMatrix5({
22052
22154
  data,
22053
22155
  size: [rows, 1]
22054
22156
  });
@@ -22118,10 +22220,10 @@ var createLsolve = /* @__PURE__ */ factory(
22118
22220
  multiplyScalar: multiplyScalar2,
22119
22221
  subtractScalar: subtractScalar2,
22120
22222
  equalScalar: equalScalar3,
22121
- DenseMatrix: DenseMatrix4
22223
+ DenseMatrix: DenseMatrix5
22122
22224
  }) => {
22123
22225
  const solveValidation = createSolveValidation({
22124
- DenseMatrix: DenseMatrix4
22226
+ DenseMatrix: DenseMatrix5
22125
22227
  });
22126
22228
  return typed3(name101, {
22127
22229
  "SparseMatrix, Array | Matrix": function(m, b) {
@@ -22159,7 +22261,7 @@ var createLsolve = /* @__PURE__ */ factory(
22159
22261
  for (let i = 0; i < rows; i++) {
22160
22262
  x2[i] = [resultAlloc.array[i]];
22161
22263
  }
22162
- return new DenseMatrix4({
22264
+ return new DenseMatrix5({
22163
22265
  data: x2,
22164
22266
  size: [rows, 1]
22165
22267
  });
@@ -22192,7 +22294,7 @@ var createLsolve = /* @__PURE__ */ factory(
22192
22294
  }
22193
22295
  x[j] = [xj];
22194
22296
  }
22195
- return new DenseMatrix4({
22297
+ return new DenseMatrix5({
22196
22298
  data: x,
22197
22299
  size: [rows, 1]
22198
22300
  });
@@ -22236,7 +22338,7 @@ var createLsolve = /* @__PURE__ */ factory(
22236
22338
  x[j] = [0];
22237
22339
  }
22238
22340
  }
22239
- return new DenseMatrix4({
22341
+ return new DenseMatrix5({
22240
22342
  data: x,
22241
22343
  size: [rows, 1]
22242
22344
  });
@@ -22265,10 +22367,10 @@ var createLsolveAll = /* @__PURE__ */ factory(
22265
22367
  multiplyScalar: multiplyScalar2,
22266
22368
  subtractScalar: subtractScalar2,
22267
22369
  equalScalar: equalScalar3,
22268
- DenseMatrix: DenseMatrix4
22370
+ DenseMatrix: DenseMatrix5
22269
22371
  }) => {
22270
22372
  const solveValidation = createSolveValidation({
22271
- DenseMatrix: DenseMatrix4
22373
+ DenseMatrix: DenseMatrix5
22272
22374
  });
22273
22375
  return typed3(name102, {
22274
22376
  "SparseMatrix, Array | Matrix": function(m, b) {
@@ -22318,7 +22420,7 @@ var createLsolveAll = /* @__PURE__ */ factory(
22318
22420
  }
22319
22421
  }
22320
22422
  return B.map(
22321
- (x) => new DenseMatrix4({
22423
+ (x) => new DenseMatrix5({
22322
22424
  data: x.map((e) => [e]),
22323
22425
  size: [rows, 1]
22324
22426
  })
@@ -22377,7 +22479,7 @@ var createLsolveAll = /* @__PURE__ */ factory(
22377
22479
  }
22378
22480
  }
22379
22481
  return B.map(
22380
- (x) => new DenseMatrix4({
22482
+ (x) => new DenseMatrix5({
22381
22483
  data: x.map((e) => [e]),
22382
22484
  size: [rows, 1]
22383
22485
  })
@@ -22443,10 +22545,10 @@ var createUsolve = /* @__PURE__ */ factory(
22443
22545
  multiplyScalar: multiplyScalar2,
22444
22546
  subtractScalar: subtractScalar2,
22445
22547
  equalScalar: equalScalar3,
22446
- DenseMatrix: DenseMatrix4
22548
+ DenseMatrix: DenseMatrix5
22447
22549
  }) => {
22448
22550
  const solveValidation = createSolveValidation({
22449
- DenseMatrix: DenseMatrix4
22551
+ DenseMatrix: DenseMatrix5
22450
22552
  });
22451
22553
  return typed3(name103, {
22452
22554
  "SparseMatrix, Array | Matrix": function(m, b) {
@@ -22484,7 +22586,7 @@ var createUsolve = /* @__PURE__ */ factory(
22484
22586
  for (let i = 0; i < rows; i++) {
22485
22587
  x2[i] = [resultAlloc.array[i]];
22486
22588
  }
22487
- return new DenseMatrix4({
22589
+ return new DenseMatrix5({
22488
22590
  data: x2,
22489
22591
  size: [rows, 1]
22490
22592
  });
@@ -22517,7 +22619,7 @@ var createUsolve = /* @__PURE__ */ factory(
22517
22619
  }
22518
22620
  x[j] = [xj];
22519
22621
  }
22520
- return new DenseMatrix4({
22622
+ return new DenseMatrix5({
22521
22623
  data: x,
22522
22624
  size: [rows, 1]
22523
22625
  });
@@ -22561,7 +22663,7 @@ var createUsolve = /* @__PURE__ */ factory(
22561
22663
  x[j] = [0];
22562
22664
  }
22563
22665
  }
22564
- return new DenseMatrix4({
22666
+ return new DenseMatrix5({
22565
22667
  data: x,
22566
22668
  size: [rows, 1]
22567
22669
  });
@@ -22590,10 +22692,10 @@ var createUsolveAll = /* @__PURE__ */ factory(
22590
22692
  multiplyScalar: multiplyScalar2,
22591
22693
  subtractScalar: subtractScalar2,
22592
22694
  equalScalar: equalScalar3,
22593
- DenseMatrix: DenseMatrix4
22695
+ DenseMatrix: DenseMatrix5
22594
22696
  }) => {
22595
22697
  const solveValidation = createSolveValidation({
22596
- DenseMatrix: DenseMatrix4
22698
+ DenseMatrix: DenseMatrix5
22597
22699
  });
22598
22700
  return typed3(name104, {
22599
22701
  "SparseMatrix, Array | Matrix": function(m, b) {
@@ -22643,7 +22745,7 @@ var createUsolveAll = /* @__PURE__ */ factory(
22643
22745
  }
22644
22746
  }
22645
22747
  return B.map(
22646
- (x) => new DenseMatrix4({
22748
+ (x) => new DenseMatrix5({
22647
22749
  data: x.map((e) => [e]),
22648
22750
  size: [rows, 1]
22649
22751
  })
@@ -22702,7 +22804,7 @@ var createUsolveAll = /* @__PURE__ */ factory(
22702
22804
  }
22703
22805
  }
22704
22806
  return B.map(
22705
- (x) => new DenseMatrix4({
22807
+ (x) => new DenseMatrix5({
22706
22808
  data: x.map((e) => [e]),
22707
22809
  size: [rows, 1]
22708
22810
  })
@@ -22920,7 +23022,7 @@ var dependencies108 = ["typed", "DenseMatrix"];
22920
23022
  var createMatAlgo12xSfs = /* @__PURE__ */ factory(
22921
23023
  name108,
22922
23024
  dependencies108,
22923
- ({ typed: typed3, DenseMatrix: DenseMatrix4 }) => {
23025
+ ({ typed: typed3, DenseMatrix: DenseMatrix5 }) => {
22924
23026
  return function matAlgo12xSfs(s, b, callback, inverse) {
22925
23027
  const avalues = s._values;
22926
23028
  const aindex = s._index;
@@ -22960,7 +23062,7 @@ var createMatAlgo12xSfs = /* @__PURE__ */ factory(
22960
23062
  }
22961
23063
  }
22962
23064
  }
22963
- return new DenseMatrix4({
23065
+ return new DenseMatrix5({
22964
23066
  data: cdata,
22965
23067
  size: [rows, columns],
22966
23068
  datatype: dt
@@ -23033,9 +23135,9 @@ var dependencies110 = [
23033
23135
  var createRound = /* @__PURE__ */ factory(
23034
23136
  name110,
23035
23137
  dependencies110,
23036
- ({ typed: typed3, config, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, BigNumber: BigNumber9, DenseMatrix: DenseMatrix4 }) => {
23138
+ ({ typed: typed3, config, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, BigNumber: BigNumber9, DenseMatrix: DenseMatrix5 }) => {
23037
23139
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
23038
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
23140
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
23039
23141
  const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
23040
23142
  function toExponent(epsilon) {
23041
23143
  return Math.abs(splitNumber(epsilon).exponent);
@@ -24444,13 +24546,13 @@ var createEqual = /* @__PURE__ */ factory(
24444
24546
  typed: typed3,
24445
24547
  matrix: matrix2,
24446
24548
  equalScalar: equalScalar3,
24447
- DenseMatrix: DenseMatrix4,
24549
+ DenseMatrix: DenseMatrix5,
24448
24550
  concat: _concat,
24449
24551
  SparseMatrix
24450
24552
  }) => {
24451
24553
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
24452
24554
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
24453
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
24555
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
24454
24556
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({ typed: typed3, matrix: matrix2 });
24455
24557
  return typed3(
24456
24558
  name124,
@@ -24840,7 +24942,7 @@ var createDotDivide = /* @__PURE__ */ factory(
24840
24942
  matrix: matrix2,
24841
24943
  equalScalar: equalScalar3,
24842
24944
  divideScalar: divideScalar2,
24843
- DenseMatrix: DenseMatrix4,
24945
+ DenseMatrix: DenseMatrix5,
24844
24946
  concat: concat3,
24845
24947
  SparseMatrix
24846
24948
  }) => {
@@ -24848,7 +24950,7 @@ var createDotDivide = /* @__PURE__ */ factory(
24848
24950
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
24849
24951
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
24850
24952
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
24851
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
24953
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
24852
24954
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
24853
24955
  typed: typed3,
24854
24956
  matrix: matrix2,
@@ -25017,9 +25119,9 @@ var createFloorNumber = /* @__PURE__ */ factory(
25017
25119
  var createFloor = /* @__PURE__ */ factory(
25018
25120
  name133,
25019
25121
  dependencies133,
25020
- ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4 }) => {
25122
+ ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5 }) => {
25021
25123
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
25022
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
25124
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
25023
25125
  const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
25024
25126
  const floorNumber = createFloorNumber({ typed: typed3, config, round: round2 });
25025
25127
  function _bigFloor(x) {
@@ -25230,7 +25332,7 @@ var dependencies135 = [
25230
25332
  var createMod = /* @__PURE__ */ factory(
25231
25333
  name135,
25232
25334
  dependencies135,
25233
- ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
25335
+ ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
25234
25336
  const floor2 = createFloor({
25235
25337
  typed: typed3,
25236
25338
  config,
@@ -25238,13 +25340,13 @@ var createMod = /* @__PURE__ */ factory(
25238
25340
  matrix: matrix2,
25239
25341
  equalScalar: equalScalar3,
25240
25342
  zeros: zeros3,
25241
- DenseMatrix: DenseMatrix4
25343
+ DenseMatrix: DenseMatrix5
25242
25344
  });
25243
25345
  const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
25244
25346
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
25245
25347
  const matAlgo05xSfSf = createMatAlgo05xSfSf({ typed: typed3, equalScalar: equalScalar3 });
25246
25348
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
25247
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
25349
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
25248
25350
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
25249
25351
  typed: typed3,
25250
25352
  matrix: matrix2,
@@ -25457,7 +25559,7 @@ var dependencies138 = ["typed", "DenseMatrix"];
25457
25559
  var createMatAlgo10xSids = /* @__PURE__ */ factory(
25458
25560
  name138,
25459
25561
  dependencies138,
25460
- ({ typed: typed3, DenseMatrix: DenseMatrix4 }) => {
25562
+ ({ typed: typed3, DenseMatrix: DenseMatrix5 }) => {
25461
25563
  return function matAlgo10xSids(s, b, callback, inverse) {
25462
25564
  const avalues = s._values;
25463
25565
  const aindex = s._index;
@@ -25497,7 +25599,7 @@ var createMatAlgo10xSids = /* @__PURE__ */ factory(
25497
25599
  }
25498
25600
  }
25499
25601
  }
25500
- return new DenseMatrix4({
25602
+ return new DenseMatrix5({
25501
25603
  data: cdata,
25502
25604
  size: [rows, columns],
25503
25605
  datatype: dt
@@ -25536,7 +25638,7 @@ var createGcd = /* @__PURE__ */ factory(
25536
25638
  equalScalar: equalScalar3,
25537
25639
  zeros: zeros3,
25538
25640
  BigNumber: BigNumber9,
25539
- DenseMatrix: DenseMatrix4,
25641
+ DenseMatrix: DenseMatrix5,
25540
25642
  concat: concat3
25541
25643
  }) => {
25542
25644
  const mod2 = createMod({
@@ -25546,12 +25648,12 @@ var createGcd = /* @__PURE__ */ factory(
25546
25648
  matrix: matrix2,
25547
25649
  equalScalar: equalScalar3,
25548
25650
  zeros: zeros3,
25549
- DenseMatrix: DenseMatrix4,
25651
+ DenseMatrix: DenseMatrix5,
25550
25652
  concat: concat3
25551
25653
  });
25552
25654
  const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
25553
25655
  const matAlgo04xSidSid = createMatAlgo04xSidSid({ typed: typed3, equalScalar: equalScalar3 });
25554
- const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
25656
+ const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
25555
25657
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
25556
25658
  typed: typed3,
25557
25659
  matrix: matrix2,
@@ -26118,9 +26220,9 @@ var createCeilNumber = /* @__PURE__ */ factory(
26118
26220
  var createCeil = /* @__PURE__ */ factory(
26119
26221
  name145,
26120
26222
  dependencies145,
26121
- ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4 }) => {
26223
+ ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5 }) => {
26122
26224
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
26123
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
26225
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
26124
26226
  const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
26125
26227
  const ceilNumber = createCeilNumber({ typed: typed3, config, round: round2 });
26126
26228
  function _bigCeil(x) {
@@ -26254,10 +26356,10 @@ var dependencies147 = ["typed", "matrix", "equalScalar", "DenseMatrix", "concat"
26254
26356
  var createBitOr = /* @__PURE__ */ factory(
26255
26357
  name147,
26256
26358
  dependencies147,
26257
- ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
26359
+ ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
26258
26360
  const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
26259
26361
  const matAlgo04xSidSid = createMatAlgo04xSidSid({ typed: typed3, equalScalar: equalScalar3 });
26260
- const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
26362
+ const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
26261
26363
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
26262
26364
  typed: typed3,
26263
26365
  matrix: matrix2,
@@ -26285,10 +26387,10 @@ var dependencies148 = ["typed", "matrix", "DenseMatrix", "concat", "SparseMatrix
26285
26387
  var createBitXor = /* @__PURE__ */ factory(
26286
26388
  name148,
26287
26389
  dependencies148,
26288
- ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
26390
+ ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
26289
26391
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
26290
26392
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
26291
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
26393
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
26292
26394
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
26293
26395
  typed: typed3,
26294
26396
  matrix: matrix2,
@@ -26429,11 +26531,11 @@ var dependencies150 = ["typed", "matrix", "equalScalar", "zeros", "DenseMatrix",
26429
26531
  var createLeftShift = /* @__PURE__ */ factory(
26430
26532
  name150,
26431
26533
  dependencies150,
26432
- ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
26534
+ ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
26433
26535
  const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
26434
26536
  const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
26435
26537
  const matAlgo08xS0Sid = createMatAlgo08xS0Sid({ typed: typed3, equalScalar: equalScalar3 });
26436
- const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
26538
+ const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
26437
26539
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
26438
26540
  const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
26439
26541
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
@@ -26500,11 +26602,11 @@ var dependencies151 = ["typed", "matrix", "equalScalar", "zeros", "DenseMatrix",
26500
26602
  var createRightArithShift = /* @__PURE__ */ factory(
26501
26603
  name151,
26502
26604
  dependencies151,
26503
- ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
26605
+ ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
26504
26606
  const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
26505
26607
  const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
26506
26608
  const matAlgo08xS0Sid = createMatAlgo08xS0Sid({ typed: typed3, equalScalar: equalScalar3 });
26507
- const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
26609
+ const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
26508
26610
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
26509
26611
  const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
26510
26612
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
@@ -26571,11 +26673,11 @@ var dependencies152 = ["typed", "matrix", "equalScalar", "zeros", "DenseMatrix",
26571
26673
  var createRightLogShift = /* @__PURE__ */ factory(
26572
26674
  name152,
26573
26675
  dependencies152,
26574
- ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
26676
+ ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
26575
26677
  const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
26576
26678
  const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
26577
26679
  const matAlgo08xS0Sid = createMatAlgo08xS0Sid({ typed: typed3, equalScalar: equalScalar3 });
26578
- const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
26680
+ const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
26579
26681
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
26580
26682
  const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
26581
26683
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
@@ -26641,10 +26743,10 @@ var dependencies153 = ["typed", "matrix", "equalScalar", "DenseMatrix", "concat"
26641
26743
  var createOr = /* @__PURE__ */ factory(
26642
26744
  name153,
26643
26745
  dependencies153,
26644
- ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
26746
+ ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
26645
26747
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
26646
26748
  const matAlgo05xSfSf = createMatAlgo05xSfSf({ typed: typed3, equalScalar: equalScalar3 });
26647
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
26749
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
26648
26750
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
26649
26751
  typed: typed3,
26650
26752
  matrix: matrix2,
@@ -26680,10 +26782,10 @@ var dependencies154 = ["typed", "matrix", "DenseMatrix", "concat", "SparseMatrix
26680
26782
  var createXor = /* @__PURE__ */ factory(
26681
26783
  name154,
26682
26784
  dependencies154,
26683
- ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
26785
+ ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
26684
26786
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
26685
26787
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
26686
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
26788
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
26687
26789
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
26688
26790
  typed: typed3,
26689
26791
  matrix: matrix2,
@@ -26870,12 +26972,12 @@ var createCompare = /* @__PURE__ */ factory(
26870
26972
  matrix: matrix2,
26871
26973
  BigNumber: BigNumber9,
26872
26974
  Fraction: Fraction5,
26873
- DenseMatrix: DenseMatrix4,
26975
+ DenseMatrix: DenseMatrix5,
26874
26976
  concat: concat3
26875
26977
  }) => {
26876
26978
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
26877
26979
  const matAlgo05xSfSf = createMatAlgo05xSfSf({ typed: typed3, equalScalar: equalScalar3 });
26878
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
26980
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
26879
26981
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
26880
26982
  typed: typed3,
26881
26983
  matrix: matrix2,
@@ -27000,10 +27102,10 @@ var dependencies159 = [
27000
27102
  var createLarger = /* @__PURE__ */ factory(
27001
27103
  name159,
27002
27104
  dependencies159,
27003
- ({ typed: typed3, config, bignumber: bignumber2, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
27105
+ ({ typed: typed3, config, bignumber: bignumber2, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
27004
27106
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
27005
27107
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
27006
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
27108
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
27007
27109
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
27008
27110
  typed: typed3,
27009
27111
  matrix: matrix2,
@@ -27058,10 +27160,10 @@ var dependencies160 = ["typed", "config", "matrix", "DenseMatrix", "concat", "Sp
27058
27160
  var createLargerEq = /* @__PURE__ */ factory(
27059
27161
  name160,
27060
27162
  dependencies160,
27061
- ({ typed: typed3, config, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
27163
+ ({ typed: typed3, config, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
27062
27164
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
27063
27165
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
27064
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
27166
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
27065
27167
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
27066
27168
  typed: typed3,
27067
27169
  matrix: matrix2,
@@ -27124,13 +27226,13 @@ var createSmaller = /* @__PURE__ */ factory(
27124
27226
  config,
27125
27227
  bignumber: bignumber2,
27126
27228
  matrix: matrix2,
27127
- DenseMatrix: DenseMatrix4,
27229
+ DenseMatrix: DenseMatrix5,
27128
27230
  concat: concat3,
27129
27231
  SparseMatrix
27130
27232
  }) => {
27131
27233
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
27132
27234
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
27133
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
27235
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
27134
27236
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
27135
27237
  typed: typed3,
27136
27238
  matrix: matrix2,
@@ -27185,10 +27287,10 @@ var dependencies162 = ["typed", "config", "matrix", "DenseMatrix", "concat", "Sp
27185
27287
  var createSmallerEq = /* @__PURE__ */ factory(
27186
27288
  name162,
27187
27289
  dependencies162,
27188
- ({ typed: typed3, config, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
27290
+ ({ typed: typed3, config, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
27189
27291
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
27190
27292
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
27191
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
27293
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
27192
27294
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
27193
27295
  typed: typed3,
27194
27296
  matrix: matrix2,
@@ -27249,13 +27351,13 @@ var createUnequal = /* @__PURE__ */ factory(
27249
27351
  config: _config,
27250
27352
  equalScalar: equalScalar3,
27251
27353
  matrix: matrix2,
27252
- DenseMatrix: DenseMatrix4,
27354
+ DenseMatrix: DenseMatrix5,
27253
27355
  concat: concat3,
27254
27356
  SparseMatrix
27255
27357
  }) => {
27256
27358
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
27257
27359
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
27258
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
27360
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
27259
27361
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
27260
27362
  typed: typed3,
27261
27363
  matrix: matrix2,
@@ -27306,12 +27408,12 @@ var dependencies164 = ["typed", "matrix", "equalScalar", "BigNumber", "DenseMatr
27306
27408
  var createAtan2 = /* @__PURE__ */ factory(
27307
27409
  name164,
27308
27410
  dependencies164,
27309
- ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, BigNumber: BigNumber9, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
27411
+ ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, BigNumber: BigNumber9, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
27310
27412
  const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
27311
27413
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
27312
27414
  const matAlgo09xS0Sf = createMatAlgo09xS0Sf({ typed: typed3, equalScalar: equalScalar3 });
27313
27415
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
27314
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
27416
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
27315
27417
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
27316
27418
  typed: typed3,
27317
27419
  matrix: matrix2,
@@ -27380,14 +27482,14 @@ var createDotPow = /* @__PURE__ */ factory(
27380
27482
  equalScalar: equalScalar3,
27381
27483
  matrix: matrix2,
27382
27484
  pow: pow2,
27383
- DenseMatrix: DenseMatrix4,
27485
+ DenseMatrix: DenseMatrix5,
27384
27486
  concat: concat3,
27385
27487
  SparseMatrix
27386
27488
  }) => {
27387
27489
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
27388
27490
  const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
27389
27491
  const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
27390
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
27492
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
27391
27493
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
27392
27494
  typed: typed3,
27393
27495
  matrix: matrix2,
@@ -27444,8 +27546,8 @@ var createFixNumber = /* @__PURE__ */ factory(
27444
27546
  var createFix = /* @__PURE__ */ factory(
27445
27547
  name167,
27446
27548
  dependencies167,
27447
- ({ typed: typed3, Complex: Complex10, matrix: matrix2, ceil: ceil2, floor: floor2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4 }) => {
27448
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
27549
+ ({ typed: typed3, Complex: Complex10, matrix: matrix2, ceil: ceil2, floor: floor2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5 }) => {
27550
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
27449
27551
  const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
27450
27552
  const fixNumber = createFixNumber({ typed: typed3, ceil: ceil2, floor: floor2 });
27451
27553
  return typed3("fix", {
@@ -30709,15 +30811,15 @@ var createSubtract = /* @__PURE__ */ factory(
30709
30811
  equalScalar: equalScalar3,
30710
30812
  subtractScalar: subtractScalar2,
30711
30813
  unaryMinus: _unaryMinus,
30712
- DenseMatrix: DenseMatrix4,
30814
+ DenseMatrix: DenseMatrix5,
30713
30815
  concat: concat3,
30714
30816
  nodeOperations: nodeOperations2
30715
30817
  }) => {
30716
30818
  const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
30717
30819
  const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
30718
30820
  const matAlgo05xSfSf = createMatAlgo05xSfSf({ typed: typed3, equalScalar: equalScalar3 });
30719
- const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
30720
- const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
30821
+ const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
30822
+ const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
30721
30823
  const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
30722
30824
  typed: typed3,
30723
30825
  matrix: matrix2,
@@ -34457,8 +34559,8 @@ var dependencies207 = ["smaller", "DenseMatrix"];
34457
34559
  var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
34458
34560
  name207,
34459
34561
  dependencies207,
34460
- ({ smaller: smaller2, DenseMatrix: DenseMatrix4 }) => {
34461
- class ImmutableDenseMatrix extends DenseMatrix4 {
34562
+ ({ smaller: smaller2, DenseMatrix: DenseMatrix5 }) => {
34563
+ class ImmutableDenseMatrix extends DenseMatrix5 {
34462
34564
  /**
34463
34565
  * Type identifier
34464
34566
  */
@@ -34498,7 +34600,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
34498
34600
  throw new Error("Invalid datatype: " + datatype);
34499
34601
  }
34500
34602
  if (isMatrix(data) || isArray(data)) {
34501
- const matrix2 = new DenseMatrix4(data, datatype);
34603
+ const matrix2 = new DenseMatrix5(data, datatype);
34502
34604
  this._data = matrix2._data;
34503
34605
  this._size = matrix2._size;
34504
34606
  this._datatype = matrix2._datatype;
@@ -34537,7 +34639,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
34537
34639
  subset(index, _replacement, _defaultValue) {
34538
34640
  switch (arguments.length) {
34539
34641
  case 1: {
34540
- const m = DenseMatrix4.prototype.subset.call(this, index);
34642
+ const m = DenseMatrix5.prototype.subset.call(this, index);
34541
34643
  if (isMatrix(m)) {
34542
34644
  const mm = m;
34543
34645
  return new ImmutableDenseMatrix({
@@ -34649,7 +34751,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
34649
34751
  if (this._min === null) {
34650
34752
  let m = null;
34651
34753
  const smallerFn = smaller2;
34652
- DenseMatrix4.prototype.forEach.call(this, function(v) {
34754
+ DenseMatrix5.prototype.forEach.call(this, function(v) {
34653
34755
  if (m === null || smallerFn(v, m)) {
34654
34756
  m = v;
34655
34757
  }
@@ -34666,7 +34768,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
34666
34768
  if (this._max === null) {
34667
34769
  let m = null;
34668
34770
  const smallerFn = smaller2;
34669
- DenseMatrix4.prototype.forEach.call(this, function(v) {
34771
+ DenseMatrix5.prototype.forEach.call(this, function(v) {
34670
34772
  if (m === null || smallerFn(m, v)) {
34671
34773
  m = v;
34672
34774
  }
@@ -34676,7 +34778,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
34676
34778
  return this._max ?? void 0;
34677
34779
  }
34678
34780
  }
34679
- Object.setPrototypeOf(ImmutableDenseMatrix.prototype, DenseMatrix4.prototype);
34781
+ Object.setPrototypeOf(ImmutableDenseMatrix.prototype, DenseMatrix5.prototype);
34680
34782
  const proto = ImmutableDenseMatrix.prototype;
34681
34783
  proto.constructor = ImmutableDenseMatrix;
34682
34784
  proto.type = "ImmutableDenseMatrix";
@@ -35056,7 +35158,7 @@ var createLup = /* @__PURE__ */ factory(
35056
35158
  larger: larger2,
35057
35159
  equalScalar: equalScalar3,
35058
35160
  unaryMinus: unaryMinus2,
35059
- DenseMatrix: DenseMatrix4,
35161
+ DenseMatrix: DenseMatrix5,
35060
35162
  SparseMatrix,
35061
35163
  Spa
35062
35164
  }) => {
@@ -35115,11 +35217,11 @@ var createLup = /* @__PURE__ */ factory(
35115
35217
  for (let i2 = 0; i2 < rows; i2++) {
35116
35218
  p2[i2] = permAlloc.array[i2];
35117
35219
  }
35118
- const l2 = new DenseMatrix4({
35220
+ const l2 = new DenseMatrix5({
35119
35221
  data: ldata2,
35120
35222
  size: [rows, n]
35121
35223
  });
35122
- const u2 = new DenseMatrix4({
35224
+ const u2 = new DenseMatrix5({
35123
35225
  data: udata2,
35124
35226
  size: [n, columns]
35125
35227
  });
@@ -35176,7 +35278,7 @@ var createLup = /* @__PURE__ */ factory(
35176
35278
  }
35177
35279
  if (j !== pi) {
35178
35280
  p[j] = [p[pi], p[pi] = p[j]][0];
35179
- DenseMatrix4._swapRows(j, pi, data);
35281
+ DenseMatrix5._swapRows(j, pi, data);
35180
35282
  }
35181
35283
  if (j < rows) {
35182
35284
  for (i = j + 1; i < rows; i++) {
@@ -35221,11 +35323,11 @@ var createLup = /* @__PURE__ */ factory(
35221
35323
  }
35222
35324
  }
35223
35325
  }
35224
- const l = new DenseMatrix4({
35326
+ const l = new DenseMatrix5({
35225
35327
  data: ldata,
35226
35328
  size: lsize
35227
35329
  });
35228
- const u = new DenseMatrix4({
35330
+ const u = new DenseMatrix5({
35229
35331
  data: udata,
35230
35332
  size: usize
35231
35333
  });
@@ -36763,7 +36865,7 @@ var dependencies227 = ["typed", "size", "subset", "compareNatural", "Index", "De
36763
36865
  var createSetCartesian = /* @__PURE__ */ factory(
36764
36866
  name227,
36765
36867
  dependencies227,
36766
- ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix4 }) => {
36868
+ ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix5 }) => {
36767
36869
  return typed3(name227, {
36768
36870
  "Array | Matrix, Array | Matrix": function(a1, a2) {
36769
36871
  let result = [];
@@ -36780,7 +36882,7 @@ var createSetCartesian = /* @__PURE__ */ factory(
36780
36882
  if (Array.isArray(a1) && Array.isArray(a2)) {
36781
36883
  return result;
36782
36884
  }
36783
- return new DenseMatrix4(result);
36885
+ return new DenseMatrix5(result);
36784
36886
  }
36785
36887
  });
36786
36888
  }
@@ -36792,7 +36894,7 @@ var dependencies228 = ["typed", "size", "subset", "compareNatural", "Index", "De
36792
36894
  var createSetDifference = /* @__PURE__ */ factory(
36793
36895
  name228,
36794
36896
  dependencies228,
36795
- ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix4 }) => {
36897
+ ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix5 }) => {
36796
36898
  return typed3(name228, {
36797
36899
  "Array | Matrix, Array | Matrix": function(a1, a2) {
36798
36900
  let result;
@@ -36821,7 +36923,7 @@ var createSetDifference = /* @__PURE__ */ factory(
36821
36923
  if (Array.isArray(a1) && Array.isArray(a2)) {
36822
36924
  return generalize(result);
36823
36925
  }
36824
- return new DenseMatrix4(generalize(result));
36926
+ return new DenseMatrix5(generalize(result));
36825
36927
  }
36826
36928
  });
36827
36929
  }
@@ -36833,7 +36935,7 @@ var dependencies229 = ["typed", "size", "subset", "compareNatural", "Index", "De
36833
36935
  var createSetDistinct = /* @__PURE__ */ factory(
36834
36936
  name229,
36835
36937
  dependencies229,
36836
- ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix4 }) => {
36938
+ ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix5 }) => {
36837
36939
  return typed3(name229, {
36838
36940
  "Array | Matrix": function(a) {
36839
36941
  let result;
@@ -36852,7 +36954,7 @@ var createSetDistinct = /* @__PURE__ */ factory(
36852
36954
  if (Array.isArray(a)) {
36853
36955
  return result;
36854
36956
  }
36855
- return new DenseMatrix4(result);
36957
+ return new DenseMatrix5(result);
36856
36958
  }
36857
36959
  });
36858
36960
  }
@@ -36864,7 +36966,7 @@ var dependencies230 = ["typed", "size", "subset", "compareNatural", "Index", "De
36864
36966
  var createSetIntersect = /* @__PURE__ */ factory(
36865
36967
  name230,
36866
36968
  dependencies230,
36867
- ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix4 }) => {
36969
+ ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix5 }) => {
36868
36970
  return typed3(name230, {
36869
36971
  "Array | Matrix, Array | Matrix": function(a1, a2) {
36870
36972
  let result;
@@ -36886,7 +36988,7 @@ var createSetIntersect = /* @__PURE__ */ factory(
36886
36988
  if (Array.isArray(a1) && Array.isArray(a2)) {
36887
36989
  return generalize(result);
36888
36990
  }
36889
- return new DenseMatrix4(generalize(result));
36991
+ return new DenseMatrix5(generalize(result));
36890
36992
  }
36891
36993
  });
36892
36994
  }
@@ -39309,9 +39411,9 @@ var dependencies240 = ["typed", "matrix", "lup", "slu", "usolve", "lsolve", "Den
39309
39411
  var createLusolve = /* @__PURE__ */ factory(
39310
39412
  name240,
39311
39413
  dependencies240,
39312
- ({ typed: typed3, matrix: matrix2, lup: lup2, slu: slu2, usolve: usolve2, lsolve: lsolve2, DenseMatrix: DenseMatrix4 }) => {
39414
+ ({ typed: typed3, matrix: matrix2, lup: lup2, slu: slu2, usolve: usolve2, lsolve: lsolve2, DenseMatrix: DenseMatrix5 }) => {
39313
39415
  const solveValidation = createSolveValidation({
39314
- DenseMatrix: DenseMatrix4
39416
+ DenseMatrix: DenseMatrix5
39315
39417
  });
39316
39418
  return typed3(name240, {
39317
39419
  "Array, Array | Matrix": function(a, b) {
@@ -39394,7 +39496,7 @@ var createLusolve = /* @__PURE__ */ factory(
39394
39496
  for (let i = 0; i < rows; i++) {
39395
39497
  x[i] = [resultAlloc.array[i]];
39396
39498
  }
39397
- return new DenseMatrix4({
39499
+ return new DenseMatrix5({
39398
39500
  data: x,
39399
39501
  size: [rows, 1]
39400
39502
  });
@@ -41856,7 +41958,7 @@ var createRotationMatrix = /* @__PURE__ */ factory(
41856
41958
  norm: norm3,
41857
41959
  BigNumber: BigNumber9,
41858
41960
  matrix: matrix2,
41859
- DenseMatrix: DenseMatrix4,
41961
+ DenseMatrix: DenseMatrix5,
41860
41962
  SparseMatrix,
41861
41963
  cos: cos2,
41862
41964
  sin: sin2
@@ -41920,7 +42022,7 @@ var createRotationMatrix = /* @__PURE__ */ factory(
41920
42022
  return new SparseMatrix(data);
41921
42023
  }
41922
42024
  if (format6 === "dense") {
41923
- return new DenseMatrix4(data);
42025
+ return new DenseMatrix5(data);
41924
42026
  }
41925
42027
  throw new TypeError(`Unknown matrix type "${format6}"`);
41926
42028
  }
@@ -42687,8 +42789,9 @@ factoryScope.add = factoryScope.addScalar;
42687
42789
  var trace2 = createTrace(factoryScope);
42688
42790
  factoryScope.trace = trace2;
42689
42791
  factoryScope.multiply = factoryScope.multiplyScalar;
42690
- var det = createDet(factoryScope);
42691
- factoryScope.det = det;
42792
+ var factoryDet = createDet(factoryScope);
42793
+ factoryScope.det = factoryDet;
42794
+ var det = ((m) => acceleratedDet(m, factoryDet));
42692
42795
  factoryScope.isInteger = (x) => typeof x === "number" && Number.isInteger(x);
42693
42796
  var reshape2 = createReshape(factoryScope);
42694
42797
  factoryScope.reshape = reshape2;
@@ -42704,8 +42807,9 @@ var resize2 = createResize(factoryScope);
42704
42807
  factoryScope.resize = resize2;
42705
42808
  var subset = createSubset(factoryScope);
42706
42809
  factoryScope.subset = subset;
42707
- var inv = createInv(factoryScope);
42708
- factoryScope.inv = inv;
42810
+ var factoryInv = createInv(factoryScope);
42811
+ factoryScope.inv = factoryInv;
42812
+ var inv = ((m) => acceleratedInv(m, factoryInv));
42709
42813
  var factory_cbrt = createCbrt(factoryScope);
42710
42814
  factoryScope.cbrt = factory_cbrt;
42711
42815
  var nthRoots = createNthRoots(factoryScope);
@@ -44950,6 +45054,24 @@ function help(search) {
44950
45054
  }
44951
45055
  return new Help(doc);
44952
45056
  }
45057
+
45058
+ // src/grad-forward.ts
45059
+ import { Dual as Dual3 } from "@danielsimonjr/mathts-core";
45060
+ function derivativeAt(fn, x0) {
45061
+ return fn(Dual3.variable(x0)).deriv;
45062
+ }
45063
+ function valueAndDerivativeAt(fn, x0) {
45064
+ const r = fn(Dual3.variable(x0));
45065
+ return { value: r.value, deriv: r.deriv };
45066
+ }
45067
+ function gradientAt(fn, x) {
45068
+ const grad = new Array(x.length);
45069
+ for (let k = 0; k < x.length; k++) {
45070
+ const duals = x.map((xi, i) => i === k ? Dual3.variable(xi) : Dual3.constant(xi));
45071
+ grad[k] = fn(duals).deriv;
45072
+ }
45073
+ return grad;
45074
+ }
44953
45075
  export {
44954
45076
  ARRAY_WORKER_THRESHOLD,
44955
45077
  CAS_BATCH_THRESHOLD,
@@ -45109,6 +45231,7 @@ export {
45109
45231
  degree,
45110
45232
  delaunayTriangulation,
45111
45233
  derivative,
45234
+ derivativeAt,
45112
45235
  det,
45113
45236
  deuteronMass,
45114
45237
  diag,
@@ -45279,6 +45402,7 @@ export {
45279
45402
  gpuMatmul,
45280
45403
  gpuScale,
45281
45404
  gpuTranspose,
45405
+ gradientAt,
45282
45406
  gradientSymbolic,
45283
45407
  graphDistance,
45284
45408
  gravitationConstant,
@@ -45678,6 +45802,7 @@ export {
45678
45802
  usolveAll,
45679
45803
  vacuumImpedance,
45680
45804
  validateClosureSource,
45805
+ valueAndDerivativeAt,
45681
45806
  variables,
45682
45807
  variance,
45683
45808
  voronoiDiagram,
@@ -0,0 +1,21 @@
1
+ /** Below this dimension the factory path is already sub-millisecond; native
2
+ * conversion overhead isn't worth it and the well-tested factory path runs. */
3
+ export declare const NATIVE_MATRIX_THRESHOLD = 8;
4
+ /** True for a square `number[][]` of side ≥ NATIVE_MATRIX_THRESHOLD. */
5
+ export declare function isLargeNumericSquare(a: unknown): a is number[][];
6
+ /**
7
+ * Determinant via native Float64Array LU: det = parity(P) · ∏ diag(U).
8
+ * Falls back to `factoryDet` for non-(large numeric square) inputs or any
9
+ * non-singular native error; a singular matrix correctly yields 0.
10
+ */
11
+ export declare function acceleratedDet<T>(a: unknown, factoryDet: (m: unknown) => T): T;
12
+ /**
13
+ * Inverse via native Float64Array LU + forward/back substitution: solve A·xⱼ = eⱼ
14
+ * for each column. Falls back to `factoryInv` for non-(large numeric square)
15
+ * inputs; on a singular matrix the native LU throws and we delegate to the
16
+ * factory, which raises the proper "Cannot calculate inverse, determinant is
17
+ * zero" error — so behaviour is identical to the factory, just faster.
18
+ * Verified against numpy: max|inv−numpy| and max|A·inv−I| ~1e-15.
19
+ */
20
+ export declare function acceleratedInv<T>(a: unknown, factoryInv: (m: unknown) => T): T;
21
+ //# sourceMappingURL=native-accel.d.ts.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"native-accel.d.ts","sourceRoot":"","sources":["../../src/matrix/native-accel.ts"],"names":[],"mappings":"AAkBA;gFACgF;AAChF,eAAO,MAAM,uBAAuB,IAAI,CAAC;AAEzC,wEAAwE;AACxE,wBAAgB,oBAAoB,CAAC,CAAC,EAAE,OAAO,GAAG,CAAC,IAAI,MAAM,EAAE,EAAE,CAQhE;AAsBD;;;;GAIG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,EAAE,UAAU,EAAE,CAAC,CAAC,EAAE,OAAO,KAAK,CAAC,GAAG,CAAC,CAY9E;AAED;;;;;;;GAOG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,EAAE,UAAU,EAAE,CAAC,CAAC,EAAE,OAAO,KAAK,CAAC,GAAG,CAAC,CAoC9E"}
@@ -1 +1 @@
1
- {"version":3,"file":"arithmetic.d.ts","sourceRoot":"","sources":["../../src/typed/arithmetic.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAKH,OAAO,EAAe,WAAW,EAAE,MAAM,gCAAgC,CAAC;AAU1E,4BAA4B;AAC5B,KAAK,GAAG,GAAG,MAAM,CAAC;AASlB;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,GAAG,wCA2Cd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCA4BnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAoEnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,MAAM,wCA6BjB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,UAAU,wCAYrB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAOpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAqBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAQd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAYjB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAIlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAYd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAad,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAiBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAoBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAgBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAKd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCA+Bd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,IAAI,wCA6Bf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAoBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAmBH,eAAO,MAAM,KAAK,wCAOhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAMlB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAMjB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAMnB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAMlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;;;;GAKG;AACH,MAAM,MAAM,gBAAgB,GAAG,UAAU,GAAG,aAAa,GAAG,QAAQ,CAAC;AA6BrE;;;;;;GAMG;AACH,eAAO,MAAM,QAAQ,wCAenB,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAMH;;;;;;;GAOG;AACH,wBAAsB,MAAM,CAC1B,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,EACV,KAAK,EAAE,GAAG,EACV,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,GACT,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,SAAS,CAAC,IAAI,EAAE,YAAY,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,GAAG,OAAO,CAAC,YAAY,CAAC,CAG/F;AAED;;GAEG;AACH,wBAAsB,MAAM,CAC1B,MAAM,EAAE,YAAY,EACpB,IAAI,EAAE,GAAG,EACT,IAAI,EAAE,GAAG,EACT,MAAM,EAAE,YAAY,GACnB,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,KAAK,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,OAAO,CAAC,YAAY,CAAC,CAGnF;AAMD;;GAEG;AACH,wBAAsB,cAAc,IAAI,OAAO,CAAC,IAAI,CAAC,CAEpD;AAED;;GAEG;AACH,wBAAsB,aAAa,IAAI,OAAO,CAAC,IAAI,CAAC,CAEnD;AAED;;GAEG;AACH,wBAAgB,iBAAiB,CAAC,YAAY,EAAE,GAAG,GAAG,OAAO,CAE5D;AAED;;GAEG;AACH,wBAAgB,cAAc,IAAI,WAAW,CAE5C;AAMD,eAAO,MAAM,eAAe;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA+E3B,CAAC"}
1
+ {"version":3,"file":"arithmetic.d.ts","sourceRoot":"","sources":["../../src/typed/arithmetic.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAKH,OAAO,EAAe,WAAW,EAAE,MAAM,gCAAgC,CAAC;AAU1E,4BAA4B;AAC5B,KAAK,GAAG,GAAG,MAAM,CAAC;AASlB;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,GAAG,wCAgDd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAiCnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAyEnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAkCjB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,UAAU,wCAYrB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAOpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAed,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAqBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAWd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAajB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAkBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAIlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAad,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAiBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAoBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAgBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAKd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCA+Bd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAgBd,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,IAAI,wCA6Bf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAoBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAmBH,eAAO,MAAM,KAAK,wCAOhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAMlB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAMjB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAMpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAMnB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAMlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;;;;GAKG;AACH,MAAM,MAAM,gBAAgB,GAAG,UAAU,GAAG,aAAa,GAAG,QAAQ,CAAC;AA6BrE;;;;;;GAMG;AACH,eAAO,MAAM,QAAQ,wCAenB,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAMH;;;;;;;GAOG;AACH,wBAAsB,MAAM,CAC1B,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,EACV,KAAK,EAAE,GAAG,EACV,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,GACT,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,SAAS,CAAC,IAAI,EAAE,YAAY,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,GAAG,OAAO,CAAC,YAAY,CAAC,CAG/F;AAED;;GAEG;AACH,wBAAsB,MAAM,CAC1B,MAAM,EAAE,YAAY,EACpB,IAAI,EAAE,GAAG,EACT,IAAI,EAAE,GAAG,EACT,MAAM,EAAE,YAAY,GACnB,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,KAAK,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,OAAO,CAAC,YAAY,CAAC,CAGnF;AAMD;;GAEG;AACH,wBAAsB,cAAc,IAAI,OAAO,CAAC,IAAI,CAAC,CAEpD;AAED;;GAEG;AACH,wBAAsB,aAAa,IAAI,OAAO,CAAC,IAAI,CAAC,CAEnD;AAED;;GAEG;AACH,wBAAgB,iBAAiB,CAAC,YAAY,EAAE,GAAG,GAAG,OAAO,CAE5D;AAED;;GAEG;AACH,wBAAgB,cAAc,IAAI,WAAW,CAE5C;AAMD,eAAO,MAAM,eAAe;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA+E3B,CAAC"}
@@ -1 +1 @@
1
- {"version":3,"file":"trigonometry.d.ts","sourceRoot":"","sources":["../../src/typed/trigonometry.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAkBH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAYd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAYd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAYd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAqBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAoBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAGhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAKf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAKf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAKf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAiBhB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAEpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAEpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAKhB,CAAC;AAMH,eAAO,MAAM,iBAAiB;;;;;;;;;;;;;;;;;;;;CA2B7B,CAAC"}
1
+ {"version":3,"file":"trigonometry.d.ts","sourceRoot":"","sources":["../../src/typed/trigonometry.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAkBH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAad,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAad,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAad,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAqBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAoBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAGhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAKf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAKf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAKf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAiBhB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAEpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAEpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAKhB,CAAC;AAMH,eAAO,MAAM,iBAAiB;;;;;;;;;;;;;;;;;;;;CA2B7B,CAAC"}
package/package.json CHANGED
@@ -1,6 +1,6 @@
1
1
  {
2
2
  "name": "@danielsimonjr/mathts-functions",
3
- "version": "0.3.0",
3
+ "version": "0.4.0",
4
4
  "description": "Mathematical functions for MathTS - arithmetic, algebra, trigonometry, statistics, and more",
5
5
  "author": "Daniel Simon Jr.",
6
6
  "license": "MIT",
@@ -33,9 +33,9 @@
33
33
  "build:prod": "tsup src/index.ts --format esm --clean --minify --treeshake"
34
34
  },
35
35
  "dependencies": {
36
- "@danielsimonjr/mathts-core": "^0.1.5",
37
- "@danielsimonjr/mathts-expression": "^0.4.0",
38
- "@danielsimonjr/mathts-matrix": "^0.1.10",
36
+ "@danielsimonjr/mathts-core": "^0.2.0",
37
+ "@danielsimonjr/mathts-expression": "^0.4.1",
38
+ "@danielsimonjr/mathts-matrix": "^0.1.11",
39
39
  "@danielsimonjr/mathts-parallel": "^0.3.0",
40
40
  "bignumber.js": "^9.1.2",
41
41
  "complex.js": "^2.2.5",