Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d2_variable_periodic_2.m @ 686:5ccf6aaf6d6b feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add D2VariablePeriodic orders 4 and 6.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 08 Feb 2018 16:44:46 -0800
parents 7368affc8f78
children bbf303c1f0cf
comparison
equal deleted inserted replaced
685:b035902869a8 686:5ccf6aaf6d6b
1 function [H, HI, D1, D2, e_l, e_r, d1_l, d1_r] = d2_variable_periodic_2(m,h) 1 function [H, HI, D1, D2, e_l, e_r, d1_l, d1_r] = d2_variable_periodic_2(m,h)
2 % m = number of unique grid points, i.e. h = L/m; 2 % m = number of unique grid points, i.e. h = L/m;
3 3
4 BP = 1; 4 if(m<3)
5 if(m<2*BP) 5 error(['Operator requires at least ' num2str(3) ' grid points']);
6 error(['Operator requires at least ' num2str(2*BP) ' grid points']);
7 end 6 end
8 7
9 % Norm 8 % Norm
10 Hv = ones(m,1); 9 Hv = ones(m,1);
11 Hv = h*Hv; 10 Hv = h*Hv;