Apply constant offset to contact region in 'polynomial flatten' to avoid discontinuity

diff --git a/hooke/plugin/vclamp.py b/hooke/plugin/vclamp.py
index 1115cde9803e032e39c2d084423c45f42f5dccb9..770a2240010d77aa65f2dbcc51e4d07b30f2ae23 100644 (file)
--- a/hooke/plugin/vclamp.py
+++ b/hooke/plugin/vclamp.py
@@ -638,7 +638,9 @@ Name of the flattening information in the `.info` dictionary.
             'degree':degree,
             'polynomial values':poly_values,
             }
-        out = def_data - mask*scipy.polyval(poly_values, dist_data)
+        out = (def_data
+               - numpy.invert(mask)*scipy.polyval(poly_values[-1:], dist_data)
+               - mask*scipy.polyval(poly_values, dist_data))
         self._set_column(hooke=hooke, params=params,
                          column_name='output deflection column',
                          values=out)