DeepAI AI Chat
Log In Sign Up

Efficient computation of the sinc matrix function for the integration of second-order differential equations

by   Lidia Aceto, et al.

This work deals with the numerical solution of systems of oscillatory second-order differential equations which often arise from the semi-discretization in space of partial differential equations. Since these differential equations exhibit pronounced or highly) oscillatory behavior, standard numerical methods are known to perform poorly. Our approach consists in directly discretizing the problem by means of Gautschi-type integrators based on sinc matrix functions. The novelty contained here is that of using a suitable rational approximation formula for the sinc matrix function to apply a rational Krylov-like approximation method with suitable choices of poles. In particular, we discuss the application of the whole strategy to a finite element discretization of the wave equation.


page 15

page 25


Exponential integrators for second-order in time partial differential equations

Two types of second-order in time partial differential equations (PDEs),...

Approximate Solutions to Second-Order Parabolic Equations: evolution systems and discretization

We study the discretization of a linear evolution partial differential e...

The Wynn identity as the long sought criterion for the choice of the optimal Padé approximant

The performed numerical analysis reveals that Wynn's identity for the co...

Characteristics-based Simulink implementation of first-order quasilinear partial differential equations

The paper deals with solving first-order quasilinear partial differentia...

A Feynman-Kac based numerical method for the exit time probability of a class of transport problems

The exit time probability, which gives the likelihood that an initial co...