hldo-sphere/rot__whitney__one__div__matrix__provider_8cc_source.html

#include "rot_whitney_one_div_matrix_provider.h"


#include <lf/uscalfe/lagr_fe.h>


#include <cmath>

#include <iomanip>

#include <iostream>


namespace projects::hldo_sphere::assemble {


Eigen::MatrixXd RotWhitneyOneDivMatrixProvider::Eval(

    const lf::mesh::Entity &entity) const {

  // Only triangles are supported

  LF_VERIFY_MSG(entity.RefEl() == lf::base::RefEl::kTria(),

                "Unsupported cell type " << entity.RefEl());


  // correct for orientation

  auto edgeOrientations = entity.RelativeOrientations();

  double s0 = lf::mesh::to_sign(edgeOrientations[0]);

  double s1 = lf::mesh::to_sign(edgeOrientations[1]);

  double s2 = lf::mesh::to_sign(edgeOrientations[2]);


  // Compute the element matrix for the product of baricentric functions

  Eigen::MatrixXd elem_mat(3, 1);

  elem_mat << -s0, -s1, -s2;


  return elem_mat;

}


}  // namespace projects::hldo_sphere::assemble

lf::base::RefEl::kTria
static constexpr RefEl kTria()
Returns the reference triangle.
Definition: ref_el.h:158

lf::mesh::Entity
Interface class representing a topological entity in a cellular complex
Definition: entity.h:39

lf::mesh::Entity::RelativeOrientations
virtual nonstd::span< const Orientation > RelativeOrientations() const =0
return span of relative orientations of sub-entities of the next higher co-dimension.

lf::mesh::Entity::RefEl
virtual base::RefEl RefEl() const =0
Describes the reference element type of this entity.

projects::hldo_sphere::assemble::RotWhitneyOneDivMatrixProvider::Eval
Eigen::MatrixXd Eval(const lf::mesh::Entity &entity) const
Compute the element matrix for some triangle of a mesh.
Definition: rot_whitney_one_div_matrix_provider.cc:11

lf::mesh::to_sign
int to_sign(Orientation o)
Definition: entity.cc:7

projects::hldo_sphere::assemble
Collection of matrix and vector element providers.
Definition: assemble.h:26