Un esquema paralelo para el cálculo del pseudoespectro de matrices de gran magnitud
Visualitza/Obre
10.1016/j.rimni.2013.10.004
Inclou dades d'ús des de 2022
Cita com:
hdl:2117/26196
Càtedra / Departament / Institut
Universitat Politècnica de Catalunya. Departament d'Arquitectura de Computadors
Tipus de documentArticle
Data publicació2015-01-01
EditorUniversitat Politècnica de Catalunya. CIMNE
Condicions d'accésAccés obert
Tots els drets reservats. Aquesta obra està protegida pels drets de propietat intel·lectual i
industrial corresponents. Sense perjudici de les exempcions legals existents, queda prohibida la seva
reproducció, distribució, comunicació pública o transformació sense l'autorització del titular dels drets
Abstract
El pseudoespectro es una gran herramienta para el estudio de sistemas dinámicos asociados a matrices no normales. En las últimas décadas su estudio y aplicación se ha intensificado, por lo que realizar su cálculo de forma eficiente resulta de especial interés para la comunidad científica. Para matrices de grandes dimensiones se han empleado diferentes métodos, entre los cuales destacan los métodos de proyección en espacios de Krylov. En este trabajo se utiliza la idea propuesta por Wright y Trefethen para aproximar el pseudoespectro de la matriz usando una proyección Hm de menor tamaño. Adicionalmente, se propone una descomposición del dominio en subregiones asignando a cada procesador una subregión. Cada procesador calcula el menor valor singular de la matriz (zI − Hm ) para todos los valores z = x + yi que representan los puntos de la subregión asignada. Se realizaron diferentes experimentos numéricos cuyos resultados fueron contrastados con los obtenidos en la bibliografía. En todos los casos estudiados el método implementado muestra una reducción del tiempo de ejecución respecto a la versión secuencial del programa desde 41x hasta 101x . 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.
CitacióOtero, 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.
ISSN0213-1315
1886-158X
1886-158X
Fitxers | Descripció | Mida | Format | Visualitza |
---|---|---|---|---|
Pseudoespectra_Postprint.pdf | Article pseudospectra | 575,5Kb | Visualitza/Obre |