Skip to content

Commit f67bfd5

Browse files
committed
use sum(abs2()) instead of dot
1 parent 14d75de commit f67bfd5

File tree

1 file changed

+2
-9
lines changed

1 file changed

+2
-9
lines changed

src/pdmat.jl

Lines changed: 2 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -68,15 +68,8 @@ end
6868

6969
### quadratic forms
7070

71-
function quad(a::PDMat, x::AbstractVector)
72-
z = a.chol.U * x
73-
dot(z, z)
74-
end
75-
76-
function invquad(a::PDMat, x::AbstractVector)
77-
z = a.chol.L \ x
78-
dot(z, z)
79-
end
71+
quad(a::PDMat, x::AbstractVector) = sum(abs2, a.chol.U * x)
72+
invquad(a::PDMat, x::AbstractVector) = sum(abs2, a.chol.L \ x)
8073

8174
"""
8275
quad!(r::AbstractArray, a::AbstractPDMat, x::StridedMatrix)

0 commit comments

Comments
 (0)