Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d4_variable_4.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_4(m,h)
2 S3_m, S_1, S_m] = d4_compatible_halfvariable_4(m,h)
3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
4 %%% 4:de ordn. SBP Finita differens %%% 3 %%% 4:de ordn. SBP Finita differens %%%
5 %%% %%% 4 %%% %%%
6 %%% H (Normen) %%% 5 %%% H (Normen) %%%
7 %%% D1=H^(-1)Q (approx f?rsta derivatan) %%% 6 %%% D1=H^(-1)Q (approx f?rsta derivatan) %%%
10 9
11 10
12 %m=20; %problemstorlek 11 %m=20; %problemstorlek
13 %h=1/(m-1); 12 %h=1/(m-1);
14 %h=1; 13 %h=1;
15 14
16 BP = 6; 15 BP = 6;
17 if(m<2*BP) 16 if(m<2*BP)
18 error(['Operator requires at least ' num2str(2*BP) ' grid points']); 17 error(['Operator requires at least ' num2str(2*BP) ' grid points']);
19 end 18 end
20 19