Function ExportReadAnsys

comp/read_ansys.cpp:1061–1365

using namespace ngcomp;

void ExportReadAnsys (py::module & m)
{
    m.def("ReadAnsysInp", [](const string & filename, double scale) -> py::dict
    {
      AnsysModel model = ParseAbaqusInp(filename);

      map<string,vector<int>> vol_regions, bnd_regions, bbnd_regions, mat_cells, element_sets;
      vector<int> id2vertex;
      auto ngmesh = BuildAnsysMesh(model, scale, vol_regions, bnd_regions,
                                   bbnd_regions, mat_cells, element_sets, id2vertex);

      auto ma = make_shared<MeshAccess>(ngmesh);
      ngmesh->updateSignal.Connect(ma.get(), [p=ma.get()] { p->UpdateBuffers(); });

      const int nvol  = ma->GetNRegions(VOL);
      const int nbnd  = ma->GetNRegions(BND);
      const int nbbnd = ma->GetNRegions(BBND);

      auto MakeRegion = [&](VorB vb, const vector<int> & nrs, int nregions)
      {
        BitArray ba(nregions);
        ba.Clear();
        for (int r : nrs) if (r >= 0 && r < nregions) ba.SetBit(r);
        return Region(ma, vb, ba);
      };

      py::dict res;
      res["mesh"] = py::cast(ma);

      // material property data
      py::dict mats;
      for (auto & [name, mat] : model.materials)
        {
          py::dict d;
          if (mat.hasE)   d["E"]   = mat.E;
          if (mat.hasNu)  d["nu"]  = mat.nu;
          if (mat.hasRho) d["rho"] = mat.rho;
          mats[py::str(name)] = d;
        }
      res["materials"] = mats;

      // "true" domains: one Region(VOL) per body/part (overlap-safe)
      py::dict domains;
      for (auto & [name, nrs] : vol_regions)
        domains[py::str(name)] = py::cast(MakeRegion(VOL, nrs, nvol));
      res["domains"] = domains;

      py::dict material_regions;
      for (auto & [name, nrs] : mat_cells)
        material_regions[py::str(name)] = py::cast(MakeRegion(VOL, nrs, nvol));
      res["material_regions"] = material_regions;

      // standalone *ELSET named selections (not bodies): element-number lists,
      // kept out of the domain partition so it stays clean
      py::dict elsets;
      for (auto & [name, els] : element_sets)
        {
          py::list l;