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]:
1.0105762539039627

In [5]:
plt.plot(sp.linearized_precertificate(x))


Out[5]:
[<matplotlib.lines.Line2D at 0x10600a210>]