Skip to content

Commit

Permalink
add comment on Nanson's formula for normals
Browse files Browse the repository at this point in the history
  • Loading branch information
jlchan committed Sep 30, 2024
1 parent fa1e004 commit c09676e
Showing 1 changed file with 3 additions and 2 deletions.
5 changes: 3 additions & 2 deletions src/geometric_functions.jl
Original file line number Diff line number Diff line change
Expand Up @@ -54,8 +54,9 @@ function geometric_factors(x, y, z, Dr, Ds, Dt, Filters=(I, I, I))
return rxJ, sxJ, txJ, ryJ, syJ, tyJ, rzJ, szJ, tzJ, J
end

# physical normals are computed via G * nhatJ, where G = matrix of J-scaled geometric terms
# Here, Vf is a face interpolation matrix which maps interpolation nodes to face nodes.
# physical outward normals are computed via Nanson's formula: G * nhatJ, where
# G = matrix of J-scaled geometric terms. Here, Vf is a face interpolation matrix
# which maps interpolation nodes to face nodes.
function compute_normals(geo::SMatrix{Dim, Dim}, Vf, nrstJ...) where {Dim}
nxyzJ = ntuple(x -> zeros(size(Vf, 1), size(first(geo), 2)), Dim)
for i = 1:Dim, j = 1:Dim
Expand Down

0 comments on commit c09676e

Please sign in to comment.