Abstract:The stability and error analysis of a second-order fast approximation are considered for the one-dimensional local and nonlocal diffusion equations in the unbounded spatial domain. We first use the conventional central difference scheme to discretize the local second-order spatial derivative operator and use an asymptotically compatible difference scheme to discretize the spatial nonlocal diffusion operator, and apply second-order backward differentiation formula (BDF2) to approximate the temporal derivative to achieve a fully discrete infinity system. To solve the resulting fully discrete systems, we develop a unified framework that is applicable to the discretization of both local and nonlocal problems. A key ingredient is to derive Dirichlet-to-Neumann (DtN)-type absorbing boundary conditions (ABCs). To do so, we apply the z-transform and solve an exterior problem using an iteration technique to derive a Dirichlet-to-Dirichlet (DtD)-type mapping as exact ABCs. After that, we use the Green formula to reformulate the DtD-type mapping equivalently as the DtN-type mapping. The resulting DtN-type mapping allows us to reduce the infinity discrete system into a finite discrete system in a truncated computational domain of interest, and also make it possible to present the stability and convergence analysis of the reduced problem under some open but reasonable assumptions. To efficiently implement the exact ABCs, we further develop a fast convolution algorithm based on approximation of the contour integral induced by the inverse z-transform. The stability and error analysis of the reduced finite discrete system based on the fast algorithm for exact ABCs are also established, and numerical examples are provided to demonstrate the effectiveness of our proposed approach.

Fast and Stably Recursive Algorithm for Computing Second Derivative of Associated Legendre Functions

Approximation of the derivatives beyond Taylor expansion

Numerical Differentiation for the Second Order Derivatives of Functions of Two Variables

Convergence of the family of the deformed Euler–Halley iterations under the Hölder condition of the second derivative

General explicit difference formulas for numerical differentiation

On the Numerical Calculation of the Roots of Special Functions Satisfying Second Order Ordinary Differential Equations

First Derivative Approximations and Applications

A new second-order midpoint approximation formula for Riemann-Liouville derivative: algorithm and its application

The Second-Order Numerical Approximation for a Modified Ericksen–Leslie Model

A general approach to computing derivatives for Hessian-based seismic inversion

Numerical Analysis of a Second-Order Algorithm for the Time-Dependent Natural Convection Problem

A Second-Order Legendre Spectral Scheme for the Long-Term Computation of Infinite Dimensional Dynamical Systems

Legendre Polynomials As A Recommended Basis for Numerical Differentiation in the Presence of Stochastic White Noise

Fourth‐order variants of Newton's method without second derivatives for solving non‐linear equations

Higher-Order Corrections to Optimisers based on Newton's Method

Discussion on Calculating High-Order Derivatives

A Fast Direct Legendre Spectral Tau Algorithm for One-Dimensional Stable Problem

An Optimal Fourth Order Derivative-Free Numerical Algorithm for Multiple Roots

A CLASS OF STIFFLY STABLE K-STEP SECOND DERIVATIVE METHODS OF ORDER 2K+1

Stability and Error Analysis for a Second-Order Fast Approximation of the Local and Nonlocal Diffusion Equations on the Real Line

A new approximate descent derivative-free algorithm for large-scale nonlinear symmetric equations