Mercurial > repos > public > sbplib

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