Add RK6 with residual viscosity update and reduce computational effort of spatial scheme - Add RK6 with residual updates - Change the D2 operator for upwind schemes to one less computationally expensive.

function prof(f)
    profile on
    try
        f();
        profile viewer
    catch e
        fprintf(2, '\n%s', getReport(e));
        profile clear
    end
end