Mercurial > repos > public > sbplib

annotate +scheme/errorSbp.m @ 815:fae41958af4f feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add support for artificial viscosity to the 1d burgers scheme. - Add support for artificial viscosity by parametrizing the discretization operator and boundary closure on the viscosity. - Add a time stepper which evaluates and updates the residual viscosity of the solution.
author Vidar Stiernstrom <vidar.stiernstrom@it.uu.se>
date Thu, 06 Sep 2018 12:43:51 +0200
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 = errorSbp(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;
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
4 err = v2 - v1;
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 e = sqrt(err'*H*err);
16bad7c459da Added abstract methods to allow for comparing of solutions. Added a bunch of error functions.
Jonatan Werpers <jonatan@werpers.com>
parents:
diff changeset
6 end