Remote Access Mathematics of Computation
Green Open Access

Mathematics of Computation

ISSN 1088-6842(online) ISSN 0025-5718(print)

Request Permissions   Purchase Content 
 

 

Finite difference weights, spectral differentiation, and superconvergence


Authors: Burhan Sadiq and Divakar Viswanath
Journal: Math. Comp. 83 (2014), 2403-2427
MSC (2010): Primary 65D05, 65D25
DOI: https://doi.org/10.1090/S0025-5718-2014-02798-1
Published electronically: January 6, 2014
MathSciNet review: 3223337
Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: Let $ z_{1},z_{2},\ldots ,z_{N}$ be a sequence of distinct grid points. A finite difference formula approximates the $ m$-th derivative $ f^{(m)}(0)$ as $ \sum w_{k}f\left (z_{k}\right )$, with $ w_{k}$ being the weight at $ z_{k}$. We derive an algorithm for finding the weights $ w_{k}$ which uses fewer arithmetic operations and less memory than the algorithm in current use (Fornberg, Mathematics of Computation, vol. 51 (1988), pp. 699-706). The algorithm we derive uses fewer arithmetic operations by a factor of $ (5m+5)/4$ in the limit of large $ N$. The optimized C++ implementation we describe is a hundred to five hundred times faster than MATLAB. The method of Fornberg is faster by a factor of five in MATLAB, however, and thus remains the most attractive option for MATLAB users.

The algorithm generalizes easily to the calculation of spectral differentiation matrices, or equivalently, finite difference weights at several different points with a fixed grid. Unlike the algorithm in current use for the calculation of spectral differentiation matrices, the algorithm we derive suffers from no numerical instability.

The order of accuracy of the finite difference formula for $ f^{(m)}(0)$ with grid points $ hz_{k}$, $ 1\leq k\leq N$, is typically $ \mathcal {O}\left (h^{N-m}\right )$. However, the most commonly used finite difference formulas have an order of accuracy that is higher than normal. For instance, the centered difference approximation $ \left (f(h)\!-\!2f(0)\!+\!f(-h)\right )$
$ /h^{2}$ to $ f''(0)$ has an order of accuracy equal to $ 2$ not $ 1$. Even unsymmetric finite difference formulas can exhibit such superconvergence or boosted order of accuracy, as shown by the explicit algebraic condition that we derive. If the grid points are real, we prove a basic result stating that the order of accuracy can never be boosted by more than $ 1$.


References [Enhancements On Off] (What's this?)


Similar Articles

Retrieve articles in Mathematics of Computation with MSC (2010): 65D05, 65D25

Retrieve articles in all journals with MSC (2010): 65D05, 65D25


Additional Information

Burhan Sadiq
Affiliation: Department of Mathematics, University of Michigan, Ann Arbor, Michigan 48109
Email: bsadiq@umich.edu

Divakar Viswanath
Affiliation: Department of Mathematics, University of Michigan, Ann Arbor, Michigan 48109
Email: divakar@umich.edu

DOI: https://doi.org/10.1090/S0025-5718-2014-02798-1
Received by editor(s): July 20, 2012
Received by editor(s) in revised form: December 22, 2012, and January 17, 2013
Published electronically: January 6, 2014
Additional Notes: The authors were supported by NSF grants DMS-0715510, DMS-1115277, and SCREMS-1026317.
Article copyright: © Copyright 2014 American Mathematical Society