Mercurial > repos > public > sbplib

comparison 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
comparison
equal deleted inserted replaced
1233:57df0bf741dc 1235:48c9a83260c8
80 if ind_delta(end)<m 80 if ind_delta(end)<m
81 ind_delta = [ind_delta; ind_delta(end)+1]; 81 ind_delta = [ind_delta; ind_delta(end)+1];
82 else 82 else
83 ind_delta = [ind_delta(1)-1; ind_delta]; 83 ind_delta = [ind_delta(1)-1; ind_delta];
84 end 84 end
85
86 index = ind_delta;
85 polynomial = (x(ind_delta)-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1))); 87 polynomial = (x(ind_delta)-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1)));
86 x_0 = (x_s-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1))); 88 x_0 = (x_s-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1)));
87 norm = fnorm(ind_delta)/h; 89 norm = fnorm(ind_delta)/h;
88 index = ind_delta;
89 90
90 % Interior 91 % Interior
91 else 92 else
93 index = ind_delta;
92 polynomial = (x(ind_delta)-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1))); 94 polynomial = (x(ind_delta)-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1)));
93 x_0 = (x_s-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1))); 95 x_0 = (x_s-x(ind_delta(1)))/(x(ind_delta(end))-x(ind_delta(1)));
94 norm = fnorm(ind_delta)/h; 96 norm = fnorm(ind_delta)/h;
95 index = ind_delta;
96 end 97 end
97 98
98 h_polynomial = polynomial(2)-polynomial(1); 99 h_polynomial = polynomial(2)-polynomial(1);
99 b = zeros(m_order+s_order,1); 100 b = zeros(m_order+s_order,1);
100 101