builtins: fix gaussian kernel for non-square

devito/builtins/initializers.py

@@ -145,9 +145,9 @@ def create_gaussian_weights(sigma, lw):
         weights = [w/w.sum() for w in (np.exp(-0.5/s**2*(np.linspace(-l, l, 2*l+1))**2)
                    for s, l in zip(sigma, lw))]
         processed = []
-        for w in weights:
+        for (w, l) in zip(weights, lw):
             temp = list(w)
-            while len(temp) < 2*max(lw)+1:
+            while len(temp) < 2*l+1:
                 temp.insert(0, 0)
                 temp.append(0)
             processed.append(np.array(temp))

devito/finite_differences/tools.py

@@ -229,7 +229,7 @@ def numeric_weights(function, deriv_order, indices, x0):
 coeff_priority = {'taylor': 1, 'standard': 1}
 
 
-def generate_indices(expr, dim, order, side=None, matvec=None, x0=None, nweights=None):
+def generate_indices(expr, dim, order, side=None, matvec=None, x0=None, nweights=0):
     """
     Indices for the finite-difference scheme.