164 {
165 double *A;
168
170
173 A = (
double* )malloc( (*n)*
LDA*
sizeof( double ) );
174 for( i=0; i<*n; i++ )
175 for( j=0; j<*n; j++ )
176 A[
LDA*i+j ]=a[ (*lda)*j+i ];
178 for( i=0; i<*n; i++ )
179 for( j=0; j<*n; j++ )
180 a[ (*lda)*j+i ]=A[
LDA*i+j ];
181 free(A);
182 }
183 else
185}
void cblas_dsyr2(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INT N, const double alpha, const double *X, const CBLAS_INT incX, const double *Y, const CBLAS_INT incY, double *A, const CBLAS_INT lda)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)