def test_compute_omega(): X = numpy.array([[1,2],[3,4],[5,6]]) M = numpy.array([[0,1],[1,0],[1,1]]) K = 2 PMF = VariationalPMF(X,M,K) PMF.compute_omega() expected_omega = [(0,1),(1,0),(2,0),(2,1)] expected_omega_I = [[1],[0],[0,1]] expected_omega_J = [[1,2],[0,2]] assert numpy.array_equal(PMF.omega,expected_omega) assert numpy.array_equal(PMF.omega_I,expected_omega_I) assert PMF.omega_J[0] == expected_omega_J[0] and PMF.omega_J[1] == expected_omega_J[1]