Mercurial > repos > public > sbplib

comparison +time/SBPInTimeSecondOrderForm.m @ 1113:47e86b5270ad feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change name of property k to dt in time.Timestepper
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 10 Apr 2019 22:40:55 +0200
parents f5e14e5986b5
children
comparison
equal deleted inserted replaced
1112:835c8fa456ec 1113:47e86b5270ad
3 A,B,C 3 A,B,C
4 M, f 4 M, f
5 5
6 n 6 n
7 t 7 t
8 k 8 dt
9 9
10 firstOrderTimeStepper 10 firstOrderTimeStepper
11 end 11 end
12 12
13 methods 13 methods
14 % Solves u_tt = Au + Bu_t + C 14 % Solves u_tt = Au + Bu_t + C
15 % A, B can either both be constants or both be function handles, 15 % A, B can either both be constants or both be function handles,
16 % They can also be omitted by setting them equal to the empty matrix. 16 % They can also be omitted by setting them equal to the empty matrix.
17 function obj = SBPInTimeSecondOrderForm(A, B, C, k, t0, v0, v0t, TYPE, order, blockSize) 17 function obj = SBPInTimeSecondOrderForm(A, B, C, dt, t0, v0, v0t, TYPE, order, blockSize)
18 default_arg('TYPE', []); 18 default_arg('TYPE', []);
19 default_arg('order', []); 19 default_arg('order', []);
20 default_arg('blockSize',[]); 20 default_arg('blockSize',[]);
21 21
22 m = length(v0); 22 m = length(v0);
37 C; 37 C;
38 ]; 38 ];
39 39
40 w0 = [v0; v0t]; 40 w0 = [v0; v0t];
41 41
42 obj.k = k; 42 obj.dt = dt;
43 obj.t = t0; 43 obj.t = t0;
44 obj.n = 0; 44 obj.n = 0;
45 45
46 obj.firstOrderTimeStepper = time.SBPInTime(obj.M, obj.f, obj.k, obj.t, w0, TYPE, order, blockSize); 46 obj.firstOrderTimeStepper = time.SBPInTime(obj.M, obj.f, obj.dt, obj.t, w0, TYPE, order, blockSize);
47 end 47 end
48 48
49 function [v,t] = getV(obj) 49 function [v,t] = getV(obj)
50 w = obj.firstOrderTimeStepper.getV(); 50 w = obj.firstOrderTimeStepper.getV();
51 v = w(1:end/2); 51 v = w(1:end/2);