Mercurial > repos > public > sbplib

changeset 648:9e5dd0d3cf60 feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix right boundary in diracDiscr
author Martin Almquist <malmquist@stanford.edu>
date Tue, 14 Nov 2017 15:35:51 -0800
parents 46c40711830f
children 1bdbe026abbc
files diracDiscr.m
diffstat 1 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/diracDiscr.m	Tue Nov 14 15:35:28 2017 -0800
+++ b/diracDiscr.m	Tue Nov 14 15:35:51 2017 -0800
@@ -24,10 +24,10 @@
 
 % Use last tot grid points
 elseif length(poss)<tot && eta(end)<eta(1)
-    index = length(x)-tot:tot;
-    pol = (x(end-tot:end)-x(end-tot))/(x(end)-x(end-tot));
-    norm = fnorm(end-tot:end)/h;
-    x_0 = (x_0in-x(end-tot))/(x(end)-x(end-tot));
+    index = length(x)-tot+1:length(x);
+    pol = (x(end-tot+1:end)-x(end-tot+1))/(x(end)-x(end-tot+1));
+    norm = fnorm(end-tot+1:end)/h;
+    x_0 = (x_0in-x(end-tot+1))/(x(end)-x(end-tot+1));
 
 % Interior
 else