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

Reply via email to