Dear all, I need the outward normal vector on a face (side) of the element to apply the boundary condition.
Is fe_face->get_normals() (always) outward or inward? The code is like: /// dim = 3; AutoPtr<QBase> qrule_face = QBase::build (GAUSS, dim -1, 2) AutoPtr<FEBase> fe_face(FEBase::build(dim, fe_type)); fe_face->attach_quadrature_rule(qrule_face.get()); const std::vector<libMesh::Point>& normal_face = fe_face->get_normals(); /// Thanks, Walter ------------------------------------------------------------------------------ Site24x7 APM Insight: Get Deep Visibility into Application Performance APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month Monitor end-to-end web transactions and take corrective actions now Troubleshoot faster and improve end-user experience. Signup Now! http://pubads.g.doubleclick.net/gampad/clk?id=272487151&iu=/4140 _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users