In this paper we consider methods for evaluating both exp(A) and exp(τA)q1 where exp(·) is the exponential function, A is a sparse skew-symmetric matrix of large dimension, q1 is a given vector, and τ is a scaling factor. The proposed method is based on two main steps: A is factorized into its tridiagonal form H by the well-known Lanczos iterative process, and then exp(A) is derived making use of an effective Schur decomposition of H. The procedure takes full advantage of the sparsity of A and of the decay behavior of exp(H). Several applications and numerical tests are also reported.

Computation of Exponentials of Large Sparse Skew Symmetric Matrices

DEL BUONO, Nicoletta;LOPEZ, Luciano;
2005

Abstract

In this paper we consider methods for evaluating both exp(A) and exp(τA)q1 where exp(·) is the exponential function, A is a sparse skew-symmetric matrix of large dimension, q1 is a given vector, and τ is a scaling factor. The proposed method is based on two main steps: A is factorized into its tridiagonal form H by the well-known Lanczos iterative process, and then exp(A) is derived making use of an effective Schur decomposition of H. The procedure takes full advantage of the sparsity of A and of the decay behavior of exp(H). Several applications and numerical tests are also reported.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: http://hdl.handle.net/11586/113313
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 27
  • ???jsp.display-item.citation.isi??? 23
social impact