LAPACK  3.10.0 LAPACK: Linear Algebra PACKage

◆ F77_dsymm()

 void F77_dsymm ( int * layout, char * rtlf, char * uplow, int * m, int * n, double * alpha, double * a, int * lda, double * b, int * ldb, double * beta, double * c, int * ldc )

Definition at line 76 of file c_dblas3.c.

78  {
79
80  double *A, *B, *C;
81  int i,j,LDA, LDB, LDC;
82  CBLAS_UPLO uplo;
83  CBLAS_SIDE side;
84
85  get_uplo_type(uplow,&uplo);
86  get_side_type(rtlf,&side);
87
88  if (*layout == TEST_ROW_MJR) {
89  if (side == CblasLeft) {
90  LDA = *m+1;
91  A = ( double* )malloc( (*m)*LDA*sizeof( double ) );
92  for( i=0; i<*m; i++ )
93  for( j=0; j<*m; j++ )
94  A[i*LDA+j]=a[j*(*lda)+i];
95  }
96  else{
97  LDA = *n+1;
98  A = ( double* )malloc( (*n)*LDA*sizeof( double ) );
99  for( i=0; i<*n; i++ )
100  for( j=0; j<*n; j++ )
101  A[i*LDA+j]=a[j*(*lda)+i];
102  }
103  LDB = *n+1;
104  B = ( double* )malloc( (*m)*LDB*sizeof( double ) );
105  for( i=0; i<*m; i++ )
106  for( j=0; j<*n; j++ )
107  B[i*LDB+j]=b[j*(*ldb)+i];
108  LDC = *n+1;
109  C = ( double* )malloc( (*m)*LDC*sizeof( double ) );
110  for( j=0; j<*n; j++ )
111  for( i=0; i<*m; i++ )
112  C[i*LDC+j]=c[j*(*ldc)+i];
113  cblas_dsymm( CblasRowMajor, side, uplo, *m, *n, *alpha, A, LDA, B, LDB,
114  *beta, C, LDC );
115  for( j=0; j<*n; j++ )
116  for( i=0; i<*m; i++ )
117  c[j*(*ldc)+i]=C[i*LDC+j];
118  free(A);
119  free(B);
120  free(C);
121  }
122  else if (*layout == TEST_COL_MJR)
123  cblas_dsymm( CblasColMajor, side, uplo, *m, *n, *alpha, a, *lda, b, *ldb,
124  *beta, c, *ldc );
125  else
126  cblas_dsymm( UNDEFINED, side, uplo, *m, *n, *alpha, a, *lda, b, *ldb,
127  *beta, c, *ldc );
128 }
#define UNDEFINED
Definition: c_dblas3.c:12
#define TEST_ROW_MJR
Definition: c_dblas3.c:11
#define TEST_COL_MJR
Definition: c_dblas3.c:10
CBLAS_UPLO
Definition: cblas.h:21
CBLAS_SIDE
Definition: cblas.h:23
@ CblasLeft
Definition: cblas.h:23
@ CblasColMajor
Definition: cblas.h:19
@ CblasRowMajor
Definition: cblas.h:19
void cblas_dsymm(CBLAS_LAYOUT layout, CBLAS_SIDE Side, CBLAS_UPLO Uplo, const CBLAS_INDEX M, const CBLAS_INDEX N, const double alpha, const double *A, const CBLAS_INDEX lda, const double *B, const CBLAS_INDEX ldb, const double beta, double *C, const CBLAS_INDEX ldc)
Definition: cblas_dsymm.c:12
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
#define LDB
Definition: example_user.c:13
#define LDA
Definition: example_user.c:12
Here is the call graph for this function:
Here is the caller graph for this function: