Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d4_variable_6.m @ 311:713b125038a3 feature/beams

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed function names in renamed files.
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 23 Sep 2016 14:59:55 +0200
parents ffa5d557942b
children 9230c056a574
comparison
equal deleted inserted replaced
310:ffa5d557942b 311:713b125038a3
1 function [H, HI, D2, D4, e_1, e_m, M4, S2_1, S2_m, S3_1,... 1 function [H, HI, D2, D4, e_1, e_m, M4, S2_1, S2_m, S3_1, S3_m, S_1, S_m] = d4_variable_6(m,h)
2 S3_m, S_1, S_m] = d4_compatible_halfvariable_6(m,h)
3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4 %%% 6:te ordn. SBP Finita differens %%% 3 %%% 6:te ordn. SBP Finita differens %%%
5 %%% operatorer med diagonal norm %%% 4 %%% operatorer med diagonal norm %%%
6 %%% Extension to variable koeff %%% 5 %%% Extension to variable koeff %%%
7 %%% %%% 6 %%% %%%
23 22
24 % Variable koefficicients are stored in vector: c, size m, 23 % Variable koefficicients are stored in vector: c, size m,
25 % with the unknown stored as c(1), c(2), ..., c_m 24 % with the unknown stored as c(1), c(2), ..., c_m
26 % x=1:h:m*h;x=x'; 25 % x=1:h:m*h;x=x';
27 % c=x.^0; 26 % c=x.^0;
28 27
29 BP = 8; 28 BP = 8;
30 if(m<2*BP) 29 if(m<2*BP)
31 error(['Operator requires at least ' num2str(2*BP) ' grid points']); 30 error(['Operator requires at least ' num2str(2*BP) ' grid points']);
32 end 31 end
33 32
41 40
42 % x1=0.70127127127127; 41 % x1=0.70127127127127;
43 42
44 43
45 % D1=(1/60*diag(ones(m-3,1),3)-9/60*diag(ones(m-2,1),2)+45/60*diag(ones(m-1,1),1)-45/60*diag(ones(m-1,1),-1)+9/60*diag(ones(m-2,1),-2)-1/60*diag(ones(m-3,1),-3)); 44 % D1=(1/60*diag(ones(m-3,1),3)-9/60*diag(ones(m-2,1),2)+45/60*diag(ones(m-1,1),1)-45/60*diag(ones(m-1,1),-1)+9/60*diag(ones(m-2,1),-2)-1/60*diag(ones(m-3,1),-3));
46 % 45 %
47 % 46 %
48 % 47 %
49 % D1(1:6,1:9)=[-21600/13649, 43200/13649*x1-7624/40947, -172800/13649*x1+ ... 48 % D1(1:6,1:9)=[-21600/13649, 43200/13649*x1-7624/40947, -172800/13649*x1+ ...
50 % 715489/81894, 259200/13649*x1-187917/13649, -172800/13649* ... 49 % 715489/81894, 259200/13649*x1-187917/13649, -172800/13649* ...
51 % x1+735635/81894, 43200/13649*x1-89387/40947, 0, 0, 0; ... 50 % x1+735635/81894, 43200/13649*x1-89387/40947, 0, 0, 0; ...
52 % -8640/12013*x1+7624/180195, 0, 86400/12013*x1-57139/12013, ... 51 % -8640/12013*x1+7624/180195, 0, 86400/12013*x1-57139/12013, ...
53 % -172800/12013*x1+745733/72078, 129600/12013*x1-91715/12013, ... 52 % -172800/12013*x1+745733/72078, 129600/12013*x1-91715/12013, ...