def test_quality_struct():
x = np.arange(-10, 10, 5)
y = np.arange(-10, 10, 5)
z = np.arange(-10, 10, 5)
x, y, z = np.meshgrid(x, y, z)
grid = pv.StructuredGrid(x, y, z)
qual = pyansys.quality(grid)
assert (qual == 1).all()
cell_center = grid.points[cells].mean(1)
# extract cells below the 0 xy plane
mask = cell_center[:, 2] < 0
cell_ind = mask.nonzero()[0]
subgrid = grid.extract_cells(cell_ind)
# advanced plotting
plotter = pv.Plotter()
plotter.add_mesh(subgrid, 'lightgrey', lighting=True, show_edges=True)
plotter.add_mesh(sphere, 'r', 'wireframe')
plotter.add_legend([[' Input Mesh ', 'r'], [' Tesselated Mesh ', 'black']])
plotter.show()
###############################################################################
# Use pyansys to compute cell quality
import pyansys
cell_qual = pyansys.quality(subgrid)
# plot quality
subgrid.plot(
scalars=cell_qual,
stitle='Quality',
cmap='bwr',
clim=[0, 1],
flip_scalars=True,
show_edges=True,
)
def test_quality():
archive = pyansys.Archive(examples.hexarchivefile)
grid = archive.grid
qual = pyansys.quality(grid)
assert (qual == 1).all()
def test_quality_type_error():
with pytest.raises(TypeError):
pyansys.quality(pv.PolyData())