98 {
99
104
108
111 LDA = *k+1;
113 for( i=0; i<*n; i++ )
114 for( j=0; j<*k; j++ ) {
117 }
118 }
119 else {
120 LDA = *n+1;
122 for( i=0; i<*k; i++ )
123 for( j=0; j<*n; j++ ) {
126 }
127 }
128
130 LDB = *n+1;
132 for( i=0; i<*k; i++ )
133 for( j=0; j<*n; j++ ) {
136 }
137 }
138 else {
139 LDB = *k+1;
141 for( i=0; i<*n; i++ )
142 for( j=0; j<*k; j++ ) {
145 }
146 }
147
148 LDC = *n+1;
150 for( j=0; j<*n; j++ )
151 for( i=0; i<*n; i++ ) {
154 }
156 B, LDB, beta, C, LDC );
157 for( j=0; j<*n; j++ )
158 for( i=0; i<*n; i++ ) {
161 }
162 free(A);
163 free(B);
164 free(C);
165 }
168 b, *ldb, beta, c, *ldc );
169 else
171 b, *ldb, beta, c, *ldc );
172}
void cblas_cgemmtr(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_TRANSPOSE TransB, 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 void *beta, void *C, const CBLAS_INT ldc)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)