Thanks John, What would be the straightforward way to append normals and tangents in a single tensor? The idea is to define them as a base so I can operate and iterate.
rgds, Renato On Thu, Apr 22, 2021 at 10:21 AM John Peterson <jwpeter...@gmail.com> wrote: > > > 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
