def runG1Calculation(self):
h = Hubbard([[0, -1], [-1, 0]], 1)
ensemble = GrandcanonicalEnsemble(h, 1, .5)
statePair = (('up', 1),('up', 1))
ensemble.calcG1([statePair])
ensemble.g1.setMesh(100, -2, 2)
ensemble.g1.setRetarded(numpy.pi)
a = ensemble.g1.getSpectralFunction(statePair)
ensemble.g1.setMatsubaraMesh(20, 1)
ensemble.g1.setMatsubara()
giw = ensemble.g1.getMatsubara(statePair)
