rdot.F

Go to the documentation of this file.

      FUNCTION rdot(N,X,INCX,Y,INCY)
      USE parkind1, ONLY : jpim, jprb
      IMPLICIT NONE
      INTEGER(KIND=JPIM) :: N
      INTEGER(KIND=JPIM) :: INCX
      INTEGER(KIND=JPIM) :: INCY
      REAL(KIND=JPRB) :: X(n*incx)
      REAL(KIND=JPRB) :: Y(n*incy)
      REAL(KIND=JPRB) :: RDOT
      INTEGER :: IX,IY,J
      REAL(KIND=JPRB) :: ZDOT
      zdot=0.0_jprb
      ix=1
      iy=1
      DO j=1,n
        zdot=zdot+x(ix)*y(iy)
        ix=ix+incx
        iy=iy+incy
      ENDDO
      rdot=zdot
      RETURN
      ENDFUNCTION rdot