Mercurial > repos > public > sbplib

comparison +scheme/Utux.m @ 367:05947fc2505c feature/hypsyst

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added different operators to Utux
author Ylva Rydin <ylva.rydin@telia.com>
date Tue, 27 Dec 2016 09:11:29 +0100
parents dbac99d2c318
children 324c927d8b1d 38c3da9675a5
comparison
equal deleted inserted replaced
366:7ada2db63268 367:05947fc2505c
15 v0 15 v0
16 end 16 end
17 17
18 18
19 methods 19 methods
20 function obj = Utux(m,xlim,order) 20 function obj = Utux(m,xlim,order,operator)
21 default_arg('a',1); 21 default_arg('a',1);
22 22
23 %Old operators 23 %Old operators
24 % [x, h] = util.get_grid(xlim{:},m); 24 % [x, h] = util.get_grid(xlim{:},m);
25 %ops = sbp.Ordinary(m,h,order); 25 %ops = sbp.Ordinary(m,h,order);
26 26
27 % ops = sbp.D1Nonequidistant(m,xlim,order); 27
28 switch operator
29 case 'NonEquidistant'
30 ops = sbp.D1Nonequidistant(m,xlim,order);
31 obj.D1 = ops.D1;
32 case 'Standard'
28 ops = sbp.D2Standard(m,xlim,order); 33 ops = sbp.D2Standard(m,xlim,order);
29 obj.D1 = ops.D1; 34 obj.D1 = ops.D1;
30 % ops = sbp.D1Upwind(m,xlim,order); 35 case 'Upwind'
31 % obj.D1 = ops.Dm; 36 ops = sbp.D1Upwind(m,xlim,order);
32 obj.x=ops.x; 37 obj.D1 = ops.Dm;
38 otherwise
39 error('Unvalid operator')
40 end
41 obj.x=ops.x;
33 42
34 43
35 obj.H = ops.H; 44 obj.H = ops.H;
36 obj.Hi = ops.HI; 45 obj.Hi = ops.HI;
37 46