Mercurial > repos > public > sbplib
comparison +time/CdiffNonlin.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 | d1f9dd55a2b0 |
| children | f5e14e5986b5 |
comparison
equal
deleted
inserted
replaced
| 580:2ce903f28193 | 816:b5e5b195da1e |
|---|---|
| 48 function [vt,t] = getVt(obj) | 48 function [vt,t] = getVt(obj) |
| 49 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) | 49 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) |
| 50 t = obj.t; | 50 t = obj.t; |
| 51 end | 51 end |
| 52 | 52 |
| 53 function state = getState(obj) | |
| 54 [v, t] = obj.getV(); | |
| 55 [vt] = obj.getVt(); | |
| 56 state = struct('v', v, 'vt', vt, 't', t, 'k', obj.k); | |
| 57 end | |
| 58 | |
| 53 function obj = step(obj) | 59 function obj = step(obj) |
| 54 D = obj.D(obj.v); | 60 D = obj.D(obj.v); |
| 55 E = obj.E(obj.v); | 61 E = obj.E(obj.v); |
| 56 S = obj.S(obj.v,obj.t); | 62 S = obj.S(obj.v,obj.t); |
| 57 | 63 |