Changeset 17891

Timestamp:

04/30/14 11:10:43 (11 years ago)

Author:

cborstad

Message:

CHG: default to gaussian kernel for radial interpolation smoothing function

File:

• issm/trunk-jpl/src/m/interp/interp.py (modified) (2 diffs)

issm/trunk-jpl/src/m/interp/interp.py

@@ -187 +187 @@
         x,y:                    data point coordinates
         data:                   data to be interpolated (same length as x,y)
-
-                function: form of radial basis function for interpolation:
-                        'multiquadric': sqrt((r/self.epsilon)**2 + 1) (default)
-                        'inverse': 1.0/sqrt((r/self.epsilon)**2 + 1)
-                        'gaussian': exp(-(r/self.epsilon)**2)
-                        'linear': r
-                        'cubic': r**3
-                        'quintic': r**5
-                        'thin_plate': r**2 * log(r)
-                epsilon: adjustable constant for scaling radial distance.  Defaults to 
+        function:       form of radial basis function for interpolation:
+                                        'multiquadric': sqrt((r/self.epsilon)**2 + 1) (default)
+                                        'inverse': 1.0/sqrt((r/self.epsilon)**2 + 1)
+                                        'gaussian': exp(-(r/self.epsilon)**2)
+                                        'linear': r
+                                        'cubic': r**3
+                                        'quintic': r**5
+                                        'thin_plate': r**2 * log(r)
+        epsilon:                adjustable constant for scaling radial distance.  Defaults to 
                                         approximate average distance between nodes.
-                smooth: float>0, adjusts the amount of smoothing applied.  Defaults to 0,
+        smooth:         float>0, adjusts the amount of smoothing applied.  Defaults to 0,
                                         such that the function always passes through nodal points.
-                z:      coordinate array if interpolating in 3 dimensions
+        z:                              coordinate array if interpolating in 3 dimensions
 
         Usage:
@@ -211 +210 @@
 
         # unpack kwargs
-        function=kwargs.pop('function','multiquadric')
+        function=kwargs.pop('function','gaussian')
         if 'function' in kwargs: del kwargs['function']
         epsilon=kwargs.pop('epsilon',None)