Skip to content

Avoid inverting matrices in EigenAnalysis #10

@juliohm

Description

@juliohm

We can be more clever and avoid inversions of S here:

function drsproj(λ, V)
Λ = Diagonal(sqrt.(λ))
S = V * inv(Λ)
S, inv(S)
end
function sdsproj(λ, V)
Λ = Diagonal(sqrt.(λ))
S = V * inv(Λ) * transpose(V)
S, inv(S)
end

We should be able to express these inverses by swapping the order of the products and by transposing the matrices instead.

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions