Double initial guess for the surface kink position.

[hooke.git] / hooke / util / fit.py

diff --git a/hooke/util/fit.py b/hooke/util/fit.py
index 13606773be4417a02d9384b7ff0bfee0648929e6..6d2e30333ac4e48c2729211c738cc73e099b7e06 100644 (file)
--- a/hooke/util/fit.py
+++ b/hooke/util/fit.py
@@ -187,6 +187,12 @@ class ModelFitter (object):
         if self._rescale == True:
             params = [p*s for p,s in zip(params, self._param_scale_factors)]
         residual = self._data - self.model(params)
+        if False:  # fit debugging
+            if not hasattr(self, '_i_'):
+                self._i_ = 0
+            self._data.tofile('data.%d' % self._i_, sep='\n')
+            self.model(params).tofile('model.%d' % self._i_, sep='\n')
+            self._i_ += 1
         if self._rescale == True:
             residual /= self._data_scale_factor
         return residual