Mercurial > repos > public > sbplib

comparison +time/Ode45.m @ 889:f5e14e5986b5 feature/timesteppers

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove getState from TimeStepper and derived classes - Currently no need for that functionality. Remove it.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 15 Nov 2018 17:29:30 -0800
parents b5e5b195da1e
children 47e86b5270ad
comparison
equal deleted inserted replaced
888:8732d6bd9890 889:f5e14e5986b5
46 function [vt,t] = getVt(obj) 46 function [vt,t] = getVt(obj)
47 vt = obj.w(end/2+1:end); 47 vt = obj.w(end/2+1:end);
48 t = obj.t; 48 t = obj.t;
49 end 49 end
50 50
51 function state = getState(obj)
52 [v, t] = obj.getV();
53 [vt] = obj.getVt();
54 state = struct('v', v, 'vt', vt, 't', t, 'k', obj.k);
55 end
56
57 function obj = step(obj) 51 function obj = step(obj)
58 [t,w] = ode45(@(t,w)(obj.F(w,t)),[obj.t obj.t+obj.k],obj.w); 52 [t,w] = ode45(@(t,w)(obj.F(w,t)),[obj.t obj.t+obj.k],obj.w);
59 53
60 obj.t = t(end); 54 obj.t = t(end);
61 obj.w = w(end,:)'; 55 obj.w = w(end,:)';