Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
def test_cylinder():
radius = 1.0
z0 = 0.0
z1 = 1.0
edge_length = 0.1
s0 = pygalmesh.Cylinder(z0, z1, radius, edge_length)
mesh = pygalmesh.generate_mesh(
s0, cell_size=0.1, edge_size=edge_length, verbose=False
)
tol = 1.0e-1
assert abs(max(mesh.points[:, 0]) - radius) < tol
assert abs(min(mesh.points[:, 0]) + radius) < tol
assert abs(max(mesh.points[:, 1]) - radius) < tol
assert abs(min(mesh.points[:, 1]) + radius) < tol
assert abs(max(mesh.points[:, 2]) - z1) < tol
assert abs(min(mesh.points[:, 2]) + z0) < tol
vol = sum(helpers.compute_volumes(mesh.points, mesh.cells["tetra"]))
ref_vol = numpy.pi * radius * radius * (z1 - z0)
assert abs(vol - ref_vol) < tol
return