    /*
     * cblas_ddot.c
     *
     * The program is a C interface to ddot.
     * It calls the fortran wrapper before calling ddot.
     *
     * Written by Keita Teranishi.  2/11/1998
     *
     */
    #include "cblas.h"
    #include "cblas_f77.h"
    double cblas_ddot( const CBLAS_INT_TYPE N, const double *X,
                          const CBLAS_INT_TYPE incX, const double *Y, const CBLAS_INT_TYPE incY)
    {
       double dot;
    #ifdef F77_INT
       F77_INT F77_N=N, F77_incX=incX, F77_incY=incY;
    #else
       #define F77_N N
       #define F77_incX incX
       #define F77_incY incY
    #endif
       F77_ddot_sub( &F77_N, X, &F77_incX, Y, &F77_incY, &dot);
       return dot;
    }