@danielsimonjr/mathts-functions 0.48.0 → 0.49.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.
- package/dist/index.d.ts +1 -0
- package/dist/index.d.ts.map +1 -1
- package/dist/index.js +168 -8
- package/dist/special/wave-functions.d.ts +123 -0
- package/dist/special/wave-functions.d.ts.map +1 -0
- package/package.json +1 -1
package/dist/index.d.ts
CHANGED
|
@@ -106,6 +106,7 @@ export { jacobiSN, jacobiCN, jacobiDN } from './special/jacobi-elliptic.js';
|
|
|
106
106
|
export { rootsLegendre } from './numeric/gauss-nodes.js';
|
|
107
107
|
export type { RootsLegendreResult } from './numeric/gauss-nodes.js';
|
|
108
108
|
export { polylog, struveH, struveL, kelvinBer, kelvinBei, barnesG } from './special/niche.js';
|
|
109
|
+
export { siegelZ, riemannSiegelZ, lerchPhi, parabolicCylinderD, coulombF, } from './special/wave-functions.js';
|
|
109
110
|
export { bfs, dfs, floydWarshall, bellmanFord, closenessCentrality, harmonicCentrality, } from './graph/traversal-centrality.js';
|
|
110
111
|
export type { BellmanFordResult } from './graph/traversal-centrality.js';
|
|
111
112
|
export { maxFlow, minCut, astar, hungarian } from './graph/optimization.js';
|
package/dist/index.d.ts.map
CHANGED
|
@@ -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;AACrC,OAAO,EAAE,MAAM,EAAE,MAAM,iBAAiB,CAAC;AAOzC,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;AAMjC,OAAO,EACL,iBAAiB,EACjB,YAAY,EACZ,UAAU,EACV,WAAW,EACX,aAAa,EACb,mBAAmB,EACnB,eAAe,EACf,eAAe,EACf,eAAe,GAChB,MAAM,yBAAyB,CAAC;AAGjC,OAAO,EAAE,YAAY,EAAE,oBAAoB,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AACnF,YAAY,EAAE,MAAM,EAAE,MAAM,mBAAmB,CAAC;AAGhD,OAAO,EACL,KAAK,EACL,KAAK,EACL,MAAM,EACN,QAAQ,EACR,QAAQ,EACR,GAAG,EACH,GAAG,EACH,MAAM,EACN,GAAG,EACH,QAAQ,EACR,QAAQ,EACR,QAAQ,EACR,UAAU,EACV,UAAU,EACV,QAAQ,EACR,SAAS,EACT,UAAU,EACV,cAAc,EACd,GAAG,EACH,SAAS,EACT,WAAW,EACX,QAAQ,EACR,SAAS,GACV,MAAM,wBAAwB,CAAC;AAChC,YAAY,EACV,gBAAgB,EAChB,qBAAqB,EACrB,oBAAoB,EACpB,cAAc,EACd,eAAe,GAChB,MAAM,wBAAwB,CAAC;AAGhC,OAAO,EACL,KAAK,EACL,OAAO,EACP,SAAS,EACT,OAAO,EACP,OAAO,EACP,MAAM,EACN,MAAM,EACN,QAAQ,GACT,MAAM,oBAAoB,CAAC;AAG5B,OAAO,EACL,cAAc,EACd,WAAW,EACX,gBAAgB,EAChB,aAAa,EACb,kBAAkB,EAClB,IAAI,EACJ,SAAS,EACT,QAAQ,EACR,aAAa,EACb,OAAO,EACP,YAAY,EACZ,SAAS,EACT,SAAS,EACT,cAAc,EACd,UAAU,EACV,UAAU,EACV,eAAe,EACf,WAAW,EACX,WAAW,EACX,gBAAgB,GACjB,MAAM,6BAA6B,CAAC;AAGrC,OAAO,EACL,KAAK,EACL,UAAU,EACV,aAAa,EACb,QAAQ,EACR,WAAW,EACX,mBAAmB,EACnB,wBAAwB,EACxB,QAAQ,GACT,MAAM,uBAAuB,CAAC;AAC/B,YAAY,EACV,WAAW,EACX,gBAAgB,EAChB,aAAa,EACb,cAAc,EACd,iBAAiB,EACjB,eAAe,GAChB,MAAM,uBAAuB,CAAC;AAI/B,OAAO,EACL,IAAI,EACJ,IAAI,EACJ,MAAM,EACN,QAAQ,EACR,SAAS,EACT,SAAS,EACT,MAAM,EACN,eAAe,EACf,cAAc,EACd,EAAE,GACH,MAAM,mBAAmB,CAAC;AAG3B,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,qBAAqB,CAAC;AAKxD,OAAO,EAAE,GAAG,EAAE,MAAM,8BAA8B,CAAC;AACnD,YAAY,EAAE,SAAS,EAAE,UAAU,EAAE,MAAM,8BAA8B,CAAC;AAC1E,OAAO,EAAE,IAAI,EAAE,MAAM,uBAAuB,CAAC;AAC7C,YAAY,EAAE,WAAW,EAAE,MAAM,uBAAuB,CAAC;AAIzD,OAAO,EAAE,eAAe,EAAE,MAAM,+BAA+B,CAAC;AAChE,YAAY,EAAE,WAAW,EAAE,sBAAsB,EAAE,MAAM,+BAA+B,CAAC;AAIzF,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,gCAAgC,CAAC;AACxE,YAAY,EAAE,aAAa,EAAE,aAAa,EAAE,aAAa,EAAE,MAAM,gCAAgC,CAAC;AAIlG,OAAO,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,qBAAqB,CAAC;AACnD,YAAY,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AAMzD,OAAO,EAAE,EAAE,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,YAAY,EAAE,kBAAkB,EAAE,MAAM,qBAAqB,CAAC;AACpG,YAAY,EACV,mBAAmB,EACnB,cAAc,EACd,aAAa,EACb,YAAY,EACZ,YAAY,GACb,MAAM,qBAAqB,CAAC;AAM7B,OAAO,EAAE,KAAK,EAAE,MAAM,oBAAoB,CAAC;AAC3C,YAAY,EAAE,kBAAkB,EAAE,YAAY,EAAE,WAAW,EAAE,MAAM,oBAAoB,CAAC;AAKxF,OAAO,EAAE,IAAI,EAAE,MAAM,mBAAmB,CAAC;AACzC,YAAY,EAAE,WAAW,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AAOjE,OAAO,EAAE,WAAW,EAAE,WAAW,EAAE,aAAa,EAAE,GAAG,EAAE,MAAM,iCAAiC,CAAC;AAC/F,YAAY,EAAE,SAAS,EAAE,MAAM,iCAAiC,CAAC;AAQjE,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,UAAU,EAAE,UAAU,EAAE,MAAM,+BAA+B,CAAC;AACzF,YAAY,EACV,YAAY,EACZ,aAAa,EACb,qBAAqB,GACtB,MAAM,+BAA+B,CAAC;AAQvC,OAAO,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,gCAAgC,CAAC;AAKnE,OAAO,EAAE,cAAc,EAAE,MAAM,8BAA8B,CAAC;AAC9D,YAAY,EAAE,qBAAqB,EAAE,oBAAoB,EAAE,MAAM,8BAA8B,CAAC;AAKhG,OAAO,EAAE,IAAI,EAAE,MAAM,4BAA4B,CAAC;AAClD,YAAY,EAAE,WAAW,EAAE,UAAU,EAAE,MAAM,4BAA4B,CAAC;AAM1E,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,YAAY,EAAE,YAAY,EAAE,MAAM,sBAAsB,CAAC;AAOzD,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,sBAAsB,CAAC;AAC/D,YAAY,EAAE,iBAAiB,EAAE,YAAY,EAAE,MAAM,sBAAsB,CAAC;AAM5E,OAAO,EAAE,mBAAmB,EAAE,MAAM,0BAA0B,CAAC;AAC/D,YAAY,EAAE,KAAK,EAAE,iBAAiB,EAAE,gBAAgB,EAAE,MAAM,0BAA0B,CAAC;AAG3F,OAAO,EAAE,aAAa,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,MAAM,uBAAuB,CAAC;AAI1E,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,YAAY,EAAE,QAAQ,EAAE,MAAM,uBAAuB,CAAC;AAC/E,YAAY,EAAE,cAAc,EAAE,cAAc,EAAE,MAAM,uBAAuB,CAAC;AAG5E,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,YAAY,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AAG9D,OAAO,EAAE,GAAG,EAAE,MAAM,aAAa,CAAC;AAClC,YAAY,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAGzD,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,UAAU,EAAE,MAAM,gCAAgC,CAAC;AAC1E,YAAY,EACV,YAAY,EACZ,wBAAwB,EACxB,2BAA2B,GAC5B,MAAM,gCAAgC,CAAC;AAGxC,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AACjE,YAAY,EACV,yBAAyB,EACzB,wBAAwB,GACzB,MAAM,6BAA6B,CAAC;AAGrC,OAAO,EAAE,UAAU,EAAE,eAAe,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAC1F,YAAY,EAAE,cAAc,EAAE,QAAQ,EAAE,MAAM,yBAAyB,CAAC;AAGxE,OAAO,EAAE,MAAM,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAInE,OAAO,EAAE,MAAM,EAAE,WAAW,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAIrE,OAAO,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAC1C,YAAY,EAAE,kBAAkB,EAAE,iBAAiB,EAAE,MAAM,aAAa,CAAC;AAIzE,OAAO,EAAE,eAAe,EAAE,YAAY,EAAE,MAAM,4BAA4B,CAAC;AAC3E,YAAY,EACV,sBAAsB,EACtB,qBAAqB,EACrB,kBAAkB,GACnB,MAAM,4BAA4B,CAAC;AAIpC,OAAO,EAAE,eAAe,EAAE,MAAM,6BAA6B,CAAC;AAC9D,YAAY,EAAE,gBAAgB,EAAE,qBAAqB,EAAE,MAAM,6BAA6B,CAAC;AAM3F,OAAO,EAAE,IAAI,EAAE,MAAM,mBAAmB,CAAC;AACzC,YAAY,EAAE,WAAW,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AACjE,YAAY,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AAK1D,OAAO,EAAE,IAAI,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AACrD,YAAY,EACV,WAAW,EACX,UAAU,EACV,iBAAiB,EACjB,gBAAgB,GACjB,MAAM,mBAAmB,CAAC;AAG3B,OAAO,EAAE,gBAAgB,EAAE,MAAM,sBAAsB,CAAC;AAGxD,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAGxF,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,QAAQ,EAAE,SAAS,EAAE,OAAO,EAAE,QAAQ,EAAE,IAAI,EAAE,MAAM,yBAAyB,CAAC;AAIpG,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,MAAM,wBAAwB,CAAC;AAIpG,OAAO,EACL,cAAc,EACd,KAAK,EACL,KAAK,EACL,MAAM,EACN,MAAM,EACN,UAAU,GACX,MAAM,2BAA2B,CAAC;AACnC,YAAY,EAAE,gBAAgB,EAAE,MAAM,2BAA2B,CAAC;AAClE,YAAY,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AAI5D,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,MAAM,sBAAsB,CAAC;AAKnE,OAAO,EACL,SAAS,EACT,UAAU,EACV,GAAG,EACH,SAAS,EACT,IAAI,EACJ,KAAK,EACL,QAAQ,GACT,MAAM,4BAA4B,CAAC;AACpC,YAAY,EACV,gBAAgB,EAChB,UAAU,EACV,WAAW,EACX,UAAU,GACX,MAAM,4BAA4B,CAAC;AAGpC,OAAO,EACL,SAAS,EACT,eAAe,EACf,KAAK,EACL,kBAAkB,EAClB,mBAAmB,EACnB,mBAAmB,EACnB,uBAAuB,EACvB,gBAAgB,EAChB,0BAA0B,GAC3B,MAAM,qBAAqB,CAAC;AAI7B,OAAO,EACL,iBAAiB,EACjB,eAAe,EACf,iBAAiB,EACjB,aAAa,EACb,aAAa,EACb,aAAa,EACb,WAAW,EACX,UAAU,EACV,SAAS,EACT,YAAY,EACZ,aAAa,EACb,QAAQ,EACR,WAAW,GACZ,MAAM,8BAA8B,CAAC;AACtC,YAAY,EAAE,gBAAgB,EAAE,MAAM,8BAA8B,CAAC;AAGrE,OAAO,EACL,oBAAoB,EACpB,iBAAiB,EACjB,qBAAqB,GACtB,MAAM,2BAA2B,CAAC;AACnC,YAAY,EAAE,MAAM,EAAE,oBAAoB,EAAE,MAAM,2BAA2B,CAAC;AAU9E,OAAO,EACL,iBAAiB,EACjB,cAAc,EACd,cAAc,EACd,QAAQ,EACR,OAAO,EACP,OAAO,EACP,WAAW,EACX,OAAO,EACP,QAAQ,GACT,MAAM,6BAA6B,CAAC;AACrC,YAAY,EACV,eAAe,EACf,cAAc,EACd,aAAa,EACb,cAAc,GACf,MAAM,6BAA6B,CAAC;AAIrC,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,6BAA6B,CAAC;AAK1E,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,OAAO,EAAE,WAAW,EAAE,MAAM,kCAAkC,CAAC;AAM7F,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,8BAA8B,CAAC;AAC5E,OAAO,EAAE,aAAa,EAAE,MAAM,0BAA0B,CAAC;AACzD,YAAY,EAAE,mBAAmB,EAAE,MAAM,0BAA0B,CAAC;AAIpE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,SAAS,EAAE,SAAS,EAAE,OAAO,EAAE,MAAM,oBAAoB,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;AACrC,OAAO,EAAE,MAAM,EAAE,MAAM,iBAAiB,CAAC;AAOzC,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;AAMjC,OAAO,EACL,iBAAiB,EACjB,YAAY,EACZ,UAAU,EACV,WAAW,EACX,aAAa,EACb,mBAAmB,EACnB,eAAe,EACf,eAAe,EACf,eAAe,GAChB,MAAM,yBAAyB,CAAC;AAGjC,OAAO,EAAE,YAAY,EAAE,oBAAoB,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AACnF,YAAY,EAAE,MAAM,EAAE,MAAM,mBAAmB,CAAC;AAGhD,OAAO,EACL,KAAK,EACL,KAAK,EACL,MAAM,EACN,QAAQ,EACR,QAAQ,EACR,GAAG,EACH,GAAG,EACH,MAAM,EACN,GAAG,EACH,QAAQ,EACR,QAAQ,EACR,QAAQ,EACR,UAAU,EACV,UAAU,EACV,QAAQ,EACR,SAAS,EACT,UAAU,EACV,cAAc,EACd,GAAG,EACH,SAAS,EACT,WAAW,EACX,QAAQ,EACR,SAAS,GACV,MAAM,wBAAwB,CAAC;AAChC,YAAY,EACV,gBAAgB,EAChB,qBAAqB,EACrB,oBAAoB,EACpB,cAAc,EACd,eAAe,GAChB,MAAM,wBAAwB,CAAC;AAGhC,OAAO,EACL,KAAK,EACL,OAAO,EACP,SAAS,EACT,OAAO,EACP,OAAO,EACP,MAAM,EACN,MAAM,EACN,QAAQ,GACT,MAAM,oBAAoB,CAAC;AAG5B,OAAO,EACL,cAAc,EACd,WAAW,EACX,gBAAgB,EAChB,aAAa,EACb,kBAAkB,EAClB,IAAI,EACJ,SAAS,EACT,QAAQ,EACR,aAAa,EACb,OAAO,EACP,YAAY,EACZ,SAAS,EACT,SAAS,EACT,cAAc,EACd,UAAU,EACV,UAAU,EACV,eAAe,EACf,WAAW,EACX,WAAW,EACX,gBAAgB,GACjB,MAAM,6BAA6B,CAAC;AAGrC,OAAO,EACL,KAAK,EACL,UAAU,EACV,aAAa,EACb,QAAQ,EACR,WAAW,EACX,mBAAmB,EACnB,wBAAwB,EACxB,QAAQ,GACT,MAAM,uBAAuB,CAAC;AAC/B,YAAY,EACV,WAAW,EACX,gBAAgB,EAChB,aAAa,EACb,cAAc,EACd,iBAAiB,EACjB,eAAe,GAChB,MAAM,uBAAuB,CAAC;AAI/B,OAAO,EACL,IAAI,EACJ,IAAI,EACJ,MAAM,EACN,QAAQ,EACR,SAAS,EACT,SAAS,EACT,MAAM,EACN,eAAe,EACf,cAAc,EACd,EAAE,GACH,MAAM,mBAAmB,CAAC;AAG3B,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,qBAAqB,CAAC;AAKxD,OAAO,EAAE,GAAG,EAAE,MAAM,8BAA8B,CAAC;AACnD,YAAY,EAAE,SAAS,EAAE,UAAU,EAAE,MAAM,8BAA8B,CAAC;AAC1E,OAAO,EAAE,IAAI,EAAE,MAAM,uBAAuB,CAAC;AAC7C,YAAY,EAAE,WAAW,EAAE,MAAM,uBAAuB,CAAC;AAIzD,OAAO,EAAE,eAAe,EAAE,MAAM,+BAA+B,CAAC;AAChE,YAAY,EAAE,WAAW,EAAE,sBAAsB,EAAE,MAAM,+BAA+B,CAAC;AAIzF,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,gCAAgC,CAAC;AACxE,YAAY,EAAE,aAAa,EAAE,aAAa,EAAE,aAAa,EAAE,MAAM,gCAAgC,CAAC;AAIlG,OAAO,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,qBAAqB,CAAC;AACnD,YAAY,EAAE,aAAa,EAAE,MAAM,qBAAqB,CAAC;AAMzD,OAAO,EAAE,EAAE,EAAE,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,YAAY,EAAE,kBAAkB,EAAE,MAAM,qBAAqB,CAAC;AACpG,YAAY,EACV,mBAAmB,EACnB,cAAc,EACd,aAAa,EACb,YAAY,EACZ,YAAY,GACb,MAAM,qBAAqB,CAAC;AAM7B,OAAO,EAAE,KAAK,EAAE,MAAM,oBAAoB,CAAC;AAC3C,YAAY,EAAE,kBAAkB,EAAE,YAAY,EAAE,WAAW,EAAE,MAAM,oBAAoB,CAAC;AAKxF,OAAO,EAAE,IAAI,EAAE,MAAM,mBAAmB,CAAC;AACzC,YAAY,EAAE,WAAW,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AAOjE,OAAO,EAAE,WAAW,EAAE,WAAW,EAAE,aAAa,EAAE,GAAG,EAAE,MAAM,iCAAiC,CAAC;AAC/F,YAAY,EAAE,SAAS,EAAE,MAAM,iCAAiC,CAAC;AAQjE,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,UAAU,EAAE,UAAU,EAAE,MAAM,+BAA+B,CAAC;AACzF,YAAY,EACV,YAAY,EACZ,aAAa,EACb,qBAAqB,GACtB,MAAM,+BAA+B,CAAC;AAQvC,OAAO,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,gCAAgC,CAAC;AAKnE,OAAO,EAAE,cAAc,EAAE,MAAM,8BAA8B,CAAC;AAC9D,YAAY,EAAE,qBAAqB,EAAE,oBAAoB,EAAE,MAAM,8BAA8B,CAAC;AAKhG,OAAO,EAAE,IAAI,EAAE,MAAM,4BAA4B,CAAC;AAClD,YAAY,EAAE,WAAW,EAAE,UAAU,EAAE,MAAM,4BAA4B,CAAC;AAM1E,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,YAAY,EAAE,YAAY,EAAE,MAAM,sBAAsB,CAAC;AAOzD,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,MAAM,sBAAsB,CAAC;AAC/D,YAAY,EAAE,iBAAiB,EAAE,YAAY,EAAE,MAAM,sBAAsB,CAAC;AAM5E,OAAO,EAAE,mBAAmB,EAAE,MAAM,0BAA0B,CAAC;AAC/D,YAAY,EAAE,KAAK,EAAE,iBAAiB,EAAE,gBAAgB,EAAE,MAAM,0BAA0B,CAAC;AAG3F,OAAO,EAAE,aAAa,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,MAAM,uBAAuB,CAAC;AAI1E,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,YAAY,EAAE,QAAQ,EAAE,MAAM,uBAAuB,CAAC;AAC/E,YAAY,EAAE,cAAc,EAAE,cAAc,EAAE,MAAM,uBAAuB,CAAC;AAG5E,OAAO,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AACzD,YAAY,EAAE,gBAAgB,EAAE,MAAM,uBAAuB,CAAC;AAG9D,OAAO,EAAE,GAAG,EAAE,MAAM,aAAa,CAAC;AAClC,YAAY,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAGzD,OAAO,EAAE,KAAK,EAAE,KAAK,EAAE,UAAU,EAAE,MAAM,gCAAgC,CAAC;AAC1E,YAAY,EACV,YAAY,EACZ,wBAAwB,EACxB,2BAA2B,GAC5B,MAAM,gCAAgC,CAAC;AAGxC,OAAO,EAAE,kBAAkB,EAAE,MAAM,6BAA6B,CAAC;AACjE,YAAY,EACV,yBAAyB,EACzB,wBAAwB,GACzB,MAAM,6BAA6B,CAAC;AAGrC,OAAO,EAAE,UAAU,EAAE,eAAe,EAAE,kBAAkB,EAAE,MAAM,yBAAyB,CAAC;AAC1F,YAAY,EAAE,cAAc,EAAE,QAAQ,EAAE,MAAM,yBAAyB,CAAC;AAGxE,OAAO,EAAE,MAAM,EAAE,kBAAkB,EAAE,MAAM,uBAAuB,CAAC;AAInE,OAAO,EAAE,MAAM,EAAE,WAAW,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAIrE,OAAO,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAC1C,YAAY,EAAE,kBAAkB,EAAE,iBAAiB,EAAE,MAAM,aAAa,CAAC;AAIzE,OAAO,EAAE,eAAe,EAAE,YAAY,EAAE,MAAM,4BAA4B,CAAC;AAC3E,YAAY,EACV,sBAAsB,EACtB,qBAAqB,EACrB,kBAAkB,GACnB,MAAM,4BAA4B,CAAC;AAIpC,OAAO,EAAE,eAAe,EAAE,MAAM,6BAA6B,CAAC;AAC9D,YAAY,EAAE,gBAAgB,EAAE,qBAAqB,EAAE,MAAM,6BAA6B,CAAC;AAM3F,OAAO,EAAE,IAAI,EAAE,MAAM,mBAAmB,CAAC;AACzC,YAAY,EAAE,WAAW,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AACjE,YAAY,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AAK1D,OAAO,EAAE,IAAI,EAAE,UAAU,EAAE,MAAM,mBAAmB,CAAC;AACrD,YAAY,EACV,WAAW,EACX,UAAU,EACV,iBAAiB,EACjB,gBAAgB,GACjB,MAAM,mBAAmB,CAAC;AAG3B,OAAO,EAAE,gBAAgB,EAAE,MAAM,sBAAsB,CAAC;AAGxD,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAGxF,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,QAAQ,EAAE,SAAS,EAAE,OAAO,EAAE,QAAQ,EAAE,IAAI,EAAE,MAAM,yBAAyB,CAAC;AAIpG,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,MAAM,wBAAwB,CAAC;AAIpG,OAAO,EACL,cAAc,EACd,KAAK,EACL,KAAK,EACL,MAAM,EACN,MAAM,EACN,UAAU,GACX,MAAM,2BAA2B,CAAC;AACnC,YAAY,EAAE,gBAAgB,EAAE,MAAM,2BAA2B,CAAC;AAClE,YAAY,EAAE,SAAS,EAAE,MAAM,4BAA4B,CAAC;AAI5D,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,OAAO,EAAE,GAAG,EAAE,MAAM,sBAAsB,CAAC;AAKnE,OAAO,EACL,SAAS,EACT,UAAU,EACV,GAAG,EACH,SAAS,EACT,IAAI,EACJ,KAAK,EACL,QAAQ,GACT,MAAM,4BAA4B,CAAC;AACpC,YAAY,EACV,gBAAgB,EAChB,UAAU,EACV,WAAW,EACX,UAAU,GACX,MAAM,4BAA4B,CAAC;AAGpC,OAAO,EACL,SAAS,EACT,eAAe,EACf,KAAK,EACL,kBAAkB,EAClB,mBAAmB,EACnB,mBAAmB,EACnB,uBAAuB,EACvB,gBAAgB,EAChB,0BAA0B,GAC3B,MAAM,qBAAqB,CAAC;AAI7B,OAAO,EACL,iBAAiB,EACjB,eAAe,EACf,iBAAiB,EACjB,aAAa,EACb,aAAa,EACb,aAAa,EACb,WAAW,EACX,UAAU,EACV,SAAS,EACT,YAAY,EACZ,aAAa,EACb,QAAQ,EACR,WAAW,GACZ,MAAM,8BAA8B,CAAC;AACtC,YAAY,EAAE,gBAAgB,EAAE,MAAM,8BAA8B,CAAC;AAGrE,OAAO,EACL,oBAAoB,EACpB,iBAAiB,EACjB,qBAAqB,GACtB,MAAM,2BAA2B,CAAC;AACnC,YAAY,EAAE,MAAM,EAAE,oBAAoB,EAAE,MAAM,2BAA2B,CAAC;AAU9E,OAAO,EACL,iBAAiB,EACjB,cAAc,EACd,cAAc,EACd,QAAQ,EACR,OAAO,EACP,OAAO,EACP,WAAW,EACX,OAAO,EACP,QAAQ,GACT,MAAM,6BAA6B,CAAC;AACrC,YAAY,EACV,eAAe,EACf,cAAc,EACd,aAAa,EACb,cAAc,GACf,MAAM,6BAA6B,CAAC;AAIrC,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,6BAA6B,CAAC;AAK1E,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,OAAO,EAAE,WAAW,EAAE,MAAM,kCAAkC,CAAC;AAM7F,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,8BAA8B,CAAC;AAC5E,OAAO,EAAE,aAAa,EAAE,MAAM,0BAA0B,CAAC;AACzD,YAAY,EAAE,mBAAmB,EAAE,MAAM,0BAA0B,CAAC;AAIpE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,OAAO,EAAE,SAAS,EAAE,SAAS,EAAE,OAAO,EAAE,MAAM,oBAAoB,CAAC;AAK9F,OAAO,EACL,OAAO,EACP,cAAc,EACd,QAAQ,EACR,kBAAkB,EAClB,QAAQ,GACT,MAAM,6BAA6B,CAAC;AAMrC,OAAO,EACL,GAAG,EACH,GAAG,EACH,aAAa,EACb,WAAW,EACX,mBAAmB,EACnB,kBAAkB,GACnB,MAAM,iCAAiC,CAAC;AACzC,YAAY,EAAE,iBAAiB,EAAE,MAAM,iCAAiC,CAAC;AAKzE,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAC5E,YAAY,EACV,aAAa,EACb,YAAY,EACZ,WAAW,EACX,eAAe,GAChB,MAAM,yBAAyB,CAAC;AAMjC,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,uBAAuB,CAAC;AAO3D,OAAO,EACL,aAAa,EACb,SAAS,EACT,kBAAkB,EAClB,cAAc,EACd,YAAY,GACb,MAAM,+BAA+B,CAAC;AAKvC,OAAO,EAAE,GAAG,EAAE,MAAM,aAAa,CAAC;AAClC,YAAY,EAAE,SAAS,EAAE,OAAO,EAAE,UAAU,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAK7E,OAAO,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AACnD,YAAY,EAAE,SAAS,EAAE,MAAM,EAAE,gBAAgB,EAAE,MAAM,gBAAgB,CAAC;AAI1E,OAAO,EAAE,UAAU,EAAE,MAAM,2BAA2B,CAAC;AACvD,YAAY,EAAE,qBAAqB,EAAE,iBAAiB,EAAE,MAAM,2BAA2B,CAAC;AAI1F,OAAO,EAAE,yBAAyB,EAAE,YAAY,EAAE,MAAM,6BAA6B,CAAC;AACtF,YAAY,EAAE,QAAQ,EAAE,SAAS,EAAE,YAAY,EAAE,OAAO,EAAE,MAAM,6BAA6B,CAAC;AAI9F,OAAO,EAAE,eAAe,EAAE,YAAY,EAAE,aAAa,EAAE,MAAM,kCAAkC,CAAC;AAChG,YAAY,EAAE,iBAAiB,EAAE,MAAM,kCAAkC,CAAC"}
|
package/dist/index.js
CHANGED
|
@@ -31747,9 +31747,9 @@ var dependencies187 = ["typed", "factorial"];
|
|
|
31747
31747
|
var createPermutations = /* @__PURE__ */ factory(
|
|
31748
31748
|
name187,
|
|
31749
31749
|
dependencies187,
|
|
31750
|
-
({ typed: typed3, factorial:
|
|
31750
|
+
({ typed: typed3, factorial: factorial6 }) => {
|
|
31751
31751
|
return typed3(name187, {
|
|
31752
|
-
"number | BigNumber":
|
|
31752
|
+
"number | BigNumber": factorial6,
|
|
31753
31753
|
"number, number": function(n, k) {
|
|
31754
31754
|
if (!isInteger(n) || n < 0) {
|
|
31755
31755
|
throw new TypeError("Positive integer value expected in function permutations");
|
|
@@ -35973,7 +35973,7 @@ var dependencies224 = ["typed", "add", "divide", "multiply", "factorial", "isInt
|
|
|
35973
35973
|
var createMultinomial = /* @__PURE__ */ factory(
|
|
35974
35974
|
name224,
|
|
35975
35975
|
dependencies224,
|
|
35976
|
-
({ typed: typed3, add: add4, divide: divide2, multiply: multiply2, factorial:
|
|
35976
|
+
({ typed: typed3, add: add4, divide: divide2, multiply: multiply2, factorial: factorial6, isInteger: isInteger3, isPositive: isPositive2 }) => {
|
|
35977
35977
|
return typed3(name224, {
|
|
35978
35978
|
"Array | Matrix": function(a) {
|
|
35979
35979
|
let sum3 = 0;
|
|
@@ -35983,9 +35983,9 @@ var createMultinomial = /* @__PURE__ */ factory(
|
|
|
35983
35983
|
throw new TypeError("Positive integer value expected in function multinomial");
|
|
35984
35984
|
}
|
|
35985
35985
|
sum3 = add4(sum3, ai);
|
|
35986
|
-
denom = multiply2(denom,
|
|
35986
|
+
denom = multiply2(denom, factorial6(ai));
|
|
35987
35987
|
});
|
|
35988
|
-
return divide2(
|
|
35988
|
+
return divide2(factorial6(sum3), denom);
|
|
35989
35989
|
}
|
|
35990
35990
|
});
|
|
35991
35991
|
}
|
|
@@ -37042,7 +37042,7 @@ var createZeta = /* @__PURE__ */ factory(
|
|
|
37042
37042
|
multiply: multiply2,
|
|
37043
37043
|
pow: pow2,
|
|
37044
37044
|
divide: divide2,
|
|
37045
|
-
factorial:
|
|
37045
|
+
factorial: factorial6,
|
|
37046
37046
|
equal: equal2,
|
|
37047
37047
|
smallerEq: smallerEq2,
|
|
37048
37048
|
isBounded: isBounded2,
|
|
@@ -37120,8 +37120,8 @@ var createZeta = /* @__PURE__ */ factory(
|
|
|
37120
37120
|
let S = k;
|
|
37121
37121
|
for (let j = k; smallerEq2(j, n); j = add4(j, 1)) {
|
|
37122
37122
|
const factor2 = divide2(
|
|
37123
|
-
multiply2(
|
|
37124
|
-
multiply2(
|
|
37123
|
+
multiply2(factorial6(add4(n, subtract3(j, 1))), pow2(4, j)),
|
|
37124
|
+
multiply2(factorial6(subtract3(n, j)), factorial6(multiply2(2, j)))
|
|
37125
37125
|
);
|
|
37126
37126
|
S = add4(S, factor2);
|
|
37127
37127
|
}
|
|
@@ -51137,6 +51137,161 @@ function barnesG(z) {
|
|
|
51137
51137
|
return Math.exp(lnG);
|
|
51138
51138
|
}
|
|
51139
51139
|
|
|
51140
|
+
// src/special/wave-functions.ts
|
|
51141
|
+
var RELTOL3 = 1e-16;
|
|
51142
|
+
function cMul3(a, b) {
|
|
51143
|
+
return { re: a.re * b.re - a.im * b.im, im: a.re * b.im + a.im * b.re };
|
|
51144
|
+
}
|
|
51145
|
+
function cInv(a) {
|
|
51146
|
+
const d = a.re * a.re + a.im * a.im;
|
|
51147
|
+
return { re: a.re / d, im: -a.im / d };
|
|
51148
|
+
}
|
|
51149
|
+
function cLog(a) {
|
|
51150
|
+
return { re: 0.5 * Math.log(a.re * a.re + a.im * a.im), im: Math.atan2(a.im, a.re) };
|
|
51151
|
+
}
|
|
51152
|
+
function cSin(a) {
|
|
51153
|
+
return { re: Math.sin(a.re) * Math.cosh(a.im), im: Math.cos(a.re) * Math.sinh(a.im) };
|
|
51154
|
+
}
|
|
51155
|
+
var LANCZOS_G = 7;
|
|
51156
|
+
var LANCZOS_C = [
|
|
51157
|
+
0.9999999999998099,
|
|
51158
|
+
676.5203681218851,
|
|
51159
|
+
-1259.1392167224028,
|
|
51160
|
+
771.3234287776531,
|
|
51161
|
+
-176.6150291621406,
|
|
51162
|
+
12.507343278686905,
|
|
51163
|
+
-0.13857109526572012,
|
|
51164
|
+
9984369578019572e-21,
|
|
51165
|
+
15056327351493116e-23
|
|
51166
|
+
];
|
|
51167
|
+
function lnGammaComplex(z) {
|
|
51168
|
+
if (z.re < 0.5) {
|
|
51169
|
+
const sinPiZ = cSin({ re: Math.PI * z.re, im: Math.PI * z.im });
|
|
51170
|
+
const lnSin = cLog(sinPiZ);
|
|
51171
|
+
const rest = lnGammaComplex({ re: 1 - z.re, im: -z.im });
|
|
51172
|
+
return { re: Math.log(Math.PI) - lnSin.re - rest.re, im: -lnSin.im - rest.im };
|
|
51173
|
+
}
|
|
51174
|
+
const x = { re: z.re - 1, im: z.im };
|
|
51175
|
+
let sum3 = { re: LANCZOS_C[0], im: 0 };
|
|
51176
|
+
for (let i = 1; i < LANCZOS_G + 2; i++) {
|
|
51177
|
+
const dr = x.re + i;
|
|
51178
|
+
const di = x.im;
|
|
51179
|
+
const denom = dr * dr + di * di;
|
|
51180
|
+
sum3 = { re: sum3.re + LANCZOS_C[i] * dr / denom, im: sum3.im - LANCZOS_C[i] * di / denom };
|
|
51181
|
+
}
|
|
51182
|
+
const t = { re: x.re + LANCZOS_G + 0.5, im: x.im };
|
|
51183
|
+
const lnT = cLog(t);
|
|
51184
|
+
const lnSum = cLog(sum3);
|
|
51185
|
+
const half = { re: x.re + 0.5, im: x.im };
|
|
51186
|
+
const p = cMul3(half, lnT);
|
|
51187
|
+
return {
|
|
51188
|
+
re: 0.5 * Math.log(2 * Math.PI) + p.re - t.re + lnSum.re,
|
|
51189
|
+
im: p.im - t.im + lnSum.im
|
|
51190
|
+
};
|
|
51191
|
+
}
|
|
51192
|
+
function factorial5(n) {
|
|
51193
|
+
let r = 1;
|
|
51194
|
+
for (let i = 2; i <= n; i++) r *= i;
|
|
51195
|
+
return r;
|
|
51196
|
+
}
|
|
51197
|
+
function riemannSiegelTheta(t) {
|
|
51198
|
+
const lg = lnGammaComplex({ re: 0.25, im: t / 2 });
|
|
51199
|
+
return lg.im - t / 2 * Math.log(Math.PI);
|
|
51200
|
+
}
|
|
51201
|
+
function borweinD(k, n) {
|
|
51202
|
+
let S = k;
|
|
51203
|
+
for (let j = k; j <= n; j++) {
|
|
51204
|
+
S += factorial5(n + j - 1) * Math.pow(4, j) / (factorial5(n - j) * factorial5(2 * j));
|
|
51205
|
+
}
|
|
51206
|
+
return n * S;
|
|
51207
|
+
}
|
|
51208
|
+
function zetaCriticalLine(t) {
|
|
51209
|
+
const n = Math.round(1.3 * 20 + 0.9 * Math.abs(t));
|
|
51210
|
+
const ln2 = Math.LN2;
|
|
51211
|
+
const mag = Math.SQRT2;
|
|
51212
|
+
const den = { re: 1 - mag * Math.cos(t * ln2), im: mag * Math.sin(t * ln2) };
|
|
51213
|
+
const d0 = borweinD(0, n);
|
|
51214
|
+
const c = cInv(cMul3({ re: d0, im: 0 }, den));
|
|
51215
|
+
let Sr = 0;
|
|
51216
|
+
let Si = 0;
|
|
51217
|
+
for (let k = 1; k <= n; k++) {
|
|
51218
|
+
const sign3 = (k - 1) % 2 === 0 ? 1 : -1;
|
|
51219
|
+
const dk = borweinD(k, n);
|
|
51220
|
+
const lnk = Math.log(k);
|
|
51221
|
+
const kpow = 1 / Math.sqrt(k);
|
|
51222
|
+
Sr += sign3 * dk * kpow * Math.cos(t * lnk);
|
|
51223
|
+
Si += sign3 * dk * kpow * -Math.sin(t * lnk);
|
|
51224
|
+
}
|
|
51225
|
+
return cMul3(c, { re: Sr, im: Si });
|
|
51226
|
+
}
|
|
51227
|
+
function siegelZ(t) {
|
|
51228
|
+
const at = Math.abs(t);
|
|
51229
|
+
const z = zetaCriticalLine(at);
|
|
51230
|
+
const th = riemannSiegelTheta(at);
|
|
51231
|
+
return Math.cos(th) * z.re - Math.sin(th) * z.im;
|
|
51232
|
+
}
|
|
51233
|
+
function riemannSiegelZ(t) {
|
|
51234
|
+
return siegelZ(t);
|
|
51235
|
+
}
|
|
51236
|
+
var LERCH_MAX_TERMS = 2e5;
|
|
51237
|
+
function lerchPhi(z, s, a) {
|
|
51238
|
+
if (Math.abs(z) >= 1) {
|
|
51239
|
+
throw new Error("lerchPhi: only |z| < 1 is supported (analytic continuation not implemented)");
|
|
51240
|
+
}
|
|
51241
|
+
if (a <= 0) {
|
|
51242
|
+
throw new Error("lerchPhi: a must be > 0");
|
|
51243
|
+
}
|
|
51244
|
+
let sum3 = 0;
|
|
51245
|
+
let zk = 1;
|
|
51246
|
+
for (let k = 0; k < LERCH_MAX_TERMS; k++) {
|
|
51247
|
+
const term = zk / Math.pow(a + k, s);
|
|
51248
|
+
sum3 += term;
|
|
51249
|
+
if (k > 1 && Math.abs(term) < Math.abs(sum3) * RELTOL3) break;
|
|
51250
|
+
zk *= z;
|
|
51251
|
+
}
|
|
51252
|
+
return sum3;
|
|
51253
|
+
}
|
|
51254
|
+
function recipGamma(y) {
|
|
51255
|
+
if (y > 0) return Math.exp(-_lgamma(y));
|
|
51256
|
+
if (y === Math.floor(y)) return 0;
|
|
51257
|
+
return Math.sin(Math.PI * y) * Math.exp(_lgamma(1 - y)) / Math.PI;
|
|
51258
|
+
}
|
|
51259
|
+
function parabolicCylinderD(nu, x) {
|
|
51260
|
+
const z = x * x / 2;
|
|
51261
|
+
const even = recipGamma((1 - nu) / 2) * hyp1f1(-nu / 2, 0.5, z);
|
|
51262
|
+
const odd = recipGamma(-nu / 2) * hyp1f1((1 - nu) / 2, 1.5, z);
|
|
51263
|
+
return Math.pow(2, nu / 2) * Math.exp(-(x * x) / 4) * (Math.sqrt(Math.PI) * even - Math.sqrt(2 * Math.PI) * x * odd);
|
|
51264
|
+
}
|
|
51265
|
+
var COULOMB_MAX_TERMS = 5e3;
|
|
51266
|
+
function coulombC(L, eta) {
|
|
51267
|
+
const lg = lnGammaComplex({ re: L + 1, im: eta });
|
|
51268
|
+
const absGamma = Math.exp(lg.re);
|
|
51269
|
+
return Math.pow(2, L) * Math.exp(-Math.PI * eta / 2) * absGamma / Math.exp(_lgamma(2 * L + 2));
|
|
51270
|
+
}
|
|
51271
|
+
function coulombF(L, eta, rho) {
|
|
51272
|
+
if (rho < 0) {
|
|
51273
|
+
throw new Error("coulombF: rho must be >= 0");
|
|
51274
|
+
}
|
|
51275
|
+
if (L < 0) {
|
|
51276
|
+
throw new Error("coulombF: L must be >= 0");
|
|
51277
|
+
}
|
|
51278
|
+
if (rho === 0) return 0;
|
|
51279
|
+
let aPrev2 = 0;
|
|
51280
|
+
let aPrev1 = 1;
|
|
51281
|
+
let sum3 = 1;
|
|
51282
|
+
let rhoK = 1;
|
|
51283
|
+
for (let k = 1; k < COULOMB_MAX_TERMS; k++) {
|
|
51284
|
+
const ak = (2 * eta * aPrev1 - aPrev2) / (k * (k + 2 * L + 1));
|
|
51285
|
+
rhoK *= rho;
|
|
51286
|
+
const term = ak * rhoK;
|
|
51287
|
+
sum3 += term;
|
|
51288
|
+
aPrev2 = aPrev1;
|
|
51289
|
+
aPrev1 = ak;
|
|
51290
|
+
if (k > 2 && term !== 0 && Math.abs(term) < Math.abs(sum3) * RELTOL3) break;
|
|
51291
|
+
}
|
|
51292
|
+
return coulombC(L, eta) * Math.pow(rho, L + 1) * sum3;
|
|
51293
|
+
}
|
|
51294
|
+
|
|
51140
51295
|
// src/graph/traversal-centrality.ts
|
|
51141
51296
|
function bfs(adj, start) {
|
|
51142
51297
|
const n = adj.length;
|
|
@@ -52525,6 +52680,7 @@ export {
|
|
|
52525
52680
|
coth,
|
|
52526
52681
|
coulomb,
|
|
52527
52682
|
coulombConstant,
|
|
52683
|
+
coulombF,
|
|
52528
52684
|
count,
|
|
52529
52685
|
cov,
|
|
52530
52686
|
createUnit,
|
|
@@ -52914,6 +53070,7 @@ export {
|
|
|
52914
53070
|
leftShift,
|
|
52915
53071
|
legendreFit,
|
|
52916
53072
|
legendreP,
|
|
53073
|
+
lerchPhi,
|
|
52917
53074
|
levenbergMarquardt,
|
|
52918
53075
|
leveneTest,
|
|
52919
53076
|
lfilter,
|
|
@@ -53056,6 +53213,7 @@ export {
|
|
|
53056
53213
|
pacf,
|
|
53057
53214
|
padeApproximant,
|
|
53058
53215
|
pageRank,
|
|
53216
|
+
parabolicCylinderD,
|
|
53059
53217
|
parallelAutoCorr,
|
|
53060
53218
|
parallelConv,
|
|
53061
53219
|
parallelFFT,
|
|
@@ -53186,6 +53344,7 @@ export {
|
|
|
53186
53344
|
rfft,
|
|
53187
53345
|
rfftfreq,
|
|
53188
53346
|
ridge,
|
|
53347
|
+
riemannSiegelZ,
|
|
53189
53348
|
rightArithShift,
|
|
53190
53349
|
rightLogShift,
|
|
53191
53350
|
risingFactorial,
|
|
@@ -53228,6 +53387,7 @@ export {
|
|
|
53228
53387
|
shapiroWilkTest,
|
|
53229
53388
|
shortestPath,
|
|
53230
53389
|
shouldParallelize,
|
|
53390
|
+
siegelZ,
|
|
53231
53391
|
sigmoid,
|
|
53232
53392
|
sign,
|
|
53233
53393
|
simplify,
|
|
@@ -0,0 +1,123 @@
|
|
|
1
|
+
/**
|
|
2
|
+
* Advanced niche special functions: the Riemann–Siegel Z-function, the Lerch
|
|
3
|
+
* transcendent, the parabolic-cylinder function D_ν, and the regular Coulomb
|
|
4
|
+
* wave function F_L.
|
|
5
|
+
*
|
|
6
|
+
* These follow the same plain-exported-function pattern as `niche.ts` and
|
|
7
|
+
* `hypergeometric.ts` — pure real-valued `number -> number` math with no
|
|
8
|
+
* typed-function array/WASM dispatch overloads. Each is oracle-pinned against
|
|
9
|
+
* `mpmath` (dps=30) in `functions/tests/gap-special-wave-oracle.test.ts`.
|
|
10
|
+
*
|
|
11
|
+
* The self-contained complex-arithmetic helpers here (complex log-gamma via
|
|
12
|
+
* Lanczos, and a critical-line Borwein ζ evaluator) exist so these functions
|
|
13
|
+
* stay pure `number`-in/`number`-out and do not need the factory `zeta`
|
|
14
|
+
* (which requires a full math instance with a Complex constructor).
|
|
15
|
+
*
|
|
16
|
+
* @packageDocumentation
|
|
17
|
+
*/
|
|
18
|
+
/** 64-bit float (default for decimals) */
|
|
19
|
+
type f64 = number;
|
|
20
|
+
/**
|
|
21
|
+
* Riemann–Siegel Z-function Z(t) = e^{iθ(t)} ζ(1/2 + it), the real-valued
|
|
22
|
+
* function on the critical line whose real zeros coincide with the imaginary
|
|
23
|
+
* parts of the non-trivial zeros of ζ. Z is even: Z(−t) = Z(t), and
|
|
24
|
+
* Z(0) = ζ(1/2) ≈ −1.4603545088.
|
|
25
|
+
*
|
|
26
|
+
* Computed by evaluating ζ(1/2 + it) via the Borwein-accelerated series and
|
|
27
|
+
* the theta function via the complex log-gamma, then taking the (real) product
|
|
28
|
+
* — i.e. Z(t) = cos θ · Re ζ − sin θ · Im ζ. This exact-theta approach is
|
|
29
|
+
* accurate for all t (the classical Riemann–Siegel asymptotic expansion, by
|
|
30
|
+
* contrast, is inaccurate for the small t of the first few zeros).
|
|
31
|
+
*
|
|
32
|
+
* Verified against `mpmath.siegelz` (dps=30): relative error < 1e-11 for
|
|
33
|
+
* O(1)-magnitude values and absolute error < 1e-9 near the first zeros
|
|
34
|
+
* (t ≈ 14.13, 21.02, 25.01), for |t| up to ~40.
|
|
35
|
+
*
|
|
36
|
+
* @param t - Height on the critical line (real; any sign)
|
|
37
|
+
* @returns Z(t)
|
|
38
|
+
*
|
|
39
|
+
* @example
|
|
40
|
+
* siegelZ(15) // ~0.7199423913
|
|
41
|
+
* siegelZ(14.134725) // ~-1.12e-7 (near the first zeta zero)
|
|
42
|
+
*/
|
|
43
|
+
export declare function siegelZ(t: f64): f64;
|
|
44
|
+
/**
|
|
45
|
+
* Alias for {@link siegelZ} — the Riemann–Siegel Z-function Z(t).
|
|
46
|
+
*/
|
|
47
|
+
export declare function riemannSiegelZ(t: f64): f64;
|
|
48
|
+
/**
|
|
49
|
+
* Lerch transcendent Φ(z, s, a) = Σ_{k≥0} z^k / (a + k)^s, evaluated by its
|
|
50
|
+
* defining power series (DLMF 25.14.1), which converges for `|z| < 1` and any
|
|
51
|
+
* real `s`, provided `a > 0` (so no denominator vanishes).
|
|
52
|
+
*
|
|
53
|
+
* Analytic continuation to `|z| ≥ 1` is out of scope and throws. The boundary
|
|
54
|
+
* case Φ(1, s, a) = ζ(s, a) is the Hurwitz zeta, not implemented here.
|
|
55
|
+
*
|
|
56
|
+
* Cross-check: the polylogarithm satisfies `Li_s(z) = z · Φ(z, s, 1)`.
|
|
57
|
+
*
|
|
58
|
+
* Verified against `mpmath.lerchphi` (dps=30) to relative error < 1e-12 across
|
|
59
|
+
* real `z ∈ (−1, 1)` and assorted `s`, `a`.
|
|
60
|
+
*
|
|
61
|
+
* @param z - Argument, must satisfy `|z| < 1`
|
|
62
|
+
* @param s - Order (any real number)
|
|
63
|
+
* @param a - Offset (real, `a > 0`)
|
|
64
|
+
* @returns Φ(z, s, a)
|
|
65
|
+
* @throws {Error} If `|z| >= 1`, or if `a + k <= 0` for some term (a <= 0)
|
|
66
|
+
*
|
|
67
|
+
* @example
|
|
68
|
+
* lerchPhi(0.5, 2, 1) // ~1.1644810529 (= 2·Li_2(0.5))
|
|
69
|
+
*/
|
|
70
|
+
export declare function lerchPhi(z: f64, s: f64, a: f64): f64;
|
|
71
|
+
/**
|
|
72
|
+
* Parabolic cylinder function D_ν(x) (Whittaker's form), via the confluent
|
|
73
|
+
* hypergeometric representation (Abramowitz & Stegun 19.3.1 / DLMF 12.4.1):
|
|
74
|
+
*
|
|
75
|
+
* D_ν(x) = 2^{ν/2} e^{−x²/4} [ √π / Γ((1−ν)/2) · M(−ν/2; 1/2; x²/2)
|
|
76
|
+
* − √(2π) x / Γ(−ν/2) · M((1−ν)/2; 3/2; x²/2) ]
|
|
77
|
+
*
|
|
78
|
+
* where M = ₁F₁ is Kummer's function (`hyp1f1`). The even/odd split makes this
|
|
79
|
+
* valid for any real order ν and any real x (D_ν is entire in x); accuracy is
|
|
80
|
+
* best for moderate |x| where the ₁F₁ series converges without heavy
|
|
81
|
+
* cancellation.
|
|
82
|
+
*
|
|
83
|
+
* Verified against `mpmath.pcfd` (dps=30) to relative error < 1e-10 on a corpus
|
|
84
|
+
* of real (ν, x) with |x| ≤ ~2.
|
|
85
|
+
*
|
|
86
|
+
* @param nu - Order (real)
|
|
87
|
+
* @param x - Argument (real)
|
|
88
|
+
* @returns D_ν(x)
|
|
89
|
+
*
|
|
90
|
+
* @example
|
|
91
|
+
* parabolicCylinderD(0, 1) // ~0.7788007831 (= e^{-1/4})
|
|
92
|
+
*/
|
|
93
|
+
export declare function parabolicCylinderD(nu: f64, x: f64): f64;
|
|
94
|
+
/**
|
|
95
|
+
* Regular Coulomb wave function F_L(η, ρ), via its ascending power series
|
|
96
|
+
* (DLMF 33.6.1–33.6.2):
|
|
97
|
+
*
|
|
98
|
+
* F_L(η, ρ) = C_L(η) ρ^{L+1} Σ_{k≥0} A_k ρ^k,
|
|
99
|
+
* A_0 = 1, A_1 = η/(L+1), A_k = (2η A_{k−1} − A_{k−2}) / (k (k + 2L + 1)).
|
|
100
|
+
*
|
|
101
|
+
* The series is entire in ρ, so this converges for all ρ ≥ 0 (accuracy is best
|
|
102
|
+
* for small-to-moderate ρ before large-argument cancellation sets in). The
|
|
103
|
+
* *irregular* companion G_L is intentionally not implemented here — it requires
|
|
104
|
+
* a numerically-delicate asymptotic/continued-fraction treatment (see the
|
|
105
|
+
* module test and CHANGELOG for the surfaced design).
|
|
106
|
+
*
|
|
107
|
+
* Verified against `mpmath.coulombf` (dps=30) to relative error < 1e-6 (in
|
|
108
|
+
* practice < 1e-9 away from zeros) for L ∈ {0,1,2,3}, |η| ≤ 5, ρ ≤ 10.
|
|
109
|
+
*
|
|
110
|
+
* Special case: F_0(0, ρ) = sin ρ.
|
|
111
|
+
*
|
|
112
|
+
* @param L - Orbital angular momentum (real, `L >= 0`)
|
|
113
|
+
* @param eta - Sommerfeld parameter (real)
|
|
114
|
+
* @param rho - Radial variable (real, `rho >= 0`)
|
|
115
|
+
* @returns F_L(η, ρ)
|
|
116
|
+
* @throws {Error} If `rho < 0` or `L < 0`
|
|
117
|
+
*
|
|
118
|
+
* @example
|
|
119
|
+
* coulombF(0, 0, 1) // ~0.8414709848 (= sin 1)
|
|
120
|
+
*/
|
|
121
|
+
export declare function coulombF(L: f64, eta: f64, rho: f64): f64;
|
|
122
|
+
export {};
|
|
123
|
+
//# sourceMappingURL=wave-functions.d.ts.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"wave-functions.d.ts","sourceRoot":"","sources":["../../src/special/wave-functions.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAKH,0CAA0C;AAC1C,KAAK,GAAG,GAAG,MAAM,CAAC;AAyIlB;;;;;;;;;;;;;;;;;;;;;;GAsBG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAKnC;AAED;;GAEG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,CAE1C;AASD;;;;;;;;;;;;;;;;;;;;;GAqBG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CAiBpD;AAcD;;;;;;;;;;;;;;;;;;;;;GAqBG;AACH,wBAAgB,kBAAkB,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,GAAG,GAAG,CASvD;AAuBD;;;;;;;;;;;;;;;;;;;;;;;;;;GA0BG;AACH,wBAAgB,QAAQ,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,CAyBxD"}
|
package/package.json
CHANGED