160 {
161 float *A;
164
166
168 LDA = *n+1;
169 A = ( float* )malloc( (*n)*LDA*sizeof( float ) );
170 for( i=0; i<*n; i++ )
171 for( j=0; j<*n; j++ )
172 A[ LDA*i+j ]=a[ (*lda)*j+i ];
174 for( i=0; i<*n; i++ )
175 for( j=0; j<*n; j++ )
176 a[ (*lda)*j+i ]=A[ LDA*i+j ];
177 free(A);
178 }
179 else
181}
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)