changeset 1955:627d94d0f333 feature/sbp_operators/laplace_curvilinear

Fix boundary iterator in normal_derivative
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 10 Feb 2025 15:38:38 +0100
parents b0915f43b122
children 496b8d3903c6
files src/SbpOperators/boundaryops/normal_derivative.jl
diffstat 1 files changed, 1 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/src/SbpOperators/boundaryops/normal_derivative.jl	Sat Feb 08 09:38:58 2025 +0100
+++ b/src/SbpOperators/boundaryops/normal_derivative.jl	Mon Feb 10 15:38:38 2025 +0100
@@ -35,7 +35,7 @@
 
     # Compute the weights for the logical derivatives
     g⁻¹ = map(inv, metric_tensor(g))
-    α = map(CartesianIndices(g⁻¹)[b_indices...]) do I # TODO: Fix iterator here
+    α = map(b_indices) do I
         gᵏⁱ = g⁻¹[I][k,:]
         gᵏᵏ = g⁻¹[I][k,k]