Related papers: On Fast Implementation of Higher Order Hermite-Fej…
This work proposes a higher-order iterative framework for solving matrix equations, inspired by the structure and functionality of neural networks. A modification of the classical Jacobi iterative method is introduced to compute…
We consider the problem of computing a Wasserstein barycenter for a set of discrete probability distributions with finite supports, which finds many applications in areas such as statistics, machine learning and image processing. When the…
A practical and efficient scheme for the higher order integration of the Landau-Lifschitz-Gilbert (LLG) equation is presented. The method is based on extrapolation of the two-step explicit midpoint rule and incorporates adaptive time step…
In this paper, superconvergence points are located for the approximation of the Riesz derivative of order $\alpha$ using classical Lobatto-type polynomials when $\alpha \in (0,1)$ and generalized Jacobi functions (GJF) for arbitrary $\alpha…
New implicit and implicit-explicit time-stepping methods for the wave equation in second-order form are described with application to two and three-dimensional problems discretized on overset grids. The implicit schemes are single step,…
We present a new formula for the Hermite multivariate interpolation problem in the framework of the Chung--Yao approach. By using the respective univariate interpolation formula, we obtain a direct and explicit solution to the classical…
As is well-known, the advantage of the high-order compact difference scheme (H-OCD) is unconditionally stable and convergent with the order $O(\tau^2+h^4)$ under the maximum norm. In this article, a new numerical gradient scheme based on…
In this paper, a class of high order numerical schemes is proposed for solving Hamilton-Jacobi (H-J) equations. This work is regarded as an extension of our previous work for nonlinear degenerate parabolic equations, see Christlieb et al.…
We determine properties of the lattice Boltzmann method for semiclassical fluids, which is based on the Boltzmann equation and the equilibrium distribution function is given either by the Bose-Einstein or the Fermi-Dirac ones. New…
Fast Multipole Methods (FMMs) based on the oscillatory Helmholtz kernel can reduce the cost of solving N-body problems arising from Boundary Integral Equations (BIEs) in acoustic or electromagnetics. However, their cost strongly increases…
The interpolation on Grassmann manifolds in the framework of parametric evolution partial differential equations is presented. Interpolation points on the Grassmann manifold are the subspaces spanned by the POD bases of the available…
We present the detailed formalism of the extremely correlated Fermi liquid theory, developed for treating the physics of the t-J model. We start from the exact Schwinger equation of motion for the Greens function for projected electrons,…
We consider regular polynomial interpolation algorithms on recursively defined sets of interpolation points which approximate global solutions of arbitrary well-posed systems of linear partial differential equations. Convergence of the…
We use the well-known observation that the solutions of Jacobi's differential equation can be represented via non-oscillatory phase and amplitude functions to develop a fast algorithm for computing multi-dimensional Jacobi polynomial…
We consider the problem of computing univariate polynomial matrices over a field that represent minimal solution bases for a general interpolation problem, some forms of which are the vector M-Pad\'e approximation problem in [Van Barel and…
An advantageous feature of piecewise constant policy timestepping for Hamilton-Jacobi-Bellman (HJB) equations is that different linear approximation schemes, and indeed different meshes, can be used for the resulting linear equations for…
Range functions are a fundamental tool for certified computations in geometric modeling, computer graphics, and robotics, but traditional range functions have only quadratic convergence order ($m=2$). For ``superior'' convergence order…
We develop a method for the rigorous estimation of Hausdorff dimensions of limit sets produced by continued fraction iterated function systems. Our method is based on the approximation of a Perron-Frobenius operator using the finite element…
We present a method to obtain weight functions associated with linear and quadratic lattices that yield discrete orthogonality with respect to a quasi-definite moment functional of the orthogonal polynomial sequences in the Askey scheme,…
This paper presents an efficient method to compute the maximum likelihood (ML) estimation of the parameters of a complex 2-D sinusoidal, with the complexity order of the FFT. The method is based on an accurate barycentric formula for…