Fix bug in SAT terms for upwind operators - Set the appropriate boundary derivatives in respect to choice of upwind discretization of D2 - Distinguish between upwind discretizations e.g upwind+ vs upwind-

function done = timeTask(fmt, varargin)
    fprintf(fmt, varargin{:});
    tStart = tic;

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