452 A = (
float* )malloc( LDA*LDA*
sizeof(
float ) );
453 AP = (
float* )malloc( (((LDA+1)*LDA)/2)*
sizeof(
float ) );
455 for( j=0, k=0; j<*n; j++ )
456 for( i=0; i<j+1; i++, k++ )
457 A[ LDA*i+j ]=ap[ k ];
458 for( i=0, k=0; i<*n; i++ )
459 for( j=i; j<*n; j++, k++ )
460 AP[ k ]=A[ LDA*i+j ];
464 for( j=0, k=0; j<*n; j++ )
465 for( i=j; i<*n; i++, k++ )
466 A[ LDA*i+j ]=ap[ k ];
467 for( i=0, k=0; i<*n; i++ )
468 for( j=0; j<i+1; j++, k++ )
469 AP[ k ]=A[ LDA*i+j ];
void get_transpose_type(char *type, CBLAS_TRANSPOSE *trans)
void get_diag_type(char *type, CBLAS_DIAG *diag)
void cblas_stpsv(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_uplo_type(char *type, CBLAS_UPLO *uplo)