Mercurial > repos > public > sbplib

comparison +scheme/errorSbp.m @ 26:ed6a704b028d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Made some changes to error functions and comparison functions before finalizing convergence script.
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 25 Sep 2015 14:54:26 +0200
parents 16bad7c459da
children
comparison
equal deleted inserted replaced
25:791decafe6e4 26:ed6a704b028d
1 function e = errorSbp(schm, v1, v2) 1 function e = errorSbp(discr, v1, v2)
2 H = schm.H; 2 % If v1 and v2 are more complex types, something like grid functions... Then we may use .getVectorFrom here!
3 H = discr.H;
3 err = v2 - v1; 4 err = v2 - v1;
4 e = sqrt(err'*H*err); 5 e = sqrt(err'*H*err);
5 end 6 end