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 done = timeTask(fmt, varargin)
    fprintf(fmt, varargin{:});
    tStart = tic;

    function done_fun()
        fprintf(' - done %fs\n', toc(tStart));
    end
    done = @done_fun;
end