Mercurial > repos > public > sbplib
comparison +rv/+time/RungekuttaExteriorRV.m @ 1020:5359a61cb4d9 feature/advectionRV
Add utility for constructing the operators used by a discretization emplying RV-stabilization
- Add constructDiffOps which creates the operatores used for a scheme with RV-stabilization
- Minor clean up in RungekuttaExteriorRV
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Fri, 14 Dec 2018 18:33:10 +0100 |
| parents | 2d7c1333bd6c |
| children | dce08a74e0ad |
comparison
equal
deleted
inserted
replaced
| 1019:f029b97dbc72 | 1020:5359a61cb4d9 |
|---|---|
| 68 function state = getState(obj) | 68 function state = getState(obj) |
| 69 state = struct('v', obj.v, 'residual', obj.residual, 'dvdt', obj.dvdt, 'Df', obj.Df, 'viscosity', obj.viscosity, 't', obj.t); | 69 state = struct('v', obj.v, 'residual', obj.residual, 'dvdt', obj.dvdt, 'Df', obj.Df, 'viscosity', obj.viscosity, 't', obj.t); |
| 70 end | 70 end |
| 71 | 71 |
| 72 function obj = step(obj) | 72 function obj = step(obj) |
| 73 % Store current time level and update v_prev | 73 % % Store current time level and update v_prev |
| 74 % numStoredStages = size(obj.v_prev,2); | 74 % numStoredStages = size(obj.v_prev,2); |
| 75 % if (numStoredStages < obj.upperBdfOrder) | 75 % if (numStoredStages < obj.upperBdfOrder) |
| 76 % obj.v_prev = [obj.v, obj.v_prev]; | 76 % obj.v_prev = [obj.v, obj.v_prev]; |
| 77 % numStoredStages = numStoredStages+1; | 77 % numStoredStages = numStoredStages+1; |
| 78 % else | 78 % else |