413 A = (
float* )malloc( LDA*LDA*
sizeof(
float ) );
414 AP = (
float* )malloc( (((LDA+1)*LDA)/2)*
sizeof(
float ) );
416 for( j=0, k=0; j<*n; j++ )
417 for( i=0; i<j+1; i++, k++ )
418 A[ LDA*i+j ]=ap[ k ];
419 for( i=0, k=0; i<*n; i++ )
420 for( j=i; j<*n; j++, k++ )
421 AP[ k ]=A[ LDA*i+j ];
424 for( j=0, k=0; j<*n; j++ )
425 for( i=j; i<*n; i++, k++ )
426 A[ LDA*i+j ]=ap[ k ];
427 for( i=0, k=0; i<*n; i++ )
428 for( j=0; j<i+1; j++, k++ )
429 AP[ k ]=A[ LDA*i+j ];
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)
void cblas_stpmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, CBLAS_TRANSPOSE TransA, CBLAS_DIAG Diag, const int N, const float *Ap, float *X, const int incX)
void get_diag_type(char *type, CBLAS_DIAG *diag)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)