changeset 1242:ff613067dec6 feature/dirac_discr

Merge heads
author Martin Almquist <malmquist@stanford.edu>
date Tue, 19 Nov 2019 17:04:32 -0800
parents d1b201fe328e (current diff) a0940c1db455 (diff)
children 4e0b88f3def1
files diracDiscr.m
diffstat 1 files changed, 1 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/diracDiscr.m	Tue Nov 19 17:01:38 2019 -0800
+++ b/diracDiscr.m	Tue Nov 19 17:04:32 2019 -0800
@@ -32,8 +32,6 @@
 % Helper function for 1D delta functions
 function ret = diracDiscr1D(x_s, x, m_order, s_order, H)
 
-    m = length(x);
-
     % Return zeros if x0 is outside grid
     if x_s < x(1) || x_s > x(end)
         ret = zeros(size(x));
@@ -44,7 +42,7 @@
         M = [];
 
         % Get interior grid spacing
-        middle = floor(m/2);
+        middle = floor(length(x)/2);
         h = x(middle+1) - x(middle); % Use middle point to allow for staggered grids.
 
         index = sourceIndices(x_s, x, tot_order, h);