From 8499dc63b98262e4614a2326f108045da6b2db39 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Sun, 9 Jul 2023 22:55:52 -0600 Subject: [PATCH] Minor doc update in manifold.cc. --- source/grid/manifold.cc | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/source/grid/manifold.cc b/source/grid/manifold.cc index 4fb7e3fb13..1be7de257a 100644 --- a/source/grid/manifold.cc +++ b/source/grid/manifold.cc @@ -922,14 +922,17 @@ FlatManifold::normal_vector( // F(xi) = sum_v vertex[v] phi_v(xi) // We get the shape functions phi_v from an object of type FE_Q(1) - // we start with the point xi=1/2, xi=(1/2,1/2), ... + // We start at the center of the cell. If the face is a line or + // square, then the center is at 0.5 or (0.5,0.5). If the face is + // a triangle, then we start at the point (1/3,1/3). const unsigned int facedim = dim - 1; Point xi; const auto face_reference_cell = face->reference_cell(); - if (face_reference_cell == ReferenceCells::get_hypercube()) + if ((dim <= 2) || + (face_reference_cell == ReferenceCells::get_hypercube())) { for (unsigned int i = 0; i < facedim; ++i) xi[i] = 1. / 2; -- 2.39.5