def test_radius_ratio_triangle():
# Create mesh and compute rations
mesh = UnitSquareMesh(MPI.comm_world, 12, 12)
ratios = MeshQuality.radius_ratios(mesh)
for c in Cells(mesh):
assert round(ratios[c] - 0.828427124746, 7) == 0
def test_radius_ratio_tetrahedron():
# Create mesh and compute ratios
mesh = UnitCubeMesh(MPI.comm_world, 14, 14, 14)
ratios = MeshQuality.radius_ratios(mesh)
for c in Cells(mesh):
assert round(ratios[c] - 0.717438935214, 7) == 0
def find_low(mesh, V, Ause):
V2dm = V.dofmap()
cq2 = MeshQuality.radius_ratios(mesh)
cq = cq2.array()
indices = np.where(cq < 0.1)[0]
dof_set = []
cell_set = []
for i in indices:
cell = Cell(mesh, i)
for v in vertices(cell):
for c in cells(v):
cell_set += [c.index()]
dof_set.extend(V2dm.cell_dofs(c.index()))
bad_set = list(set(dof_set))
bad_cells = list(set(cell_set))
# print('BAD CELLS=', bad_cells)
# print(len(bad_cells))
# print('BAD DOFS=', bad_dofs)
# print(len(bad_dofs))
# check redundancy
re_dofs = []
for d1 in bad_set:
for d2 in bad_set:
if Ause[d1, d2] != 0:
re_dofs.append(d2)
bad_dofs = list(set(re_dofs))
return bad_dofs
