The pseudospectra is a powerful tool to study the behavior of dynamic systems associated to non-normalmatrices. Studies and applications have increased in the last decades, thus, its efficient computation hasbecome of interest for the scientific community. In the large scale setting, different approaches have beenproposed, some of them based on projection on Krylov subspaces. In this work we use the idea proposed byWright and Trefethen to approximate the pseudospectra of a matrix A using a projection Hmof smallersize. Additionally, we propose a domain decomposition of the interest region into subregions whichare assigned to a set of processors. Each processor calculates the minimal singular values of matrices(zI - Hm) where z = x + yi represents a point of the corresponding subregion. We conduct a numericalexperimentation comparing the results with those on the literature of the topic. In all cases the proposedscheme shows a reduction in CPU time with respect to the sequential version, achieving from 41x to 101x.
CitationOtero, B.; Astudillo, R.; Castillo, Z. Un esquema paralelo para el cálculo del pseudoespectro de matrices de gran magnitud. "Revista internacional de métodos numéricos para cálculo y diseño en ingeniería", 01 Gener 2015, vol. 31, núm. 1, p. 8-12.
All rights reserved. This work is protected by the corresponding intellectual and industrial property rights. Without prejudice to any existing legal exemptions, reproduction, distribution, public communication or transformation of this work are prohibited without permission of the copyright holder. If you wish to make any use of the work not provided for in the law, please contact: email@example.com