Merge pull request #119 from ElOceanografo/patch-1

[Breaking] Make inverting a PDSparseMat throw an error
JuliaStats · Sep 19, 2020 · d50e35a · d50e35a
2 parents 9e1d882 + eca5218
commit d50e35a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/pdsparsemat.jl b/src/pdsparsemat.jl
@@ -51,7 +51,7 @@ end
 
 ### Algebra
 
-Base.inv(a::PDSparseMat{T}) where {T<:Real} = PDMat( a\eye(T,a.dim) )
+Base.inv(a::PDSparseMat{T}) where {T<:Real} = PDMat(inv(a.mat))
 LinearAlgebra.logdet(a::PDSparseMat) = logdet(a.chol)
 
 ### whiten and unwhiten