Levin methods for highly oscillatory integrals with singularities

12/20/2019
by   Yinkun Wang, et al.
0

In this paper, new Levin methods are presented for calculating oscillatory integrals with algebraic and/or logarithmic singularities. To avoid singularity, the technique of singularity separation is applied and then the singular ODE occurring in classic Levin methods is converted into two kinds of non-singular ODEs. The solutions of one can be obtained explicitly, while those of the other can be solved efficiently by collocation methods. The proposed methods can attach arbitrarily high asymptotic orders and also enjoy superalgebraic convergence with respect to the number of collocation points. Several numerical experiments are presented to validate the efficiency of the proposed methods.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
05/12/2021

Variable stepsize SDIMSIMs for ordinary differential equations

Second derivative general linear methods (SGLMs) have been already imple...
research
08/02/2023

Subgradient Langevin Methods for Sampling from Non-smooth Potentials

This paper is concerned with sampling from probability distributions π o...
research
08/02/2022

Solving singular generalized eigenvalue problems. Part II: projection and augmentation

Generalized eigenvalue problems involving a singular pencil may be very ...
research
08/20/2022

A Modified Trapezoidal Rule for a Class of Weakly Singular Integrals in n Dimensions

In this paper we propose and analyze a general arbitrarily high-order mo...
research
10/04/2022

Efficient implicit solvers for models of neuronal networks

We introduce economical versions of standard implicit ODE solvers that a...
research
07/10/2016

A New Bengali Readability Score

In this paper we have proposed methods to analyze the readability of Ben...
research
10/18/2021

Model Order Estimation for A Sum of Complex Exponentials

In this paper, we present a new method for estimating the number of term...

Please sign up or login with your details

Forgot password? Click here to reset