On Moore‐Penrose Pseudoinverse Computation for Stiffness Matrices Resulting from Higher Order Approximation

Abstract

Computing the pseudoinverse of a matrix is an essential component of many computational methods. It arises in statistics, graphics, robotics, numerical modeling, and many more areas. Therefore, it is desirable to select reliable algorithms that can perform this operation efficiently and robustly. A demanding benchmark test for the pseudoinverse computation was introduced. The stiffness matrices for higher order approximation turned out to be such tough problems and therefore can serve as good benchmarks for algorithms of the pseudoinverse computation. It was found out that only one algorithm, out of five known from literature, enabled us to obtain acceptable results for the pseudoinverse of the proposed benchmark test.

Impulse Delay in the Cardiac Conduction System

“…Therefore, each matrix here is presupposed Moore-Penrose invertible. Details on Moore-Penrose pseudoinverse of matrices may be found in [17,18] and in profuse literatures. The beauty of non-square systems is their less amenability to modelling errors [19].…”

“…consisting of triple matrices P, Q, R assumed to be positive definite, and where in the present case A 0 and A 1 are Moore-Penrose invertible [18,19,20]). Therefore:…”