An efficient bound for the condition number of the matrix exponential

A new bound for the condition number of the matrix exponential is presented. Using the bound, we propose an efficient approximation to the condition number, denoted by kappa(g) (S, X), that avoids the computation of the Frechet derivative of the matrix exponential that underlies condition number estimation in the existing algorithms. We exploit the identity e(X) = (e(X/2))(2s) for a nonnegative integer s with the properties of the Frechet derivative operator to obtain the bound. Our cost analysis reveals that considerable computational savings are possible since estimating the condition number by the existing algorithms requires several invocation of the Frechet derivative of the matrix exponential whose single invocation costs as twice as the cost of the matrix exponential itself. The bound and hence kappa(g)(s, X) only involve Frechet derivative of a monomial of degree 2(s), which can be computed exactly in 2s matrix multiplications. We propose two versions of the scaling and squaring algorithm that implement kappa(g)(s, X). Our numerical experiments show that kappa(g)(s, X) captures the behavior of the condition number and moreover outperforms the condition number in the estimation of relative forward errors for a wide range of problems. (C) 2015 The Author. Production and hosting by Elsevier B.V. on behalf of Taibah University.