Skip to content
This repository was archived by the owner on Dec 16, 2024. It is now read-only.

Commit b8ca474

Browse files
authored
Exploit the first law of Hermitian matrices
Don't know why this one went unnoticed for so long.
1 parent 8e65c8c commit b8ca474

File tree

1 file changed

+1
-14
lines changed

1 file changed

+1
-14
lines changed

src/cholmod.jl

Lines changed: 1 addition & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -1583,20 +1583,7 @@ function \(A::RealHermSymComplexHermF64SSL, B::StridedVecOrMatInclAdjAndTrans)
15831583
end
15841584
end
15851585
end
1586-
function \(adjA::Adjoint{<:Any,<:RealHermSymComplexHermF64SSL}, B::StridedVecOrMatInclAdjAndTrans)
1587-
A = adjA.parent
1588-
F = cholesky(A; check = false)
1589-
if issuccess(F)
1590-
return \(adjoint(F), B)
1591-
else
1592-
ldlt!(F, A; check = false)
1593-
if issuccess(F)
1594-
return \(adjoint(F), B)
1595-
else
1596-
return \(adjoint(lu(SparseMatrixCSC{eltype(A), SuiteSparse_long}(A))), B)
1597-
end
1598-
end
1599-
end
1586+
\(adjA::Adjoint{<:Any,<:RealHermSymComplexHermF64SSL}, B::StridedVecOrMatInclAdjAndTrans) = parent(adjA)\B
16001587

16011588
## Other convenience methods
16021589
function diag(F::Factor{Tv}) where Tv

0 commit comments

Comments
 (0)