def test_mesh_utils():
mesh1 = utils.load_example_mesh("unit_disc")
mesh2 = utils.load_example_mesh("10x10_plane")
mesh2.vertices += np.array([[0, 0, 20]])
mesh3 = utils.combine_meshes((mesh1, mesh2))
assert len(mesh3.vertices) == len(mesh1.vertices) + len(mesh2.vertices)
planemesh = load_example_mesh("10x10_plane_hires")
# Specify coil plane geometry
center_offset = np.array([0, 0, 0])
standoff = np.array([0, 3, 0])
# Create coil plane pairs
coil_plus = trimesh.Trimesh(planemesh.vertices + center_offset + standoff,
planemesh.faces,
process=False)
coil_minus = trimesh.Trimesh(planemesh.vertices + center_offset - standoff,
planemesh.faces,
process=False)
joined_planes = combine_meshes((coil_plus, coil_minus))
# Create mesh class object
coil = MeshConductor(mesh_obj=joined_planes,
fix_normals=True,
basis_name="suh",
N_suh=100)
#%%
# Set up target and stray field points
# Here, the target points are on a volumetric grid within a sphere
center = np.array([0, 0, 0])
sidelength = 1.5
planemesh = load_example_mesh("10x10_plane_hires")
# Specify coil plane geometry
center_offset = np.array([0, 0, 0])
standoff = np.array([0, 4, 0])
# Create coil plane pairs
coil_plus = trimesh.Trimesh(planemesh.vertices + center_offset + standoff,
planemesh.faces,
process=False)
coil_minus = trimesh.Trimesh(planemesh.vertices + center_offset - standoff,
planemesh.faces,
process=False)
mesh1 = combine_meshes((coil_plus, coil_minus))
mesh2 = mesh1.copy()
mesh2.apply_scale(1.4)
coil = MeshConductor(mesh_obj=mesh1, basis_name="inner", N_sph=4)
shieldcoil = MeshConductor(mesh_obj=mesh2, basis_name="inner", N_sph=4)
#%%
# Plot geometry
f = mlab.figure(None,
bgcolor=(1, 1, 1),
fgcolor=(0.5, 0.5, 0.5),
size=(800, 800))
coil.plot_mesh(opacity=0.2, figure=f)
shieldcoil.plot_mesh(opacity=0.2, figure=f)
