def test_connectivity_loop(): vertices = np.array([(0.0, 1.0), (1.0, 0.0)]) element_to_vertex = np.array([[0, 1], [1, 0]]) m = from_vertices_and_etov(vertices, element_to_vertex) assert(m.is_neighbor(0, 1, 'left')) assert(m.is_neighbor(0, 1, 'right')) assert(m.is_neighbor(1, 0, 'left')) assert(m.is_neighbor(1, 0, 'right'))
def test_element_widths(): vertices = np.array([(0.0, 0.0), (1.0, 0.0), (3.0, 0.0)]) etov = np.array([(0, 1), (1, 2)]) m = from_vertices_and_etov(vertices, etov) assert(m.elements[0].length == 1.0) assert(m.elements[1].length == 2.0)
def test_misorientation_correct_harder(): v = np.array([[0, 1], [0, 2], [0, 3], [0, 4]]) e = np.array([[0, 1],[2, 1],[2, 3]]) mesh = from_vertices_and_etov(v, e, flip = True) mesh.check_misorientation()