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

◆ F77_cher2k()

void F77_cher2k ( 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,
float * beta,
CBLAS_TEST_COMPLEX * c,
CBLAS_INT * ldc )

Definition at line 421 of file c_cblas3.c.

428 {
429 CBLAS_INT i,j,LDA,LDB,LDC;
430 CBLAS_TEST_COMPLEX *A, *B, *C;
431 CBLAS_UPLO uplo;
432 CBLAS_TRANSPOSE trans;
433
434 get_uplo_type(uplow,&uplo);
435 get_transpose_type(transp,&trans);
436
437 if (*layout == TEST_ROW_MJR) {
438 if (trans == CblasNoTrans) {
439 LDA = *k+1;
440 LDB = *k+1;
441 A=(CBLAS_TEST_COMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_COMPLEX ));
442 B=(CBLAS_TEST_COMPLEX* )malloc((*n)*LDB*sizeof(CBLAS_TEST_COMPLEX ));
443 for( i=0; i<*n; i++ )
444 for( j=0; j<*k; j++ ) {
445 A[i*LDA+j].real=a[j*(*lda)+i].real;
446 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
447 B[i*LDB+j].real=b[j*(*ldb)+i].real;
448 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
449 }
450 }
451 else {
452 LDA = *n+1;
453 LDB = *n+1;
454 A=(CBLAS_TEST_COMPLEX* )malloc( LDA*(*k)*sizeof(CBLAS_TEST_COMPLEX ) );
455 B=(CBLAS_TEST_COMPLEX* )malloc( LDB*(*k)*sizeof(CBLAS_TEST_COMPLEX ) );
456 for( i=0; i<*k; i++ )
457 for( j=0; j<*n; j++ ){
458 A[i*LDA+j].real=a[j*(*lda)+i].real;
459 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
460 B[i*LDB+j].real=b[j*(*ldb)+i].real;
461 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
462 }
463 }
464 LDC = *n+1;
465 C=(CBLAS_TEST_COMPLEX* )malloc( (*n)*LDC*sizeof(CBLAS_TEST_COMPLEX ) );
466 for( i=0; i<*n; i++ )
467 for( j=0; j<*n; j++ ) {
468 C[i*LDC+j].real=c[j*(*ldc)+i].real;
469 C[i*LDC+j].imag=c[j*(*ldc)+i].imag;
470 }
471 cblas_cher2k(CblasRowMajor, uplo, trans, *n, *k, alpha, A, LDA,
472 B, LDB, *beta, C, LDC );
473 for( j=0; j<*n; j++ )
474 for( i=0; i<*n; i++ ) {
475 c[j*(*ldc)+i].real=C[i*LDC+j].real;
476 c[j*(*ldc)+i].imag=C[i*LDC+j].imag;
477 }
478 free(A);
479 free(B);
480 free(C);
481 }
482 else if (*layout == TEST_COL_MJR)
483 cblas_cher2k(CblasColMajor, uplo, trans, *n, *k, alpha, a, *lda,
484 b, *ldb, *beta, c, *ldc );
485 else
486 cblas_cher2k(UNDEFINED, uplo, trans, *n, *k, alpha, a, *lda,
487 b, *ldb, *beta, c, *ldc );
488}
#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
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)
CBLAS_UPLO
Definition cblas.h:41
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: