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_cone():
base_radius = 1.0
height = 2.0
edge_size = 0.1
s0 = pygalmesh.Cone(base_radius, height, edge_size)
mesh = pygalmesh.generate_mesh(
s0, cell_size=0.1, edge_size=edge_size, verbose=False
)
tol = 2.0e-1
assert abs(max(mesh.points[:, 0]) - base_radius) < tol
assert abs(min(mesh.points[:, 0]) + base_radius) < tol
assert abs(max(mesh.points[:, 1]) - base_radius) < tol
assert abs(min(mesh.points[:, 1]) + base_radius) < tol
assert abs(max(mesh.points[:, 2]) - height) < tol
assert abs(min(mesh.points[:, 2]) + 0.0) < tol
vol = sum(helpers.compute_volumes(mesh.points, mesh.cells["tetra"]))
ref_vol = numpy.pi * base_radius * base_radius / 3.0 * height
assert abs(vol - ref_vol) < tol
return