def logpdf_joint(self, x, y): return gu.logsumexp([np.log(.25) + norm.logpdf(x, loc=mx, scale=self.noise) + norm.logpdf(y, loc=my, scale=self.noise) for (mx,my) in zip(self.mx, self.my)])