Sources Esope | Cast3M

sdot
C SDOT      SOURCE    CHAT      05/01/13    03:13:04     5004
      FUNCTION SDOT(N,V,M,W,L)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C     PRODUIT DE DEUX VECTEURS
      DIMENSION V(1),W(1)
      P=0.D0
      I1=1
      J1=1
      IF(M.EQ.1.AND.L.EQ.1) GO TO 2
      DO 1 I=1,N
      P = P + V(I1) * W(J1)
      I1=I1+M
      J1=J1+L
   1  CONTINUE
      SDOT = P
      RETURN
  2   DO 3 I = 1,N
   3  P = P + V(I) * W(I)
      SDOT = P
      RETURN
      END