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_)
