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 = av_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 = av_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 = av_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)