Function TetraDubinerP

src/Numerical/Functions.cpp:128–140

}

double TetraDubinerP(const std::array<unsigned, 3>& i, const std::array<double, 3>& xi) {
  const double rNum = 2.0 * xi[0] - 1.0 + xi[1] + xi[2];
  const double sNum = 2.0 * xi[1] - 1.0 + xi[2];
  const double t = 2.0 * xi[2] - 1.0;
  const double sigmatheta = 1.0 - xi[1] - xi[2];
  const double theta = 1.0 - xi[2];

  const double ti = SingularityFreeJacobiP(i[0], 0, 0, rNum, sigmatheta);
  const double tij = SingularityFreeJacobiP(i[1], 2 * i[0] + 1, 0, sNum, theta);
  const double tijk = SingularityFreeJacobiP(i[2], 2 * i[0] + 2 * i[1] + 2, 0, t, 1.0);

  return ti * tij * tijk;
}

std::array<double, 3> gradTetraDubinerP(const std::array<unsigned, 3>& i,
                                        const std::array<double, 3>& xi) {