High order fast algorithm for the Caputo fractional derivative

Abstract

In the paper, we present a high order fast algorithm with almost optimum memory for the Caputo fractional derivative, which can be expressed as a convolution of u'(t) with the kernel (tn-t)-α. In the fast algorithm, the interval [0,tn-1] is split into nonuniform subintervals. The number of the subintervals is in the order of n at the n-th time step. The fractional kernel function is approximated by a polynomial function of K-th degree with a uniform absolute error on each subinterval. We save K+1 integrals on each subinterval, which can be written as a convolution of u'(t) with a polynomial base function. As compared with the direct method, the proposed fast algorithm reduces the storage requirement and computational cost from O(n) to O((K+1) n) at the n-th time step. We prove that the convergence rate of the fast algorithm is the same as the direct method even a high order direct method is considered. The convergence rate and efficiency of the fast algorithm are illustrated via several numerical examples.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…