def dist(a, b):
    return sum((i-j)**2 for i, j in zip(a, b))
def forward(self, x):
        mean_x = self.mean_module(x)
        covar_x = self.covar_module(x)
        return MultivariateNormal(mean_x, covar_x)