From: Wolfgang Bangerth Date: Mon, 10 Jul 2023 04:55:52 +0000 (-0600) Subject: Minor doc update in manifold.cc. X-Git-Tag: relicensing~711^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=8499dc63b98262e4614a2326f108045da6b2db39;p=dealii.git Minor doc update in manifold.cc. --- 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;