Pass element ID to Python ST/BC
Why?
Fast lookup of certain properties in the Python script.
Non-functional requirements
- existing Python scripts must not break
- which element ID to pass? Bulk or boundary mesh element ID? Both?
Alternatives
Handwritten lookup in the Python script. E.g. using VTK's vtkCellLocator
locator = vtkCellLocator()
locator.SetDataSet(bulk_mesh)
locator.BuildLocator()
...
locator.FindCell(...) # probably locator.FindCell((x,y,z), 1e-16)
See https://vtk.org/doc/nightly/html/classvtkCellLocator.html#a7400cb179b18aed1ea21e5a31ca421e1