Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/normal_derivative.jl @ 1955:627d94d0f333 feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix boundary iterator in normal_derivative
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 10 Feb 2025 15:38:38 +0100
parents 835b1dcee38e
children
comparison
equal deleted inserted replaced
1954:b0915f43b122 1955:627d94d0f333
33 k = grid_id(boundary) 33 k = grid_id(boundary)
34 b_indices = boundary_indices(g, boundary) 34 b_indices = boundary_indices(g, boundary)
35 35
36 # Compute the weights for the logical derivatives 36 # Compute the weights for the logical derivatives
37 g⁻¹ = map(inv, metric_tensor(g)) 37 g⁻¹ = map(inv, metric_tensor(g))
38 α = map(CartesianIndices(g⁻¹)[b_indices...]) do I # TODO: Fix iterator here 38 α = map(b_indices) do I
39 gᵏⁱ = g⁻¹[I][k,:] 39 gᵏⁱ = g⁻¹[I][k,:]
40 gᵏᵏ = g⁻¹[I][k,k] 40 gᵏᵏ = g⁻¹[I][k,k]
41 41
42 gᵏⁱ./sqrt(gᵏᵏ) 42 gᵏⁱ./sqrt(gᵏᵏ)
43 end 43 end