def test_plot_eigenvalues_03(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 200).reshape(50, 4)
weights = np.ones((50, 1)) / 50
ss = ActiveSubspaces()
ss.compute(gradients=gradients, weights=weights, nboot=200)
with assert_plot_figures_added():
ss.plot_eigenvalues(n_evals=3, figsize=(7, 7), title='Eigenvalues')
def test_compute_05(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 60).reshape(15, 4)
outputs = np.random.uniform(0, 5, 15)
ss = ActiveSubspaces()
ss.compute(inputs=inputs, outputs=outputs, method='local', nboot=150)
true_evals = np.array([13.794711, 11.102377, 3.467318, 1.116324])
np.testing.assert_array_almost_equal(true_evals, ss.evals)
def test_compute_bootstrap_ranges_01(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 60).reshape(30, 2)
weights = np.ones((30, 1)) / 30
ss = ActiveSubspaces()
ss.compute(gradients=gradients, weights=weights, nboot=100)
true_bounds_evals = np.array([[0.3000497, 0.59008536],
[0.17398718, 0.40959827]])
np.testing.assert_array_almost_equal(true_bounds_evals, ss.evals_br)
def test_plot_sufficient_summary_02(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 200).reshape(50, 4)
weights = np.ones((50, 1)) / 50
ss = ActiveSubspaces()
ss.compute(gradients=gradients, weights=weights, nboot=200)
ss.partition(3)
with self.assertRaises(ValueError):
ss.plot_sufficient_summary(10, 10)
def test_compute_bootstrap_ranges_03(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 180).reshape(30, 3, 2)
weights = np.ones((30, 1)) / 30
ss = ActiveSubspaces()
ss.compute(gradients=gradients, weights=weights, nboot=100)
true_bounds_evals = np.array([[0.99330673, 1.62694823],
[0.65987633, 1.11751475]])
np.testing.assert_array_almost_equal(true_bounds_evals, ss.evals_br)
def test_forward_01(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 60).reshape(15, 4)
outputs = np.random.uniform(0, 5, 15)
ss = ActiveSubspaces()
ss.compute(inputs=inputs, outputs=outputs, method='local', nboot=250)
ss.partition(2)
active = ss.forward(np.random.uniform(-1, 1, 8).reshape(2, 4))[0]
true_active = np.array([[-0.004748, 0.331107], [-0.949099, 0.347534]])
np.testing.assert_array_almost_equal(true_active, active)
def test_plot_eigenvectors_03(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 200).reshape(50, 4)
weights = np.ones((50, 1)) / 50
ss = ActiveSubspaces()
ss.compute(gradients=gradients, weights=weights, nboot=200)
with assert_plot_figures_added():
ss.plot_eigenvectors(n_evects=2,
figsize=(5, 8),
labels=[r'$x$', r'$y$', r'$r$', r'$z$'])
def test_forward_05(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 60).reshape(15, 4)
outputs = np.random.uniform(0, 5, 15)
ss = ActiveSubspaces()
ss.compute(inputs=inputs, outputs=outputs, method='local', nboot=250)
ss.partition(2)
inactive = ss.forward(np.random.uniform(-1, 1, 8).reshape(2, 4))[1]
true_inactive = np.array([[-1.03574242, -0.04662904],
[-0.49850367, -0.37146678]])
np.testing.assert_array_almost_equal(true_inactive, inactive)
def test_forward_04(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 60).reshape(15, 4)
outputs = np.random.uniform(0, 5, 45).reshape(15, 3)
ss = ActiveSubspaces()
ss.compute(inputs=inputs, outputs=outputs, method='local', nboot=250)
ss.partition(2)
active = ss.forward(np.random.uniform(-1, 1, 8).reshape(2, 4))[0]
true_active = np.array([[0.15284753, 0.67109407],
[0.69006622, -0.4165206]])
np.testing.assert_array_almost_equal(true_active, active)
def test_compute_06(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 60).reshape(15, 4)
outputs = np.random.uniform(0, 5, 15)
ss = ActiveSubspaces()
ss.compute(inputs=inputs, outputs=outputs, method='local', nboot=200)
true_evects = np.array([[-0.164383, 0.717021, -0.237246, -0.634486],
[-0.885808, 0.177628, 0.004112, 0.428691],
[0.255722, 0.558199, 0.734083, 0.290071],
[0.350612, 0.377813, -0.636254, 0.574029]])
np.testing.assert_array_almost_equal(true_evects, ss.evects)
def test_plot_sufficient_summary_03(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 200).reshape(50, 4)
weights = np.ones((50, 1)) / 50
ss = ActiveSubspaces()
ss.compute(gradients=gradients, weights=weights, nboot=200)
ss.partition(2)
with assert_plot_figures_added():
ss.plot_sufficient_summary(
np.random.uniform(-1, 1, 100).reshape(25, 4),
np.random.uniform(-1, 1, 25).reshape(-1, 1))
def test_compute_02(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 60).reshape(15, 4)
weights = np.ones((15, 1)) / 15
ss = ActiveSubspaces()
ss.compute(gradients=gradients,
weights=weights,
method='exact',
nboot=150)
true_evals = np.array([0.571596, 0.465819, 0.272198, 0.175012])
np.testing.assert_array_almost_equal(true_evals, ss.evals)
def test_compute_bootstrap_ranges_04(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 180).reshape(30, 3, 2)
weights = np.ones((30, 1)) / 30
ss = ActiveSubspaces()
ss.compute(gradients=gradients,
weights=weights,
method='exact',
nboot=100)
true_bounds_subspace = np.array([[0.00109331, 0.30254992, 0.90447872]])
np.testing.assert_array_almost_equal(true_bounds_subspace, ss.subs_br)
def test_compute_bootstrap_ranges_02(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 60).reshape(30, 2)
weights = np.ones((30, 1)) / 30
ss = ActiveSubspaces()
ss.compute(gradients=gradients,
weights=weights,
method='exact',
nboot=100)
true_bounds_subspace = np.array([[0.00261813, 0.58863862, 0.99998352]])
np.testing.assert_array_almost_equal(true_bounds_subspace, ss.subs_br)
def test_compute_07(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 180).reshape(15, 3, 4)
weights = np.ones((15, 1)) / 15
ss = ActiveSubspaces()
ss.compute(gradients=gradients,
weights=weights,
method='exact',
nboot=150)
true_evals = np.array([1.32606312, 1.20519582, 0.94811868, 0.68505712])
np.testing.assert_array_almost_equal(true_evals, ss.evals)
def test_backward_02(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 80).reshape(16, 5)
outputs = np.random.uniform(-1, 3, 16)
ss = ActiveSubspaces()
ss.compute(inputs=inputs, outputs=outputs, method='local', nboot=250)
ss.partition(2)
new_inputs = np.random.uniform(-1, 1, 15).reshape(3, 5)
active = ss.forward(new_inputs)[0]
new_inputs = ss.backward(reduced_inputs=active, n_points=500)[0]
np.testing.assert_array_almost_equal(
np.kron(active, np.ones((500, 1))), new_inputs.dot(ss.W1))
def test_compute_10(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 60).reshape(15, 4)
outputs = np.random.uniform(0, 5, 45).reshape(15, 3)
weights = np.ones((15, 1)) / 15
ss = ActiveSubspaces()
ss.compute(inputs=inputs,
outputs=outputs,
weights=weights,
method='local',
nboot=150)
true_evals = np.array(
[84.08055975, 25.87980349, 9.61982202, 8.29248646])
np.testing.assert_array_almost_equal(true_evals, ss.evals)
def test_compute_03(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 60).reshape(15, 4)
weights = np.ones((15, 1)) / 15
ss = ActiveSubspaces()
ss.compute(gradients=gradients,
weights=weights,
method='exact',
nboot=200)
true_evects = np.array([[-0.019091, -0.408566, 0.861223, -0.301669],
[-0.767799, 0.199069, 0.268823, 0.546434],
[-0.463451, -0.758442, -0.427696, -0.164486],
[-0.441965, 0.467131, -0.055723, -0.763774]])
np.testing.assert_array_almost_equal(true_evects, ss.evects)
def test_forward_03(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 60).reshape(15, 4)
outputs = np.random.uniform(0, 5, 45).reshape(15, 3)
ss = ActiveSubspaces()
ss.compute(inputs=inputs,
outputs=outputs,
method='local',
nboot=250,
metric=np.diag(np.ones(3)))
ss.partition(2)
new_inputs = np.random.uniform(-1, 1, 8).reshape(2, 4)
active, inactive = ss.forward(new_inputs)
reconstructed_inputs = active.dot(ss.W1.T) + inactive.dot(ss.W2.T)
np.testing.assert_array_almost_equal(new_inputs, reconstructed_inputs)
def test_compute_08(self):
np.random.seed(42)
gradients = np.random.uniform(-1, 1, 180).reshape(15, 3, 4)
weights = np.ones((15, 1)) / 15
ss = ActiveSubspaces()
ss.compute(gradients=gradients,
weights=weights,
method='exact',
nboot=150)
true_evects = np.array(
[[0.67237041, 0.49917148, 0.50889687, 0.1994238],
[0.20398894, -0.66183856, 0.09970486, 0.71443486],
[-0.52895262, -0.11348076, 0.83802337, -0.07104981],
[0.47593663, -0.54764923, 0.16970489, -0.66690696]])
np.testing.assert_array_almost_equal(true_evects, ss.evects)
def test_hit_and_run_inactive_01(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 60).reshape(15, 4)
outputs = np.random.uniform(0, 5, 15)
ss = ActiveSubspaces()
ss.compute(inputs=inputs, outputs=outputs, method='local', nboot=250)
ss.partition(1)
new_inputs = np.random.uniform(-1, 1, 8).reshape(2, 4)
active = ss.forward(new_inputs)[0]
inactive_swap = np.array([
ss._hit_and_run_inactive(reduced_input=red_inp, n_points=1)
for red_inp in active
])
inactive_inputs = np.swapaxes(inactive_swap, 1, 2)
new_inputs = ss._rotate_x(reduced_inputs=active,
inactive_inputs=inactive_inputs)[0]
np.testing.assert_array_almost_equal(active, new_inputs.dot(ss.W1))
def test_compute_11(self):
np.random.seed(42)
inputs = np.random.uniform(-1, 1, 60).reshape(15, 4)
outputs = np.random.uniform(0, 5, 45).reshape(15, 3)
weights = np.ones((15, 1)) / 15
metric = np.diag(2 * np.ones(3))
ss = ActiveSubspaces()
ss.compute(inputs=inputs,
outputs=outputs,
weights=weights,
method='local',
nboot=150,
metric=metric)
true_evects = np.array(
[[0.75159386, 0.34814972, 0.5093598, 0.23334745],
[-0.44047075, 0.89306499, -0.00249322, 0.09172904],
[0.35007644, 0.2587911, -0.30548878, -0.84684725],
[0.34429445, 0.11938954, -0.8045017, 0.46902504]])
np.testing.assert_array_almost_equal(true_evects, ss.evects)
def test_compute_01(self):
ss = ActiveSubspaces()
with self.assertRaises(ValueError):
ss.compute()
