void accel_multMatrixVec(double* m1,double* v, int ndata,int npol, double* out){
int m,n,k;
double alpha=1.0,beta=0;
int inc=1;
m=npol;
n=ndata;
/*
* An important note about this. FORTRAN effectively transposes all matricies becaue
* the memory is ordered differently. Therefore we specify "T" to do a transpose
*
*
* M.Keith 2013.
*/
F77_dgemv("T",&m,&n,&alpha ,m1,&m,v,&inc,&beta,out,&inc);
/*
m=ndata;
n=npol;
k=ndata;
F77_dgemm("N","T",&m,&n,&k,&alpha,m1,&m,m2,&n,&beta,out,&m);
*/
}
void cblas_dgemv(enum CBLAS_ORDER order,
enum CBLAS_TRANSPOSE TransA, f77_int M, f77_int N,
double alpha, const double *A, f77_int lda,
const double *X, f77_int incX, double beta,
double *Y, f77_int incY)
{
char TA;
#ifdef F77_CHAR
F77_CHAR F77_TA;
#else
#define F77_TA &TA
#endif
#ifdef F77_INT
F77_INT F77_M=M, F77_N=N, F77_lda=lda, F77_incX=incX, F77_incY=incY;
#else
#define F77_M M
#define F77_N N
#define F77_lda lda
#define F77_incX incX
#define F77_incY incY
#endif
extern int CBLAS_CallFromC;
extern int RowMajorStrg;
RowMajorStrg = 0;
CBLAS_CallFromC = 1;
if (order == CblasColMajor)
{
if (TransA == CblasNoTrans) TA = 'N';
else if (TransA == CblasTrans) TA = 'T';
else if (TransA == CblasConjTrans) TA = 'C';
else
{
cblas_xerbla(2, "cblas_dgemv","Illegal TransA setting, %d\n", TransA);
CBLAS_CallFromC = 0;
RowMajorStrg = 0;
return;
}
#ifdef F77_CHAR
F77_TA = C2F_CHAR(&TA);
#endif
F77_dgemv(F77_TA, &F77_M, &F77_N, &alpha, A, &F77_lda, X, &F77_incX,
&beta, Y, &F77_incY);
}
else if (order == CblasRowMajor)
{
RowMajorStrg = 1;
if (TransA == CblasNoTrans) TA = 'T';
else if (TransA == CblasTrans) TA = 'N';
else if (TransA == CblasConjTrans) TA = 'N';
else
{
cblas_xerbla(2, "cblas_dgemv","Illegal TransA setting, %d\n", TransA);
CBLAS_CallFromC = 0;
RowMajorStrg = 0;
return;
}
#ifdef F77_CHAR
F77_TA = C2F_CHAR(&TA);
#endif
F77_dgemv(F77_TA, &F77_N, &F77_M, &alpha, A, &F77_lda, X,
&F77_incX, &beta, Y, &F77_incY);
}
else cblas_xerbla(1, "cblas_dgemv", "Illegal Order setting, %d\n", order);
CBLAS_CallFromC = 0;
RowMajorStrg = 0;
return;
}
