Function MapBernstein

comp/webgui.cpp:174–196

}

vector<string> MapBernstein(FlatTensor<3, double> input, ngfem::ELEMENT_TYPE eltype,
                            int order, int n_components) {
  auto [ne, nip, comps] = input.Shape();
  Tensor<3, double> output(nip, ne, comps);
  Tensor<3, float> output_f(nip, ne, comps);
  auto trafo = GetIBernsteinBasis(eltype, order);

  for (size_t i = 0; i < ne; i++)
    output(STAR, i, STAR) = trafo * input(i, STAR, STAR);

  // convert double to float
  auto n = output.GetTotalSize();
  for (auto i : Range(n)) output_f.Data()[i] = output.Data()[i];

  vector<string> ret;
  for (auto i : Range(nip)) {
    auto vals = output_f(i, STAR, STAR);
    ret.push_back(base64_encode(FlatArray<unsigned char>(
        vals.Height() * vals.Width() * sizeof(vals(0, 0)),
        reinterpret_cast<unsigned char *>(vals.Data()))));
  }
  return ret;
}

vector<string> GenerateWireframeData( shared_ptr<MeshAccess> ma,
                                      shared_ptr<CoefficientFunction> cf,