Benchmarks of Cuda-Based GMRES Solver for Toeplitz and Hankel Matrices and Applications to Topology Optimization of Photonic Components

Iu B. Minin, S. A. Matveev, M. V. Fedorov, I. E. Zacharov, S. G. Rykovanov

Research output: Contribution to journalArticlepeer-review

Abstract

Generalized Minimal Residual Method (GMRES) was benchmarked on many types of GPUs for solving linear systems based on dense and sparse matrices. However, there are still no GMRES implementation benchmarks on Tesla V100 compared to GTX 1080 Ti ones or even for Toeplitz-like matrices. The introduced software consists of a Python module and a C++ library which enable to manage streams for concurrent computations of separated linear systems on a GPU (and GPUs). The GMRES solver is parallelized for running on a NVIDIA GPGPU accelerator. The parallelization efficiency is explored when GMRES is applied to solve (Helmholtz equation) linear systems based on the use of Green’s Function Integral Equation Method (GFIEM) for computing electric field distribution in the design domain. The proposed implementation shew the maximal speedup of 55 (t¯=0.017s) and of 125 (t¯=0.77s) for 1024 × 1024 (on GTX 1080 Ti) and 8192 × 8192 (on Tesla V100) dense Toeplitz matrices generated from GFIEM. 1024 × 1024 resolution provides accuracy 6.1% that can be acceptable according to testing and demonstrating on gradient computations and topology optimization. We open up possibilities for robust topology optimization of passive photonic integrated components. That has the advantage, e. g., of faster and more accurate designing photonic components on a PC without a supercomputer.

Original languageEnglish
Pages (from-to)438-452
Number of pages15
JournalComputational Mathematics and Modeling
Volume32
Issue number4
DOIs
Publication statusPublished - Oct 2021

Keywords

  • CUDA C++
  • GMRES
  • GPU
  • Helmholtz equation
  • parallelization
  • performance
  • Python module
  • Toeplitz-like matrix
  • topology optimization

Fingerprint

Dive into the research topics of 'Benchmarks of Cuda-Based GMRES Solver for Toeplitz and Hankel Matrices and Applications to Topology Optimization of Photonic Components'. Together they form a unique fingerprint.

Cite this