dimtag -> dataset mapping

[chat-nti-fm]

Mechanism that could solve many related issues:

• solve issues possibly treated by functional approach bgem.gmsh Functional approach #62

• similar to Field application #64, but more detailed design

• region setting issue ObjectSet issues #51

• will resolve Region marking for resulting mesh. #34

• Mapa dimtag -> dataset_ID, propaguje se během bool operací

• Snad možné časem implmentovat mapování i pro boundary objekty,
  to by řešilo problémy s nastavování BC apod. Nutno z bool operací vymáčkonout mapování pro
  všechny zasažené dimtagy všech nižších dimenzí (tj. Explicitně do operace přidat hranice)

• Při vytvoření objektu, možno nastavit jeho dataset;
  použít pro fractures, nastavit ID pukliny; jinak default dataset

• Pro existující bgem objekt (výsledek bol operací) možno přidat/editovat datasety jeho dimtagů:

  1. Identifikace všech dataset ID pro BGEM ShapeSet
  2. Pro každý dataset ID kopie + merge s novým datasetem
  3. Aplikace map starých na nové datasety do globální mapy dimtagů -> datasety

• Po síťování, konstrukce mapy element -> shape -> datasetID -> dataset
  Our mesh wrapper class, extend one that exists, but need better internal structure and reader,
  possibly based on meshio

• Extraktor(dataset_var, regions=ALL): dict: element -> el_dataset[‘dataset_var’]

Možný poblém: mnoho operací s Python dict, pomalé. Solve later.
Místo dict, použít postupně rostoucí numpy arrays mapující [0: max shape ID] na dataset.

Allows:

• Postpone region mapping completely after meshing (kepp region var in all datasets
• Set mesh step through datasets, again handled more transprently.
• Propagate geometry related data to the mesh.
• Write datasets into separate file, allowing passing their info to the solver possibly through a field