Mercurial > repos > public > sbplib
annotate +scheme/error2d.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 | ed6a704b028d |
children |
rev | line source |
---|---|
26
ed6a704b028d
Made some changes to error functions and comparison functions before finalizing convergence script.
Jonatan Werpers <jonatan@werpers.com>
parents:
15
diff
changeset
|
1 function e = error2d(discr, v1, v2) |
ed6a704b028d
Made some changes to error functions and comparison functions before finalizing convergence script.
Jonatan Werpers <jonatan@werpers.com>
parents:
15
diff
changeset
|
2 % If v1 and v2 are more complex types, something like grid functions... Then we may use .getVectorFrom here! |
ed6a704b028d
Made some changes to error functions and comparison functions before finalizing convergence script.
Jonatan Werpers <jonatan@werpers.com>
parents:
15
diff
changeset
|
3 h = discr.h; |
ed6a704b028d
Made some changes to error functions and comparison functions before finalizing convergence script.
Jonatan Werpers <jonatan@werpers.com>
parents:
15
diff
changeset
|
4 e = sqrt(h.^2*sum((v1-v2).^2)); |
15
16bad7c459da
Added abstract methods to allow for comparing of solutions. Added a bunch of error functions.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff
changeset
|
5 end |