def test_basis2(): # test orthonormality in Aitchison inner product b = basis(4) for i in range(len(b)): ii = inner(b[i], b[i]) assert ii == pytest.approx(1.0) for j in range(i + 1, len(b)): ij = inner(b[i], b[j]) assert ij == pytest.approx(0.0)
def test_basis(): b = np.array([[1 / 1., -1, 0, 0, 0], [1 / 2., 1 / 2., -1, 0, 0], [1 / 3., 1 / 3., 1 / 3., -1, 0], [1 / 4., 1 / 4., 1 / 4., 1 / 4., -1]]) b *= np.sqrt([i / (i + 1) for i in range(1, 5)])[:, np.newaxis] b = closure(exp2(b)) b_ = basis(4) assert b.shape == b_.shape assert np.allclose(b, b_)
def test_basis2(): # test orthonormality in Aitchison inner product b = basis(4) for i in range(len(b)): ii = inner(b[i], b[i]) assert_almost_equal(ii, 1.0) for j in range(i + 1, len(b)): ij = inner(b[i], b[j]) assert_almost_equal(ij, 0.0)
def test_basis(): b = np.array([[1 / 1., -1, 0, 0, 0], [1 / 2., 1 / 2., -1, 0, 0], [1 / 3., 1 / 3., 1 / 3., -1, 0], [1 / 4., 1 / 4., 1 / 4., 1 / 4., -1]]) b *= np.sqrt([i / (i + 1) for i in range(1, 5)])[:, np.newaxis] b = closure(exp2(b)) b_ = basis(4) assert_equal_shape(b, b_) assert_almost_equal_array(b, b_)
def test_basis2(): # test orthonormality in Aitchison inner product b = basis(4) for i in range(len(b)): ii = inner(b[i], b[i]) assert_almost_equal(ii, 1.0) for j in range(i+1, len(b)): ij = inner(b[i], b[j]) assert_almost_equal(ij, 0.0)
def test_basis(): b = np.array([[1/1., -1, 0, 0, 0], [1/2., 1/2., -1, 0, 0], [1/3., 1/3., 1/3., -1, 0], [1/4., 1/4., 1/4., 1/4., -1]]) b *= np.sqrt([i/(i+1) for i in range(1, 5)])[:, np.newaxis] b = closure(exp2(b)) b_ = basis(4) assert_equal_shape(b, b_) assert_almost_equal_array(b, b_)