search

SeisSol / PUMGen

Mesh generation for SeisSol
BSD 3-Clause "New" or "Revised" License
7 stars 3 forks source link

Boundary conditions incorrectly set in bypass-apf branch #64

Closed Thomas-Ulrich closed 9 months ago

Thomas-Ulrich commented 1 year ago

I'm generating a mesh with https://github.com/SeisSol/Examples/blob/master/tpv5/generating_mirror_mesh.sh And vizualizing the BC with: python ~/SeisSol/Meshing/vizualizeBoundaryConditions/vizualizeBoundaryConditions.py tpv5_f200_halfdom.xdmf 0 I get with the bypass-apf branch image and the mesh64 branch image

davschneller commented 1 year ago

Ok, I see where the problem may lie... I forgot a logical "or" operator (i.e. each cell gets only the boundary condition of one face set)—after that we seem to get the full boundary conditions again, at least for the test cases I had let run. The branch is updated.