Mercurial > repos > public > sbplib

annotate +scheme/errorSbp.m @ 1036:8a9393084b30 feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change argument order to the "correct" order, i.e providing diffOp specific parameters before the opSet.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Fri, 18 Jan 2019 08:58:26 +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 = 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