Skip to content

A mesh utility to get element rotation matrices

wenqing requested to merge wenqing/ogs:get_rotation_matrix into master

For all elements of a mesh, the presented function firstly checks whether there are inclined elements. If there is no inclined element at all in the mesh, it returns an empty vector. Otherwise, the rotation matrix, which is used to map the local vector to the global system, is computed for each inclined element and stored to the corresponding entry of a vector for rotation matrix.

The function returns a vector of element rotation matrix pointers, the usage of which can be found in !3632 (merged).

The size of the rotation matrix of an inclined element is n_s \times n_e, with n_s the dimension of the space occupied by the mesh, and n_e the element dimension.

  1. Feature description was added to the changelog
  2. Tests covering your feature were added? Yes, with unit tests
  3. Any new feature or behavior change was documented? Yes.
Edited by Dmitry Yu. Naumov

Merge request reports