Function GenerateVolumeData

comp/webgui.cpp:278–318

}

vector<string> GenerateVolumeData( shared_ptr<MeshAccess> ma,
                                    shared_ptr<CoefficientFunction> cf,
                                    int order) {
  LocalHeapMem<100000> lh("webgui_volume_data");
  Array<double> volume_data;
  auto comps = cf->Dimension();
  for (auto elnr : Range(ma->GetNElements(3))) {
    HeapReset hr(lh);

    auto el = ma->GetElement({VOL, 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);
    volume_data.Append(
        FlatArray<double>(values.AsVector().Size(), values.Data()));
  }
  auto nip = GetElementPoints(ET_TET, order).Size();
  auto nel = volume_data.Size() / nip /
             comps;  // might be different from ma->GetNElements(3), because
                     // hexes/prisms are divided into tets
  FlatTensor<3, double> data(nel, nip, comps, volume_data.Data());
  Tensor<3, double> output(nip, nel, comps);

  for (size_t i = 0; i < nel; i++)
    output(STAR, i, STAR) = data(i, STAR, STAR);

  Tensor<3, float> output_f(nip, nel, comps);
  auto n = output.GetTotalSize();
  for (auto i : Range(n)) output_f.Data()[i] = output.Data()[i];

  std::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;
}

unique_ptr<WebguiData> GenerateWebguiData(shared_ptr<MeshAccess> ma,
                                          shared_ptr<CoefficientFunction> cf,