http://www.cnr.it/ontology/cnr/individuo/prodotto/ID61352

A Sparse Nonsymmetric Eigensolver for Distributed Memory Architectures (Articolo in rivista)

Type

Label

A Sparse Nonsymmetric Eigensolver for Distributed Memory Architectures (Articolo in rivista) (literal)

Anno

Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#doi

Alternative label

Mario Rosario Guarracino (1), Francesca Perla (2), Paolo Zanetti (2) (2008)
A Sparse Nonsymmetric Eigensolver for Distributed Memory Architectures
in International journal of parallel, emergent and distributed systems (Print); TAYLOR & FRANCIS LTD, London (Regno Unito)
(literal)

Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#autori

Pagina inizio

Pagina fine

Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#url

Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#numeroVolume

Rivista

Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#pagineTotali

Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#numeroFascicolo

Note

Http://www.cnr.it/ontology/cnr/pubblicazioni.owl#affiliazioni

(1) Institute for High Performance Computing and Networking, Italian National Research Council, Via P. Castellino 111, 80131, Naples, Italy (2) University of Naples Parthenope, Via Medina 40, 80133, Naples, Italy (literal)

Titolo

A Sparse Nonsymmetric Eigensolver for Distributed Memory Architectures (literal)

Abstract

In this work, we propose an efficient parallel implementation of the nonsymmetric block Lanczos algorithm for the computation of few extreme eigenvalues, and corresponding eigenvectors, of real nonhermitian matrices for distributed memory multicomputers. The reorganisation of the block Lanczos algorithm implemented allows to exploit a coarse-grained parallelism and to harness the computational power of the target architectures. The computational kernels of the algorithm are matrix-matrix multiplications, with dense and sparse factors, QR factorisation and singular value decomposition. To reduce the total amount of communication involved in the matrix-matrix multiplication with a sparse factor, we substitute each matrix appearing in the algorithm with its transpose. Then, we develop an efficient parallelisation of the matrix-matrix multiplication when the second factor is sparse. Some other linear algebra operations are performed using ScaLAPACK library. The parallel eigensolver has been tested on a cluster of PCs. All reported results show the proposed algorithm is efficient on the target architectures for problems of adequate dimension. (literal)

Editore