vix.ing · top · new · best · stats · spec

A Fast Time Stepping Method for Evaluating Fractional Integrals

2010/01/01 by Jing-Rebecca Li, Jing‐Rebecca Li · 9 citations
Mathematics · #Fractional Differential Equations Solutions #Mathematical functions and polynomials #Numerical methods for differential equations

paper · doi:10.1137/080736533

openalex publication_date 2010/01/01 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/30

Abstract

We evaluate the fractional integral Iα[f](t)=(1)/(Γ(α))∫0t(t-τ)α-1 f(τ) dτ, 0<α<1, at time steps t=Δ t,2Δ t,…,NΔ t by making use of the integral representation of the convolution kernel tα-1=(1)/(Γ(1-α))∫0e-ξ t ξ dξ. We construct an efficient Q-point quadrature of this integral representation and use it as a part of a fast time stepping method. The new method has algorithmic complexity O(NQ) and storage requirement O(Q). The number of quadrature nodes Q is independent of N and grows like O((-logε-logΔ t)2), where ε is the quadrature error tolerance and Δ t is the size of the time step. The (possible) singularity of f near τ=0 is taken into account. This new method is particularly well-suited for long time simulations.

Citations

Cited by

Related