It seems that the only way to add two vectors is to use VecAXPY(x,1,y); likewise, to subtract one must call VecAXPY(x,-1,y). Are these routines optimized to recognize these special cases to avoid the extraneous scalar multiplications?