Mercurial > repos > public > sbplib
comparison +time/CdiffTimeDep.m @ 425:e56dbd9e4196 feature/grids
Merge feature/beams
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 07 Feb 2017 16:09:02 +0100 |
| parents | 151ab2b5a686 |
| children | 8894e9c49e40 |
comparison
equal
deleted
inserted
replaced
| 423:a2cb0d4f4a02 | 425:e56dbd9e4196 |
|---|---|
| 1 classdef CdiffTimeDep < time.Timestepper | |
| 2 properties | |
| 3 D | |
| 4 E | |
| 5 S | |
| 6 k | |
| 7 t | |
| 8 v | |
| 9 v_prev | |
| 10 n | |
| 11 end | |
| 12 | |
| 13 | |
| 14 methods | |
| 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, | |
| 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) | |
| 19 function obj = CdiffTimeDep(D, E, S, k, t0, n0, v, v_prev) | |
| 20 m = length(v); | |
| 21 default_arg('E', @(t)sparse(m,m)); | |
| 22 default_arg('S', @(t)sparse(m,1)); | |
| 23 | |
| 24 obj.D = D; | |
| 25 obj.E = E; | |
| 26 obj.S = S; | |
| 27 | |
| 28 obj.k = k; | |
| 29 obj.t = t0; | |
| 30 obj.n = n0; | |
| 31 obj.v = v; | |
| 32 obj.v_prev = v_prev; | |
| 33 end | |
| 34 | |
| 35 function [v,t] = getV(obj) | |
| 36 v = obj.v; | |
| 37 t = obj.t; | |
| 38 end | |
| 39 | |
| 40 function [vt,t] = getVt(obj) | |
| 41 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) | |
| 42 t = obj.t; | |
| 43 end | |
| 44 | |
| 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)); | |
| 47 obj.t = obj.t + obj.k; | |
| 48 obj.n = obj.n + 1; | |
| 49 end | |
| 50 end | |
| 51 end |