232 {
234 float *A, *B;
239
244
247 LDA = *m+1;
248 A = ( float* )malloc( (*m)*LDA*sizeof( float ) );
249 for( i=0; i<*m; i++ )
250 for( j=0; j<*m; j++ )
251 A[i*LDA+j]=a[j*(*lda)+i];
252 }
253 else{
254 LDA = *n+1;
255 A = ( float* )malloc( (*n)*LDA*sizeof( float ) );
256 for( i=0; i<*n; i++ )
257 for( j=0; j<*n; j++ )
258 A[i*LDA+j]=a[j*(*lda)+i];
259 }
260 LDB = *n+1;
261 B = ( float* )malloc( (*m)*LDB*sizeof( float ) );
262 for( i=0; i<*m; i++ )
263 for( j=0; j<*n; j++ )
264 B[i*LDB+j]=b[j*(*ldb)+i];
266 A, LDA, B, LDB );
267 for( j=0; j<*n; j++ )
268 for( i=0; i<*m; i++ )
269 b[j*(*ldb)+i]=B[i*LDB+j];
270 free(A);
271 free(B);
272 }
275 a, *lda, b, *ldb);
276 else
278 a, *lda, b, *ldb);
279}
void cblas_strmm(CBLAS_LAYOUT layout, CBLAS_SIDE Side, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const CBLAS_INT M, const CBLAS_INT N, const float alpha, const float *A, const CBLAS_INT lda, float *B, const CBLAS_INT ldb)
void get_diag_type(char *type, CBLAS_DIAG *diag)
void get_side_type(char *type, CBLAS_SIDE *side)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)