de/d50/dvvdot_8f_source.html

00001       SUBROUTINE DVVDOT( N, DOT, X, INCX, Y, INCY )
00002 *
00003 *  -- PBLAS auxiliary routine (version 2.0) --
00004 *     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
00005 *     and University of California, Berkeley.
00006 *     April 1, 1998
00007 *
00008 *     .. Scalar Arguments ..
00009       INTEGER            INCX, INCY, N
00010       DOUBLE PRECISION   DOT
00011 *     ..
00012 *     .. Array Arguments ..
00013       DOUBLE PRECISION   X( * ), Y( * )
00014 *     ..
00015 *
00016 *  Purpose
00017 *  =======
00018 *
00019 *  DVVDOT computes the following dot product:
00020 *
00021 *     dot = dot + x**T * y,
00022 *
00023 *  where x and y are n vectors.
00024 *
00025 *  Arguments
00026 *  =========
00027 *
00028 *  N       (input) INTEGER
00029 *          On entry, N  specifies the length of the  vectors x and y.  N
00030 *          must be at least zero.
00031 *
00032 *  DOT     (input/output) DOUBLE PRECISION
00033 *          On exit, DOT is updated with the dot product of the vectors x
00034 *          and y.
00035 *
00036 *  X       (input) DOUBLE PRECISION array of dimension at least
00037 *          ( 1 + ( n - 1 )*abs( INCX ) ). Before entry,  the incremented
00038 *          array X must contain the vector x.
00039 *
00040 *  INCX    (input) INTEGER
00041 *          On entry, INCX specifies the increment for the elements of X.
00042 *          INCX must not be zero.
00043 *
00044 *  Y       (input) DOUBLE PRECISION array of dimension at least
00045 *          ( 1 + ( n - 1 )*abs( INCY ) ). Before entry,  the incremented
00046 *          array Y must contain the vector y.
00047 *
00048 *  INCY    (input) INTEGER
00049 *          On entry, INCY specifies the increment for the elements of Y.
00050 *          INCY must not be zero.
00051 *
00052 *  -- Written on April 1, 1998 by
00053 *     Antoine Petitet, University  of  Tennessee, Knoxville 37996, USA.
00054 *
00055 *  =====================================================================
00056 *
00057 *     .. External Functions ..
00058       DOUBLE PRECISION   DDOT
00059       EXTERNAL           DDOT
00060 *     ..
00061 *     .. Executable Statements ..
00062 *
00063       DOT = DOT + DDOT( N, X, INCX, Y, INCY )
00064 *
00065       RETURN
00066 *
00067 *     End of DVVDOT
00068 *
00069       END