Related papers: Rapid numerical solutions for the Mukhanov-Sasaki …
The numerical waveforms for the extreme mass-ratio inspirals (EMRIs) require a huge amount of homogeneous solutions of the Teukolsky equation in the frequency domain. The calculation accuracy and efficiency of the homogeneous solutions are…
A fast method for the computation of layer potentials that arise in acoustic scattering is introduced. The principal idea is to split the singular kernel into a smooth and a local part. The potential due to the smooth part is computed…
This paper introduces an efficient algorithm for finding the dominant generalized eigenvectors of a pair of symmetric matrices. Combining tools from approximation theory and convex optimization, we develop a simple scalable algorithm with…
In this work we discuss the possibility to reduce the computational complexity of modal methods, i.e. methods based on eigenmodes expansion, from the third power to the second power of the number of eigenmodes. The proposed approach is…
The simulation of turbulence in the boundary region of a tokamak is crucial for understanding and optimizing the performance of fusion reactors. In this work, the use of low-rank linear algebra techniques is shown to enhance the efficiency…
We introduce a direct method allowing to solve numerically inverse type problems for linear hyperbolic equations. We first consider the reconstruction of the full solution of the wave equation posed in $\Omega\times (0,T)$ - $\Omega$ a…
We propose a novel numerical homogenization method based on the edge multiscale approach for solving indefinite time-harmonic Maxwell equations in heterogeneous media with large wavenumber. Numerical methods for these equations in…
We present a novel procedure to solve the Schr\"odinger equation, which in optics is the paraxial wave equation, with an initial multisingular vortex Gaussian beam. This initial condition has a number of singularities in a plane transversal…
In this work, a class of non-linear weakly singular fractional integro-differential equations is considered, and we first prove existence, uniqueness, and smoothness properties of the solution under certain assumptions on the given data. We…
We introduce a new class of numerical methods for solving McKean-Vlasov stochastic differential equations, which are relevant in the context of distribution-dependent or mean-field models, under super-linear growth conditions for both the…
There are two usual computational methods for linear (waves and instabilities) problem: eigenvalue (dispersion relation) solver and initial value solver. In fact, we can introduce an idea of the combination of them, i.e., we keep time…
A numerical method is proposed for computing time-periodic and relative time-periodic solutions in dissipative wave systems. In such solutions, the temporal period, and possibly other additional internal parameters such as the propagation…
The numerical algorithm of the inverse quantum scattering is developed. This algorithm is based on the Marchenko theory, and includes three steps. The first one is the algebraic Pade approximation of the unitary S-matrix, what is realized…
The matrix Numerov method provides an efficient framework for solving the time-independent Schr\"odinger equation as a matrix eigenvalue problem. However, for singular potentials such as the Coulomb interaction, the expected fourth-order…
We present a novel numerical routine (oscode) with a C++ and Python interface for the efficient solution of one-dimensional, second-order, ordinary differential equations with rapidly oscillating solutions. The method is based on a…
This paper studies the numerical methods to approximate the solutions for a sort of McKean-Vlasov neutral stochastic differential delay equations (MV-NSDDEs) that the growth of the drift coefficients is super-linear. First, We obtain that…
We propose a spectral collocation method to approximate the exact boundary control of the wave equation in a square domain. The idea is to introduce a suitable approximate control problem that we solve in the finite-dimensional space of…
Finding reliably and efficiently the spectrum of the resonant states of an optical system under varying parameters of the medium surrounding it is a technologically important task, primarily due to various sensing applications.…
Global spectral methods offer the potential to compute solutions of partial differential equations numerically to very high accuracy. In this work, we develop a novel global spectral method for linear partial differential equations on cubes…
We propose two improvements to the well-known power series method for confined one-dimensional quantum-mechanical problems. They consist of the addition of a variational step were the energy plays the role of a variational parameter. We…