Mercurial > repos > public > sbplib

comparison +rv/+time/RungekuttaRvBdf.m @ 1195:a4c00628a39d feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add higher order approximations to BDFDerivative
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 07 Aug 2019 13:27:36 +0200
parents d02e5b8a0b24
children
comparison
equal deleted inserted replaced
1194:bd5383809917 1195:a4c00628a39d
27 obj.v = v0; 27 obj.v = v0;
28 obj.n = 0; 28 obj.n = 0;
29 obj.RV = RV; 29 obj.RV = RV;
30 obj.lowerBdfOrder = bdfOrders.lowerBdfOrder; 30 obj.lowerBdfOrder = bdfOrders.lowerBdfOrder;
31 obj.upperBdfOrder = bdfOrders.upperBdfOrder; 31 obj.upperBdfOrder = bdfOrders.upperBdfOrder;
32 assert((obj.lowerBdfOrder >= 1) && (obj.upperBdfOrder <= 6)); 32 assert((obj.lowerBdfOrder >= 1) && (obj.upperBdfOrder <= 9));
33 obj.v_prev = []; 33 obj.v_prev = [];
34 obj.DvDt = rv.time.BdfDerivative(); 34 obj.DvDt = rv.time.BdfDerivative();
35 35
36 if (rkOrder == 4) % Use specialized RK4 scheme 36 if (rkOrder == 4) % Use specialized RK4 scheme
37 obj.rkScheme = @time.rk.rungekutta_4; 37 obj.rkScheme = @time.rk.rungekutta_4;