def METHOD_NAME(self, x, y):
    """
    Calculates the elementwise determinant of the log Jacobian, dx/dy
    """
    if self.training:
        var = torch.var(y, dim=0, keepdim=True)
    else:
        # NOTE: You wouldn't typically run this function in eval mode, but included for gradient tests
        var = self.moving_variance
    return -self.constrained_gamma.log() + 0.5 * torch.log(var + self.epsilon)