diff --git a/GPy/models/sparse_GP.py b/GPy/models/sparse_GP.py
index b816e684..7da36835 100644
--- a/GPy/models/sparse_GP.py
+++ b/GPy/models/sparse_GP.py
@@ -185,11 +185,12 @@ class sparse_GP(GP):
         self.likelihood._set_params(p[self.Z.size+self.kern.Nparam:])
         self._compute_kernel_matrices()
         if self.auto_scale_factor:
-            if self.likelihood.is_heteroscedastic:
-                self.scale_factor = max(100.,(self.psi2_beta_scaled.sum(0).max()))
-                print self.scale_factor
-            else:
-                self.scale_factor = np.sqrt(self.psi2.sum(0).mean()*self.likelihood.precision)
+            self.scale_factor = np.sqrt(self.psi2.sum(0).mean()*self.likelihood.precision)
+        #if self.auto_scale_factor:
+        #    if self.likelihood.is_heteroscedastic:
+        #        self.scale_factor = max(1,np.sqrt(self.psi2_beta_scaled.sum(0).mean()))
+        #    else:
+        #        self.scale_factor = np.sqrt(self.psi2.sum(0).mean()*self.likelihood.precision)
         self._computations()
 
     def _get_params(self):