64 {
65 double *A;
70
74
76 LDA = *n+1;
77 A = ( double* )malloc( (*n)*LDA*sizeof( double ) );
78 for( i=0; i<*n; i++ )
79 for( j=0; j<*n; j++ )
80 A[ LDA*i+j ]=a[ (*lda)*j+i ];
82 free(A);
83 }
86 else {
88 }
89}
void cblas_dtrmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT N, const double *A, const CBLAS_INT lda, double *X, const CBLAS_INT incX)
void get_diag_type(char *type, CBLAS_DIAG *diag)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)