Mercurial > repos > public > sbplib_julia

comparison src/SbpOperators/boundaryops/normal_derivative.jl @ 1786:1f42944d4a72 feature/sbp_operators/laplace_curvilinear

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix errors introduced by merge
author Jonatan Werpers <jonatan@werpers.com>
date Mon, 16 Sep 2024 11:03:37 +0200
parents 29b96fc75bee
children 835b1dcee38e
comparison
equal deleted inserted replaced
1785:96f8cad255b4 1786:1f42944d4a72
49 49
50 50
51 # Assemble difference operator 51 # Assemble difference operator
52 mapreduce(+,1:ndims(g)) do i 52 mapreduce(+,1:ndims(g)) do i
53 if i == k 53 if i == k
54 ∂_ξᵢ = normal_derivative(logicalgrid(g), stencil_set, boundary) 54 ∂_ξᵢ = normal_derivative(logical_grid(g), stencil_set, boundary)
55 else 55 else
56 e = boundary_restriction(logicalgrid(g), stencil_set, boundary) 56 e = boundary_restriction(logical_grid(g), stencil_set, boundary)
57 ∂_ξᵢ = σ ∘ e ∘ first_derivative(logicalgrid(g), stencil_set, i) 57 ∂_ξᵢ = σ ∘ e ∘ first_derivative(logical_grid(g), stencil_set, i)
58 end 58 end
59 59
60 αᵢ = componentview(α,i) 60 αᵢ = componentview(α,i)
61 DiagonalTensor(αᵢ) ∘ ∂_ξᵢ 61 DiagonalTensor(αᵢ) ∘ ∂_ξᵢ
62 end 62 end