 |
LAPACK 3.3.1
Linear Algebra PACKage
|
|
LAPACK 3.3.1
Linear Algebra PACKage
|
00001 REAL FUNCTION SDOT(N,SX,INCX,SY,INCY) 00002 * .. Scalar Arguments .. 00003 INTEGER INCX,INCY,N 00004 * .. 00005 * .. Array Arguments .. 00006 REAL SX(*),SY(*) 00007 * .. 00008 * 00009 * Purpose 00010 * ======= 00011 * 00012 * SDOT forms the dot product of two vectors. 00013 * uses unrolled loops for increments equal to one. 00014 * 00015 * Further Details 00016 * =============== 00017 * 00018 * jack dongarra, linpack, 3/11/78. 00019 * modified 12/3/93, array(1) declarations changed to array(*) 00020 * 00021 * ===================================================================== 00022 * 00023 * .. Local Scalars .. 00024 REAL STEMP 00025 INTEGER I,IX,IY,M,MP1 00026 * .. 00027 * .. Intrinsic Functions .. 00028 INTRINSIC MOD 00029 * .. 00030 STEMP = 0.0e0 00031 SDOT = 0.0e0 00032 IF (N.LE.0) RETURN 00033 IF (INCX.EQ.1 .AND. INCY.EQ.1) THEN 00034 * 00035 * code for both increments equal to 1 00036 * 00037 * 00038 * clean-up loop 00039 * 00040 M = MOD(N,5) 00041 IF (M.NE.0) THEN 00042 DO I = 1,M 00043 STEMP = STEMP + SX(I)*SY(I) 00044 END DO 00045 IF (N.LT.5) THEN 00046 SDOT=STEMP 00047 RETURN 00048 END IF 00049 END IF 00050 MP1 = M + 1 00051 DO I = MP1,N,5 00052 STEMP = STEMP + SX(I)*SY(I) + SX(I+1)*SY(I+1) + 00053 $ SX(I+2)*SY(I+2) + SX(I+3)*SY(I+3) + SX(I+4)*SY(I+4) 00054 END DO 00055 ELSE 00056 * 00057 * code for unequal increments or equal increments 00058 * not equal to 1 00059 * 00060 IX = 1 00061 IY = 1 00062 IF (INCX.LT.0) IX = (-N+1)*INCX + 1 00063 IF (INCY.LT.0) IY = (-N+1)*INCY + 1 00064 DO I = 1,N 00065 STEMP = STEMP + SX(IX)*SY(IY) 00066 IX = IX + INCX 00067 IY = IY + INCY 00068 END DO 00069 END IF 00070 SDOT = STEMP 00071 RETURN 00072 END
1.7.3 
