On error estimates for Galerkin spectral discretizations of parabolic problems with nonsmooth initial data

We analyze the Legendre and Chebyshev spectral Galerkin semidiscretizations of a one dimensional homogeneous parabolic problem with nonconstant coefficients. We present error estimates for both smooth and nonsmooth data. In the Chebyshev case a limit in the order of approximation is established. On the contrary, in the Legendre case we find an arbitrary high order of convegence.

     

The numerical solution of integral-algebraic equations of index 1 by polynomial spline collocation methods

In this paper, we study polynomial spline collocation methods applied to a particular class of integral-algebraic equations of Volterra type. We analyse mixed systems of second and first kind integral equations. Global convergence and local superconvergence results are established.

     

On the resolution of relative Thue equations

An efficient algorithm is given for the resolution of relative Thue equations. The essential improvement is the application of an appropriate version of Wildanger's enumeration procedure based on the ellipsoid method of Fincke and Pohst.

Recently relative Thue equations have gained an important application, e.g., in computing power integral bases in algebraic number fields. The presented methods can surely be used to speed up those algorithms.

The method is illustrated by numerical examples.

     

Analysis of third-order methods for secular equations

Third-order numerical methods are analyzed for secular equations. These equations arise in several matrix problems and numerical linear algebra applications. A closer look at an existing method shows that it can be considered as a classical method for an equivalent problem. This not only leads to other third-order methods, it also provides the means for a unifying convergence analysis of these methods and for their comparisons. Finally, we consider approximated versions of the aforementioned methods.

     

Spectral discretization of Darcy equations coupled with Stokes equations by vorticity–velocity–pressure formulations

We propose to make the numerical analysis of a model coupling the Darcy equations in a porous medium with the Stokes equations in the cracks. The coupling is provided by a pressure continuity on the interface. We describe a discretization by spectral element methods. We derive a priori optimal error estimates and we present some numerical experiments which confirm the results of the analysis.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1628–1651, 2017     

Jacobi rational approximation and spectral method for differential equations of degenerate type

We introduce an orthogonal system on the half line, induced by Jacobi polynomials. Some results on the Jacobi rational approximation are established, which play important roles in designing and analyzing the Jacobi rational spectral method for various differential equations, with the coefficients degenerating at certain points and growing up at infinity. The Jacobi rational spectral method is proposed for a model problem appearing frequently in finance. Its convergence is proved. Numerical results demonstrate the efficiency of this new approach.

     

Continuous collocation approximations to solutions of first kind Volterra equations

In this paper we give necessary and sufficient conditions for convergence of continuous collocation approximations of solutions of first kind Volterra integral equations. The results close some longstanding gaps in the theory of polynomial spline collocation methods for such equations. The convergence analysis is based on a Runge-Kutta or ODE approach.

     

Multi-quadratic mappings in Banach spaces

We prove the stability of multi-quadratic functional equations in Banach spaces.

     

Computing periodic solutions of linear differential-algebraic equations by waveform relaxation

We propose an algorithm, which is based on the waveform relaxation (WR) approach, to compute the periodic solutions of a linear system described by differential-algebraic equations. For this kind of two-point boundary problems, we derive an analytic expression of the spectral set for the periodic WR operator. We show that the periodic WR algorithm is convergent if the supremum value of the spectral radii for a series of matrices derived from the system is less than 1. Numerical examples, where discrete waveforms are computed with a backward-difference formula, further illustrate the correctness of the theoretical work in this paper.

     

Nonlinear wave equations and singular solutions

We construct solutions to nonlinear wave equations that are singular along a prescribed noncharacteristic hypersurface, which is the graph of a function satisfying not the Eikonal but another partial differential equation of the first order. The method of Fuchsian reduction is employed.

     

The equations determining intermediate integrals for Monge-Ampère PDE

In this note we will find the differential equations determining the intermediate integrals for Monge-Ampère equations in an arbitrary number of variables.

     

On Spectral Method for Volterra Functional Integro-Differential Equations of Neutral Type

The main purpose of this work is to provide a numerical method for the solution of Volterra functional integro-differential equations of neutral type based on a spectral approach. We analyze the convergence properties of the spectral method to approximate smooth solutions of Volterra functional integro-differential equations of neutral type. It is shown that for the neutral integro-differential equations, the spectral methods yield an exponential order of convergence.     

Asymptotic analysis for linear difference equations

We are concerned with asymptotic analysis for linear difference equations in a locally convex space. First we introduce the profile operator, which plays a central role in analyzing the asymptotic behaviors of the solutions. Then factorial asymptotic expansions for the solutions are given quite explicitly. Finally we obtain Gevrey estimates for the solutions. In a forthcoming paper we will develop the theory of cohomology groups for recurrence relations. The main results in this paper lay analytic foundations of such an algebraic theory, while they are of intrinsic interest in the theory of finite differences.

     

A minimum energy problem and Dirichlet spaces

We analyze a minimum energy problem for a discrete electrostatic model in the complex plane and discuss some applications. A natural characteristic distinguishing the state of minimum energy from other equilibrium states is established. It enables us to gain insight into the structure of positive trigonometric polynomials and Dirichlet spaces associated with finitely atomic measures. We also derive a related family of linear second order differential equations with polynomial solutions.

     

A globally convergent BFGS method for nonlinear monotone equations without any merit functions

Since 1965, there has been significant progress in the theoretical study on quasi-Newton methods for solving nonlinear equations, especially in the local convergence analysis. However, the study on global convergence of quasi-Newton methods is relatively fewer, especially for the BFGS method. To ensure global convergence, some merit function such as the squared norm merit function is typically used. In this paper, we propose an algorithm for solving nonlinear monotone equations, which combines the BFGS method and the hyperplane projection method. We also prove that the proposed BFGS method converges globally if the equation is monotone and Lipschitz continuous without differentiability requirement on the equation, which makes it possible to solve some nonsmooth equations. An attractive property of the proposed method is that its global convergence is independent of any merit function.We also report some numerical results to show efficiency of the proposed method.

     

Spectral methods for pantograph-type differential and integral equations with multiple delays

We analyze the convergence properties of the spectral method when used to approximate smooth solutions of delay differential or integral equations with two or more vanishing delays. It is shown that for the pantograph-type functional equations the spectral methods yield the familiar exponential order of convergence. Various numerical examples are used to illustrate these results.      

Exponential averaging for Hamiltonian evolution equations

We derive estimates on the magnitude of non-adiabatic interaction between a Hamiltonian partial differential equation and a high-frequency nonlinear oscillator. Assuming spatial analyticity of the initial conditions, we show that the dynamics can be transformed to the uncoupled dynamics of an infinite-dimensional Hamiltonian system and an anharmonic oscillator, up to coupling terms which are exponentially small in a certain power of the frequency of the oscillator. The result is derived from an abstract averaging theorem for infinite-dimensional analytic evolution equations in Gevrey spaces. Refining upon a similar result by Neishtadt for analytic ordinary differential equations, the temporal estimate crucially depends on the spatial regularity of the initial condition. The result shows to what extent the strong resonances between rapid forcing and highly oscillatory spatial modes can be suppressed by the choice of sufficiently smooth initial data. An application is provided by a system of nonlinear Schrödinger equations, coupled to a rapidly forcing single mode, representing small-scale oscillations. We provide an example showing that the estimates for partial differential equations we derive here are necessarily different from those in the context of ordinary differential equations.

     

Legendre spectral Galerkin method for second-kind Volterra integral equations

The Legendre spectral Galerkin method for the Volterra integral equations of the second kind is proposed in this paper. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors (in the L ₂ norm) will decay exponentially provided that the kernel function and the source function are sufficiently smooth. Numerical examples are given to illustrate the theoretical results.      

Evolution Galerkin methods for hyperbolic systems in two space dimensions

The subject of the paper is the analysis of three new evolution Galerkin schemes for a system of hyperbolic equations, and particularly for the wave equation system. The aim is to construct methods which take into account all of the infinitely many directions of propagation of bicharacteristics. The main idea of the evolution Galerkin methods is the following: the initial function is evolved using the characteristic cone and then projected onto a finite element space. A numerical comparison is given of the new methods with already existing methods, both those based on the use of bicharacteristics as well as commonly used finite difference and finite volume methods. We discuss the stability properties of the schemes and derive error estimates.

     

Spectral residual method without gradient information for solving large-scale nonlinear systems of equations

A fully derivative-free spectral residual method for solving large-scale nonlinear systems of equations is presented. It uses in a systematic way the residual vector as a search direction, a spectral steplength that produces a nonmonotone process and a globalization strategy that allows for this nonmonotone behavior. The global convergence analysis of the combined scheme is presented. An extensive set of numerical experiments that indicate that the new combination is competitive and frequently better than well-known Newton-Krylov methods for large-scale problems is also presented.