Mercurial > repos > public > sbplib

comparison +time/Cdiff.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 8894e9c49e40
children 7a5e770974ed
comparison
equal deleted inserted replaced
888:8732d6bd9890 889:f5e14e5986b5
41 function [vt,t] = getVt(obj) 41 function [vt,t] = getVt(obj)
42 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u)) 42 vt = (obj.v-obj.v_prev)/obj.k; % Could be improved using u_tt = f(u))
43 t = obj.t; 43 t = obj.t;
44 end 44 end
45 45
46 function state = getState(obj)
47 state = struct('v', obj.v, 'v_prev', obj.v_prev, 't', obj.t, 'k', obj.k);
48 end
49
50 function obj = step(obj) 46 function obj = step(obj)
51 [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.k, obj.D, obj.E, obj.S); 47 [obj.v, obj.v_prev] = time.cdiff.cdiff(obj.v, obj.v_prev, obj.k, obj.D, obj.E, obj.S);
52 obj.t = obj.t + obj.k; 48 obj.t = obj.t + obj.k;
53 obj.n = obj.n + 1; 49 obj.n = obj.n + 1;
54 end 50 end