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

◆ F77_chemm()

void F77_chemm ( CBLAS_INT * layout,
char * rtlf,
char * uplow,
CBLAS_INT * m,
CBLAS_INT * n,
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 175 of file c_cblas3.c.

182 {
183
184 CBLAS_TEST_COMPLEX *A, *B, *C;
185 CBLAS_INT i,j,LDA, LDB, LDC;
186 CBLAS_UPLO uplo;
187 CBLAS_SIDE side;
188
189 get_uplo_type(uplow,&uplo);
190 get_side_type(rtlf,&side);
191
192 if (*layout == TEST_ROW_MJR) {
193 if (side == CblasLeft) {
194 LDA = *m+1;
195 A= (CBLAS_TEST_COMPLEX* )malloc((*m)*LDA*sizeof(CBLAS_TEST_COMPLEX));
196 for( i=0; i<*m; i++ )
197 for( j=0; j<*m; j++ ) {
198 A[i*LDA+j].real=a[j*(*lda)+i].real;
199 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
200 }
201 }
202 else{
203 LDA = *n+1;
204 A=(CBLAS_TEST_COMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_COMPLEX ) );
205 for( i=0; i<*n; i++ )
206 for( j=0; j<*n; j++ ) {
207 A[i*LDA+j].real=a[j*(*lda)+i].real;
208 A[i*LDA+j].imag=a[j*(*lda)+i].imag;
209 }
210 }
211 LDB = *n+1;
212 B=(CBLAS_TEST_COMPLEX* )malloc( (*m)*LDB*sizeof(CBLAS_TEST_COMPLEX ) );
213 for( i=0; i<*m; i++ )
214 for( j=0; j<*n; j++ ) {
215 B[i*LDB+j].real=b[j*(*ldb)+i].real;
216 B[i*LDB+j].imag=b[j*(*ldb)+i].imag;
217 }
218 LDC = *n+1;
219 C=(CBLAS_TEST_COMPLEX* )malloc((*m)*LDC*sizeof(CBLAS_TEST_COMPLEX ) );
220 for( j=0; j<*n; j++ )
221 for( i=0; i<*m; i++ ) {
222 C[i*LDC+j].real=c[j*(*ldc)+i].real;
223 C[i*LDC+j].imag=c[j*(*ldc)+i].imag;
224 }
225 cblas_chemm( CblasRowMajor, side, uplo, *m, *n, alpha, A, LDA, B, LDB,
226 beta, C, LDC );
227 for( j=0; j<*n; j++ )
228 for( i=0; i<*m; i++ ) {
229 c[j*(*ldc)+i].real=C[i*LDC+j].real;
230 c[j*(*ldc)+i].imag=C[i*LDC+j].imag;
231 }
232 free(A);
233 free(B);
234 free(C);
235 }
236 else if (*layout == TEST_COL_MJR)
237 cblas_chemm( CblasColMajor, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
238 beta, c, *ldc );
239 else
240 cblas_chemm( UNDEFINED, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
241 beta, c, *ldc );
242}
#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
CBLAS_SIDE
Definition cblas.h:43
@ CblasLeft
Definition cblas.h:43
@ CblasColMajor
Definition cblas.h:39
@ CblasRowMajor
Definition cblas.h:39
void cblas_chemm(CBLAS_LAYOUT layout, CBLAS_SIDE Side, CBLAS_UPLO Uplo, const CBLAS_INT M, const CBLAS_INT N, 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)
Definition cblas_chemm.c:12
#define CBLAS_INT
Definition cblas.h:24
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
Here is the call graph for this function: