Hello, I found that multiplying a matrix by its inverse does not yield an identity matrix. Please see below a MWE: matrix = [5.93419 4.8067; 4.8067 5.93419] invmatrix = inv(matrix) matrix*invmatrix The last line generates the following output: 2×2 Matrix{Float64}: 1.0 4.44089e-16 4.4…

It’s nearly symmetric. Because of the imprecision inherent in floating point numbers , it’s best to use approximate comparisons like this: julia> varc ≈ varc' true The approximate equality symbol can be typed in the REPL with \approx<TAB>, or you can use the equivalent function call: isapprox(varc,…

Matrix multiplied by inverse matrix does not reproduce identity matrix

General Usage

rafael.guerra August 3, 2021, 12:58am 2

One solution is posted here, for display purposes.

1 Like

Topic		Replies	Views
Inverse of a symmetric matrix is not symmetric? General Usage	16	3442	April 4, 2018
The LinearAlgebra.jl package's `inv()` function does not preserve matrix symmetry General Usage	4	335	March 6, 2021
Inverse of symmetric Matrix does not give proper output Performance question	3	459	May 26, 2020
Is there a bug in large matrix inverse? General Usage question	6	1834	May 20, 2020
Why same calculations on matrices yield different results? General Usage question	4	396	March 28, 2022

Matrix multiplied by inverse matrix does not reproduce identity matrix

Related topics