Mercurial > repos > public > sbplib

comparison +time/SBPInTimeSecondOrderForm.m @ 816:b5e5b195da1e feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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 c92d2f8319c2
children f5e14e5986b5
comparison
equal deleted inserted replaced
580:2ce903f28193 816:b5e5b195da1e
56 w = obj.firstOrderTimeStepper.getV(); 56 w = obj.firstOrderTimeStepper.getV();
57 vt = w(end/2+1:end); 57 vt = w(end/2+1:end);
58 t = obj.t; 58 t = obj.t;
59 end 59 end
60 60
61 function state = getState(obj)
62 [v, t] = obj.getV();
63 [vt] = obj.getVt();
64 state = struct('v', v, 'vt', vt, 't', t, 'k', obj.k);
65 end
66
61 function obj = step(obj) 67 function obj = step(obj)
62 obj.firstOrderTimeStepper.step(); 68 obj.firstOrderTimeStepper.step();
63 obj.t = obj.firstOrderTimeStepper.t; 69 obj.t = obj.firstOrderTimeStepper.t;
64 obj.n = obj.firstOrderTimeStepper.n; 70 obj.n = obj.firstOrderTimeStepper.n;
65 end 71 end