LutzGross / esys-escript.github.io

Other
29 stars 13 forks source link

finley: Dirac Delta Functions and GMSH tags #38

Open aellery opened 2 years ago

aellery commented 2 years ago

It is possible to tag points in GMSH.

ReadGmsh("meshDirac.msh", 3, optimize=True) does read the tags and their locations.

It is not possible to use a DiracDeltaFunction unless ReadGmsh contains the dirac points and dirac tags explicitly.

domain2=ReadGmsh("meshDirac.msh", 3, optimize=True, diracPoints = electrodes, diracTags = diracTags)

Attached is a small gmsh file with tags for electrodes and small python code that uses the mesh.

An optimal solution wold allow Dirac tags and points read from the .msh that does not need the tags and points to be reentered in the script.

files.zip