Mercurial > repos > public > sbplib
comparison +time/CdiffTimeDep.m @ 886:8894e9c49e40 feature/timesteppers
Merge with default for latest changes
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Thu, 15 Nov 2018 16:36:21 -0800 |
| parents | +time/Cdiff.m@b5e5b195da1e +time/Cdiff.m@151ab2b5a686 |
| children | f5e14e5986b5 |
comparison
equal
deleted
inserted
replaced
| 816:b5e5b195da1e | 886:8894e9c49e40 |
|---|---|
| 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 state = getState(obj) | |
| 46 state = struct('v', obj.v, 'v_prev', obj.v_prev, 't', obj.t, 'k', obj.k); | |
| 47 end | |
| 48 | |
| 49 function obj = step(obj) | |
| 50 [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)); | |
| 51 obj.t = obj.t + obj.k; | |
| 52 obj.n = obj.n + 1; | |
| 53 end | |
| 54 end | |
| 55 end |