Mercurial > repos > public > sbplib

comparison +sbp/D1Nonequidistant.m @ 259:a94bb3d92aba operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added empty borrowing property to nonequidstant class.
author Martin Almquist <martin.almquist@it.uu.se>
date Thu, 08 Sep 2016 13:44:13 +0200
parents a8d89688ba01
children 6009f2712d13
comparison
equal deleted inserted replaced
258:a8d89688ba01 259:a94bb3d92aba
7 e_1 % Left boundary operator 7 e_1 % Left boundary operator
8 e_m % Right boundary operator 8 e_m % Right boundary operator
9 m % Number of grid points. 9 m % Number of grid points.
10 h % Step size 10 h % Step size
11 x % grid 11 x % grid
12 borrowing % Struct with borrowing limits for different norm matrices
12 end 13 end
13 14
14 methods 15 methods
15 function obj = D1Nonequidistant(m,L,order,option) 16 function obj = D1Nonequidistant(m,L,order,option)
16 17
17 default_arg('option','Accurate'); 18 default_arg('option','Accurate');
18 % 'Accurate' operators are optimized for accuracy 19 % 'Accurate' operators are optimized for accuracy
19 % 'Minimal' operators have the smallest possible boundary 20 % 'Minimal' operators have the smallest possible boundary
20 % closure 21 % closure
21 22
22 switch option 23 switch option
23 24
24 case {'Accurate','accurate','A'} 25 case {'Accurate','accurate','A'}
25 26
61 obj.e_m(m) = 1; 62 obj.e_m(m) = 1;
62 63
63 obj.HI = inv(obj.H); 64 obj.HI = inv(obj.H);
64 obj.Q = obj.H*obj.D1 - obj.e_m*obj.e_m' + obj.e_0*obj.e_0'; 65 obj.Q = obj.H*obj.D1 - obj.e_m*obj.e_m' + obj.e_0*obj.e_0';
65 66
67 obj.borrowing = [];
68
66 end 69 end
67 end 70 end
68 71
69 methods (Static) 72 methods (Static)
70 function lambda = smallestGrid(obj) 73 function lambda = smallestGrid(obj)