An FPGA-based analysis of trade-offs in the presence of ill-conditioning and different precision levels in computations

Several areas, such as physical and health sciences, require the use of matrices as fundamental tools for solving various problems. Matrices are used in real-life contexts, such as control, automation, and optimization, wherein results are expected to improve with increase of computational precision. However, special attention should be paid to ill-conditioned matrices, which can produce unstable systems; an inadequate handling of precision might worsen results since the solution found for data with errors might be too far from the one for data without errors besides increasing other costs in hardware resources and critical paths. In this paper, we make a wake-up call, using 2 × 2 matrices to show how ill-conditioning and precision can affect system design (resources, cost, etc.). We first demonstrate some examples of real-life problems where ill-conditioning is present in matrices obtained from the discretization of the operational equations (ill-posed in the sense of Hadamard) that model these problems. If these matrices are not handled appropriately (i.e., if ill-conditioning is not considered), large errors can result in the computed solutions to the systems of equations in the presence of errors. Furthermore, we illustrate the generated effect in the calculation of the inverse of an ill-conditioned matrix when its elements are approximated by truncation. We present two case studies to illustrate the effects on calculation errors caused by increasing or reducing precision to s digits. To illustrate the costs, we implemented the adjoint matrix inversion algorithm on different field-programmable gate arrays (FPGAs), namely, Spartan-7, Artix-7, Kintex-7, and Virtex-7, using the full-unrolling hardware technique. The implemented architecture is useful for analyzing trade-offs when precision is increased; this also helps analyze performance, efficiency, and energy consumption. By means of a detailed description of the trade-offs among these metrics, concerning precision and ill-conditioning, we conclude that the need for resources seems to grow not linearly when precision is increased. We also conclude that, if error is to be reduced below a certain threshold, it is necessary to determine an optimal precision point. Otherwise, the system becomes more sensitive to measurement errors and a better alternative would be to choose precision carefully, and/or to apply regularization or preconditioning methods, which would also reduce the resources required.

[1]  Gumersindo Verdú,et al.  System Matrix Analysis for Computed Tomography Imaging , 2015, PloS one.

[2]  Paul Adrien Maurice Dirac,et al.  The lorentz transformation and absolute time , 1953 .

[3]  M. M. Morín-Castillo,et al.  Prototipo de Silla de Ruedas Dirigida Usando Parpadeos , 2018 .

[4]  D. Potyagaylo,et al.  Solving the Inverse Problem of Electrocardiography on the Endocardium Using a Single Layer Source , 2019, Front. Physiol..

[5]  Jesse V Jokerst,et al.  Comparison of Deconvolution Filters for Photoacoustic Tomography , 2016, PloS one.

[6]  Simon C. Potter,et al.  A Genome-Wide Association Search for Type 2 Diabetes Genes in African Americans , 2012, PLoS ONE.

[7]  Ramon Casanova,et al.  Classification of Structural MRI Images in Alzheimer's Disease from the Perspective of Ill-Posed Problems , 2012, PloS one.

[8]  Joseph R. Cavallaro,et al.  FPGA Implementation of Matrix Inversion Using QRD-RLS Algorithm , 2005, Conference Record of the Thirty-Ninth Asilomar Conference onSignals, Systems and Computers, 2005..

[9]  Carlos H. Llanos,et al.  FPGA implementation of large-scale matrix inversion using single, double and custom floating-point precision , 2012, 2012 VIII Southern Conference on Programmable Logic.

[10]  Joakim Munkhammar,et al.  Polynomial probability distribution estimation using the method of moments , 2017, PloS one.

[11]  V. Elamaran,et al.  An approach to design a matrix inversion hardware module using FPGA , 2014, 2014 International Conference on Control, Instrumentation, Communication and Computational Technologies (ICCICCT).

[12]  Ulrik Plesner Jacobsen,et al.  Interpretation of Appearance: The Effect of Facial Features on First Impressions and Personality , 2014, PloS one.

[13]  Malabika Sarker,et al.  Effective maternal, newborn and child health programming among Rohingya refugees in Cox’s Bazar, Bangladesh: Implementation challenges and potential solutions , 2020, PloS one.

[14]  Zhenghui Hu,et al.  Exploiting Magnetic Resonance Angiography Imaging Improves Model Estimation of BOLD Signal , 2012, PloS one.

[15]  R. Feynman,et al.  The Theory of a general quantum system interacting with a linear dissipative system , 1963 .

[16]  Suheel Abdullah,et al.  Accelerated sparsity based reconstruction of compressively sensed multichannel EEG signals , 2020, PloS one.

[17]  A. Kirsch An Introduction to the Mathematical Theory of Inverse Problems , 1996, Applied Mathematical Sciences.

[18]  Martin D. Smith,et al.  The Fishery Performance Indicators: A Management Tool for Triple Bottom Line Outcomes , 2015, PloS one.

[19]  Noriyuki Kushida Condition Number Estimation of Preconditioned Matrices , 2015, PloS one.

[20]  Bin Xing,et al.  MW151 Inhibited IL-1β Levels after Traumatic Brain Injury with No Effect on Microglia Physiological Responses , 2016, PloS one.

[21]  Vladimir E. Zuev,et al.  Inverse problems of lidar sensing of the atmosphere , 1983 .

[22]  Fozia Hanif Khan,et al.  A Space Efficient Flexible Pivot Selection Approach to Evaluate Determinant and Inverse of a Matrix , 2014, PloS one.

[23]  Jose Jacobo Oliveros Oliveros,et al.  STABLE IDENTIFICATION OF SOURCES LOCATED ON SEPARATION INTERFACES OF TWO DIFFERENT HOMOGENEOUS MEDIA , 2019, Advances in Differential Equations and Control Processes.

[24]  J. J. C. Mones,et al.  Stable numerical solution of the Cauchy problem for the Laplace equation in irregular annular regions , 2017 .