496 {
501
504
507 LDA = *k+1;
508 LDB = *k+1;
511 for( i=0; i<*n; i++ )
512 for( j=0; j<*k; j++ ) {
517 }
518 }
519 else {
520 LDA = *n+1;
521 LDB = *n+1;
524 for( i=0; i<*k; i++ )
525 for( j=0; j<*n; j++ ){
530 }
531 }
532 LDC = *n+1;
534 for( i=0; i<*n; i++ )
535 for( j=0; j<*n; j++ ) {
538 }
540 B, LDB, beta, C, LDC );
541 for( j=0; j<*n; j++ )
542 for( i=0; i<*n; i++ ) {
545 }
546 free(A);
547 free(B);
548 free(C);
549 }
552 b, *ldb, beta, c, *ldc );
553 else
555 b, *ldb, beta, c, *ldc );
556}
void cblas_csyr2k(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE Trans, const CBLAS_INT N, const CBLAS_INT K, const void *alpha, const void *A, const CBLAS_INT lda, const void *B, const CBLAS_INT ldb, const void *beta, void *C, const CBLAS_INT ldc)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)