Mercurial > repos > public > sbplib
comparison +time/Cdiff.m @ 816:b5e5b195da1e feature/timesteppers
Add getState to timesteppers, returning the relevant state of the timestepper
- Add getState which returns the 'state' of the specialized timestepper.
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Mon, 10 Sep 2018 16:19:16 +0200 |
| parents | 484b48e95c83 |
| children | 8894e9c49e40 |
comparison
equal
deleted
inserted
replaced
| 580:2ce903f28193 | 816:b5e5b195da1e |
|---|---|
| 43 function [vt,t] = getVt(obj) | 43 function [vt,t] = getVt(obj) |
| 44 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) | 44 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) |
| 45 t = obj.t; | 45 t = obj.t; |
| 46 end | 46 end |
| 47 | 47 |
| 48 function state = getState(obj) | |
| 49 state = struct('v', obj.v, 'v_prev', obj.v_prev, 't', obj.t, 'k', obj.k); | |
| 50 end | |
| 51 | |
| 48 function obj = step(obj) | 52 function obj = step(obj) |
| 49 [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.k, obj.D, obj.E, obj.S); | 53 [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.k, obj.D, obj.E, obj.S); |
| 50 obj.t = obj.t + obj.k; | 54 obj.t = obj.t + obj.k; |
| 51 obj.n = obj.n + 1; | 55 obj.n = obj.n + 1; |
| 52 end | 56 end |