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
Symbolics Writing element (matrix) inverse General Usage	0	142	October 19, 2023
Symbolics.jl: Why doesn't a matrix multiplied by its inverse simplify to the identity? Specific Domains symbolics	3	1102	February 20, 2023
Inverting a matrix of symbolic equations General Usage linearalgebra , symbolics	6	468	April 6, 2022
LinearAlgebra.inv! not exported, in-place inversion Performance linearalgebra	5	237	February 1, 2024
Symbolics - Identity matrices for beginners General Usage linearalgebra , symbolics	5	1539	August 31, 2021

Matrix multiplied by inverse matrix does not reproduce identity matrix

Related topics