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

function r = rowVector(v)
    if  size(v,1) == 1
        r = v
    elseif size(v,2) == 1
        r = v';
    else
        error('v is not a matrix');
    end
end