428 {
433
436
439 LDA = *k+1;
440 LDB = *k+1;
443 for( i=0; i<*n; i++ )
444 for( j=0; j<*k; j++ ) {
449 }
450 }
451 else {
452 LDA = *n+1;
453 LDB = *n+1;
456 for( i=0; i<*k; i++ )
457 for( j=0; j<*n; j++ ){
462 }
463 }
464 LDC = *n+1;
466 for( i=0; i<*n; i++ )
467 for( j=0; j<*n; j++ ) {
470 }
472 B, LDB, *beta, C, LDC );
473 for( j=0; j<*n; j++ )
474 for( i=0; i<*n; i++ ) {
477 }
478 free(A);
479 free(B);
480 free(C);
481 }
484 b, *ldb, *beta, c, *ldc );
485 else
487 b, *ldb, *beta, c, *ldc );
488}
void cblas_cher2k(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 float 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)