Skip to content

Commit c7478a3

Browse files
committed
avoid Adjoint wrapper for chol_lower(::Matrix)
1 parent 45a17ae commit c7478a3

File tree

1 file changed

+2
-1
lines changed

1 file changed

+2
-1
lines changed

src/chol.jl

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -5,7 +5,8 @@ CholType{T,S<:AbstractMatrix} = Cholesky{T,S}
55
chol_lower(a::CholType) = a.uplo === 'L' ? a.L : a.U'
66
chol_upper(a::CholType) = a.uplo === 'U' ? a.U : a.L'
77

8-
chol_lower(a::Matrix) = chol_lower(cholesky(a))
8+
# For a dense Matrix, the following allows us to avoid the Adjoint wrapper:
9+
chol_lower(a::Matrix) = cholesky(Hermitian(a, :L)).L
910

1011
if HAVE_CHOLMOD
1112
CholTypeSparse{T} = SuiteSparse.CHOLMOD.Factor{T}

0 commit comments

Comments
 (0)