In [1]:
from __future__ import division
import numpy as np
from locore import LowRank
In [2]:
n0 = 10
n = n0 * n0
q = int(0.8 * n)
Phi = np.random.randn(q,n)
sp = LowRank(Phi)
In [3]:
r = 2
x = dot(np.random.randn(n0, r), np.random.randn(r, n0))
In [4]:
sp.ic(x)
Out[4]:
In [5]:
plt.plot(sp.linearized_precertificate(x))
Out[5]: