Mercurial > repos > public > sbplib
diff +rv/ResidualViscosity.m @ 1014:e547794a9407 feature/advectionRV
Add boot-strapping to RungeKuttaExteriorRV
- Higher order BDF approximations are successively used as increasing number of time levels are obtained.
| author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
|---|---|
| date | Thu, 06 Dec 2018 11:30:47 +0100 |
| parents | 1e437c9e5132 |
| children | 9b7fcd5e4480 |
line wrap: on
line diff
--- a/+rv/ResidualViscosity.m Wed Dec 05 15:04:44 2018 +0100 +++ b/+rv/ResidualViscosity.m Thu Dec 06 11:30:47 2018 +0100 @@ -40,6 +40,13 @@ obj.viscosity = min(obj.Cmax*obj.h*abs(obj.waveSpeed(v)), obj.Cres*obj.h^2*abs(obj.residual)/norm(v-mean(v),inf)); end + function smoothendVector = smoothen(obj, vector) + smoothendVector = vector; + for i = 2:length(vector)-1 + smoothendVector(i) = (1/6)*(vector(i-1) + 4*vector(i) + vector(i+1)); + end + end + function [residual, u_t, grad_f] = getResidual(obj) residual = obj.residual; u_t = obj.u_t;