A new efficient and accurate spline algorithm for the matrix exponential computation

A new efficient and accurate spline algorithm for the matrix exponential computation, Emilio Defez, Javier Ibáñez, Jorge Sastre, Jesús Peinado, Pedro Alonso. Journal of Computational and Applied Mathematics, Volume 337, pp. 354-365. August 2018. Preprint, Matlab code expmspl.m.

Abstract: In this work an accurate and efficient method based on matrix splines for computing matrix exponential is given. An algorithm and a MATLAB implementation have been developed and compared with the state-of-the-art algorithms for computing the matrix exponential. We also developed a parallel implementation for large scale problems. This implementation allowed us to get a much better performance when working with this kind of problems.

Modelling acoustics on the Poincaré half-plane

Modelling acoustics on the Poincaré half-plane. Michael M. Tung. Journal of Computational and Applied Mathematics DOI10.1016/j.cam.2017.10.037

Abstract: Novel advances in the field of metamaterial research have permitted the engineering of devices with extraordinary characteristics. Here, we explore the possibilities in transformation acoustics to implement a model for the simulation of acoustic wave propagation on the Poincaré half-plane-the simplest model possessing hyperbolic geometry and also of considerable historical interest. We start off from a variational principle on the given spacetime manifold to find the design description of the model in the laboratory. After examining some significant geometrical and physical properties of the Poincaré half-plane model, we derive a general formal solution for its acoustic wave propagation. A numerical example for the evolution of the acoustic potential on a rectangular region of the Poincaré half-plane concludes this discussion.

A new type of Hermite matrix polynomial series

A new type of Hermite matrix polynomial series (E. DEFEZ, M. TUNG), Quaestiones Mathematicae, 41(2), 2018, 205–212.

Abstract: Conventional Hermite polynomials emerge in a great diversity of applications
in mathematical physics, engineering, and related fields. However, in physical
systems with higher degrees of freedom it will be of practical interest to extend
the scalar Hermite functions to their matrix analogue. This work introduces various
new generating functions for Hermite matrix polynomials and examines existence and
convergence of their associated series expansion by using Mehlerís formula for the
general matrix case. Moreover, we derive interesting new relations for even- and
odd-power summation in the generating-function expansion containing Hermite matrix
polynomials. Some new results for the scalar case are also presented.

 

Two algorithms for computing the matrix cosine based on new Hermite approximations

 

Efficient evaluation of matrix polynomials

Efficient evaluation of matrix polynomials, Jorge Sastre, Linear Algebra Applications, Vol. 539, pp. 229-250, Feb. 2018. Preprint.

Abstract: This paper presents a new family of methods for evaluating matrix polynomials more efficiently than the state-of-the-art Paterson–Stockmeyer method. Examples of the application of the methods to the Taylor polynomial approximation of matrix functions like the matrix exponential and matrix cosine are given. Their efficiency is compared with that of the best existing evaluation schemes for general polynomial and rational approximations, and also with a recent method based on mixed rational and polynomial approximants. For many years, the Paterson–Stockmeyer method has been considered the most efficient general method for the evaluation of matrix polynomials. In this paper we show that this statement is no longer true. Moreover, for many years rational approximations have been considered more efficient than polynomial approximations, although recently it has been shown that often this is not the case in the computation of the matrix exponential and matrix cosine. In this paper we show that in fact polynomial approximations provide a higher order of approximation than the state-of-the-art computational methods for rational approximations for the same cost in terms of matrix products. For an early unpublished version of this word submitted on Feb. 18, 2016 to Appl. Math. Comput. see AMC-S-16-00951.pdf.