def test_rs_bnd_2d_int(self):
np.random.seed(42)
X0 = np.random.uniform(-1.0,1.0,size=(50,3))
f0 = np.random.normal(size=(50,1))
df0 = np.random.normal(size=(50,3))
sub = ss.Subspaces()
sub.compute(df=df0)
sub.partition(2)
avd = dom.BoundedActiveVariableDomain(sub)
avm = dom.BoundedActiveVariableMap(avd)
asrs = asm.ActiveSubspaceResponseSurface(avm)
asrs.train_with_data(X0, f0)
I = asi.av_integrate(asrs, avm, 10)
def test_rs_bnd_int(self):
data = helper.load_test_npz('test_rs_1.npz')
X, f, df = data['X'], data['f'], data['df']
sub = ss.Subspaces()
sub.compute(df)
avd = dom.BoundedActiveVariableDomain(sub)
avm = dom.BoundedActiveVariableMap(avd)
asrs = asm.ActiveSubspaceResponseSurface(avm)
asrs.train_with_data(X, f)
np.random.seed(43)
I = asi.av_integrate(asrs, avm, 10)
if self.writeData:
np.savez('data/test_int_0_2',I=I)
data_test = helper.load_test_npz('test_int_0_2.npz')
np.testing.assert_almost_equal(I, data_test['I'])
print '\n'
print 'rs bnd: {:6.4f}'.format(I)
