On the Computation of Square Roots and Inverse Square Roots of Gram Matrices for Surface Integral Equations in Electromagnetics

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Chen, Rui, Merlini, Adrien, Andriulli, Francesco P.
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866916738699886592
author Chen, Rui
Merlini, Adrien
Andriulli, Francesco P.
author_facet Chen, Rui
Merlini, Adrien
Andriulli, Francesco P.
contents Surface integral equations (SIEs)-based boundary element methods are widely used for analyzing electromagnetic scattering scenarii. However, after discretization of SIEs, the spectrum and eigenvectors of the boundary element matrices are not usually representative of the spectrum and eigenfunctions of the underlying surface integral operators, which can be problematic for methods that rely heavily on spectral properties. To address this issue, we delineate some efficient algorithms that allow for the computation of matrix square roots and inverse square roots of the Gram matrices corresponding to the discretization scheme, which can be used for revealing the spectrum of standard electromagnetic integral operators. The algorithms, which are based on properly chosen expansions of the square root and inverse square root functions, are quite effective when applied to several of the most relevant Gram matrices used for boundary element discretizations in electromagnetics. Tables containing different sets of expansion coefficients are provided along with comparative numerical experiments that evidence advantages and disadvantages of the different approaches. In addition, to demonstrate the spectrum-revealing properties of the proposed techniques, they are applied to the discretization of the problem of scattering by a sphere for which the analytic spectrum is known.
format Preprint
id arxiv_https___arxiv_org_abs_2405_11167
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle On the Computation of Square Roots and Inverse Square Roots of Gram Matrices for Surface Integral Equations in Electromagnetics
Chen, Rui
Merlini, Adrien
Andriulli, Francesco P.
Numerical Analysis
Surface integral equations (SIEs)-based boundary element methods are widely used for analyzing electromagnetic scattering scenarii. However, after discretization of SIEs, the spectrum and eigenvectors of the boundary element matrices are not usually representative of the spectrum and eigenfunctions of the underlying surface integral operators, which can be problematic for methods that rely heavily on spectral properties. To address this issue, we delineate some efficient algorithms that allow for the computation of matrix square roots and inverse square roots of the Gram matrices corresponding to the discretization scheme, which can be used for revealing the spectrum of standard electromagnetic integral operators. The algorithms, which are based on properly chosen expansions of the square root and inverse square root functions, are quite effective when applied to several of the most relevant Gram matrices used for boundary element discretizations in electromagnetics. Tables containing different sets of expansion coefficients are provided along with comparative numerical experiments that evidence advantages and disadvantages of the different approaches. In addition, to demonstrate the spectrum-revealing properties of the proposed techniques, they are applied to the discretization of the problem of scattering by a sphere for which the analytic spectrum is known.
title On the Computation of Square Roots and Inverse Square Roots of Gram Matrices for Surface Integral Equations in Electromagnetics
topic Numerical Analysis
url https://arxiv.org/abs/2405.11167