Allow for the residual normalization function to return a vector. Also change the default normalization from normalizing on norm(u-mean(u),inf) to norm(u/2), since this proved to be better for 1d burgers

% Assert that array A has the size s.
function assert_size(A,s)
    warning('Use assertSize() instead!')
    assertSize(A,s);
end