def test_bread_cp():
"""
Test CANDECOMP/PARFAC decomposition using bread dataset.
"""
X, meta = load_bread()
assert_raises(ValueError, cp, X)
U1 = cp(X, 2, init_type="hosvd")
U2 = _cp3(X, 2, tol=1E-4, max_iter=500, init_type="hosvd")
for n, i in enumerate(U1):
assert_almost_equal(U1[n], U2[n])
def test_bread_cp():
"""
Test CANDECOMP/PARFAC decomposition using bread dataset.
"""
X, meta = load_bread()
assert_raises(ValueError, cp, X)
U1 = cp(X, 2, init_type="hosvd")
U2 = _cp3(X, 2, tol=1E-4, max_iter=500, init_type="hosvd")
for n, i in enumerate(U1):
assert_almost_equal(U1[n], U2[n])
"""An example of using tensorlib and scipy to decompose a builtin dataset."""
from tensorlib.datasets import load_bread
from tensorlib.decomposition import cp
from scipy import linalg
import numpy as np
import matplotlib.pyplot as plt
X, meta = load_bread()
U0, U1, U2 = cp(X, n_components=2, random_state=1999)
X0_flat = X.reshape(X.shape[0], -1)
U, S, V = linalg.svd(X0_flat, full_matrices=False)
svd_proj = V[:, :2]
t1 = np.dot(X0_flat.T, U0).T
t2 = np.dot(X0_flat.T, svd_proj).T
plt.figure()
plt.title("Tensor decomposition of bread data")
plt.scatter(t1[0, :], t1[1, :], color="darkred")
plt.figure()
plt.title("SVD of bread data")
plt.scatter(t2[0, :], t2[1, :], color="steelblue")
plt.show()
def test_fetch_bread():
raise SkipTest
load_bread()
