Please check xAXPBY in the Reference BLAS for the following conditions:
N>0,INCY<0,xA=0, andxB/=0.
Then, if I get it correctly, xSCAL should be called,
CALL xSCAL(N, xB, xY, INCY).
However, xSCAL returns early when the increment is non-positive, i.e., does not scale xY by xB.
This does not look like the intended operation of xAXPBY.
Thanks!
(As a side note, there are two typos:
End of ZAXBPY and End of CAXBPY instead of ZAXPBY and CAXPBY, respectively.)
Please check
xAXPBYin the Reference BLAS for the following conditions:N>0,INCY<0,xA=0, andxB/=0.Then, if I get it correctly,
xSCALshould be called,CALL xSCAL(N, xB, xY, INCY).However,
xSCALreturns early when the increment is non-positive, i.e., does not scalexYbyxB.This does not look like the intended operation of
xAXPBY.Thanks!
(As a side note, there are two typos:
End of ZAXBPYandEnd of CAXBPYinstead ofZAXPBYandCAXPBY, respectively.)