class TestMlfa:
def setUp(self):
self.mlfa = MlFactorAnalysis(1,20,{})
self.cyy = np.array([[1,0.664], [0.664,1]])
self.lambdas = np.array([[0.7,0.2]])
self.taus = np.diag([0.06,0.06])
def tearDown(self):
print "TEARDOWN"
def test_objective(self):
res = self.mlfa.objective(self.cyy, self.lambdas, self.taus)
assert res < 0
def test_woodbury(self):
res = woodbury(self.lambdas, self.taus)
true = linalg.inv(self.taus + np.transpose(self.lambdas).dot(self.lambdas))
assert_almost_equal(np.sum(res), np.sum(true), 10)
def setUp(self):
self.mlfa = MlFactorAnalysis(1,20,{})
self.cyy = np.array([[1,0.664], [0.664,1]])
self.lambdas = np.array([[0.7,0.2]])
self.taus = np.diag([0.06,0.06])
