A GPU parallel randomized CUR compression method for the Method of Moments

In this work, we propose a GPU parallel implementation of the randomized CUR (or Pseudo Skeleton) Approximation to compress the -matrices of linear systems that arise in the discretization of integral equations modeling electromagnetic scattering problems. This compression method is highly paralleli...

Descripción completa

Detalles Bibliográficos
Autores: López Menchón, Héctor|||0000-0001-5776-1423; Heldring, Alexander|||0000-0003-2011-2096; Úbeda Farré, Eduard|||0000-0001-6759-0445; Rius Casals, Juan Manuel|||0000-0003-0606-5422
Tipo de documento: Artículo
Fecha de publicación:2023
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Lenguaje:Inglés
OAI Identifier:oai:upcommons.upc.edu:2117/391415
Acceso en línea:https://hdl.handle.net/2117/391415
https://dx.doi.org/10.1016/j.cpc.2023.108696
Nivel de acceso:Acceso abierto
Palabra clave:Estimation theory
Electric field integral equation
Graphics processing unit
Low-rank approximation
Method of Moments
H-matrices
Randomized
Estimació, Teoria de l'
Àrees temàtiques de la UPC::Enginyeria de la telecomunicació
Descripción
Sumario:In this work, we propose a GPU parallel implementation of the randomized CUR (or Pseudo Skeleton) Approximation to compress the -matrices of linear systems that arise in the discretization of integral equations modeling electromagnetic scattering problems. This compression method is highly parallelizable, in contrast with other similar methods such as the Adaptive Cross Approximation. It involves dense linear algebra computations that can be efficiently implemented on a GPU device. Besides, a stochastic convergence criterion is introduced to minimize the communication between the host and the device. Testing the code with standard cases shows the efficiency and accuracy of the method.