Go to the source code of this file.
|
double precision function | ddot (N, DX, INCX, DY, INCY) |
| DDOT More...
|
|
double precision function ddot |
( |
integer |
N, |
|
|
double precision, dimension(*) |
DX, |
|
|
integer |
INCX, |
|
|
double precision, dimension(*) |
DY, |
|
|
integer |
INCY |
|
) |
| |
DDOT
- Purpose:
DDOT forms the dot product of two vectors.
uses unrolled loops for increments equal to one.
- Author
- Univ. of Tennessee
-
Univ. of California Berkeley
-
Univ. of Colorado Denver
-
NAG Ltd.
- Date
- November 2011
- Further Details:
jack dongarra, linpack, 3/11/78.
modified 12/3/93, array(1) declarations changed to array(*)
Definition at line 53 of file ddot.f.
63 DOUBLE PRECISION dx(*),dy(*)
69 DOUBLE PRECISION dtemp
78 IF (incx.EQ.1 .AND. incy.EQ.1)
THEN
88 dtemp = dtemp + dx(i)*dy(i)
97 dtemp = dtemp + dx(i)*dy(i) + dx(i+1)*dy(i+1) +
98 $ dx(i+2)*dy(i+2) + dx(i+3)*dy(i+3) + dx(i+4)*dy(i+4)
107 IF (incx.LT.0) ix = (-n+1)*incx + 1
108 IF (incy.LT.0) iy = (-n+1)*incy + 1
110 dtemp = dtemp + dx(ix)*dy(iy)
double precision function ddot(N, DX, INCX, DY, INCY)
DDOT