@danielsimonjr/mathts-functions 0.43.3 → 0.44.0

package/dist/index.d.ts

@@ -87,6 +87,7 @@ export { rfft, irfft, fftshift, ifftshift, fftfreq, rfftfreq, fftn } from './sig
 export { cheby1, cheby2, ellip, sosfilt, zpk2sos, bilinear, buttord } from './signal/iir-design.js';
 export { firwinBandpass, firls, remez, savgol, wiener, deconvolve, } from './signal/fir-smoothing.js';
 export type { DeconvolveResult } from './signal/fir-smoothing.js';
+export type { RemezType } from './signal/remez-exchange.js';
 export { idwt, wavedec, waverec, cwt } from './signal/wavelets.js';
 export { findPeaks, peakWidths, csd, coherence, stft, istft, decimate, } from './signal/spectral-peaks.js';
 export type { FindPeaksOptions, CsdOptions, StftOptions, StftResult, } from './signal/spectral-peaks.js';

package/dist/index.d.ts.map

@@ -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,MAAM,qBAAqB,CAAC;AAClE,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;AAOxF,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;AAIlE,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;AAM9F,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"}
+{"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,MAAM,qBAAqB,CAAC;AAClE,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;AAOxF,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;AAM9F,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"}

package/dist/index.js

@@ -20444,7 +20444,7 @@ var createLgamma = /* @__PURE__ */ factory(
       }
     });
     function lgammaComplex2(n) {
-      const TWOPI = 6.283185307179586;
+      const TWOPI2 = 6.283185307179586;
       const LOGPI = 1.1447298858494002;
       const REFLECTION = 0.1;
       if (n.isNaN()) {
@@ -20454,7 +20454,7 @@ var createLgamma = /* @__PURE__ */ factory(
       } else if (n.re >= SMALL_RE || Math.abs(n.im) >= SMALL_IM) {
         return lgammaStirling(n);
       } else if (n.re <= REFLECTION) {
-        const tmp = copysign(TWOPI, n.im) * Math.floor(0.5 * n.re + 0.25);
+        const tmp = copysign(TWOPI2, n.im) * Math.floor(0.5 * n.re + 0.25);
         const a = n.mul(Math.PI).sin().log();
         const b = lgammaComplex2(new Complex13(1 - n.re, -n.im));
         return new Complex13(LOGPI, tmp).sub(a).sub(b);
@@ -48860,20 +48860,156 @@ function bilinear(bIn, aIn, fs) {
   const a0 = aDesc[0];
   return { b: bDesc.map((v) => v / a0), a: aDesc.map((v) => v / a0) };
 }
+function bandStopObj(wp, ind, passb, stopb, gpass, gstop) {
+  const passbC = [passb[0], passb[1]];
+  passbC[ind] = wp;
+  const nat0 = stopb[0] * (passbC[0] - passbC[1]) / (stopb[0] * stopb[0] - passbC[0] * passbC[1]);
+  const nat1 = stopb[1] * (passbC[0] - passbC[1]) / (stopb[1] * stopb[1] - passbC[0] * passbC[1]);
+  const nat = Math.min(Math.abs(nat0), Math.abs(nat1));
+  const GSTOP = Math.pow(10, 0.1 * Math.abs(gstop));
+  const GPASS = Math.pow(10, 0.1 * Math.abs(gpass));
+  return Math.log10((GSTOP - 1) / (GPASS - 1)) / (2 * Math.log10(nat));
+}
+function fminbound(f, x1, x2) {
+  const xatol = 1e-5;
+  const maxfun = 500;
+  const sqrtEps = Math.sqrt(22e-17);
+  const goldenMean = 0.5 * (3 - Math.sqrt(5));
+  const sign3 = (v) => v > 0 ? 1 : v < 0 ? -1 : 0;
+  let a = x1;
+  let b = x2;
+  let fulc = a + goldenMean * (b - a);
+  let nfc = fulc;
+  let xf = fulc;
+  let rat = 0;
+  let e = 0;
+  let x = xf;
+  let fx = f(x);
+  let num2 = 1;
+  let fu = Infinity;
+  let ffulc = fx;
+  let fnfc = fx;
+  let xm = 0.5 * (a + b);
+  let tol1 = sqrtEps * Math.abs(xf) + xatol / 3;
+  let tol2 = 2 * tol1;
+  while (Math.abs(xf - xm) > tol2 - 0.5 * (b - a)) {
+    let golden = true;
+    if (Math.abs(e) > tol1) {
+      golden = false;
+      let r = (xf - nfc) * (fx - ffulc);
+      let q = (xf - fulc) * (fx - fnfc);
+      let p = (xf - fulc) * q - (xf - nfc) * r;
+      q = 2 * (q - r);
+      if (q > 0) p = -p;
+      q = Math.abs(q);
+      r = e;
+      e = rat;
+      if (Math.abs(p) < Math.abs(0.5 * q * r) && p > q * (a - xf) && p < q * (b - xf)) {
+        rat = p / q;
+        x = xf + rat;
+        if (x - a < tol2 || b - x < tol2) {
+          const si2 = sign3(xm - xf) + (xm - xf === 0 ? 1 : 0);
+          rat = tol1 * si2;
+        }
+      } else {
+        golden = true;
+      }
+    }
+    if (golden) {
+      e = xf >= xm ? a - xf : b - xf;
+      rat = goldenMean * e;
+    }
+    const si = sign3(rat) + (rat === 0 ? 1 : 0);
+    x = xf + si * Math.max(Math.abs(rat), tol1);
+    fu = f(x);
+    num2 += 1;
+    if (fu <= fx) {
+      if (x >= xf) a = xf;
+      else b = xf;
+      fulc = nfc;
+      ffulc = fnfc;
+      nfc = xf;
+      fnfc = fx;
+      xf = x;
+      fx = fu;
+    } else {
+      if (x < xf) a = x;
+      else b = x;
+      if (fu <= fnfc || nfc === xf) {
+        fulc = nfc;
+        ffulc = fnfc;
+        nfc = x;
+        fnfc = fu;
+      } else if (fu <= ffulc || fulc === xf || fulc === nfc) {
+        fulc = x;
+        ffulc = fu;
+      }
+    }
+    xm = 0.5 * (a + b);
+    tol1 = sqrtEps * Math.abs(xf) + xatol / 3;
+    tol2 = 2 * tol1;
+    if (num2 >= maxfun) break;
+  }
+  return xf;
+}
 function buttord(wp, ws, gpass, gstop) {
-  if (!(wp > 0 && wp < 1)) throw new Error("buttord: wp must satisfy 0 < wp < 1");
-  if (!(ws > 0 && ws < 1)) throw new Error("buttord: ws must satisfy 0 < ws < 1");
-  if (wp === ws) throw new Error("buttord: wp and ws must differ");
-  const isLow = wp < ws;
-  const passb = Math.tan(Math.PI * wp / 2);
-  const stopb = Math.tan(Math.PI * ws / 2);
-  const nat = isLow ? stopb / passb : passb / stopb;
+  const wpArr = typeof wp === "number" ? [wp] : [wp[0], wp[1]];
+  const wsArr = typeof ws === "number" ? [ws] : [ws[0], ws[1]];
+  if (wpArr.length !== wsArr.length)
+    throw new Error("buttord: wp and ws must both be scalars or both be [low, high] pairs");
+  for (const v of [...wpArr, ...wsArr])
+    if (!(v > 0 && v < 1)) throw new Error("buttord: all edge frequencies must satisfy 0 < f < 1");
+  let filterType = 2 * (wpArr.length - 1) + 1;
+  if (wpArr[0] >= wsArr[0]) filterType += 1;
+  const passb = wpArr.map((f) => Math.tan(Math.PI * f / 2));
+  const stopb = wsArr.map((f) => Math.tan(Math.PI * f / 2));
   const GSTOP = Math.pow(10, 0.1 * Math.abs(gstop));
   const GPASS = Math.pow(10, 0.1 * Math.abs(gpass));
+  let nat;
+  if (filterType === 1) {
+    nat = stopb[0] / passb[0];
+  } else if (filterType === 2) {
+    nat = passb[0] / stopb[0];
+  } else if (filterType === 3) {
+    const sb = [stopb[0], stopb[1]];
+    const pb0 = [passb[0], passb[1]];
+    const wp0 = fminbound(
+      (v) => bandStopObj(v, 0, pb0, sb, gpass, gstop),
+      passb[0],
+      stopb[0] - 1e-12
+    );
+    const wp1 = fminbound(
+      (v) => bandStopObj(v, 1, pb0, sb, gpass, gstop),
+      stopb[1] + 1e-12,
+      passb[1]
+    );
+    passb[0] = wp0;
+    passb[1] = wp1;
+    const nat0 = stopb[0] * (passb[0] - passb[1]) / (stopb[0] * stopb[0] - passb[0] * passb[1]);
+    const nat1 = stopb[1] * (passb[0] - passb[1]) / (stopb[1] * stopb[1] - passb[0] * passb[1]);
+    nat = Math.min(Math.abs(nat0), Math.abs(nat1));
+  } else {
+    const nat0 = (stopb[0] * stopb[0] - passb[0] * passb[1]) / (stopb[0] * (passb[0] - passb[1]));
+    const nat1 = (stopb[1] * stopb[1] - passb[0] * passb[1]) / (stopb[1] * (passb[0] - passb[1]));
+    nat = Math.min(Math.abs(nat0), Math.abs(nat1));
+  }
   const N = Math.ceil(Math.log10((GSTOP - 1) / (GPASS - 1)) / (2 * Math.log10(nat)));
   const W0 = Math.pow(GPASS - 1, -1 / (2 * N));
-  const WN = isLow ? W0 * passb : passb / W0;
-  return { N, Wn: Math.atan(WN) * 2 / Math.PI };
+  const toDigital = (WN) => Math.atan(WN) * 2 / Math.PI;
+  if (filterType === 1) return { N, Wn: toDigital(W0 * passb[0]) };
+  if (filterType === 2) return { N, Wn: toDigital(passb[0] / W0) };
+  if (filterType === 3) {
+    const discr = Math.sqrt(Math.pow(passb[1] - passb[0], 2) + 4 * W0 * W0 * passb[0] * passb[1]);
+    const wn0 = (passb[1] - passb[0] + discr) / (2 * W0);
+    const wn1 = (passb[1] - passb[0] - discr) / (2 * W0);
+    const sorted2 = [Math.abs(wn0), Math.abs(wn1)].sort((x, y) => x - y);
+    return { N, Wn: [toDigital(sorted2[0]), toDigital(sorted2[1])] };
+  }
+  const w0s = [-W0, W0].map(
+    (s) => -s * ((passb[1] - passb[0]) / 2) + Math.sqrt(s * s / 4 * Math.pow(passb[1] - passb[0], 2) + passb[0] * passb[1])
+  );
+  const sorted = w0s.map(Math.abs).sort((x, y) => x - y);
+  return { N, Wn: [toDigital(sorted[0]), toDigital(sorted[1])] };
 }
 function toComplex(x) {
   return x instanceof Complex12 ? x : new Complex12(x, 0);
@@ -48978,6 +49114,551 @@ function convDirect(x, h, mode2 = "full") {
   }
 }
 
+// src/signal/remez-exchange.ts
+var PI = Math.PI;
+var TWOPI = 2 * Math.PI;
+function lagrangeInterp2(k, n, m, x) {
+  let retval = 1;
+  const q = x[k];
+  for (let l = 1; l <= m; l++) {
+    for (let j = l; j <= n; j += m) {
+      if (j !== k) retval *= 2 * (q - x[j]);
+    }
+  }
+  return 1 / retval;
+}
+function freqEval(k, n, grid, x, y, ad) {
+  let d = 0;
+  let p = 0;
+  const xf = Math.cos(TWOPI * grid[k]);
+  for (let j = 1; j <= n; j++) {
+    const c = ad[j] / (xf - x[j]);
+    d += c;
+    p += c * y[j];
+  }
+  return p / d;
+}
+function eff(freq, fx, lband, jtype) {
+  return jtype !== 2 ? fx[lband] : fx[lband] * freq;
+}
+function wate(freq, fx, wtx, lband, jtype) {
+  if (jtype !== 2) return wtx[lband];
+  if (fx[lband] >= 1e-4) return wtx[lband] / freq;
+  return wtx[lband];
+}
+function remezCore(des, grid, edge, wt, ngrid, nbands, iext, alpha, nfcns, itrmax) {
+  const size2 = ngrid + 2 * nfcns + 8;
+  const a = new Array(size2).fill(0);
+  const p = new Array(size2).fill(0);
+  const q = new Array(size2).fill(0);
+  const ad = new Array(size2).fill(0);
+  const x = new Array(size2).fill(0);
+  const y = new Array(size2).fill(0);
+  let devl = -1;
+  const nz = nfcns + 1;
+  const nzz = nfcns + 2;
+  let niter = 0;
+  let dev = 0;
+  let comp = 0;
+  let ynz = 0;
+  let y1 = 0;
+  let jchnge = 0;
+  let k1 = 0;
+  let knz = 0;
+  let klow = 0;
+  let nut = 0;
+  let nut1 = 0;
+  let nu = 0;
+  let j = 0;
+  let l = 0;
+  let kup = 0;
+  let luck = 0;
+  let err = 0;
+  let state = 100;
+  let done = false;
+  let failed = false;
+  outer: for (; ; ) {
+    switch (state) {
+      case 100: {
+        iext[nzz] = ngrid + 1;
+        ++niter;
+        if (niter > itrmax) {
+          done = true;
+          break outer;
+        }
+        for (j = 1; j <= nz; j++) x[j] = Math.cos(grid[iext[j]] * TWOPI);
+        const jet = Math.floor((nfcns - 1) / 15) + 1;
+        for (j = 1; j <= nz; j++) ad[j] = lagrangeInterp2(j, nz, jet, x);
+        let dnum = 0;
+        let dden2 = 0;
+        let ksign = 1;
+        for (j = 1; j <= nz; j++) {
+          l = iext[j];
+          dnum += ad[j] * des[l];
+          dden2 += ksign * ad[j] / wt[l];
+          ksign = -ksign;
+        }
+        dev = dnum / dden2;
+        nu = 1;
+        if (dev > 0) nu = -1;
+        dev = -nu * dev;
+        ksign = nu;
+        for (j = 1; j <= nz; j++) {
+          l = iext[j];
+          y[j] = des[l] + ksign * dev / wt[l];
+          ksign = -ksign;
+        }
+        if (dev <= devl) {
+          failed = true;
+          break outer;
+        }
+        devl = dev;
+        jchnge = 0;
+        k1 = iext[1];
+        knz = iext[nz];
+        klow = 0;
+        nut = -nu;
+        j = 1;
+        state = 200;
+        continue;
+      }
+      case 200: {
+        if (j === nzz) ynz = comp;
+        if (j >= nzz) {
+          state = 300;
+          continue;
+        }
+        kup = iext[j + 1];
+        l = iext[j] + 1;
+        nut = -nut;
+        if (j === 2) y1 = comp;
+        comp = dev;
+        if (l >= kup) {
+          state = 220;
+          continue;
+        }
+        err = (freqEval(l, nz, grid, x, y, ad) - des[l]) * wt[l];
+        if (nut * err - comp <= 0) {
+          state = 220;
+          continue;
+        }
+        comp = nut * err;
+        state = 210;
+        continue;
+      }
+      case 210: {
+        if (++l >= kup) {
+          state = 215;
+          continue;
+        }
+        err = (freqEval(l, nz, grid, x, y, ad) - des[l]) * wt[l];
+        if (nut * err - comp <= 0) {
+          state = 215;
+          continue;
+        }
+        comp = nut * err;
+        state = 210;
+        continue;
+      }
+      case 215: {
+        iext[j++] = l - 1;
+        klow = l - 1;
+        ++jchnge;
+        state = 200;
+        continue;
+      }
+      case 220: {
+        --l;
+        state = 225;
+        continue;
+      }
+      case 225: {
+        if (--l <= klow) {
+          state = 250;
+          continue;
+        }
+        err = (freqEval(l, nz, grid, x, y, ad) - des[l]) * wt[l];
+        if (nut * err - comp > 0) {
+          state = 230;
+          continue;
+        }
+        if (jchnge <= 0) {
+          state = 225;
+          continue;
+        }
+        state = 260;
+        continue;
+      }
+      case 230: {
+        comp = nut * err;
+        state = 235;
+        continue;
+      }
+      case 235: {
+        if (--l <= klow) {
+          state = 240;
+          continue;
+        }
+        err = (freqEval(l, nz, grid, x, y, ad) - des[l]) * wt[l];
+        if (nut * err - comp <= 0) {
+          state = 240;
+          continue;
+        }
+        comp = nut * err;
+        state = 235;
+        continue;
+      }
+      case 240: {
+        klow = iext[j];
+        iext[j] = l + 1;
+        ++j;
+        ++jchnge;
+        state = 200;
+        continue;
+      }
+      case 250: {
+        l = iext[j] + 1;
+        if (jchnge > 0) {
+          state = 215;
+          continue;
+        }
+        state = 255;
+        continue;
+      }
+      case 255: {
+        if (++l >= kup) {
+          state = 260;
+          continue;
+        }
+        err = (freqEval(l, nz, grid, x, y, ad) - des[l]) * wt[l];
+        if (nut * err - comp <= 0) {
+          state = 255;
+          continue;
+        }
+        comp = nut * err;
+        state = 210;
+        continue;
+      }
+      case 260: {
+        klow = iext[j++];
+        state = 200;
+        continue;
+      }
+      case 300: {
+        if (j > nzz) {
+          state = 320;
+          continue;
+        }
+        if (k1 > iext[1]) k1 = iext[1];
+        if (knz < iext[nz]) knz = iext[nz];
+        nut1 = nut;
+        nut = -nu;
+        l = 0;
+        kup = k1;
+        comp = ynz * 1.00001;
+        luck = 1;
+        state = 310;
+        continue;
+      }
+      case 310: {
+        if (++l >= kup) {
+          state = 315;
+          continue;
+        }
+        err = (freqEval(l, nz, grid, x, y, ad) - des[l]) * wt[l];
+        if (nut * err - comp <= 0) {
+          state = 310;
+          continue;
+        }
+        comp = nut * err;
+        j = nzz;
+        state = 210;
+        continue;
+      }
+      case 315: {
+        luck = 6;
+        state = 325;
+        continue;
+      }
+      case 320: {
+        if (luck > 9) {
+          state = 350;
+          continue;
+        }
+        if (comp > y1) y1 = comp;
+        k1 = iext[nzz];
+        state = 325;
+        continue;
+      }
+      case 325: {
+        l = ngrid + 1;
+        klow = knz;
+        nut = -nut1;
+        comp = y1 * 1.00001;
+        state = 330;
+        continue;
+      }
+      case 330: {
+        if (--l <= klow) {
+          state = 340;
+          continue;
+        }
+        err = (freqEval(l, nz, grid, x, y, ad) - des[l]) * wt[l];
+        if (nut * err - comp <= 0) {
+          state = 330;
+          continue;
+        }
+        j = nzz;
+        comp = nut * err;
+        luck = luck + 10;
+        state = 235;
+        continue;
+      }
+      case 340: {
+        if (luck === 6) {
+          state = 370;
+          continue;
+        }
+        for (j = 1; j <= nfcns; j++) iext[nzz - j] = iext[nz - j];
+        iext[1] = k1;
+        state = 100;
+        continue;
+      }
+      case 350: {
+        const kn = iext[nzz];
+        for (j = 1; j <= nfcns; j++) iext[j] = iext[j + 1];
+        iext[nz] = kn;
+        state = 100;
+        continue;
+      }
+      case 370: {
+        if (jchnge > 0) {
+          state = 100;
+          continue;
+        }
+        done = true;
+        break outer;
+      }
+      default:
+        break outer;
+    }
+  }
+  if (failed) return { status: -1, dev, niter };
+  if (!done) return { status: -1, dev, niter };
+  const nm1 = nfcns - 1;
+  const fsh = 1e-6;
+  const gtemp = grid[1];
+  x[nzz] = -2;
+  const cn = 2 * nfcns - 1;
+  const delf = 1 / cn;
+  l = 1;
+  let kkk = 0;
+  if (edge[1] === 0 && edge[2 * nbands] === 0.5) kkk = 1;
+  if (nfcns <= 3) kkk = 1;
+  let aa = 0;
+  let bb = 0;
+  if (kkk !== 1) {
+    const dtemp = Math.cos(TWOPI * grid[1]);
+    const dnum = Math.cos(TWOPI * grid[ngrid]);
+    aa = 2 / (dtemp - dnum);
+    bb = -(dtemp + dnum) / (dtemp - dnum);
+  }
+  for (j = 1; j <= nfcns; j++) {
+    let ft = (j - 1) * delf;
+    let xt = Math.cos(TWOPI * ft);
+    if (kkk !== 1) {
+      xt = (xt - bb) / aa;
+      ft = Math.acos(xt) / TWOPI;
+    }
+    let assigned = false;
+    for (; ; ) {
+      const xe = x[l];
+      if (xt > xe) {
+        if (xt - xe < fsh) {
+          a[j] = y[l];
+          assigned = true;
+          break;
+        }
+        grid[1] = ft;
+        a[j] = freqEval(1, nz, grid, x, y, ad);
+        assigned = true;
+        break;
+      }
+      if (xe - xt < fsh) {
+        a[j] = y[l];
+        assigned = true;
+        break;
+      }
+      ++l;
+    }
+    if (assigned && l > 1) l = l - 1;
+  }
+  grid[1] = gtemp;
+  const dden = TWOPI / cn;
+  for (j = 1; j <= nfcns; j++) {
+    let dtemp = 0;
+    const dnum = (j - 1) * dden;
+    if (nm1 >= 1) {
+      for (let k = 1; k <= nm1; k++) dtemp += a[k + 1] * Math.cos(dnum * k);
+    }
+    alpha[j] = 2 * dtemp + a[1];
+  }
+  for (j = 2; j <= nfcns; j++) alpha[j] *= 2 / cn;
+  alpha[1] /= cn;
+  if (kkk !== 1) {
+    p[1] = 2 * alpha[nfcns] * bb + alpha[nm1];
+    p[2] = 2 * aa * alpha[nfcns];
+    q[1] = alpha[nfcns - 2] - alpha[nfcns];
+    for (j = 2; j <= nm1; j++) {
+      if (j >= nm1) {
+        aa *= 0.5;
+        bb *= 0.5;
+      }
+      p[j + 1] = 0;
+      for (let k = 1; k <= j; k++) {
+        a[k] = p[k];
+        p[k] = 2 * bb * a[k];
+      }
+      p[2] += a[1] * 2 * aa;
+      const jm1 = j - 1;
+      for (let k = 1; k <= jm1; k++) p[k] += q[k] + aa * a[k + 1];
+      const jp1 = j + 1;
+      for (let k = 3; k <= jp1; k++) p[k] += aa * a[k - 1];
+      if (j !== nm1) {
+        for (let k = 1; k <= j; k++) q[k] = -a[k];
+        q[1] += alpha[nfcns - 1 - j];
+      }
+    }
+    for (j = 1; j <= nfcns; j++) alpha[j] = p[j];
+  }
+  if (nfcns <= 3) {
+    alpha[nfcns + 1] = 0;
+    alpha[nfcns + 2] = 0;
+  }
+  return { status: 0, dev, niter };
+}
+function remezExchange(numtaps, bands, desired, weight, type = "bandpass", maxiter = 25, gridDensity = 16) {
+  if (!Number.isInteger(numtaps) || numtaps < 1)
+    throw new Error("remez: numtaps must be a positive integer");
+  if (bands.length < 2 || bands.length % 2 !== 0)
+    throw new Error("remez: bands must be a non-empty list of band-edge pairs");
+  const numbands = bands.length / 2;
+  if (desired.length !== numbands)
+    throw new Error("remez: desired must have length bands.length / 2 (one value per band)");
+  const w = weight ?? new Array(numbands).fill(1);
+  if (w.length !== numbands)
+    throw new Error("remez: weight must have length bands.length / 2 (one value per band)");
+  for (let i = 0; i < bands.length; i++) {
+    if (!(bands[i] >= 0 && bands[i] <= 0.5))
+      throw new Error("remez: band edges must be in [0, 0.5] (fs = 1, 0.5 = Nyquist)");
+    if (i > 0 && bands[i] < bands[i - 1])
+      throw new Error("remez: bands must be monotonically non-decreasing");
+  }
+  const jtype = type === "bandpass" ? 1 : type === "differentiator" ? 2 : 3;
+  const edge = [0, ...bands];
+  const fx = [0, ...desired];
+  const wtx = [0, ...w];
+  const dimsize = Math.ceil(numtaps / 2 + 2);
+  const wrksize = gridDensity * dimsize;
+  const neg = jtype === 1 ? 0 : 1;
+  const nodd = numtaps % 2;
+  let nfcns = Math.floor(numtaps / 2);
+  if (nodd === 1 && neg === 0) nfcns += 1;
+  const des = new Array(wrksize + 2).fill(0);
+  const grid = new Array(wrksize + 2).fill(0);
+  const wt = new Array(wrksize + 2).fill(0);
+  const alpha = new Array(dimsize + 3).fill(0);
+  const iext = new Array(dimsize + 3).fill(0);
+  let delf = gridDensity * nfcns;
+  delf = 0.5 / delf;
+  grid[1] = edge[1];
+  if (neg !== 0 && edge[1] < delf) grid[1] = delf;
+  let j = 1;
+  let l = 1;
+  let lband = 1;
+  for (; ; ) {
+    const fup = edge[l + 1];
+    for (; ; ) {
+      const temp2 = grid[j];
+      des[j] = eff(temp2, fx, lband, jtype);
+      wt[j] = wate(temp2, fx, wtx, lband, jtype);
+      if (++j > wrksize) throw new Error("remez: grid too dense / too many points");
+      grid[j] = temp2 + delf;
+      if (grid[j] > fup) break;
+    }
+    grid[j - 1] = fup;
+    des[j - 1] = eff(fup, fx, lband, jtype);
+    wt[j - 1] = wate(fup, fx, wtx, lband, jtype);
+    ++lband;
+    l += 2;
+    if (lband > numbands) break;
+    grid[j] = edge[l];
+  }
+  let ngrid = j - 1;
+  if (neg === nodd && grid[ngrid] > 0.5 - delf) --ngrid;
+  if (ngrid < nfcns + 1)
+    throw new Error("remez: band too narrow for the requested numtaps (grid underflow)");
+  if (neg <= 0) {
+    if (nodd !== 1) {
+      for (j = 1; j <= ngrid; j++) {
+        const change = Math.cos(PI * grid[j]);
+        des[j] = des[j] / change;
+        wt[j] = wt[j] * change;
+      }
+    }
+  } else {
+    if (nodd !== 1) {
+      for (j = 1; j <= ngrid; j++) {
+        const change = Math.sin(PI * grid[j]);
+        des[j] = des[j] / change;
+        wt[j] = wt[j] * change;
+      }
+    } else {
+      for (j = 1; j <= ngrid; j++) {
+        const change = Math.sin(TWOPI * grid[j]);
+        des[j] = des[j] / change;
+        wt[j] = wt[j] * change;
+      }
+    }
+  }
+  const temp = (ngrid - 1) / nfcns;
+  for (j = 1; j <= nfcns; j++) iext[j] = Math.floor((j - 1) * temp) + 1;
+  iext[nfcns + 1] = ngrid;
+  const res = remezCore(des, grid, edge, wt, ngrid, numbands, iext, alpha, nfcns, maxiter);
+  if (res.status < 0) throw new Error("remez: Remez exchange failed to converge");
+  const h = new Array(numtaps + 1).fill(0);
+  const nz = nfcns + 1;
+  const nm1 = nfcns - 1;
+  if (neg <= 0) {
+    if (nodd !== 0) {
+      for (j = 1; j <= nm1; j++) h[j] = 0.5 * alpha[nz - j];
+      h[nfcns] = alpha[1];
+    } else {
+      h[1] = 0.25 * alpha[nfcns];
+      for (j = 2; j <= nm1; j++) h[j] = 0.25 * (alpha[nz - j] + alpha[nfcns + 2 - j]);
+      h[nfcns] = 0.5 * alpha[1] + 0.25 * alpha[2];
+    }
+  } else {
+    if (nodd !== 0) {
+      h[1] = 0.25 * alpha[nfcns];
+      h[2] = 0.25 * alpha[nm1];
+      for (j = 3; j <= nm1; j++) h[j] = 0.25 * (alpha[nz - j] - alpha[nfcns + 3 - j]);
+      h[nfcns] = 0.5 * alpha[1] - 0.25 * alpha[3];
+      h[nz] = 0;
+    } else {
+      h[1] = 0.25 * alpha[nfcns];
+      for (j = 2; j <= nm1; j++) h[j] = 0.25 * (alpha[nz - j] - alpha[nfcns + 2 - j]);
+      h[nfcns] = 0.5 * alpha[1] - 0.25 * alpha[2];
+    }
+  }
+  for (j = 1; j <= nfcns; j++) {
+    const k = numtaps + 1 - j;
+    h[k] = neg === 0 ? h[j] : -h[j];
+  }
+  if (neg === 1 && nodd === 1) h[nz] = 0;
+  return h.slice(1);
+}
+
 // src/signal/fir-smoothing.ts
 var arr9 = (x) => Array.isArray(x) ? x : Array.from(x);
 function solveLinearSystem(Ain, bin2) {
@@ -49179,11 +49860,6 @@ function solveWeightedBasis(basis, pts, weights) {
   }
   return solveLinearSystem(A, b);
 }
-function evalBasis(basis, coeffs, w) {
-  let s = 0;
-  for (let k = 0; k < basis.nBasis; k++) s += coeffs[k] * basis.phi(k, w);
-  return s;
-}
 function firls(numtaps, bands, desired) {
   validateBandsDesired(numtaps, bands, desired);
   const basis = makeBasis(numtaps);
@@ -49195,22 +49871,8 @@ function firls(numtaps, bands, desired) {
   );
   return coeffsToTaps(numtaps, basis, coeffs);
 }
-function remez(numtaps, bands, desired) {
-  validateBandsDesired(numtaps, bands, desired);
-  const basis = makeBasis(numtaps);
-  const pts = sampleBands(bands, desired);
-  let weights = pts.map(() => 1);
-  let coeffs = solveWeightedBasis(basis, pts, weights);
-  const ITERATIONS = 15;
-  for (let iter = 1; iter < ITERATIONS; iter++) {
-    const errs = pts.map((pt) => Math.abs(evalBasis(basis, coeffs, pt.w) - pt.D));
-    const maxErr = Math.max(...errs, 1e-12);
-    const reweighted = weights.map((wgt, i) => wgt * Math.pow(errs[i] / maxErr + 1e-6, 0.5));
-    const meanW = reweighted.reduce((a, b) => a + b, 0) / reweighted.length;
-    weights = reweighted.map((wgt) => wgt / meanW);
-    coeffs = solveWeightedBasis(basis, pts, weights);
-  }
-  return coeffsToTaps(numtaps, basis, coeffs);
+function remez(numtaps, bands, desired, weight, type = "bandpass") {
+  return remezExchange(numtaps, bands, desired, weight, type);
 }
 
 // src/signal/wavelets.ts

package/dist/numeric/control-equations.d.ts

@@ -11,15 +11,27 @@
  *   - {@link care} — continuous algebraic Riccati equation
  *     `AᵀX + XA − X B R⁻¹ Bᵀ X + Q = 0`, solved via the **matrix sign
  *     function** applied to the Hamiltonian `H = [[A, −BR⁻¹Bᵀ], [−Q, −Aᵀ]]`.
- *     This avoids both complex eigenvectors (this codebase's `eig` only
- *     returns real eigenvector columns — see `matrix-functions.ts` — so a
- *     Hamiltonian-eigenvector approach isn't usable here) and the need for a
- *     stabilizing initial gain that a Kleinman/Newton iteration would
- *     require (the natural starting point `X₀ = 0` is not stabilizing for
- *     e.g. the classic double-integrator `A = [[0,1],[0,0]]`, which is only
+ *     This needs no stabilizing initial gain that a Kleinman/Newton iteration
+ *     would require (the natural starting point `X₀ = 0` is not stabilizing
+ *     for e.g. the classic double-integrator `A = [[0,1],[0,0]]`, which is only
  *     marginally stable). The sign-function Newton iteration converges
  *     unconditionally (given the standard stabilizability/detectability
  *     assumptions) directly from `H` itself.
+ *
+ *     The classical alternative is the **Hamiltonian-eigenvector** construction
+ *     (eigendecompose `H`, take the stable invariant subspace `U=[U₁;U₂]`,
+ *     `X = U₂U₁⁻¹`). Since `matrix`'s `eig` gained complex eigenvectors
+ *     (`vectorsIm`, 2026-07-16) that path IS now implementable — but it was
+ *     prototyped and **measured against the sign function** (2026-07-18) and
+ *     found strictly **less** accurate on every corpus case, with the gap
+ *     widening on ill-conditioned Hamiltonians (near-uncontrollable /
+ *     lightly-damped: Riccati residual ~1e-10..1e-13 for the eigenvector path
+ *     vs ~1e-13..1e-14 for the sign function). Eigenvector-basis subspace
+ *     extraction is the numerically inferior classical method (ill-conditioned
+ *     eigenvectors near defectiveness), which is why LAPACK/scipy use ordered
+ *     Schur, not eigenvectors. The self-contained sign-function method is
+ *     therefore retained. (`funm`/`cosm`/`sinm` likewise use eigen*values*,
+ *     not eigenvectors — see `matrix-functions.ts`.)
  *   - {@link dare} — discrete algebraic Riccati equation
  *     `AᵀXA − X − AᵀXB(R + BᵀXB)⁻¹BᵀXA + Q = 0`, solved via the
  *     **structure-preserving doubling algorithm** (SDA), the discrete-time

package/dist/numeric/control-equations.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"control-equations.d.ts","sourceRoot":"","sources":["../../src/numeric/control-equations.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AAuJH;;;;;;;;;;GAUG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CAgB9D;AASD;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CA6C3F;AAMD;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CA+B3F"}
+{"version":3,"file":"control-equations.d.ts","sourceRoot":"","sources":["../../src/numeric/control-equations.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA8CG;AAuJH;;;;;;;;;;GAUG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CAgB9D;AASD;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CA6C3F;AAMD;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,EAAE,CAAC,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,EAAE,CA+B3F"}

package/dist/signal/fir-smoothing.d.ts

@@ -1,3 +1,4 @@
+import { type RemezType } from './remez-exchange.js';
 type Vec = readonly number[] | Float64Array;
 /** Polynomial/deconvolution division result: `signal = conv(divisor, quotient) + remainder`. */
 export interface DeconvolveResult {
@@ -49,14 +50,17 @@ export declare function wiener(x: Vec, mysize?: number): number[];
  */
 export declare function firls(numtaps: number, bands: readonly number[], desired: readonly number[]): number[];
 /**
- * Parks-McClellan-style equiripple FIR design. **This is NOT the exact
- * Remez-exchange algorithm** — it is Lawson's algorithm, an iteratively
- * reweighted least squares (IRLS) scheme that approximates the minimax
- * (Chebyshev) solution by repeatedly re-solving the `firls` normal equations
- * with per-sample weights pushed up wherever the previous iteration's error
- * was largest, converging toward (but not guaranteed to reach) a true
- * equiripple response. Documented as approximate per the Task 3 spec.
+ * Optimal equiripple FIR design (`scipy.signal.remez`) via the exact
+ * Parks-McClellan / Remez exchange algorithm (delegates to `remezExchange`).
+ *
+ * **Convention note:** unlike `firls` above (which follows `scipy.signal.firls`
+ * with `1 = Nyquist` and one `desired` value per band *edge*), `remez` follows
+ * `scipy.signal.remez` exactly: band edges are normalized to `[0, 0.5]`
+ * (`fs = 1`, `0.5 = Nyquist`), and `desired`/`weight` carry one value per band
+ * (length `bands.length / 2`). `type` is `'bandpass'` (symmetric, the default),
+ * `'differentiator'`, or `'hilbert'` (antisymmetric). Coefficients match
+ * `scipy.signal.remez` to machine precision.
  */
-export declare function remez(numtaps: number, bands: readonly number[], desired: readonly number[]): number[];
+export declare function remez(numtaps: number, bands: readonly number[], desired: readonly number[], weight?: readonly number[], type?: RemezType): number[];
 export {};
 //# sourceMappingURL=fir-smoothing.d.ts.map

package/dist/signal/fir-smoothing.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"fir-smoothing.d.ts","sourceRoot":"","sources":["../../src/signal/fir-smoothing.ts"],"names":[],"mappings":"AAUA,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAG5C,gGAAgG;AAChG,MAAM,WAAW,gBAAgB;IAC/B,QAAQ,EAAE,MAAM,EAAE,CAAC;IACnB,SAAS,EAAE,MAAM,EAAE,CAAC;CACrB;AAsCD;;;;;;GAMG;AACH,wBAAgB,cAAc,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,MAAM,EAAE,CAa5F;AA+BD;;;;;;;;GAQG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,YAAY,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,MAAM,EAAE,CA4BhF;AAID;;;;;GAKG;AACH,wBAAgB,UAAU,CAAC,MAAM,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,GAAG,gBAAgB,CAiBtE;AAID;;;;;;GAMG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,MAAM,SAAI,GAAG,MAAM,EAAE,CAsBnD;AAoHD;;;;;;;;;GASG;AACH,wBAAgB,KAAK,CACnB,OAAO,EAAE,MAAM,EACf,KAAK,EAAE,SAAS,MAAM,EAAE,EACxB,OAAO,EAAE,SAAS,MAAM,EAAE,GACzB,MAAM,EAAE,CAUV;AAED;;;;;;;;GAQG;AACH,wBAAgB,KAAK,CACnB,OAAO,EAAE,MAAM,EACf,KAAK,EAAE,SAAS,MAAM,EAAE,EACxB,OAAO,EAAE,SAAS,MAAM,EAAE,GACzB,MAAM,EAAE,CAkBV"}
+{"version":3,"file":"fir-smoothing.d.ts","sourceRoot":"","sources":["../../src/signal/fir-smoothing.ts"],"names":[],"mappings":"AASA,OAAO,EAAiB,KAAK,SAAS,EAAE,MAAM,qBAAqB,CAAC;AAEpE,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAG5C,gGAAgG;AAChG,MAAM,WAAW,gBAAgB;IAC/B,QAAQ,EAAE,MAAM,EAAE,CAAC;IACnB,SAAS,EAAE,MAAM,EAAE,CAAC;CACrB;AAsCD;;;;;;GAMG;AACH,wBAAgB,cAAc,CAAC,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,MAAM,EAAE,CAa5F;AA+BD;;;;;;;;GAQG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,YAAY,EAAE,MAAM,EAAE,SAAS,EAAE,MAAM,GAAG,MAAM,EAAE,CA4BhF;AAID;;;;;GAKG;AACH,wBAAgB,UAAU,CAAC,MAAM,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,GAAG,gBAAgB,CAiBtE;AAID;;;;;;GAMG;AACH,wBAAgB,MAAM,CAAC,CAAC,EAAE,GAAG,EAAE,MAAM,SAAI,GAAG,MAAM,EAAE,CAsBnD;AA8GD;;;;;;;;;GASG;AACH,wBAAgB,KAAK,CACnB,OAAO,EAAE,MAAM,EACf,KAAK,EAAE,SAAS,MAAM,EAAE,EACxB,OAAO,EAAE,SAAS,MAAM,EAAE,GACzB,MAAM,EAAE,CAUV;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,KAAK,CACnB,OAAO,EAAE,MAAM,EACf,KAAK,EAAE,SAAS,MAAM,EAAE,EACxB,OAAO,EAAE,SAAS,MAAM,EAAE,EAC1B,MAAM,CAAC,EAAE,SAAS,MAAM,EAAE,EAC1B,IAAI,GAAE,SAAsB,GAC3B,MAAM,EAAE,CAEV"}

package/dist/signal/iir-design.d.ts

@@ -63,17 +63,26 @@ export declare function bilinear(bIn: Vec, aIn: Vec, fs: number): {
     a: number[];
 };
 /**
- * Minimum Butterworth filter order (and the corresponding `-3dB`-ish natural
- * frequency `Wn`) meeting a passband/stopband spec, for a digital lowpass
- * (`wp < ws`) or highpass (`wp > ws`) filter. `wp`/`ws` are normalized to
- * Nyquist (0..1); `gpass`/`gstop` are in dB. Matches
- * `scipy.signal.buttord(wp, ws, gpass, gstop)` (scalar/digital case only —
- * the bandpass/bandstop array form is not implemented).
+ * Minimum Butterworth filter order and the corresponding natural frequency
+ * `Wn` (the "3 dB" frequency) meeting a passband/stopband spec, for a digital
+ * lowpass, highpass, bandpass, or bandstop filter. Frequencies are normalized
+ * to Nyquist (`0..1`); `gpass`/`gstop` are in dB. For lowpass/highpass, `wp`
+ * and `ws` are scalars and `Wn` is a scalar; for bandpass/bandstop, `wp` and
+ * `ws` are `[low, high]` pairs and `Wn` is the `[low, high]` natural-frequency
+ * array. Matches `scipy.signal.buttord(wp, ws, gpass, gstop)`.
+ *
+ * Filter type follows scipy: scalar `wp < ws` → lowpass, `wp > ws` → highpass;
+ * array with `wp[0] > ws[0]` (passband inside stopband) → bandpass, `wp[0] <
+ * ws[0]` (stopband inside passband) → bandstop.
  */
 export declare function buttord(wp: number, ws: number, gpass: number, gstop: number): {
     N: number;
     Wn: number;
 };
+export declare function buttord(wp: readonly [number, number], ws: readonly [number, number], gpass: number, gstop: number): {
+    N: number;
+    Wn: [number, number];
+};
 type ZpkRoot = number | Complex;
 /**
  * Group zeros/poles/gain into cascaded second-order sections:

package/dist/signal/iir-design.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"iir-design.d.ts","sourceRoot":"","sources":["../../src/signal/iir-design.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;GAiBG;AACH,OAAO,EAAE,OAAO,EAAE,MAAM,4BAA4B,CAAC;AACrD,OAAO,EAOL,KAAK,WAAW,EACjB,MAAM,2BAA2B,CAAC;AAKnC,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAoB5C;;;;GAIG;AACH,wBAAgB,MAAM,CACpB,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GAAG,SAAS,MAAM,EAAE,EAC9B,KAAK,GAAE,WAAmB,GACzB;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CAI9B;AAgCD;;;;;GAKG;AACH,wBAAgB,MAAM,CACpB,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GAAG,SAAS,MAAM,EAAE,EAC9B,KAAK,GAAE,WAAmB,GACzB;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CAI9B;AAkJD;;;;;;;;GAQG;AACH,wBAAgB,KAAK,CACnB,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GAAG,SAAS,MAAM,EAAE,EAC9B,KAAK,GAAE,WAAmB,GACzB;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CAK9B;AA4BD;;;;;;GAMG;AACH,wBAAgB,QAAQ,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,MAAM,GAAG;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CA8BrF;AAID;;;;;;;GAOG;AACH,wBAAgB,OAAO,CACrB,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,KAAK,EAAE,MAAM,EACb,KAAK,EAAE,MAAM,GACZ;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,EAAE,EAAE,MAAM,CAAA;CAAE,CAc3B;AAID,KAAK,OAAO,GAAG,MAAM,GAAG,OAAO,CAAC;AAoDhC;;;;;;;GAOG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,SAAS,OAAO,EAAE,EAAE,CAAC,EAAE,SAAS,OAAO,EAAE,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,EAAE,EAAE,CAmB3F;AAED;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,GAAG,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CAS7E"}
+{"version":3,"file":"iir-design.d.ts","sourceRoot":"","sources":["../../src/signal/iir-design.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;GAiBG;AACH,OAAO,EAAE,OAAO,EAAE,MAAM,4BAA4B,CAAC;AACrD,OAAO,EAOL,KAAK,WAAW,EACjB,MAAM,2BAA2B,CAAC;AAKnC,KAAK,GAAG,GAAG,SAAS,MAAM,EAAE,GAAG,YAAY,CAAC;AAoB5C;;;;GAIG;AACH,wBAAgB,MAAM,CACpB,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GAAG,SAAS,MAAM,EAAE,EAC9B,KAAK,GAAE,WAAmB,GACzB;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CAI9B;AAgCD;;;;;GAKG;AACH,wBAAgB,MAAM,CACpB,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GAAG,SAAS,MAAM,EAAE,EAC9B,KAAK,GAAE,WAAmB,GACzB;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CAI9B;AAkJD;;;;;;;;GAQG;AACH,wBAAgB,KAAK,CACnB,CAAC,EAAE,MAAM,EACT,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,GAAG,SAAS,MAAM,EAAE,EAC9B,KAAK,GAAE,WAAmB,GACzB;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CAK9B;AA4BD;;;;;;GAMG;AACH,wBAAgB,QAAQ,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,MAAM,GAAG;IAAE,CAAC,EAAE,MAAM,EAAE,CAAC;IAAC,CAAC,EAAE,MAAM,EAAE,CAAA;CAAE,CA8BrF;AAyHD;;;;;;;;;;;;GAYG;AACH,wBAAgB,OAAO,CACrB,EAAE,EAAE,MAAM,EACV,EAAE,EAAE,MAAM,EACV,KAAK,EAAE,MAAM,EACb,KAAK,EAAE,MAAM,GACZ;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,EAAE,EAAE,MAAM,CAAA;CAAE,CAAC;AAC7B,wBAAgB,OAAO,CACrB,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,EAC7B,EAAE,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,EAC7B,KAAK,EAAE,MAAM,EACb,KAAK,EAAE,MAAM,GACZ;IAAE,CAAC,EAAE,MAAM,CAAC;IAAC,EAAE,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,CAAA;CAAE,CAAC;AAkFvC,KAAK,OAAO,GAAG,MAAM,GAAG,OAAO,CAAC;AAoDhC;;;;;;;GAOG;AACH,wBAAgB,OAAO,CAAC,CAAC,EAAE,SAAS,OAAO,EAAE,EAAE,CAAC,EAAE,SAAS,OAAO,EAAE,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,EAAE,EAAE,CAmB3F;AAED;;;;GAIG;AACH,wBAAgB,OAAO,CAAC,GAAG,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,EAAE,GAAG,GAAG,MAAM,EAAE,CAS7E"}

package/dist/signal/remez-exchange.d.ts

@@ -0,0 +1,28 @@
+/**
+ * Exact Parks-McClellan (Remez exchange) optimal equiripple FIR design.
+ *
+ * This is a faithful TypeScript port of the classic McClellan-Parks-Rabiner
+ * FORTRAN program (via Erik Kvaleberg's `remez.c`), the same code SciPy wraps
+ * in `scipy.signal.remez` (`_sigtoolsmodule.cc`: `pre_remez` + `remez`). It
+ * replaces the earlier Lawson-IRLS approximation. Conventions match SciPy
+ * exactly: band edges are normalized to `[0, 0.5]` (`fs = 1`, `0.5 = Nyquist`),
+ * `desired`/`weight` carry one value per band (length `bands.length / 2`), and
+ * `type` is `'bandpass'` (symmetric), `'differentiator'`, or `'hilbert'`
+ * (antisymmetric). Coefficients are pinned against `scipy.signal.remez` in
+ * `functions/tests/remez-pm.test.ts`.
+ *
+ * The port preserves the original's 1-based array indexing and control flow
+ * (the extremal-search `goto` block is expressed as a small state machine) so
+ * that it reproduces SciPy's discrete-grid extremal selection bit-for-bit.
+ */
+export type RemezType = 'bandpass' | 'differentiator' | 'hilbert';
+/**
+ * Optimal equiripple FIR design (`scipy.signal.remez`) via the exact
+ * Parks-McClellan / Remez exchange algorithm. `bands` is a flat list of band
+ * edges normalized to `[0, 0.5]` (`fs = 1`, `0.5 = Nyquist`); `desired` and
+ * `weight` carry one value per band (`bands.length / 2` entries). `type`
+ * selects the linear-phase family (`'bandpass'` symmetric; `'differentiator'`
+ * and `'hilbert'` antisymmetric). Returns `numtaps` FIR coefficients.
+ */
+export declare function remezExchange(numtaps: number, bands: readonly number[], desired: readonly number[], weight?: readonly number[], type?: RemezType, maxiter?: number, gridDensity?: number): number[];
+//# sourceMappingURL=remez-exchange.d.ts.map

package/dist/signal/remez-exchange.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"remez-exchange.d.ts","sourceRoot":"","sources":["../../src/signal/remez-exchange.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAKH,MAAM,MAAM,SAAS,GAAG,UAAU,GAAG,gBAAgB,GAAG,SAAS,CAAC;AAgelE;;;;;;;GAOG;AACH,wBAAgB,aAAa,CAC3B,OAAO,EAAE,MAAM,EACf,KAAK,EAAE,SAAS,MAAM,EAAE,EACxB,OAAO,EAAE,SAAS,MAAM,EAAE,EAC1B,MAAM,CAAC,EAAE,SAAS,MAAM,EAAE,EAC1B,IAAI,GAAE,SAAsB,EAC5B,OAAO,SAAK,EACZ,WAAW,SAAK,GACf,MAAM,EAAE,CA4IV"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@danielsimonjr/mathts-functions",
-  "version": "0.43.3",
+  "version": "0.44.0",
   "description": "Mathematical functions for MathTS - arithmetic, algebra, trigonometry, statistics, and more",
   "author": "Daniel Simon Jr.",
   "license": "MIT",