From cfe9c40aeca2b696cf5f057aab6212201c1dffc9 Mon Sep 17 00:00:00 2001
From: Michiel Cottaar <michiel.cottaar@ndcn.ox.ac.uk>
Date: Mon, 12 Feb 2024 10:24:19 +0000
Subject: [PATCH] Set default value for qvec
---
src/overlapping/abstract.jl | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/overlapping/abstract.jl b/src/overlapping/abstract.jl
index c41dd6f..689586e 100644
--- a/src/overlapping/abstract.jl
+++ b/src/overlapping/abstract.jl
@@ -168,7 +168,7 @@ function qvec(ao::AbstractOverlapping, index1::Union{Nothing, Integer}, index2::
if (index1 isa Number) && (index1 == index2)
return zeros(3)
end
- sum(qvec.(get_parts(ao, index1, index2)))
+ sum(qvec.(get_parts(ao, index1, index2)); init=zero(SVector{3, Float64}))
end
qvec(ao::AbstractOverlapping) = qvec(ao, nothing, nothing)
--
GitLab