Mercurial > repos > public > sbplib
annotate convergence.m @ 1012:1e437c9e5132 feature/advectionRV
Create residual viscosity package +rv and generalize the ResidualViscosity class
- Generalize residual viscosity, by passing user-defined flux and calculating the time derivative outside of the update.
- Create separate RungekuttaRV specifically using interior RV updates
- Separate the artifical dissipation operator from the scheme AdvectionRV1D so that the same scheme can be reused for creating the diff op used by the ResidualViscosity class
author | Vidar Stiernström <vidar.stiernstrom@it.uu.se> |
---|---|
date | Wed, 05 Dec 2018 13:44:10 +0100 |
parents | ad1a41d135a1 |
children |
rev | line source |
---|---|
23
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
1 function q = convergence(e,h) |
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
2 for i = 1:length(e)-1 |
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
3 q(i) = log(e(i)/e(i+1))/log(h(i+1)/h(i)); |
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
4 end |
ad1a41d135a1
Added function to calculate convergence given errors. Added function to print convergence table (not finished.)
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
5 end |