LAPACK 3.12.1
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches

◆ F77_ztrsm()

void F77_ztrsm ( CBLAS_INT * layout,
char * rtlf,
char * uplow,
char * transp,
char * diagn,
CBLAS_INT * m,
CBLAS_INT * n,
CBLAS_TEST_ZOMPLEX * alpha,
CBLAS_TEST_ZOMPLEX * a,
CBLAS_INT * lda,
CBLAS_TEST_ZOMPLEX * b,
CBLAS_INT * ldb )

Definition at line 622 of file c_zblas3.c.

628 {
629 CBLAS_INT i,j,LDA,LDB;
630 CBLAS_TEST_ZOMPLEX *A, *B;
631 CBLAS_SIDE side;
632 CBLAS_DIAG diag;
633 CBLAS_UPLO uplo;
634 CBLAS_TRANSPOSE trans;
635
636 get_uplo_type(uplow,&uplo);
637 get_transpose_type(transp,&trans);
638 get_diag_type(diagn,&diag);
639 get_side_type(rtlf,&side);
640
641 if (*layout == TEST_ROW_MJR) {
642 if (side == CblasLeft) {
643 LDA = *m+1;
644 A=(CBLAS_TEST_ZOMPLEX* )malloc( (*m)*LDA*sizeof(CBLAS_TEST_ZOMPLEX ) );
645 for( i=0; i<*m; i++ )
646 for( j=0; j<*m; j++ ) {
647 A[i*LDA+j].real=a[j*(*lda)+i].real;
648 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
649 }
650 }
651 else{
652 LDA = *n+1;
653 A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
654 for( i=0; i<*n; i++ )
655 for( j=0; j<*n; j++ ) {
656 A[i*LDA+j].real=a[j*(*lda)+i].real;
657 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
658 }
659 }
660 LDB = *n+1;
661 B=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDB*sizeof(CBLAS_TEST_ZOMPLEX));
662 for( i=0; i<*m; i++ )
663 for( j=0; j<*n; j++ ) {
664 B[i*LDB+j].real=b[j*(*ldb)+i].real;
665 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
666 }
667 cblas_ztrsm(CblasRowMajor, side, uplo, trans, diag, *m, *n, alpha,
668 A, LDA, B, LDB );
669 for( j=0; j<*n; j++ )
670 for( i=0; i<*m; i++ ) {
671 b[j*(*ldb)+i].real=B[i*LDB+j].real;
672 b[j*(*ldb)+i].imag=B[i*LDB+j].imag;
673 }
674 free(A);
675 free(B);
676 }
677 else if (*layout == TEST_COL_MJR)
678 cblas_ztrsm(CblasColMajor, side, uplo, trans, diag, *m, *n, alpha,
679 a, *lda, b, *ldb);
680 else
681 cblas_ztrsm(UNDEFINED, side, uplo, trans, diag, *m, *n, alpha,
682 a, *lda, b, *ldb);
683}
#define UNDEFINED
Definition c_zblas3.c:13
#define TEST_ROW_MJR
Definition c_zblas3.c:12
#define TEST_COL_MJR
Definition c_zblas3.c:11
CBLAS_UPLO
Definition cblas.h:41
CBLAS_TRANSPOSE
Definition cblas.h:40
CBLAS_SIDE
Definition cblas.h:43
@ CblasLeft
Definition cblas.h:43
@ CblasColMajor
Definition cblas.h:39
@ CblasRowMajor
Definition cblas.h:39
CBLAS_DIAG
Definition cblas.h:42
#define CBLAS_INT
Definition cblas.h:24
void cblas_ztrsm(CBLAS_LAYOUT layout, CBLAS_SIDE Side, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT M, const CBLAS_INT N, const void *alpha, const void *A, const CBLAS_INT lda, void *B, const CBLAS_INT ldb)
Definition cblas_ztrsm.c:12
void get_diag_type(char *type, CBLAS_DIAG *diag)
Definition auxiliary.c:25
void get_side_type(char *type, CBLAS_SIDE *side)
Definition auxiliary.c:32
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: