Mercurial > repos > public > sbplib

diff diracDiscr.m @ 1235:48c9a83260c8 feature/dirac_discr

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make all if-branches have the same order of variable assignments
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 19 Nov 2019 23:51:32 +0100
parents 52d774e69b1f
children 3722c2579818
line wrap: on
line diff
--- a/diracDiscr.m	Tue Nov 19 13:55:43 2019 -0800
+++ b/diracDiscr.m	Tue Nov 19 23:51:32 2019 +0100
@@ -82,17 +82,18 @@
             else
                 ind_delta = [ind_delta(1)-1; ind_delta];
             end
+            
+            index = ind_delta;
             polynomial = (x(ind_delta)-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1)));
             x_0 = (x_s-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1)));
             norm = fnorm(ind_delta)/h;
-            index = ind_delta;
 
         % Interior
         else
+            index = ind_delta;
             polynomial = (x(ind_delta)-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1)));
             x_0 = (x_s-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1)));
             norm = fnorm(ind_delta)/h;
-            index = ind_delta;
         end
 
         h_polynomial = polynomial(2)-polynomial(1);