Related papers: A Functional Package for Automatic Solution of Ord…
We announce some Python classes for numerical solution of partial differential equations, or boundary value problems of ordinary differential equations. These classes are built on routines in \texttt{numpy} and \texttt{scipy.sparse.linalg}…
We develop a spectral method to solve the heat equation in a closed cylinder, achieving a quasi-optimal $\mathcal{O}(N\log N)$ complexity and high-order, spectral accuracy. The algorithm relies on a Chebyshev--Chebyshev--Fourier (CCF)…
High-quality ordinary differential equation (ODE) solver libraries have a long history, going back to the 1970s. Over the past several years we have implemented, on top of the PETSc linear and nonlinear solver package, a new…
The Obreshkov method is a single-step multi-derivative method used in the numerical solution of differential equations and has been used in recent years in efficient circuit simulation. It has been shown that it can be made of arbitrary…
We propose a spectral solver for the Poisson equation on a square domain, achieving optimal complexity through the ultraspherical spectral method and the alternating direction implicit (ADI) method. Compared with the state-of-the-art…
For a function that is analytic on and around an interval, Chebyshev polynomial interpolation provides spectral convergence. However, if the function has a singularity close to the interval, the rate of convergence is near one. In these…
The advection-diffusion and wave equations are the fundamental equations governing any physical law and therefore arise in many areas of physics and astrophysics. For complex problems and geometries, only numerical simulations can give…
Moved by the need for rigorous and reliable numerical tools for the analysis of peridynamic materials, the authors propose a model able to capture the dispersive features of nonlocal soliton-like solutions obtained by a peridynamic…
This paper describes the algorithms, features and implementation of PyDEC, a Python library for computations related to the discretization of exterior calculus. PyDEC facilitates inquiry into both physical problems on manifolds as well as…
We present a Chebyshev collocation method for linear ODE and DDE problems. We first give a posteriori estimates for the accuracy of the approximate solution of a scalar ODE initial value problem. Examples of the success of the estimate are…
A spectral element method (SEM) is developed to solve polarized radiative transfer in multidimensional participating medium. The angular discretization is based on the discrete-ordinates approach, and the spatial discretization is conducted…
The numerical solution methods for partial differential equation (PDE) solution allow obtaining a discrete field that converges towards the solution if the method is applied to the correct problem. Nevertheless, the numerical methods…
There exist several methods for computing exact solutions of algebraic differential equations. Most of the methods, however, do not ensure existence and uniqueness of the solutions and might fail after several steps, or are restricted to…
In this study linear and nonlinear higher order singularly perturbed problems are examined by a numerical approach, the differential quadrature method. Here, the main idea is using Chebyshev polynomials to acquire the weighting coefficient…
The efficient solution of discretisations of coupled systems of partial differential equations (PDEs) is at the core of much of numerical simulation. Significant effort has been expended on scalable algorithms to precondition Krylov…
The purpose of this study is to utilize the Chebyshev spectral method neural network(CSNN) model to solve differential equations. This approach employs a single-layer neural network wherein Chebyshev spectral methods are used to construct…
This paper presents a numerical method for variable coefficient elliptic PDEs with mostly smooth solutions on two dimensional domains. The PDE is discretized via a multi-domain spectral collocation method of high local order (order 30 and…
The aim of the present work is to introduce a method based on Chebyshev polynomials for the numerical solution of a system of Cauchy type singular integral equations of the first kind on a finite segment. Moreover, an estimation error is…
Operator-splitting methods are widespread in the numerical solution of differential equations, especially the initial-value problems in ordinary differential equations that arise from a method-of-lines discretization of partial differential…
A new approach for integration of the initial value problem for ordinary differential equations is suggested. The algorithm is based on approximation of the solution by a system of functions that contains orthogonal exponential polynomials.