Mesh.save subdomains and boundaries
See original GitHub issueShould Mesh.save
save the .subdomains
and .boundaries
attributes? In a way that they might be reread by Mesh.load
.
If I/O is to continue to rely on meshio.write
and meshio.read
(as it should), I think that saving boundaries requires saving cells of type bnd_type
(e.g. 'line'
for MeshQuad
); currently skfem.Mesh.save
only saves cells of type Mesh.meshio_type
.
(I’m currently parsing a polyMesh
from OpenFOAM and constructing an skfem.Mesh
. My parser is very slow so I want to be able to save the mesh while experimenting with the solver.)
Perhaps instead for the moment I should just construct a meshio.Mesh
and meshio.write
that, leaving scikit-fem out of it, but this might be something to think about for other cases, e.g. when the mesh has been refined or adapted inside scikit-fem.
Issue Analytics
- State:
- Created 4 years ago
- Comments:44 (44 by maintainers)
Top GitHub Comments
I’m planning to use products of prime numbers to distinguish the overlaps.