def test_power(): # 1D x = np.array([1 / 3.0, 2 / 3.0]) y = 3 px_ = np.array([1 / 9.0, 8 / 9.0]) px = power(x, y) assert px_.shape == px.shape assert np.allclose(px_, px) x = np.array([1 / 4.0, 3 / 4.0]) y = 3 px_ = np.array([1 / 28.0, 27 / 28.0]) px = power(x, y) assert px_.shape == px.shape assert np.allclose(px_, px) # 2D with multiple rows x = np.array([[1 / 3.0, 2 / 3.0], [1 / 4.0, 3 / 4.0]]) y = [3, 3] px_ = np.array([[1 / 9.0, 8 / 9.0], [1 / 28.0, 27 / 28.0]]) px = power(x, y) assert px_.shape == px.shape assert np.allclose(px_, px) # 2D with one row x = np.array([[1 / 3.0, 2 / 3.0]]) y = [3] px_ = np.array([[1 / 9.0, 8 / 9.0]]) px = power(x, y) assert px_.shape == px.shape assert np.allclose(px_, px)
def test_power(): # 1D x = np.array([1/3.0, 2/3.0]) y = 3 px_ = np.array([1/9.0, 8/9.0]) px = power(x, y) assert_true(px_.shape == px.shape) assert_true(allclose(px_, px)) x = np.array([1/4.0, 3/4.0]) y = 3 px_ = np.array([1/28.0, 27/28.0]) px = power(x, y) assert_true(px_.shape == px.shape) assert_true(allclose(px_, px)) # 2D with multiple rows x = np.array([[1/3.0, 2/3.0], [1/4.0, 3/4.0]]) y = [3, 3] px_ = np.array([[1/9.0, 8/9.0], [1/28.0, 27/28.0]]) px = power(x, y) assert_true(px_.shape == px.shape) assert_true(allclose(px_, px)) # 2D with one row x = np.array([[1/3.0, 2/3.0]]) y = [3] px_ = np.array([[1/9.0, 8/9.0]]) px = power(x, y) assert_true(px_.shape == px.shape) assert_true(allclose(px_, px))
def test_power(): # 1D x = np.array([1 / 3, 2 / 3]) y = 3 px_ = np.array([1 / 9, 8 / 9]) px = power(x, y) assert px_.shape == px.shape assert np.allclose(px_, px) x = np.array([1 / 4, 3 / 4]) y = 3 px_ = np.array([1 / 28, 27 / 28]) px = power(x, y) assert px_.shape == px.shape assert np.allclose(px_, px) # 2D with multiple rows x = np.array([[1 / 3, 2 / 3], [1 / 4, 3 / 4]]) y = [3, 3] px_ = np.array([[1 / 9, 8 / 9], [1 / 28, 27 / 28]]) px = power(x, y) assert px_.shape == px.shape assert np.allclose(px_, px) # 2D with one row x = np.array([[1 / 3, 2 / 3]]) y = [3] px_ = np.array([[1 / 9, 8 / 9]]) px = power(x, y) assert px_.shape == px.shape assert np.allclose(px_, px)
def test_clr_prop(): x1 = np.array([1, 2, 3, 4]) x2 = np.array([3, 5, 1, 1]) a1, a2 = 0.3, 0.9 z1 = clr(perturbation(power(x1, a1), power(x2, a2))) z2 = a1 * clr(x1) + a2 * clr(x2) assert np.allclose(z1, z2)
def test_clr_prop(): x1 = np.array([1, 2, 3, 4]) x2 = np.array([3, 5, 1, 1]) a1, a2 = 0.3, 0.9 z1 = clr(perturbation(power(x1, a1), power(x2, a2))) z2 = a1 * clr(x1) + a2 * clr(x2) assert_almost_equal_array(z1, z2)