140 {
141 float *A;
144
146
149 A = (
float* )malloc( (*n)*
LDA*
sizeof( float ) );
150 for( i=0; i<*n; i++ )
151 for( j=0; j<*n; j++ )
152 A[
LDA*i+j ]=a[ (*lda)*j+i ];
154 for( i=0; i<*n; i++ )
155 for( j=0; j<*n; j++ )
156 a[ (*lda)*j+i ]=A[
LDA*i+j ];
157 free(A);
158 }
159 else
161}
void cblas_ssyr(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const float alpha, const float *X, const CBLAS_INT incX, float *A, const CBLAS_INT lda)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)