def test_2D_1(): # 2-D mass balance triangle A = sp.array([[-1., 0], [0, -1], [1, 1]]) b = sp.array([[0., 0, 1]]).T vs = artools.con2vert(A, b) vs_ref = sp.array([[0., 1], [0, 0], [1, 0]]) assert (artools.same_rows(vs, vs_ref) is True)
def test_1(self): # 2-D mass balance triangle A = sp.array([[-1., 0], [0, -1], [1, 1]]) b = sp.array([[0., 0, 1]]).T vs = con2vert(A, b) vs_ref = sp.array([[0., 1], [0, 0], [1, 0]]) assert (sameRows(vs, vs_ref) is True)
def test_3(self): # mass balance triangle cut at y = 0.5 A = sp.array([[-1., 0], [0, -1], [1, 1], [0, 1]]) b = sp.array([[0., 0, 1, 0.5]]).T vs = con2vert(A, b) vs_ref = sp.array([[0, 0.5], [0, 0], [1, 0], [0.5, 0.5]]) assert (sameRows(vs, vs_ref) is True)
def test_2(self): # 2-D unit square A = sp.array([[-1., 0], [0, -1], [1, 0], [0, 1]]) b = sp.array([[0., 0, 1, 1]]).T vs = con2vert(A, b) vs_ref = sp.array([[0., 1], [0, 0], [1, 0], [1, 1]]) assert (sameRows(vs, vs_ref) is True)
def test_2D_3(): # mass blanace triangle cut at y = 0.5 A = sp.array([[-1., 0], [0, -1], [1, 1], [0, 1]]) b = sp.array([[0., 0, 1, 0.5]]).T vs = artools.con2vert(A, b) vs_ref = sp.array([[0, 0.5], [0, 0], [1, 0], [0.5, 0.5]]) assert (artools.same_rows(vs, vs_ref) is True)
def test_2D_2(): # 2-D unit square A = sp.array([[-1., 0], [0, -1], [1, 0], [0, 1]]) b = sp.array([[0., 0, 1, 1]]).T vs = artools.con2vert(A, b) vs_ref = sp.array([[0., 1], [0, 0], [1, 0], [1, 1]]) assert (artools.same_rows(vs, vs_ref) is True)
def test_interior_1(self): """ This set of constraints produces an initial guess that is not an interior point of the region. """ A = sp.array([[1.0, 0.0], [1.0, 1.0], [-1.0, 1.0], [-3.0, -1.0], [0.0, -1.0]]) b = sp.array([[1.0, 1.0, 0.0, 0.0, 0.0]]).T vs_ref = sp.array([[0.5, 0.5], [0.0, 0.0], [1.0, 0.0]]) vs = con2vert(A, b) assert (sameRows(vs, vs_ref) is True)