*
* $Id: vdotn2.F,v 1.1.1.1 1996/02/15 17:53:07 mclareni Exp $
*
* $Log: vdotn2.F,v $
* Revision 1.1.1.1  1996/02/15 17:53:07  mclareni
* Kernlib
*
*
#include "sys/CERNLIB_machine.h"
#include "pilot.h"
@PROCESS DIRECTIVE ('*VDIR:') OPT(3) VECTOR
      FUNCTION VDOTN2 (X,Y,N)
C
C CERN PROGLIB# F121    VDOTN2          .VERSION KERNIBM  2.27  890622
C ORIG. 20/06/72
C 28/02/89  M.Roethlisberger            Optimization/Vectorization
C
      DIMENSION X(*),Y(*)
      REAL*8    XX8, XY8, YY8
      PARAMETER (LVMIN =     30)
C
      IF (N.GE.LVMIN) THEN
C*VDIR: PREFER VECTOR

         XX = 0.
         XY = 0.
         YY = 0.
C
         DO 90 J=1,N
            XX = XX + X(J)*X(J)
            XY = XY + X(J)*Y(J)
   90       YY = YY + Y(J)*Y(J)
C
         VDOTN2= XY**2/(XX*YY)

      ELSE IF (N.GT.0) THEN
C
C*VDIR: PREFER SCALAR

         XX8 = 0.
         XY8 = 0.
         YY8 = 0.
C
         DO 95 J=1,N
            XX8 = XX8 + DBLE(X(J))*DBLE(X(J))
            XY8 = XY8 + DBLE(X(J))*DBLE(Y(J))
   95       YY8 = YY8 + DBLE(Y(J))*DBLE(Y(J))
C

         XX=XX8
         XY=XY8
         YY=YY8
         VDOTN2= XY**2/(XX*YY)

      ELSE

         VDOTN2 = 0.0

      ENDIF
      END
#ifdef CERNLIB_TCGEN_VDOTN2
#undef CERNLIB_TCGEN_VDOTN2
#endif