Mercurial > repos > public > sbplib
comparison +time/ExplicitRungeKuttaDiscreteData.m @ 737:521076a6ffa5 feature/d1_staggered
Update Runge Kutta for discrete data to return stage rates if needed.
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Wed, 25 Apr 2018 16:58:45 -0700 |
parents | 9c3c180c7ed3 |
children | 44c46bd6913a |
comparison
equal
deleted
inserted
replaced
690:9c3c180c7ed3 | 737:521076a6ffa5 |
---|---|
43 obj.K = sparse(obj.m, obj.s); | 43 obj.K = sparse(obj.m, obj.s); |
44 obj.U = sparse(obj.m, obj.s); | 44 obj.U = sparse(obj.m, obj.s); |
45 | 45 |
46 end | 46 end |
47 | 47 |
48 function [v,t,U,T] = getV(obj) | 48 function [v,t,U,T,K] = getV(obj) |
49 v = obj.v; | 49 v = obj.v; |
50 t = obj.t; | 50 t = obj.t; |
51 U = obj.U; % Stage approximations in previous time step. | 51 U = obj.U; % Stage approximations in previous time step. |
52 T = obj.T; % Stage times in previous time step. | 52 T = obj.T; % Stage times in previous time step. |
53 K = obj.K; % Stage rates in previous time step. | |
53 end | 54 end |
54 | 55 |
55 function [a,b,c,s] = getTableau(obj) | 56 function [a,b,c,s] = getTableau(obj) |
56 a = obj.a; | 57 a = obj.a; |
57 b = obj.b; | 58 b = obj.b; |