On May 15, 2018, at 10:13 AM, John Peterson <jwpeter...@gmail.com<mailto:jwpeter...@gmail.com>> wrote:
On Tue, May 15, 2018 at 5:37 AM, Griffith, Boyce Eugene <boy...@email.unc.edu<mailto:boy...@email.unc.edu>> wrote: Folks -- Is there an easy way in the library to get normal vectors on surface meshes (e.g. 2D elements in 3D space, 1D elements in 2D space, etc)? It seems like most of the support for this assumes that these are only needed for doing surface integrals on volumetric meshes, but I am sure I am overlooking something. The "simplest" way is to build an FE object and only pre-request the normals by calling: const std::vector<Point> & face_normals = fe_face->get_normals(); I may be making a blunder here, but this doesn't seem to work if the dimension of the quadrature rule is equal to the dimension of the mesh --- the normals do not appear to be computed. (an example can be seen in adaptivity_ex4). The normals could of course be computed using only geometric information, but this is not currently implemented. -- John ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users