Mercurial > repos > public > sbplib

annotate +scheme/errorSbp.m @ 998:2b1b944deae1 feature/getBoundaryOp

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add getBoundaryOperator to all 1d schemes. Did not add getBoundaryQuadrature because it doesnt make sense in 1d (?)
author Martin Almquist <malmquist@stanford.edu>
date Sat, 12 Jan 2019 13:35:19 -0800
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