Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d2_noneq_variable_8.m @ 1330:855871e0b852 feature/D2_boundary_opt

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add size checks to the operators
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 17 Feb 2022 18:55:11 +0100
parents c2d716c4f1ed
children 8e9df030a0a5
comparison
equal deleted inserted replaced
1329:7df63b17e078 1330:855871e0b852
8 8
9 % BP: Number of boundary points 9 % BP: Number of boundary points
10 % order: Accuracy of interior stencil 10 % order: Accuracy of interior stencil
11 BP = 8; 11 BP = 8;
12 order = 8; 12 order = 8;
13 if(N<2*BP)
14 error(['Operator requires at least ' num2str(2*BP) ' grid points']);
15 end
13 16
14 %%%% Norm matrix %%%%%%%% 17 %%%% Norm matrix %%%%%%%%
15 P = zeros(BP, 1); 18 P = zeros(BP, 1);
16 P0 = 1.0758368078310e-01; 19 P0 = 1.0758368078310e-01;
17 P1 = 6.1909685107891e-01; 20 P1 = 6.1909685107891e-01;