1 /*
2  * cblas_ddot.c
3  *
4  * The program is a C interface to ddot.
5  * It calls the fortran wrapper before calling ddot.
6  *
7  * Written by Keita Teranishi.  2/11/1998
8  *
9  */
10 #include "cblas.h"
11 #include "cblas_f77.h"
cblas_ddot(const CBLAS_INDEX N,const double * X,const CBLAS_INDEX incX,const double * Y,const CBLAS_INDEX incY)12 double cblas_ddot( const CBLAS_INDEX N, const double *X,
13                       const CBLAS_INDEX incX, const double *Y, const CBLAS_INDEX incY)
14 {
15    double dot;
16 #ifdef F77_INT
17    F77_INT F77_N=N, F77_incX=incX, F77_incY=incY;
18 #else
19    #define F77_N N
20    #define F77_incX incX
21    #define F77_incY incY
22 #endif
23    F77_ddot_sub( &F77_N, X, &F77_incX, Y, &F77_incY, &dot);
24    return dot;
25 }
26