142 A = (
float* )malloc( (*n)*LDA*
sizeof( float ) );
143 for( i=0; i<*n; i++ )
144 for( j=0; j<*k; j++ )
145 A[i*LDA+j]=a[j*(*lda)+i];
149 A = (
float* )malloc( (*k)*LDA*
sizeof( float ) );
150 for( i=0; i<*k; i++ )
151 for( j=0; j<*n; j++ )
152 A[i*LDA+j]=a[j*(*lda)+i];
155 C = (
float* )malloc( (*n)*LDC*
sizeof( float ) );
156 for( i=0; i<*n; i++ )
157 for( j=0; j<*n; j++ )
158 C[i*LDC+j]=c[j*(*ldc)+i];
161 for( j=0; j<*n; j++ )
162 for( i=0; i<*n; i++ )
163 c[j*(*ldc)+i]=C[i*LDC+j];
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)
void cblas_ssyrk(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE Trans, const int N, const int K, const float alpha, const float *A, const int lda, const float beta, float *C, const int ldc)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)