Mercurial > repos > public > sbplib

annotate +scheme/errorSbp.m @ 104:ffd9e68f63cc feature/arclen-param

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Splined in the other direction! Works great. Runs faster. Removed ugly fix outside [0,1] in spline(). Camelcased numerical_derivative().
author Martin Almquist <martin.almquist@it.uu.se>
date Mon, 07 Dec 2015 18:46:59 +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