def test_zonotope_0(self): data = helper.load_test_npz('test_spec_decomp_1.npz') df0 = data['df'] sub = ss.Subspaces() sub.compute(df0) sub.partition(2) np.random.seed(42) bavd = dom.BoundedActiveVariableDomain(sub) vertices = bavd.vertY plt.zonotope_2d_plot(vertices)
def test_zonotope_1(self): data = helper.load_test_npz('test_spec_decomp_1.npz') df0 = data['df'] sub = ss.Subspaces() sub.compute(df0) sub.partition(2) np.random.seed(42) bavd = dom.BoundedActiveVariableDomain(sub) bavm = dom.BoundedActiveVariableMap(bavd) Y = as_design(bavm, 8, NMC=1)[0] vertices = bavd.vertY plt.zonotope_2d_plot(vertices, design=Y)
def test_zonotope_3(self): data = helper.load_test_npz('test_spec_decomp_1.npz') df0 = data['df'] sub = ss.Subspaces() sub.compute(df0) sub.partition(2) np.random.seed(42) bavd = dom.BoundedActiveVariableDomain(sub) bavm = dom.BoundedActiveVariableMap(bavd) Y = as_design(bavm, 8, NMC=1)[0] vertices = bavd.vertY Yp, Yw = av_quadrature_rule(bavm, 8) plt.zonotope_2d_plot(vertices, design=Y, y=Yp, f=Yw, out_label='quadrature rule')
def test_zonotope_2(self): data = helper.load_test_npz('test_spec_decomp_1.npz') df0 = data['df'] sub = ss.Subspaces() sub.compute(df0) sub.partition(2) np.random.seed(42) bavd = dom.BoundedActiveVariableDomain(sub) bavm = dom.BoundedActiveVariableMap(bavd) Y = as_design(bavm, 8, NMC=1)[0] vertices = bavd.vertY Xp = np.random.uniform(-1.0,1.0,size=(20, sub.W1.shape[0])) Yp = np.dot(Xp, sub.W1) fp = np.sum(Yp, axis=1) plt.zonotope_2d_plot(vertices, design=Y, y=Yp, f=fp)