There are a bunch of cases where the alpha are multiplied by -1 before being passed to VecMAXPY() maybe a 1/3 to 1/2 the uses in KSP
There are other cases (in TS) where the alpha are multipled by h before being passed to VecMAXPY()
Is this change simply to make the calling code simplier and not require the scaling of the alpha initially? Note that to implement the new API efficiently essentially one would scale the alpha on entry and scale them back on exit. No more efficient then today and less efficient when they are not used later and do not need to be scaled back. Is there another way to implement this efficiently?