dsp-collection 0.2.5 → 0.2.7

package/math/MutableComplex.d.ts

@@ -1,22 +1,24 @@
-import Complex from "./Complex.js";
-export default class MutableComplex extends Complex {
-    re: number;
-    im: number;
-    constructor(re?: number, im?: number);
-    static fromComplex(x: Complex): MutableComplex;
-    static expj(arg: number): MutableComplex;
-    static fromPolar(abs: number, arg: number): MutableComplex;
-    set(x: Complex): void;
-    setReIm(re: number, im?: number): void;
-    setExpj(arg: number): void;
-    addRealTo(x: number): void;
-    addTo(x: Complex): void;
-    subRealFrom(x: number): void;
-    subFrom(x: Complex): void;
-    mulByReal(x: number): void;
-    mulBy(x: Complex): void;
-    divByReal(x: number): void;
-    divBy(x: Complex): void;
-    setMul(re1: number, im1: number, re2: number, im2: number): void;
-    setDiv(re1: number, im1: number, re2: number, im2: number): void;
-}
+import Complex from "./Complex.ts";
+export default class MutableComplex extends Complex {
+    re: number;
+    im: number;
+    constructor(re?: number, im?: number);
+    static fromComplex(x: Complex): MutableComplex;
+    static expj(arg: number): MutableComplex;
+    static fromPolar(abs: number, arg: number): MutableComplex;
+    set(x: Complex): void;
+    setReIm(re: number, im?: number): void;
+    setExpj(arg: number): void;
+    setPolar(abs: number, arg: number): void;
+    addRealTo(x: number): void;
+    addTo(x: Complex): void;
+    subRealFrom(x: number): void;
+    subFrom(x: Complex): void;
+    mulByReal(x: number): void;
+    mulBy(x: Complex): void;
+    divByReal(x: number): void;
+    divBy(x: Complex): void;
+    setMul(re1: number, im1: number, re2: number, im2: number): void;
+    setDiv(re1: number, im1: number, re2: number, im2: number): void;
+}
+//# sourceMappingURL=MutableComplex.d.ts.map

package/math/MutableComplex.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"MutableComplex.d.ts","sourceRoot":"","sources":["../../src/math/MutableComplex.ts"],"names":[],"mappings":"AAAA,OAAO,OAAO,MAAM,cAAc,CAAC;AAUnC,MAAM,CAAC,OAAO,OAAO,cAAe,SAAQ,OAAO;IAKhC,EAAE,EAAE,MAAM,CAAC;IAGX,EAAE,EAAE,MAAM,CAAC;gBAYP,EAAE,GAAE,MAAU,EAAE,EAAE,GAAE,MAAU;WAQpC,WAAW,CAAE,CAAC,EAAE,OAAO,GAAI,cAAc;WAOhC,IAAI,CAAE,GAAG,EAAE,MAAM,GAAI,cAAc;WAOnC,SAAS,CAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAI,cAAc;IAQrE,GAAG,CAAE,CAAC,EAAE,OAAO;IAOf,OAAO,CAAE,EAAE,EAAE,MAAM,EAAE,EAAE,SAAI;IAO3B,OAAO,CAAE,GAAG,EAAE,MAAM;IAOpB,QAAQ,CAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM;IASlC,SAAS,CAAE,CAAC,EAAE,MAAM;IAMpB,KAAK,CAAE,CAAC,EAAE,OAAO;IAOjB,WAAW,CAAE,CAAC,EAAE,MAAM;IAMtB,OAAO,CAAE,CAAC,EAAE,OAAO;IAOnB,SAAS,CAAE,CAAC,EAAE,MAAM;IAOpB,KAAK,CAAE,CAAC,EAAE,OAAO;IAMjB,SAAS,CAAE,CAAC,EAAE,MAAM;IAOpB,KAAK,CAAE,CAAC,EAAE,OAAO;IAQjB,MAAM,CAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM;IAO1D,MAAM,CAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM;CAKhE"}

package/math/MutableComplex.js

@@ -1,65 +1,69 @@
-import Complex from "./Complex.js";
-export default class MutableComplex extends Complex {
-    constructor(re = 0, im = 0) {
-        super(re, im);
-    }
-    static fromComplex(x) {
-        return new MutableComplex(x.re, x.im);
-    }
-    static expj(arg) {
-        return new MutableComplex(Math.cos(arg), Math.sin(arg));
-    }
-    static fromPolar(abs, arg) {
-        return new MutableComplex(abs * Math.cos(arg), abs * Math.sin(arg));
-    }
-    set(x) {
-        this.re = x.re;
-        this.im = x.im;
-    }
-    setReIm(re, im = 0) {
-        this.re = re;
-        this.im = im;
-    }
-    setExpj(arg) {
-        this.re = Math.cos(arg);
-        this.im = Math.sin(arg);
-    }
-    addRealTo(x) {
-        this.re += x;
-    }
-    addTo(x) {
-        this.re += x.re;
-        this.im += x.im;
-    }
-    subRealFrom(x) {
-        this.re -= x;
-    }
-    subFrom(x) {
-        this.re -= x.re;
-        this.im -= x.im;
-    }
-    mulByReal(x) {
-        this.re *= x;
-        this.im *= x;
-    }
-    mulBy(x) {
-        this.setMul(this.re, this.im, x.re, x.im);
-    }
-    divByReal(x) {
-        this.re /= x;
-        this.im /= x;
-    }
-    divBy(x) {
-        this.setDiv(this.re, this.im, x.re, x.im);
-    }
-    setMul(re1, im1, re2, im2) {
-        this.re = re1 * re2 - im1 * im2;
-        this.im = re1 * im2 + im1 * re2;
-    }
-    setDiv(re1, im1, re2, im2) {
-        const m = re1 * re1 + im1 * im1;
-        this.re = (re1 * re2 + im1 * im2) / m;
-        this.im = (im1 * re2 - re1 * im2) / m;
-    }
-}
+import Complex from "./Complex.js";
+export default class MutableComplex extends Complex {
+    constructor(re = 0, im = 0) {
+        super(re, im);
+    }
+    static fromComplex(x) {
+        return new MutableComplex(x.re, x.im);
+    }
+    static expj(arg) {
+        return new MutableComplex(Math.cos(arg), Math.sin(arg));
+    }
+    static fromPolar(abs, arg) {
+        return new MutableComplex(abs * Math.cos(arg), abs * Math.sin(arg));
+    }
+    set(x) {
+        this.re = x.re;
+        this.im = x.im;
+    }
+    setReIm(re, im = 0) {
+        this.re = re;
+        this.im = im;
+    }
+    setExpj(arg) {
+        this.re = Math.cos(arg);
+        this.im = Math.sin(arg);
+    }
+    setPolar(abs, arg) {
+        this.re = abs * Math.cos(arg);
+        this.im = abs * Math.sin(arg);
+    }
+    addRealTo(x) {
+        this.re += x;
+    }
+    addTo(x) {
+        this.re += x.re;
+        this.im += x.im;
+    }
+    subRealFrom(x) {
+        this.re -= x;
+    }
+    subFrom(x) {
+        this.re -= x.re;
+        this.im -= x.im;
+    }
+    mulByReal(x) {
+        this.re *= x;
+        this.im *= x;
+    }
+    mulBy(x) {
+        this.setMul(this.re, this.im, x.re, x.im);
+    }
+    divByReal(x) {
+        this.re /= x;
+        this.im /= x;
+    }
+    divBy(x) {
+        this.setDiv(this.re, this.im, x.re, x.im);
+    }
+    setMul(re1, im1, re2, im2) {
+        this.re = re1 * re2 - im1 * im2;
+        this.im = re1 * im2 + im1 * re2;
+    }
+    setDiv(re1, im1, re2, im2) {
+        const m = re1 * re1 + im1 * im1;
+        this.re = (re1 * re2 + im1 * im2) / m;
+        this.im = (im1 * re2 - re1 * im2) / m;
+    }
+}
 //# sourceMappingURL=MutableComplex.js.map

package/math/MutableComplex.js.map

@@ -1 +1 @@
-{"version":3,"file":"MutableComplex.js","sourceRoot":"","sources":["../../src/math/MutableComplex.ts"],"names":[],"mappings":"AAAA,OAAO,OAAO,MAAM,cAAc,CAAC;AAUnC,MAAM,CAAC,OAAO,OAAO,cAAe,SAAQ,OAAO;IAoBhD,YAAoB,KAAa,CAAC,EAAE,KAAa,CAAC;QAC/C,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAAC,CAAC;IAOZ,MAAM,CAAC,WAAW,CAAE,CAAU;QAClC,OAAO,IAAI,cAAc,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAMpC,MAAM,CAAU,IAAI,CAAE,GAAW;QACrC,OAAO,IAAI,cAAc,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAAC,CAAC;IAMtD,MAAM,CAAU,SAAS,CAAE,GAAW,EAAE,GAAW;QACvD,OAAO,IAAI,cAAc,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAAC,CAAC;IAOlE,GAAG,CAAE,CAAU;QACnB,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;QACf,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAKb,OAAO,CAAE,EAAU,EAAE,EAAE,GAAG,CAAC;QAC/B,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IAAC,CAAC;IAKX,OAAO,CAAE,GAAW;QACxB,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACxB,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;IAOtB,SAAS,CAAE,CAAS;QACxB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;IAAC,CAAC;IAKX,KAAK,CAAE,CAAU;QACrB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC;QAChB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAKd,WAAW,CAAE,CAAS;QAC1B,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;IAAC,CAAC;IAKX,OAAO,CAAE,CAAU;QACvB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC;QAChB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAKd,SAAS,CAAE,CAAS;QACxB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;QACb,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;IAAC,CAAC;IAKX,KAAK,CAAE,CAAU;QACrB,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAKxC,SAAS,CAAE,CAAS;QACxB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;QACb,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;IAAC,CAAC;IAKX,KAAK,CAAE,CAAU;QACrB,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAOxC,MAAM,CAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW;QAC9D,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAChC,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;IAAC,CAAC;IAK9B,MAAM,CAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW;QAC9D,MAAM,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAChC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;QACtC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;CAE1C"}
+{"version":3,"file":"MutableComplex.js","sourceRoot":"","sources":["../../src/math/MutableComplex.ts"],"names":[],"mappings":"AAAA,OAAO,OAAO,MAAM,cAAc,CAAC;AAUnC,MAAM,CAAC,OAAO,OAAO,cAAe,SAAQ,OAAO;IAoBhD,YAAoB,KAAa,CAAC,EAAE,KAAa,CAAC;QAC/C,KAAK,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAAC,CAAC;IAOZ,MAAM,CAAC,WAAW,CAAE,CAAU;QAClC,OAAO,IAAI,cAAc,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAMpC,MAAM,CAAU,IAAI,CAAE,GAAW;QACrC,OAAO,IAAI,cAAc,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAAC,CAAC;IAMtD,MAAM,CAAU,SAAS,CAAE,GAAW,EAAE,GAAW;QACvD,OAAO,IAAI,cAAc,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IAAC,CAAC;IAOlE,GAAG,CAAE,CAAU;QACnB,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;QACf,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAKb,OAAO,CAAE,EAAU,EAAE,EAAE,GAAG,CAAC;QAC/B,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IAAC,CAAC;IAKX,OAAO,CAAE,GAAW;QACxB,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACxB,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;IAKtB,QAAQ,CAAE,GAAW,EAAE,GAAW;QACtC,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC9B,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;IAO5B,SAAS,CAAE,CAAS;QACxB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;IAAC,CAAC;IAKX,KAAK,CAAE,CAAU;QACrB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC;QAChB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAKd,WAAW,CAAE,CAAS;QAC1B,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;IAAC,CAAC;IAKX,OAAO,CAAE,CAAU;QACvB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC;QAChB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC;IAAC,CAAC;IAKd,SAAS,CAAE,CAAS;QACxB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;QACb,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;IAAC,CAAC;IAKX,KAAK,CAAE,CAAU;QACrB,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAKxC,SAAS,CAAE,CAAS;QACxB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;QACb,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;IAAC,CAAC;IAKX,KAAK,CAAE,CAAU;QACrB,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC;IAAC,CAAC;IAOxC,MAAM,CAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW;QAC9D,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAChC,IAAI,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;IAAC,CAAC;IAK9B,MAAM,CAAE,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,GAAW;QAC9D,MAAM,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;QAChC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;QACtC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;CAE1C"}

package/math/NumApprox.d.ts

@@ -1,3 +1,4 @@
-export declare function argMax_scanArith(f: (x: number) => number, argLo: number, argHi: number, argIncr: number): number;
-export declare function argMax_scanGeom(f: (x: number) => number, argLo: number, argHi: number, argFactor: number): number;
-export declare function argMax_goldenSectionSearch(f: (x: number) => number, argLo: number, argHi: number, tolerance: number): number;
+export declare function argMax_scanArith(f: (x: number) => number, argLo: number, argHi: number, argIncr: number): number;
+export declare function argMax_scanGeom(f: (x: number) => number, argLo: number, argHi: number, argFactor: number): number;
+export declare function argMax_goldenSectionSearch(f: (x: number) => number, argLo: number, argHi: number, tolerance: number): number;
+//# sourceMappingURL=NumApprox.d.ts.map

package/math/NumApprox.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"NumApprox.d.ts","sourceRoot":"","sources":["../../src/math/NumApprox.ts"],"names":[],"mappings":"AAiBA,wBAAgB,gBAAgB,CAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,GAAI,MAAM,CAYhG;AAiBnB,wBAAgB,eAAe,CAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAI,MAAM,CAYjG;AAiBnB,wBAAgB,0BAA0B,CAAE,CAAC,EAAE,CAAC,CAAC,EAAE,MAAM,KAAK,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAI,MAAM,CA2BxE"}

package/math/NumApprox.js

@@ -1,68 +1,68 @@
-export function argMax_scanArith(f, argLo, argHi, argIncr) {
-    let maxArg = NaN;
-    let maxVal = -Infinity;
-    let lastArg = NaN;
-    for (let x = argLo; x < argHi; x += argIncr) {
-        const v = f(x);
-        if (v > maxVal) {
-            maxArg = x;
-            maxVal = v;
-        }
-        lastArg = x;
-    }
-    if (lastArg != argHi && argHi > argLo && f(argHi) > maxVal) {
-        maxArg = argHi;
-    }
-    return maxArg;
-}
-export function argMax_scanGeom(f, argLo, argHi, argFactor) {
-    let maxArg = NaN;
-    let maxVal = -Infinity;
-    let lastArg = NaN;
-    for (let x = argLo; x < argHi; x *= argFactor) {
-        const v = f(x);
-        if (v > maxVal) {
-            maxArg = x;
-            maxVal = v;
-        }
-        lastArg = x;
-    }
-    if (lastArg != argHi && argHi > argLo && f(argHi) > maxVal) {
-        maxArg = argHi;
-    }
-    return maxArg;
-}
-export function argMax_goldenSectionSearch(f, argLo, argHi, tolerance) {
-    const invPhi = (Math.sqrt(5) - 1) / 2;
-    const invPhi2 = 1 - invPhi;
-    let lo = argLo;
-    let hi = argHi;
-    let width = hi - lo;
-    let x1 = lo + width * invPhi2;
-    let x2 = lo + width * invPhi;
-    let y1 = f(x1);
-    let y2 = f(x2);
-    if (width <= tolerance) {
-        return (lo + hi) / 2;
-    }
-    const n = Math.ceil(Math.log(tolerance / width) / Math.log(invPhi));
-    for (let i = 0; i < n; i++) {
-        width *= invPhi;
-        if (y1 >= y2) {
-            hi = x2;
-            x2 = x1;
-            y2 = y1;
-            x1 = lo + width * invPhi2;
-            y1 = f(x1);
-        }
-        else {
-            lo = x1;
-            x1 = x2;
-            y1 = y2;
-            x2 = lo + width * invPhi;
-            y2 = f(x2);
-        }
-    }
-    return (y1 >= y2) ? (lo + x2) / 2 : (x1 + hi) / 2;
-}
+export function argMax_scanArith(f, argLo, argHi, argIncr) {
+    let maxArg = NaN;
+    let maxVal = -Infinity;
+    let lastArg = NaN;
+    for (let x = argLo; x < argHi; x += argIncr) {
+        const v = f(x);
+        if (v > maxVal) {
+            maxArg = x;
+            maxVal = v;
+        }
+        lastArg = x;
+    }
+    if (lastArg != argHi && argHi > argLo && f(argHi) > maxVal) {
+        maxArg = argHi;
+    }
+    return maxArg;
+}
+export function argMax_scanGeom(f, argLo, argHi, argFactor) {
+    let maxArg = NaN;
+    let maxVal = -Infinity;
+    let lastArg = NaN;
+    for (let x = argLo; x < argHi; x *= argFactor) {
+        const v = f(x);
+        if (v > maxVal) {
+            maxArg = x;
+            maxVal = v;
+        }
+        lastArg = x;
+    }
+    if (lastArg != argHi && argHi > argLo && f(argHi) > maxVal) {
+        maxArg = argHi;
+    }
+    return maxArg;
+}
+export function argMax_goldenSectionSearch(f, argLo, argHi, tolerance) {
+    const invPhi = (Math.sqrt(5) - 1) / 2;
+    const invPhi2 = 1 - invPhi;
+    let lo = argLo;
+    let hi = argHi;
+    let width = hi - lo;
+    let x1 = lo + width * invPhi2;
+    let x2 = lo + width * invPhi;
+    let y1 = f(x1);
+    let y2 = f(x2);
+    if (width <= tolerance) {
+        return (lo + hi) / 2;
+    }
+    const n = Math.ceil(Math.log(tolerance / width) / Math.log(invPhi));
+    for (let i = 0; i < n; i++) {
+        width *= invPhi;
+        if (y1 >= y2) {
+            hi = x2;
+            x2 = x1;
+            y2 = y1;
+            x1 = lo + width * invPhi2;
+            y1 = f(x1);
+        }
+        else {
+            lo = x1;
+            x1 = x2;
+            y1 = y2;
+            x2 = lo + width * invPhi;
+            y2 = f(x2);
+        }
+    }
+    return (y1 >= y2) ? (lo + x2) / 2 : (x1 + hi) / 2;
+}
 //# sourceMappingURL=NumApprox.js.map

package/math/NumApprox.js.map

@@ -1 +1 @@
-{"version":3,"file":"NumApprox.js","sourceRoot":"","sources":["../../src/math/NumApprox.ts"],"names":[],"mappings":"AAiBA,MAAM,UAAU,gBAAgB,CAAE,CAAwB,EAAE,KAAa,EAAE,KAAa,EAAE,OAAe;IACtG,IAAI,MAAM,GAAW,GAAG,CAAC;IACzB,IAAI,MAAM,GAAW,CAAC,QAAQ,CAAC;IAC/B,IAAI,OAAO,GAAW,GAAG,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,IAAI,OAAO,EAAE;QAC1C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,IAAI,CAAC,GAAG,MAAM,EAAE;YACb,MAAM,GAAG,CAAC,CAAC;YACX,MAAM,GAAG,CAAC,CAAC;SAAE;QAChB,OAAO,GAAG,CAAC,CAAC;KAAE;IACjB,IAAI,OAAO,IAAI,KAAK,IAAI,KAAK,GAAG,KAAK,IAAI,CAAC,CAAC,KAAK,CAAC,GAAG,MAAM,EAAE;QACzD,MAAM,GAAG,KAAK,CAAC;KAAE;IACpB,OAAO,MAAM,CAAC;AAAC,CAAC;AAiBnB,MAAM,UAAU,eAAe,CAAE,CAAwB,EAAE,KAAa,EAAE,KAAa,EAAE,SAAiB;IACvG,IAAI,MAAM,GAAW,GAAG,CAAC;IACzB,IAAI,MAAM,GAAW,CAAC,QAAQ,CAAC;IAC/B,IAAI,OAAO,GAAW,GAAG,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,IAAI,SAAS,EAAE;QAC5C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,IAAI,CAAC,GAAG,MAAM,EAAE;YACb,MAAM,GAAG,CAAC,CAAC;YACX,MAAM,GAAG,CAAC,CAAC;SAAE;QAChB,OAAO,GAAG,CAAC,CAAC;KAAE;IACjB,IAAI,OAAO,IAAI,KAAK,IAAI,KAAK,GAAG,KAAK,IAAI,CAAC,CAAC,KAAK,CAAC,GAAG,MAAM,EAAE;QACzD,MAAM,GAAG,KAAK,CAAC;KAAE;IACpB,OAAO,MAAM,CAAC;AAAC,CAAC;AAiBnB,MAAM,UAAU,0BAA0B,CAAE,CAAwB,EAAE,KAAa,EAAE,KAAa,EAAE,SAAiB;IAClH,MAAM,MAAM,GAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvC,MAAM,OAAO,GAAG,CAAC,GAAG,MAAM,CAAC;IAC3B,IAAI,EAAE,GAAG,KAAK,CAAC;IACf,IAAI,EAAE,GAAG,KAAK,CAAC;IACf,IAAI,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC;IACpB,IAAI,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,OAAO,CAAC;IAC9B,IAAI,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,MAAM,CAAC;IAC7B,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;IACf,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;IACf,IAAI,KAAK,IAAI,SAAS,EAAE;QACrB,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;KAAE;IAC1B,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC;IACpE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QACzB,KAAK,IAAI,MAAM,CAAC;QAChB,IAAI,EAAE,IAAI,EAAE,EAAE;YACX,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,OAAO,CAAC;YAC1B,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;SAAE;aACV;YACH,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,MAAM,CAAC;YACzB,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;SAAE;KAAC;IACpB,OAAO,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;AAAC,CAAC"}
+{"version":3,"file":"NumApprox.js","sourceRoot":"","sources":["../../src/math/NumApprox.ts"],"names":[],"mappings":"AAiBA,MAAM,UAAU,gBAAgB,CAAE,CAAwB,EAAE,KAAa,EAAE,KAAa,EAAE,OAAe;IACtG,IAAI,MAAM,GAAW,GAAG,CAAC;IACzB,IAAI,MAAM,GAAW,CAAC,QAAQ,CAAC;IAC/B,IAAI,OAAO,GAAW,GAAG,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,IAAI,OAAO,EAAE,CAAC;QAC3C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC;YACd,MAAM,GAAG,CAAC,CAAC;YACX,MAAM,GAAG,CAAC,CAAC;QAAC,CAAC;QAChB,OAAO,GAAG,CAAC,CAAC;IAAC,CAAC;IACjB,IAAI,OAAO,IAAI,KAAK,IAAI,KAAK,GAAG,KAAK,IAAI,CAAC,CAAC,KAAK,CAAC,GAAG,MAAM,EAAE,CAAC;QAC1D,MAAM,GAAG,KAAK,CAAC;IAAC,CAAC;IACpB,OAAO,MAAM,CAAC;AAAC,CAAC;AAiBnB,MAAM,UAAU,eAAe,CAAE,CAAwB,EAAE,KAAa,EAAE,KAAa,EAAE,SAAiB;IACvG,IAAI,MAAM,GAAW,GAAG,CAAC;IACzB,IAAI,MAAM,GAAW,CAAC,QAAQ,CAAC;IAC/B,IAAI,OAAO,GAAW,GAAG,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,IAAI,SAAS,EAAE,CAAC;QAC7C,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACf,IAAI,CAAC,GAAG,MAAM,EAAE,CAAC;YACd,MAAM,GAAG,CAAC,CAAC;YACX,MAAM,GAAG,CAAC,CAAC;QAAC,CAAC;QAChB,OAAO,GAAG,CAAC,CAAC;IAAC,CAAC;IACjB,IAAI,OAAO,IAAI,KAAK,IAAI,KAAK,GAAG,KAAK,IAAI,CAAC,CAAC,KAAK,CAAC,GAAG,MAAM,EAAE,CAAC;QAC1D,MAAM,GAAG,KAAK,CAAC;IAAC,CAAC;IACpB,OAAO,MAAM,CAAC;AAAC,CAAC;AAiBnB,MAAM,UAAU,0BAA0B,CAAE,CAAwB,EAAE,KAAa,EAAE,KAAa,EAAE,SAAiB;IAClH,MAAM,MAAM,GAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvC,MAAM,OAAO,GAAG,CAAC,GAAG,MAAM,CAAC;IAC3B,IAAI,EAAE,GAAG,KAAK,CAAC;IACf,IAAI,EAAE,GAAG,KAAK,CAAC;IACf,IAAI,KAAK,GAAG,EAAE,GAAG,EAAE,CAAC;IACpB,IAAI,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,OAAO,CAAC;IAC9B,IAAI,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,MAAM,CAAC;IAC7B,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;IACf,IAAI,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;IACf,IAAI,KAAK,IAAI,SAAS,EAAE,CAAC;QACtB,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC;IAC1B,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC;IACpE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC1B,KAAK,IAAI,MAAM,CAAC;QAChB,IAAI,EAAE,IAAI,EAAE,EAAE,CAAC;YACZ,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,OAAO,CAAC;YAC1B,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAAC,CAAC;aACV,CAAC;YACJ,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,GAAG,EAAE,GAAG,KAAK,GAAG,MAAM,CAAC;YACzB,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC;QAAC,CAAC;IAAA,CAAC;IACpB,OAAO,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;AAAC,CAAC"}

package/math/PolyReal.d.ts

@@ -1,13 +1,14 @@
-import Complex from "./Complex.js";
-export declare function evaluateReal(a: ArrayLike<number>, x: number): number;
-export declare function evaluateComplex(a: ArrayLike<number>, x: Complex): Complex;
-export declare function expand(zeros: ArrayLike<number>): Float64Array;
-export declare function compareEqual(a1: ArrayLike<number>, a2: ArrayLike<number>, eps?: number): boolean;
-export declare function add(a1: ArrayLike<number>, a2: ArrayLike<number>, eps?: number): Float64Array;
-export declare function multiply(a1: ArrayLike<number>, a2: ArrayLike<number>, eps?: number): Float64Array;
-export declare function divide(a1r: ArrayLike<number>, a2r: ArrayLike<number>, eps?: number): Float64Array[];
-export declare function gcd(a1: ArrayLike<number>, a2: ArrayLike<number>, eps?: number): Float64Array;
-export declare function evaluateFractionComplex(f: ArrayLike<number>[], x: Complex): Complex;
-export declare function addFractions(f1: ArrayLike<number>[], f2: ArrayLike<number>[], eps?: number): Float64Array[];
-export declare function multiplyFractions(f1: ArrayLike<number>[], f2: ArrayLike<number>[], eps?: number): Float64Array[];
-export declare function normalizeFraction(f: ArrayLike<number>[], eps?: number): Float64Array[];
+import Complex from "./Complex.ts";
+export declare function evaluateReal(a: ArrayLike<number>, x: number): number;
+export declare function evaluateComplex(a: ArrayLike<number>, x: Complex): Complex;
+export declare function expand(zeros: ArrayLike<number>): Float64Array;
+export declare function compareEqual(a1: ArrayLike<number>, a2: ArrayLike<number>, eps?: number): boolean;
+export declare function add(a1: ArrayLike<number>, a2: ArrayLike<number>, eps?: number): Float64Array;
+export declare function multiply(a1: ArrayLike<number>, a2: ArrayLike<number>, eps?: number): Float64Array;
+export declare function divide(a1r: ArrayLike<number>, a2r: ArrayLike<number>, eps?: number): Float64Array[];
+export declare function gcd(a1: ArrayLike<number>, a2: ArrayLike<number>, eps?: number): Float64Array;
+export declare function evaluateFractionComplex(f: ArrayLike<number>[], x: Complex): Complex;
+export declare function addFractions(f1: ArrayLike<number>[], f2: ArrayLike<number>[], eps?: number): Float64Array[];
+export declare function multiplyFractions(f1: ArrayLike<number>[], f2: ArrayLike<number>[], eps?: number): Float64Array[];
+export declare function normalizeFraction(f: ArrayLike<number>[], eps?: number): Float64Array[];
+//# sourceMappingURL=PolyReal.d.ts.map

package/math/PolyReal.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"PolyReal.d.ts","sourceRoot":"","sources":["../../src/math/PolyReal.ts"],"names":[],"mappings":"AAgBA,OAAO,OAAO,MAAM,cAAc,CAAC;AAanC,wBAAgB,YAAY,CAAE,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,CAAC,EAAE,MAAM,GAAI,MAAM,CAQzD;AAYd,wBAAgB,eAAe,CAAE,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,CAAC,EAAE,OAAO,GAAI,OAAO,CAQ9D;AAKd,wBAAgB,MAAM,CAAE,KAAK,EAAE,SAAS,CAAC,MAAM,CAAC,GAAI,YAAY,CAQlD;AAKd,wBAAgB,YAAY,CAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,GAAG,SAAI,GAAI,OAAO,CAS7E;AAKjB,wBAAgB,GAAG,CAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,GAAG,SAAI,GAAI,YAAY,CAShE;AAK1B,wBAAgB,QAAQ,CAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,GAAG,SAAI,GAAI,YAAY,CAgBrE;AAM1B,wBAAgB,MAAM,CAAE,GAAG,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,GAAG,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,GAAG,SAAI,GAAI,YAAY,EAAE,CAwB/D;AAKlC,wBAAgB,GAAG,CAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,GAAG,SAAI,GAAI,YAAY,CAe1E;AA+ChB,wBAAgB,uBAAuB,CAAE,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,CAAC,EAAE,OAAO,GAAI,OAAO,CAG/D;AAKvB,wBAAgB,YAAY,CAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,GAAG,SAAI,GAAI,YAAY,EAAE,CAkB/E;AAK1B,wBAAgB,iBAAiB,CAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,GAAG,SAAI,GAAI,YAAY,EAAE,CAGpF;AAK1B,wBAAgB,iBAAiB,CAAE,CAAC,EAAE,SAAS,CAAC,MAAM,CAAC,EAAE,EAAE,GAAG,SAAI,GAAI,YAAY,EAAE,CAQ1D"}