Mercurial > repos > public > sbplib
diff +sbp/+implementations/d2_4.m @ 382:8b8672134be8 feature/beams
Edit for clarity.
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Mon, 02 Jan 2017 09:12:05 +0100 |
| parents | f7ac3cd6eeaa |
| children |
line wrap: on
line diff
--- a/+sbp/+implementations/d2_4.m Mon Jan 02 09:07:49 2017 +0100 +++ b/+sbp/+implementations/d2_4.m Mon Jan 02 09:12:05 2017 +0100 @@ -1,5 +1,5 @@ function [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = d2_4(m,h) - + BP = 4; if(m<2*BP) error(['Operator requires at least ' num2str(2*BP) ' grid points']); @@ -20,7 +20,12 @@ Q=spdiags([e -8*e 0*e 8*e -e], -2:2, m, m)/12; %Q=(-1/12*diag(ones(m-2,1),2)+8/12*diag(ones(m-1,1),1)-8/12*diag(ones(m-1,1),-1)+1/12*diag(ones(m-2,1),-2)); - Q_U = [0 0.59e2 / 0.96e2 -0.1e1 / 0.12e2 -0.1e1 / 0.32e2; -0.59e2 / 0.96e2 0 0.59e2 / 0.96e2 0; 0.1e1 / 0.12e2 -0.59e2 / 0.96e2 0 0.59e2 / 0.96e2; 0.1e1 / 0.32e2 0 -0.59e2 / 0.96e2 0;]; + Q_U = [ + 0 0.59e2/0.96e2 -0.1e1/0.12e2 -0.1e1/0.32e2; + -0.59e2/0.96e2 0 0.59e2/0.96e2 0; + 0.1e1/0.12e2 -0.59e2/0.96e2 0 0.59e2/0.96e2; + 0.1e1/0.32e2 0 -0.59e2/0.96e2 0; + ]; Q(1:4,1:4)=Q_U; Q(m-3:m,m-3:m)=rot90( -Q_U(1:4,1:4) ,2 ); @@ -30,13 +35,18 @@ %M=-(-1/12*diag(ones(m-2,1),2)+16/12*diag(ones(m-1,1),1)+16/12*diag(ones(m-1,1),-1)-1/12*diag(ones(m-2,1),-2)-30/12*diag(ones(m,1),0)); - M_U=[0.9e1 / 0.8e1 -0.59e2 / 0.48e2 0.1e1 / 0.12e2 0.1e1 / 0.48e2; -0.59e2 / 0.48e2 0.59e2 / 0.24e2 -0.59e2 / 0.48e2 0; 0.1e1 / 0.12e2 -0.59e2 / 0.48e2 0.55e2 / 0.24e2 -0.59e2 / 0.48e2; 0.1e1 / 0.48e2 0 -0.59e2 / 0.48e2 0.59e2 / 0.24e2;]; + M_U=[ + 0.9e1/0.8e1 -0.59e2/0.48e2 0.1e1/0.12e2 0.1e1/0.48e2; + -0.59e2/0.48e2 0.59e2/0.24e2 -0.59e2/0.48e2 0; + 0.1e1/0.12e2 -0.59e2/0.48e2 0.55e2/0.24e2 -0.59e2/0.48e2; + 0.1e1/0.48e2 0 -0.59e2/0.48e2 0.59e2/0.24e2; + ]; M(1:4,1:4)=M_U; M(m-3:m,m-3:m)=rot90( M_U ,2 ); M=M/h; - S_U=[-0.11e2 / 0.6e1 3 -0.3e1 / 0.2e1 0.1e1 / 0.3e1;]/h; + S_U=[-0.11e2/0.6e1 3 -0.3e1/0.2e1 0.1e1/0.3e1;]/h; S_1=sparse(sparse(1,m)); S_1(1:4)=S_U; S_m=sparse(sparse(1,m));