252 {
253
258
261
264 LDA = *m+1;
266 for( i=0; i<*m; i++ )
267 for( j=0; j<*m; j++ )
268 A[i*LDA+j]=a[j*(*lda)+i];
269 }
270 else{
271 LDA = *n+1;
273 for( i=0; i<*n; i++ )
274 for( j=0; j<*n; j++ )
275 A[i*LDA+j]=a[j*(*lda)+i];
276 }
277 LDB = *n+1;
279 for( i=0; i<*m; i++ )
280 for( j=0; j<*n; j++ )
281 B[i*LDB+j]=b[j*(*ldb)+i];
282 LDC = *n+1;
284 for( j=0; j<*n; j++ )
285 for( i=0; i<*m; i++ )
286 C[i*LDC+j]=c[j*(*ldc)+i];
287 cblas_zsymm(
CblasRowMajor, side, uplo, *m, *n, alpha, A, LDA, B, LDB,
288 beta, C, LDC );
289 for( j=0; j<*n; j++ )
290 for( i=0; i<*m; i++ )
291 c[j*(*ldc)+i]=C[i*LDC+j];
292 free(A);
293 free(B);
294 free(C);
295 }
297 cblas_zsymm(
CblasColMajor, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
298 beta, c, *ldc );
299 else
300 cblas_zsymm(
UNDEFINED, side, uplo, *m, *n, alpha, a, *lda, b, *ldb,
301 beta, c, *ldc );
302}
void cblas_zsymm(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)
void get_side_type(char *type, CBLAS_SIDE *side)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)