A parallel generator of non-Hermitian matrices computed from given spectra
Résumé
Iterative linear algebra methods to solve linear systems and eigenvalue problems with non-Hermitian matrices are important for both the simulation arising from diverse scientific fields and the applications related to big data, machine learning and artificial intelligence. The spectral property of these matrices has impacts on the convergence of these solvers. Moreover, with the increase of the size of applications, iterative methods are implemented in parallel on clusters. Analysis of their behaviors with non-Hermitian matrices on supercomputers is so complex that we need to generate large-scale matrices with different given spectra for benchmarking. These test matrices should be non-Hermitian and non-trivial, with high dimension. This paper highlights a scalable matrix generator that constructs large sparse matrices using the user-defined spectrum, and the eigenvalues of generated matrices are ensured to be the same as the predefined spectrum. This generator is implemented on CPUs and multi-GPUs platforms, with good strong and weak scaling performance on several supercomputers. We also propose a method to verify its ability to guarantee the given spectra. Finally, we give an example to evaluate the numerical properties and parallel performance of iterative methods using this matrix generator.
Origine : Fichiers produits par l'(les) auteur(s)
Loading...