mirror of
https://github.com/NGSolve/netgen.git
synced 2024-12-26 13:50:33 +05:00
23 lines
629 B
Python
23 lines
629 B
Python
from netgen.meshing import *
|
|
|
|
def ReadViaMeshIO(filename):
|
|
import meshio
|
|
import numpy as np
|
|
|
|
# print ("reading via meshio:", filename)
|
|
|
|
m = meshio.read(filename)
|
|
pts = m.points
|
|
|
|
mesh = Mesh(dim=pts.shape[1])
|
|
# mesh.AddPoints ( np.asarray(pts, dtype=np.float64) ) # needed for correct little/big endian
|
|
mesh.AddPoints ( pts )
|
|
|
|
fd = mesh.Add (FaceDescriptor(bc=1))
|
|
for cb in m.cells:
|
|
# mesh.AddElements(dim=cb.dim, index=1, data=np.asarray(cb.data, dtype=np.int32), base=0)
|
|
mesh.AddElements(dim=cb.dim, index=1, data=cb.data, base=0)
|
|
|
|
return mesh
|
|
|