Convert ev to a generated ntuple

JuliaGeometry · Oct 1, 2024 · c2e6233 · c2e6233
1 parent 1587085
commit c2e6233
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/differentiation.jl b/src/differentiation.jl
@@ -23,8 +23,8 @@ function jacobian(
     # Get the partial derivative along the n'th axis via finite difference
     #   approximation, where ts is the current parametric position
     function ∂ₙr(ts, n)
-        εv = zeros(T, length(ts))
-        εv[n] = ε
+        # Construct a tuple with ε in the n'th element and zeros otherwise
+        εv = ntuple(i -> i == n ? T(ε) : zero(T), length(ts))
 
         # Select orientation of finite-diff
         if ts[n] < T(0.01)