sdsdot(l) - Linux man page
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