177 {
179 float *A, *B, *C;
182
185
188 LDA = *k+1;
189 LDB = *k+1;
190 A = ( float* )malloc( (*n)*LDA*sizeof( float ) );
191 B = ( float* )malloc( (*n)*LDB*sizeof( float ) );
192 for( i=0; i<*n; i++ )
193 for( j=0; j<*k; j++ ) {
194 A[i*LDA+j]=a[j*(*lda)+i];
195 B[i*LDB+j]=b[j*(*ldb)+i];
196 }
197 }
198 else {
199 LDA = *n+1;
200 LDB = *n+1;
201 A = ( float* )malloc( LDA*(*k)*sizeof( float ) );
202 B = ( float* )malloc( LDB*(*k)*sizeof( float ) );
203 for( i=0; i<*k; i++ )
204 for( j=0; j<*n; j++ ){
205 A[i*LDA+j]=a[j*(*lda)+i];
206 B[i*LDB+j]=b[j*(*ldb)+i];
207 }
208 }
209 LDC = *n+1;
210 C = ( float* )malloc( (*n)*LDC*sizeof( float ) );
211 for( i=0; i<*n; i++ )
212 for( j=0; j<*n; j++ )
213 C[i*LDC+j]=c[j*(*ldc)+i];
215 B, LDB, *beta, C, LDC );
216 for( j=0; j<*n; j++ )
217 for( i=0; i<*n; i++ )
218 c[j*(*ldc)+i]=C[i*LDC+j];
219 free(A);
220 free(B);
221 free(C);
222 }
225 b, *ldb, *beta, c, *ldc );
226 else
228 b, *ldb, *beta, c, *ldc );
229}
void cblas_ssyr2k(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE Trans, const CBLAS_INT N, const CBLAS_INT K, const float alpha, const float *A, const CBLAS_INT lda, const float *B, const CBLAS_INT ldb, const float beta, float *C, const CBLAS_INT ldc)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)