On Wed, Apr 21, 2021 at 10:39 AM Renato Poli <rebp...@gmail.com> wrote:
> Hello, > > I am willing to make a coordinate rotation on a Tet face to calculate > normal and tangent displacements between discontinuous elements. > > As I understand, I need a transformation matrix gathering normal and > tangents to define a local coordinate system: > T_ij = [ nx_i ns_i nt_i ] > > I can see all the calculations are there, but I am not sure what is the > most elegant/usual way to get such matrix in libmesh. > > Any suggestion? Any documentation pointer would be helpful! > See FEAbstract::get_normals() and get_tangents() in include/fe/fe_abstract.h. These call the corresponding functions on FEMap. Similarly to the other "get_" interfaces on FE, the values are computed at quadrature points, so you would need to reinit a "side" FE after pre-requesting these values, similarly to what is done in systems_of_equations_ex4, for example. -- John _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
