#include <gsl/gsl_math.h>
#include <gsl/gsl_cblas.h>
#include "cblas.h"
#include "error_cblas_l2.h"

void
cblas_dspmv (const enum CBLAS_ORDER order, const enum CBLAS_UPLO Uplo,
             const int N, const double alpha, const double *Ap,
             const double *X, const int incX, const double beta, double *Y,
             const int incY)
{
#define BASE double
#include "source_spmv.h"
#undef BASE
}