Reliable operations on oscillatory functions

From MaRDI portal



Abstract: Approximate p-point Leibniz derivation formulas as well as interpolatory Simpson quadrature sums adapted to oscillatory functions are discussed. Both theoretical considerations and numerical evidence concerning the dependence of the discretization errors on the frequency parameter of the oscillatory functions show that the accuracy gain of the present formulas over those based on the exponential fitting approach [L. Ixaru, "Computer Physics Communications", 105 (1997) 1--19] is overwhelming.


To differentiate or integrate a function of the form \(F(x)=f(x)w(\omega x+\delta)\) where \(w\) is \(sin\) or \(cos\) or \(sinh\) or \(cosh\), a method is proposed where only the \(f\)-part of \(F\) is approximated, while the \(w\)-factor is treated exactly. For the derivatives, the Leibniz formula is used to express \(F^{(n)}\) in terms of derivatives of \(f\) and the latter are approximated by classical finite differences. For the integrals, a second degree interpolating polynomial is used for the \(f\)-part and the resulting approximation is integrated exactly. These techniques give much better results than by the method of exponential fitting [cf. \textit{L. Ixaru}, Comput. Phys. Commun. 105, No. 1, 1-19 (1997; Zbl 0930.65150)]. Moreover the error estimates do not depend on \(\omega\).











This page was built for publication: Reliable operations on oscillatory functions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1973552)