Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d4_compatible_4.m @ 266:bfa130b7abf6 operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added error message for too few grid points to all implementation files.
author Martin Almquist <martin.almquist@it.uu.se>
date Fri, 09 Sep 2016 11:03:13 +0200
parents 6009f2712d13
children f7ac3cd6eeaa
comparison
equal deleted inserted replaced
265:3e0432b3cb68 266:bfa130b7abf6
23 % I annat fall anv?nd de "traditionella" som har noggrannare 23 % I annat fall anv?nd de "traditionella" som har noggrannare
24 % randsplutningar f?r D1 och D2 24 % randsplutningar f?r D1 och D2
25 25
26 % Vi b?rjar med normen. Notera att alla SBP operatorer delar samma norm, 26 % Vi b?rjar med normen. Notera att alla SBP operatorer delar samma norm,
27 % vilket ?r n?dv?ndigt f?r stabilitet 27 % vilket ?r n?dv?ndigt f?r stabilitet
28
29 BP = 6;
30 if(m<2*BP)
31 error(['Operator requires at least ' num2str(2*BP) ' grid points']);
32 end
28 33
29 H=diag(ones(m,1),0); 34 H=diag(ones(m,1),0);
30 H_U=[0.3e1 / 0.11e2 0 0 0 0 0; 0 0.2125516311e10 / 0.1311004640e10 0 0 0 0; 0 0 0.278735189e9 / 0.1966506960e10 0 0 0; 0 0 0 0.285925927e9 / 0.163875580e9 0 0; 0 0 0 0 0.1284335339e10 / 0.1966506960e10 0; 0 0 0 0 0 0.4194024163e10 / 0.3933013920e10;]; 35 H_U=[0.3e1 / 0.11e2 0 0 0 0 0; 0 0.2125516311e10 / 0.1311004640e10 0 0 0 0; 0 0 0.278735189e9 / 0.1966506960e10 0 0 0; 0 0 0 0.285925927e9 / 0.163875580e9 0 0; 0 0 0 0 0.1284335339e10 / 0.1966506960e10 0; 0 0 0 0 0 0.4194024163e10 / 0.3933013920e10;];
31 H(1:6,1:6)=H_U; 36 H(1:6,1:6)=H_U;
32 H(m-5:m,m-5:m)=fliplr(flipud(H_U)); 37 H(m-5:m,m-5:m)=fliplr(flipud(H_U));