Closed krober10nd closed 3 years ago
import meshio import SeismicMesh as sm cube1 = sm.Cube((-0.5, 0.5, -0.5, 0.5, -0.5, 0.5)) cube2 = sm.Cube((-0.25, 0.75, -0.25, 0.75, -0.75, 0.75)) domain = sm.Union([cube1, cube2],smoothness=0.25) # domain.show(samples=1e5) points, cells = sm.generate_mesh(domain=domain, edge_length=0.10) meshio.write_points_cells("testing.vtk", points, [("tetra", cells)])