LAPACK  3.6.1
LAPACK: Linear Algebra PACKage
void F77_zsymm ( int *  layout,
char *  rtlf,
char *  uplow,
int *  m,
int *  n,
CBLAS_TEST_ZOMPLEX alpha,
CBLAS_TEST_ZOMPLEX a,
int *  lda,
CBLAS_TEST_ZOMPLEX b,
int *  ldb,
CBLAS_TEST_ZOMPLEX beta,
CBLAS_TEST_ZOMPLEX c,
int *  ldc 
)

Definition at line 154 of file c_zblas3.c.

157  {
158 
159  CBLAS_TEST_ZOMPLEX *A, *B, *C;
160  int i,j,LDA, LDB, LDC;
161  CBLAS_UPLO uplo;
162  CBLAS_SIDE side;
163 
164  get_uplo_type(uplow,&uplo);
165  get_side_type(rtlf,&side);
166 
167  if (*layout == TEST_ROW_MJR) {
168  if (side == CblasLeft) {
169  LDA = *m+1;
170  A=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDA*sizeof(CBLAS_TEST_ZOMPLEX));
171  for( i=0; i<*m; i++ )
172  for( j=0; j<*m; j++ )
173  A[i*LDA+j]=a[j*(*lda)+i];
174  }
175  else{
176  LDA = *n+1;
177  A=(CBLAS_TEST_ZOMPLEX* )malloc((*n)*LDA*sizeof(CBLAS_TEST_ZOMPLEX ) );
178  for( i=0; i<*n; i++ )
179  for( j=0; j<*n; j++ )
180  A[i*LDA+j]=a[j*(*lda)+i];
181  }
182  LDB = *n+1;
183  B=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDB*sizeof(CBLAS_TEST_ZOMPLEX ));
184  for( i=0; i<*m; i++ )
185  for( j=0; j<*n; j++ )
186  B[i*LDB+j]=b[j*(*ldb)+i];
187  LDC = *n+1;
188  C=(CBLAS_TEST_ZOMPLEX* )malloc((*m)*LDC*sizeof(CBLAS_TEST_ZOMPLEX));
189  for( j=0; j<*n; j++ )
190  for( i=0; i<*m; i++ )
191  C[i*LDC+j]=c[j*(*ldc)+i];
192  cblas_zsymm( CblasRowMajor, side, uplo, *m, *n, alpha, A, LDA, B, LDB,
193  beta, C, LDC );
194  for( j=0; j<*n; j++ )
195  for( i=0; i<*m; i++ )
196  c[j*(*ldc)+i]=C[i*LDC+j];
197  free(A);
198  free(B);
199  free(C);
200  }
201  else if (*layout == TEST_COL_MJR)
202  cblas_zsymm( CblasColMajor, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
203  beta, c, *ldc );
204  else
205  cblas_zsymm( UNDEFINED, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
206  beta, c, *ldc );
207 }
void get_side_type(char *type, CBLAS_SIDE *side)
Definition: auxiliary.c:32
#define LDB
Definition: example_user.c:13
CBLAS_SIDE
Definition: cblas.h:23
#define UNDEFINED
Definition: c_zblas3.c:12
#define TEST_COL_MJR
Definition: c_zblas3.c:10
#define LDA
Definition: example_user.c:12
CBLAS_UPLO
Definition: cblas.h:21
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
Definition: auxiliary.c:18
void cblas_zsymm(CBLAS_LAYOUT layout, CBLAS_SIDE Side, CBLAS_UPLO Uplo, const int M, const int N, const void *alpha, const void *A, const int lda, const void *B, const int ldb, const void *beta, void *C, const int ldc)
Definition: cblas_zsymm.c:12
#define TEST_ROW_MJR
Definition: c_zblas3.c:11

Here is the call graph for this function: