Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d4_variable_6.m @ 329:bf801c3709be feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bug fixes in operators.
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 28 Sep 2016 16:31:06 +0200
parents b19e142fcae1
children 43d02533bea3
comparison
equal deleted inserted replaced
328:31d6698c1edf 329:bf801c3709be
125 M(5,10)=M(10,5); 125 M(5,10)=M(10,5);
126 M(m-4,m-9)=M(m-9,m-4); 126 M(m-4,m-9)=M(m-9,m-4);
127 127
128 M=M/h; 128 M=M/h;
129 129
130 D2=HI*(-M-diag(c)*e_1*S_1'+diag(c)*e_m*S_m'); 130 D2 = HI*(-M - c(1)*e_l*d1_l' + c(m)*e_r*d1_r');
131 end 131 end
132 D2 = @D2_fun; 132 D2 = @D2_fun;
133 133
134 % Fourth derivative, 1th order accurate at first 8 boundary points (still 134 % Fourth derivative, 1th order accurate at first 8 boundary points (still
135 % yield 5th order convergence if stable: for example u_tt=-u_xxxx 135 % yield 5th order convergence if stable: for example u_tt=-u_xxxx