From c09676eb50951547fc511de71c70caf2d57b2c51 Mon Sep 17 00:00:00 2001 From: Jesse Chan Date: Mon, 30 Sep 2024 09:24:57 -0500 Subject: [PATCH] add comment on Nanson's formula for normals --- src/geometric_functions.jl | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/geometric_functions.jl b/src/geometric_functions.jl index 00e8a251..4b238af1 100644 --- a/src/geometric_functions.jl +++ b/src/geometric_functions.jl @@ -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