17 {
18
19 float *A;
22
25 LDA = *n+1;
26 A = ( float* )malloc( (*m)*LDA*sizeof( float ) );
27 for( i=0; i<*m; i++ )
28 for( j=0; j<*n; j++ )
29 A[ LDA*i+j ]=a[ (*lda)*j+i ];
31 *m, *n, *alpha, A, LDA, x, *incx, *beta, y, *incy );
32 free(A);
33 }
36 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy );
37 else
39 *m, *n, *alpha, a, *lda, x, *incx, *beta, y, *incy );
40}
void cblas_sgemv(const CBLAS_LAYOUT layout, const CBLAS_TRANSPOSE TransA, const CBLAS_INT M, const CBLAS_INT N, const float alpha, const float *A, const CBLAS_INT lda, const float *X, const CBLAS_INT incX, const float beta, float *Y, const CBLAS_INT incY)
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)