rpl/lapack/blas/dsdot.f - diff

Return to dsdot.f CVS log

Up to [local] / rpl / lapack / blas

-version 1.4, 2010/12/21 13:51:25
+version 1.5, 2011/07/22 07:38:01
  Line 72
  *     ..
        DSDOT = 0.0D0
        IF (N.LE.0) RETURN
-       IF (INCX.EQ.INCY .AND. INCX.GT.0) GO TO 20
+       IF (INCX.EQ.INCY .AND. INCX.GT.0) THEN
  *
- *     Code for unequal or nonpositive increments.
+ *     Code for equal, positive, non-unit increments.
  *
-       KX = 1
+          NS = N*INCX
-       KY = 1
+          DO I = 1,NS,INCX
-       IF (INCX.LT.0) KX = 1 + (1-N)*INCX
+             DSDOT = DSDOT + DBLE(SX(I))*DBLE(SY(I))
-       IF (INCY.LT.0) KY = 1 + (1-N)*INCY
+          END DO
-       DO 10 I = 1,N
+       ELSE
-           DSDOT = DSDOT + DBLE(SX(KX))*DBLE(SY(KY))
-           KX = KX + INCX
-           KY = KY + INCY
-CONTINUE
-       RETURN
  *
- *     Code for equal, positive, non-unit increments.
+ *     Code for unequal or nonpositive increments.
  *
-NS = N*INCX
+          KX = 1
-       DO 30 I = 1,NS,INCX
+          KY = 1
-           DSDOT = DSDOT + DBLE(SX(I))*DBLE(SY(I))
+          IF (INCX.LT.0) KX = 1 + (1-N)*INCX
-CONTINUE
+          IF (INCY.LT.0) KY = 1 + (1-N)*INCY
+          DO I = 1,N
+             DSDOT = DSDOT + DBLE(SX(KX))*DBLE(SY(KY))
+             KX = KX + INCX
+             KY = KY + INCY
+          END DO
+       END IF
        RETURN
        END