Improving the accuracy of the adaptive cross approximation with a convergence criterion based on random sampling

The accuracy of the Adaptive Cross Approximation (ACA) algorithm, a popular method for the compression of low rank matrix blocks in Method of Moment computations, is sometimes seriously compromised by unpredictable errors in the convergence criterion. This paper proposes an alternative criterion, ba...

Descripción completa

Detalles Bibliográficos
Autores: Heldring, Alexander|||0000-0003-2011-2096, Úbeda Farré, Eduard|||0000-0001-6759-0445, Rius Casals, Juan Manuel|||0000-0003-0606-5422
Tipo de recurso: artículo
Fecha de publicación:2020
País:España
Institución:Universitat Politècnica de Catalunya (UPC)
Repositorio:UPCommons. Portal del coneixement obert de la UPC
Idioma:inglés
OAI Identifier:oai:upcommons.upc.edu:2117/337061
Acceso en línea:https://hdl.handle.net/2117/337061
https://dx.doi.org/10.1109/TAP.2020.3010857
Access Level:acceso abierto
Palabra clave:Computer algorithms
Signal processing
Adaptive cross approximation (ACA)
Computational electromagnetics
Fast solvers
Method of moments
Algorismes computacionals
Tractament del senyal
Àrees temàtiques de la UPC::Enginyeria de la telecomunicació::Processament del senyal
Descripción
Sumario:The accuracy of the Adaptive Cross Approximation (ACA) algorithm, a popular method for the compression of low rank matrix blocks in Method of Moment computations, is sometimes seriously compromised by unpredictable errors in the convergence criterion. This paper proposes an alternative criterion, based on global sampling of the error in the elements of the ACA compressed matrix. The sampling error depends on the size of the sample but also on the population distribution of the error, which makes it difficult to control the error independently of the underlying problem. However, as argued and demonstrated in the paper, the distribution of the error converges to the same unique probability distribution function for all low rank matrices. Complementing the sampling criterion with a simple mechanism to detect this convergence, we arrive at a criterion that controls the error irrespective of the underlying problem. As a practical example the RCS of a moderate size metallic ogive is computed to illustrate the merits of the proposed criterion. The proposed algorithm may also be useful in other methods that approximate low-rank matrices by interpolation of a reduced set of its elements.