Function GenerateSurfaceData

comp/webgui.cpp:250–276

}

vector<string> GenerateSurfaceData( shared_ptr<MeshAccess> ma,
                                    shared_ptr<CoefficientFunction> cf,
                                    int order) {
  // generate surface data
  LocalHeapMem<100000> lh("webgui_surface_data");
  auto vb = ma->GetDimension() == 3 ? BND : VOL;
  Array<double> surface_data;
  auto comps = cf->Dimension();
  for (auto elnr : Range(ma->GetNElements(2))) {
    HeapReset hr(lh);

    auto el = ma->GetElement({vb, elnr});
    auto &ir = GetElementPoints(el.GetType(), order);
    auto &trafo = ma->GetTrafo(el, lh);
    BaseMappedIntegrationRule &mir = trafo(ir, lh);
    FlatMatrix<> values(ir.Size(), comps, lh);
    cf->Evaluate(mir, values);
    surface_data.Append(
        FlatArray<double>(values.AsVector().Size(), values.Data()));
  }
  auto nip = GetElementPoints(ET_TRIG, order).Size();
  auto nel = surface_data.Size() / nip /
             comps;  // might be different from ma->GetNElements(2), because
                     // quads are divided into 2 trigs
  FlatTensor<3, double> data(nel, nip, comps, surface_data.Data());
  return MapBernstein(data, ET_TRIG, order, cf->Dimension());
}

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