Mercurial > repos > public > sbplib

comparison +time/CdiffTimeDep.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
1 classdef CdiffTimeDep < time.Timestepper 1 classdef CdiffTimeDep < time.Timestepper
2 properties 2 properties
3 D 3 D
4 E 4 E
5 S 5 S
6 k 6 dt
7 t 7 t
8 v 8 v
9 v_prev 9 v_prev
10 n 10 n
11 end 11 end
13 13
14 methods 14 methods
15 % Solves u_tt = Du + E(t)u_t + S(t) 15 % Solves u_tt = Du + E(t)u_t + S(t)
16 % D, E, S can either all be constants or all be function handles, 16 % D, E, S can either all be constants or all be function handles,
17 % They can also be omitted by setting them equal to the empty matrix. 17 % They can also be omitted by setting them equal to the empty matrix.
18 % CdiffTimeDep(D, E, S, k, t0, n0, v, v_prev) 18 % CdiffTimeDep(D, E, S, dt, t0, n0, v, v_prev)
19 function obj = CdiffTimeDep(D, E, S, k, t0, n0, v, v_prev) 19 function obj = CdiffTimeDep(D, E, S, dt, t0, n0, v, v_prev)
20 m = length(v); 20 m = length(v);
21 default_arg('E', @(t)sparse(m,m)); 21 default_arg('E', @(t)sparse(m,m));
22 default_arg('S', @(t)sparse(m,1)); 22 default_arg('S', @(t)sparse(m,1));
23 23
24 obj.D = D; 24 obj.D = D;
25 obj.E = E; 25 obj.E = E;
26 obj.S = S; 26 obj.S = S;
27 27
28 obj.k = k; 28 obj.dt = dt;
29 obj.t = t0; 29 obj.t = t0;
30 obj.n = n0; 30 obj.n = n0;
31 obj.v = v; 31 obj.v = v;
32 obj.v_prev = v_prev; 32 obj.v_prev = v_prev;
33 end 33 end
36 v = obj.v; 36 v = obj.v;
37 t = obj.t; 37 t = obj.t;
38 end 38 end
39 39
40 function [vt,t] = getVt(obj) 40 function [vt,t] = getVt(obj)
41 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) 41 vt = (obj.v-obj.v_prev)/obj.dt; % Could be improved using u_tt = f(u))
42 t = obj.t; 42 t = obj.t;
43 end 43 end
44 44
45 function obj = step(obj) 45 function obj = step(obj)
46 [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.k, obj.D, obj.E(obj.t), obj.S(obj.t)); 46 [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.dt, obj.D, obj.E(obj.t), obj.S(obj.t));
47 obj.t = obj.t + obj.k; 47 obj.t = obj.t + obj.dt;
48 obj.n = obj.n + 1; 48 obj.n = obj.n + 1;
49 end 49 end
50 end 50 end
51 end 51 end