Mercurial > repos > public > sbplib_julia
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]
