LAPACK 3.12.1
LAPACK: Linear Algebra PACKage
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages

◆ F77_csyr2k()

void F77_csyr2k ( CBLAS_INT * layout,
char * uplow,
char * transp,
CBLAS_INT * n,
CBLAS_INT * k,
CBLAS_TEST_COMPLEX * alpha,
CBLAS_TEST_COMPLEX * a,
CBLAS_INT * lda,
CBLAS_TEST_COMPLEX * b,
CBLAS_INT * ldb,
CBLAS_TEST_COMPLEX * beta,
CBLAS_TEST_COMPLEX * c,
CBLAS_INT * ldc )

Definition at line 489 of file c_cblas3.c.

496 {
497 CBLAS_INT i,j,LDA,LDB,LDC;
498 CBLAS_TEST_COMPLEX *A, *B, *C;
499 CBLAS_UPLO uplo;
500 CBLAS_TRANSPOSE trans;
501
502 get_uplo_type(uplow,&uplo);
503 get_transpose_type(transp,&trans);
504
505 if (*layout == TEST_ROW_MJR) {
506 if (trans == CblasNoTrans) {
507 LDA = *k+1;
508 LDB = *k+1;
509 A=(CBLAS_TEST_COMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_COMPLEX));
510 B=(CBLAS_TEST_COMPLEX* )malloc((*n)*LDB*sizeof(CBLAS_TEST_COMPLEX));
511 for( i=0; i<*n; i++ )
512 for( j=0; j<*k; j++ ) {
513 A[i*LDA+j].real=a[j*(*lda)+i].real;
514 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
515 B[i*LDB+j].real=b[j*(*ldb)+i].real;
516 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
517 }
518 }
519 else {
520 LDA = *n+1;
521 LDB = *n+1;
522 A=(CBLAS_TEST_COMPLEX* )malloc(LDA*(*k)*sizeof(CBLAS_TEST_COMPLEX));
523 B=(CBLAS_TEST_COMPLEX* )malloc(LDB*(*k)*sizeof(CBLAS_TEST_COMPLEX));
524 for( i=0; i<*k; i++ )
525 for( j=0; j<*n; j++ ){
526 A[i*LDA+j].real=a[j*(*lda)+i].real;
527 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
528 B[i*LDB+j].real=b[j*(*ldb)+i].real;
529 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
530 }
531 }
532 LDC = *n+1;
533 C=(CBLAS_TEST_COMPLEX* )malloc( (*n)*LDC*sizeof(CBLAS_TEST_COMPLEX));
534 for( i=0; i<*n; i++ )
535 for( j=0; j<*n; j++ ) {
536 C[i*LDC+j].real=c[j*(*ldc)+i].real;
537 C[i*LDC+j].imag=c[j*(*ldc)+i].imag;
538 }
539 cblas_csyr2k(CblasRowMajor, uplo, trans, *n, *k, alpha, A, LDA,
540 B, LDB, beta, C, LDC );
541 for( j=0; j<*n; j++ )
542 for( i=0; i<*n; i++ ) {
543 c[j*(*ldc)+i].real=C[i*LDC+j].real;
544 c[j*(*ldc)+i].imag=C[i*LDC+j].imag;
545 }
546 free(A);
547 free(B);
548 free(C);
549 }
550 else if (*layout == TEST_COL_MJR)
551 cblas_csyr2k(CblasColMajor, uplo, trans, *n, *k, alpha, a, *lda,
552 b, *ldb, beta, c, *ldc );
553 else
554 cblas_csyr2k(UNDEFINED, uplo, trans, *n, *k, alpha, a, *lda,
555 b, *ldb, beta, c, *ldc );
556}
#define UNDEFINED
Definition c_cblas3.c:12
#define TEST_ROW_MJR
Definition c_cblas3.c:11
#define TEST_COL_MJR
Definition c_cblas3.c:10
CBLAS_UPLO
Definition cblas.h:41
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)
CBLAS_TRANSPOSE
Definition cblas.h:40
@ CblasNoTrans
Definition cblas.h:40
@ CblasColMajor
Definition cblas.h:39
@ CblasRowMajor
Definition cblas.h:39
#define CBLAS_INT
Definition cblas.h:24
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
Definition auxiliary.c:18
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)
Definition auxiliary.c:8
Here is the call graph for this function: