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
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
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