Numerical Analysis of Matrix Functions
Description of the Project|
Reports and Papers|
This 3 year EPSRC funded
project is concerned with functions f(A) of square matrices
A. The importance of matrix functions lies in the diverse
roles they play in the solution of problems in science and engineering.
For example, certain differential equations can be solved explicitly
in terms of matrix functions:
dy/dt = Ay
y(t) = eAt .
The corresponding inverse problems, arising in system identification,
can be solved by computing the matrix logarithm.
The aim of this project is to
further our understanding of the theory, computation and application of
matrix functions. New and improved algorithmic techniques will be developed
for the computation of both general and particular functions,
and made available in software. We also plan to investigate the problems
quantum chromodynamics (QCD) computations with the
aim of developing improved theory and methods.
In developing theory and algorithms it is natural to treat
general f as well as those particular f for which
special techniques can be applied.
Some of the research proposed applies to dense matrices and some to sparse
However, we emphasize that most sparse f(A)
the computation of f(B) for a dense (and much smaller) B.
Therefore work for the dense case has immediate payoffs
for the sparse case, too.
of publications concerning the theory and computation of matrix functions.
Matlab M-files for the Schur-Parlett algorithm for
computing matrix functions are available in the tar file
The toolbox is distributed under
this MIT license.
This tar file contains a MEX-file
swap.c that allows the algorithm to use
the LAPACK routine xTREXC.
Prof. Nick Higham
is the Richardson Professor of Applied Mathematics
at the University of Manchester.
Dr. Philip Davies
is a Research Associate in the Department of Mathematics at the
University of Manchester.
Physics LATtice gauge theory)
preprint archive server.
PCMF: Parallel Computation of
Accelerator Center (SLAC),
The Worldwide Web Virtual
Library: High-Energy Physics