Research Article
Approximating the Inverse of a Square Matrix with Application in Computation of the Moore-Penrose Inverse
Id = SparseArrayi_, i_} -> 1.}, {n, n; | V = DiagonalMatrix@SparseArray1/NormalDiagonalA; | DoV1 = SparseArrayV; | V2 = ChopA.V1; V3 = 3 Id + V2.(−3 Id + V2); | V4 = SparseArrayV2.V3; | V = Chop−(1/4) V1.V3.SparseArray−13 Id | + V4.(15 Id + V4.(−7 Id + V4)); | PrintV;Li = NNormId − V.A, 1; | Print"The residual norm is:" | Columni}, Frame -> All, FrameStyle -> DirectiveBlue | ColumnLi, Frame -> All, FrameStyle -> DirectiveBlue; | , {i, 1; // AbsoluteTiming |
|