def test_cuboids_union(): c0 = pygalmesh.Cuboid([0, 0, -0.5], [3, 3, 0.5]) c1 = pygalmesh.Cuboid([1, 1, -2], [2, 2, 2]) u = pygalmesh.Union([c0, c1]) pygalmesh.generate_mesh(u, 'out.mesh', cell_size=0.2, edge_size=0.2, verbose=False) vertices, cells, _, _, _ = meshio.read('out.mesh') # filter the vertices that belong to cells verts = vertices[numpy.unique(cells['tetra'])] tol = 1.0e-2 assert abs(max(verts[:, 0]) - 3.0) < tol assert abs(min(verts[:, 0]) - 0.0) < tol assert abs(max(verts[:, 1]) - 3.0) < tol assert abs(min(verts[:, 1]) - 0.0) < tol assert abs(max(verts[:, 2]) - 2.0) < tol assert abs(min(verts[:, 2]) + 2.0) < tol vol = sum(compute_volumes(vertices, cells['tetra'])) assert abs(vol - 12.0) < 0.1 return
def test_balls_union(): radius = 1.0 displacement = 0.5 s0 = pygalmesh.Ball([displacement, 0, 0], radius) s1 = pygalmesh.Ball([-displacement, 0, 0], radius) u = pygalmesh.Union([s0, s1]) a = numpy.sqrt(radius**2 - displacement**2) edge_size = 0.1 n = int(2 * numpy.pi * a / edge_size) circ = [[ 0.0, a * numpy.cos(i * 2 * numpy.pi / n), a * numpy.sin(i * 2 * numpy.pi / n) ] for i in range(n)] circ.append(circ[0]) pygalmesh.generate_mesh( u, "out.mesh", feature_edges=[circ], cell_size=0.15, edge_size=edge_size, verbose=False, ) mesh = meshio.read("out.mesh") assert abs(max(mesh.points[:, 0]) - (radius + displacement)) < 0.02 assert abs(min(mesh.points[:, 0]) + (radius + displacement)) < 0.02 assert abs(max(mesh.points[:, 1]) - radius) < 0.02 assert abs(min(mesh.points[:, 1]) + radius) < 0.02 assert abs(max(mesh.points[:, 2]) - radius) < 0.02 assert abs(min(mesh.points[:, 2]) + radius) < 0.02 vol = sum(compute_volumes(mesh.points, mesh.cells["tetra"])) h = radius - displacement ref_vol = 2 * (4.0 / 3.0 * numpy.pi * radius**3 - h * numpy.pi / 6.0 * (3 * a**2 + h**2)) assert abs(vol - ref_vol) < 0.1 return
def test_cuboids_union(): c0 = pygalmesh.Cuboid([0, 0, -0.5], [3, 3, 0.5]) c1 = pygalmesh.Cuboid([1, 1, -2], [2, 2, 2]) u = pygalmesh.Union([c0, c1]) mesh = pygalmesh.generate_mesh(u, cell_size=0.2, edge_size=0.2, verbose=False) # filter the vertices that belong to cells verts = mesh.points[numpy.unique(mesh.cells["tetra"])] tol = 1.0e-2 assert abs(max(verts[:, 0]) - 3.0) < tol assert abs(min(verts[:, 0]) - 0.0) < tol assert abs(max(verts[:, 1]) - 3.0) < tol assert abs(min(verts[:, 1]) - 0.0) < tol assert abs(max(verts[:, 2]) - 2.0) < tol assert abs(min(verts[:, 2]) + 2.0) < tol vol = sum(helpers.compute_volumes(mesh.points, mesh.cells["tetra"])) assert abs(vol - 12.0) < 0.1 return