13 {
14
15 double *A;
18
21 LDA = *n+1;
22 A = ( double* )malloc( (*m)*LDA*sizeof( double ) );
23 for( i=0; i<*m; i++ )
24 for( j=0; j<*n; j++ )
25 A[ LDA*i+j ]=a[ (*lda)*j+i ];
27 *m, *n, *alpha, A, LDA, x, *incx, *beta, y, *incy );
28 free(A);
29 }
32 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy );
33 else
35 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy );
36}
void cblas_dgemv(CBLAS_LAYOUT layout, CBLAS_TRANSPOSE TransA, const CBLAS_INT M, const CBLAS_INT N, const double alpha, const double *A, const CBLAS_INT lda, const double *X, const CBLAS_INT incX, const double beta, double *Y, const CBLAS_INT incY)
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)