Hdiv.cpp File Reference

Implementation of H-curl base. More...

Go to the source code of this file.

Functions
static double	get_ksi (const double n0, const double n1)
static FTensor::Tensor1< double, 3 >	get_diff_ksi (const double *diffN, const int n0, const int n1)

Detailed Description

Implementation of H-curl base.

Based on Hierarchic Finite Element Bases on Unstructured Tetrahedral Meshes, by Mark Ainsworth and Joe Coyle Shape functions for MBTRI/MBTET and HCurl space

Definition in file Hdiv.cpp.

Function Documentation

get_diff_ksi()

static FTensor::Tensor1< double, 3 > get_diff_ksi ( const double *  diffN, const int  n0, const int  n1  )

inlinestatic

Definition at line 70 of file Hdiv.cpp.

                                                     {
  static_assert(
      AinsworthHdivOptions::edge_coordinate ==
              AinsworthHdivOptions::EdgeCoordinate::PAPER ||
          AinsworthHdivOptions::edge_coordinate ==
              AinsworthHdivOptions::EdgeCoordinate::ENDPOINT_O ||
          AinsworthHdivOptions::edge_coordinate ==
              AinsworthHdivOptions::EdgeCoordinate::ENDPOINT_I,
      "Unsupported Ainsworth H(div) edge coordinate");
 
  FTensor::Tensor1<double, 3> t_diff_ksi;
  for (int dd = 0; dd != 3; ++dd) {
    if constexpr (AinsworthHdivOptions::edge_coordinate ==
                  AinsworthHdivOptions::EdgeCoordinate::PAPER) {
      t_diff_ksi(dd) = diffN[3 * n1 + dd] - diffN[3 * n0 + dd];
    } else if constexpr (AinsworthHdivOptions::edge_coordinate ==
                         AinsworthHdivOptions::EdgeCoordinate::ENDPOINT_O) {
      t_diff_ksi(dd) = -2 * diffN[3 * n0 + dd];
    } else if constexpr (AinsworthHdivOptions::edge_coordinate ==
                         AinsworthHdivOptions::EdgeCoordinate::ENDPOINT_I) {
      t_diff_ksi(dd) = 2 * diffN[3 * n1 + dd];
    }
  }
  return t_diff_ksi;
}

get_ksi()

static double get_ksi ( const double  n0, const double  n1  )

inlinestatic

Definition at line 47 of file Hdiv.cpp.

                                                               {
  static_assert(
      AinsworthHdivOptions::edge_coordinate ==
              AinsworthHdivOptions::EdgeCoordinate::PAPER ||
          AinsworthHdivOptions::edge_coordinate ==
              AinsworthHdivOptions::EdgeCoordinate::ENDPOINT_O ||
          AinsworthHdivOptions::edge_coordinate ==
              AinsworthHdivOptions::EdgeCoordinate::ENDPOINT_I,
      "Unsupported Ainsworth H(div) edge coordinate");
 
  if constexpr (AinsworthHdivOptions::edge_coordinate ==
                AinsworthHdivOptions::EdgeCoordinate::PAPER) {
    return n1 - n0;
  } else if constexpr (AinsworthHdivOptions::edge_coordinate ==
                       AinsworthHdivOptions::EdgeCoordinate::ENDPOINT_O) {
    return 1 - 2 * n0;
  } else if constexpr (AinsworthHdivOptions::edge_coordinate ==
                       AinsworthHdivOptions::EdgeCoordinate::ENDPOINT_I) {
    return 2 * n1 - 1;
  }
  return 0;
}