Spectral-Analysis-UPB / PyZEAL

Project dealing with the numerical calculation of zeros, poles and residues of holomorphic and meromorphic functions. It aspires to be a PYthon ZEAL (ZEros of AnaLytic functions, a Fortran90 package) successor.
https://pyzeal.readthedocs.io/en/latest/
GNU General Public License v3.0
4 stars 1 forks source link

Improve efficiency of quadrature based estimation #34

Open PhilippSchuette opened 1 year ago

PhilippSchuette commented 1 year ago

With the current status of argument estimation there is exactly one estimator capable of higher moments: QuadratureEstimator. Right now, this estimator re-calculates function values for every concrete moment order requested. This is quite inefficient as function values f'(z) / f(z) could quite easily be reused between different moment orders.

Implementing this reuse is mostly a matter of storing previously calculated function arrays in internal caches without too much additional overhead in the case where only the zeroth moment is required by some concrete FinderAlgorithm.

Note that this is quite distinct from the internal caching used in SummationEstimator because we retrieve function values on the exact complex line segment currently considered!