sdsdot

SDSDOT(1)                        BLAS routine                        SDSDOT(1)



NAME
SYNOPSIS
       REAL FUNCTION SDSDOT(N,SB,SX,INCX,SY,INCY)

           REAL      SB

           INTEGER   INCX,INCY,N

           REAL      SX(*),SY(*)

           DOUBLE    PRECISION DSDOT

           INTEGER   I,KX,KY,NS

           INTRINSIC DBLE

           DSDOT     = SB

           IF        (N.LE.0) GO TO 30

           IF        (INCX.EQ.INCY .AND. INCX.GT.0) GO TO 40

           KX        = 1

           KY        = 1

           IF        (INCX.LT.0) KX = 1 + (1-N)*INCX

           IF        (INCY.LT.0) KY = 1 + (1-N)*INCY

           DO        10 I = 1,N

           DSDOT     = DSDOT + DBLE(SX(KX))*DBLE(SY(KY))

           KX        = KX + INCX

           KY        = KY + INCY

           10        CONTINUE

           30        SDSDOT = DSDOT

           RETURN

           40        NS = N*INCX

           DO        50 I = 1,NS,INCX

           DSDOT     = DSDOT + DBLE(SX(I))*DBLE(SY(I))

           50        CONTINUE

           SDSDOT    = DSDOT

           RETURN

           END

PURPOSE
BLAS routine                     November 2008                       SDSDOT(1)