Function BernsteinTrig

comp/webgui.cpp:42–45

}

double BernsteinTrig(double x, double y, int i, int j, int n) {
  return factorial(n) / factorial(i) / factorial(j) / factorial(n - i - j) *
         pow(x, i) * pow(y, j) * pow(1 - x - y, n - i - j);
}

Matrix<double> GetIBernsteinBasis(ngfem::ELEMENT_TYPE etype, int order) {
  std::map<std::pair<ngfem::ELEMENT_TYPE, int>, Matrix<double>> cache;