In [1]:

    

cd ..


    

    




/Users/dsuess/Code/CS Algorithms

In [2]:

    

import mpnum as mp
import numpy as np

from itertools import islice
from csalgs.tt.altmin import AltminEstimator
from mpnum.special import inner_prod_mps


    

In [3]:

    

DIM = 4
RANK = 10
SITES = 4
NR_MEASUREMENTS = 1 * SITES * DIM * RANK**2

X_true = mp.random_mpa(SITES, DIM, RANK, normalized=True)
A = [mp.random_mpa(SITES, DIM, 1, normalized=True) 
     for _ in range(NR_MEASUREMENTS)]
y = [inner_prod_mps(a, X_true) for a in A]


    

In [4]:

    

reconstructions = islice(AltminEstimator(A, y, RANK), 20)
pl.plot(np.log([mp.normdist(X_true, X_sharp) 
                for X_sharp in reconstructions]))


    

    
Out[4]:





[<matplotlib.lines.Line2D at 0x10f0b5cc0>]






    

In [ ]: