Related papers: Interpolatory quadrature rules for oscillatory int…
We introduce quasi-Monte Carlo rules for the numerical integration of functions $f$ defined on $[0,1]^s$, $s \ge 1$, which satisfy the following properties: the Fourier-, Fourier cosine- or Walsh coefficients of $f$ are absolutely summable…
The modeling of electric machines and power transformers typically involves systems of nonlinear magnetostatics or -quasistatics, and their efficient and accurate simulation is required for the reliable design, control, and optimization of…
In this paper we propose methods for computing Fresnel integrals based on truncated trapezium rule approximations to integrals on the real line, these trapezium rules modified to take into account poles of the integrand near the real axis.…
Importance sampling (IS) and numerical integration methods are usually employed for approximating moments of complicated target distributions. In its basic procedure, the IS methodology randomly draws samples from a proposal distribution…
This work is an extension of previous work by Alazah et al. [M. Alazah, S. N. Chandler-Wilde, and S. La Porte, Numerische Mathematik, 128(4):635-661, 2014]. We split the computation of the Fresnel Integrals into 3 cases: a truncated Taylor…
In this paper, we propose a numerical method of computing an integral whose integrand is a slowly decaying oscillatory function. In the proposed method, we consider a complex analytic function in the upper-half complex plane, which is…
Quasinormal modes characterize the final stage of a black hole merger. In this regime, spacetime curvature is high, these modes can be used to probe potential corrections to general relativity. In this paper, we utilize the effective field…
We develop a stochastic approach to time-dependent DFT with optimally-tuned range-separated hybrids containing non-local exchange, for calculating optical spectra. The attractive electron-hole interaction, which leads to the formation of…
In phase space, we analytically obtain the characteristic functions (CFs) of a forced harmonic oscillator [Talkner et al., Phys. Rev. E, 75, 050102 (2007)], a time-dependent mass and frequency harmonic oscillator [Deffner and Lutz, Phys.…
The essentially non-oscillatory (ENO) method is an efficient high order numerical method for solving hyperbolic conservation laws designed to reduce the Gibbs oscillations, if existent, by adaptively choosing the local stencil for the…
In this paper, we introduce an algorithm that provides approximate solutions to semi-linear ordinary differential equations with highly oscillatory solutions, which, after an appropriate change of variables, can be rewritten as…
To investigate the possibility whether electron-phonon coupling can enhance orbital fluctuations in iron-based superconductors, we develop an ab initio method to construct the effective low-energy models including the phonon-related terms.…
The space-based detector LISA may observe gravitational waves from the early inspiral of stellar-mass black hole binaries, some of which could have significant eccentricity. Current gravitational waveform templates are only valid for small…
In this paper, He's frequency-amplitude formulation with some choice of location points that improve accuracy is applied to determine the periodic solution for the nonlinear oscillations of a punctual charge in the electric field of charged…
We introduce regular series expansion for weakly- and moderately-correlated fermionic systems, based on Fluctuating Local Field approach. The method relies on the explicit account of leading fluctuating mode(s) and is therefore suitable for…
In this paper, we consider a zero-order stochastic oracle model of estimating definite integrals. In this model, integral estimation methods may query an oracle function for a fixed number of noisy values of the integrand function and use…
Analytic modeling of gravitational waves from inspiraling eccentric binaries poses an interesting mathematical challenge. When constructing analytic waveforms in the frequency domain, one has to contend with the fact that the phase of the…
In order to approximate the Riemann--Stieltjes integral $\int_a^b {f\left( t \right)dg\left( t \right)}$ by $2$--point Gaussian quadrature rule, we introduce the quadrature rule \begin{align*} \int_{ - 1}^1 {f\left( t \right)dg\left( t…
Gaussian quadrature rules are a classical tool for the numerical approximation of integrals with smooth integrands and positive weight functions. We derive and expicitly list asymptotic expressions for the points and weights of Gaussian…
Exceptional points (EPs) are degeneracies of non-Hermitian systems, where both eigenvalues and eigenvectors coalesce. Classical and quantum systems exhibiting high-order EPs have recently been identified as fundamental building blocks for…