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_torus():
major_radius = 1.0
minor_radius = 0.5
s0 = pygalmesh.Torus(major_radius, minor_radius)
mesh = pygalmesh.generate_mesh(s0, cell_size=0.1, verbose=False)
tol = 1.0e-2
radii_sum = major_radius + minor_radius
assert abs(max(mesh.points[:, 0]) - radii_sum) < tol
assert abs(min(mesh.points[:, 0]) + radii_sum) < tol
assert abs(max(mesh.points[:, 1]) - radii_sum) < tol
assert abs(min(mesh.points[:, 1]) + radii_sum) < tol
assert abs(max(mesh.points[:, 2]) - minor_radius) < tol
assert abs(min(mesh.points[:, 2]) + minor_radius) < tol
vol = sum(helpers.compute_volumes(mesh.points, mesh.cells["tetra"]))
ref_vol = (numpy.pi * minor_radius * minor_radius) * (2 * numpy.pi * major_radius)
assert abs(vol - ref_vol) < 1.0e-1
return