372 A = (
float* )malloc( LDA*LDA*
sizeof(
float ) );
373 AP = (
float* )malloc( (((LDA+1)*LDA)/2)*
sizeof(
float ) );
375 for( j=0, k=0; j<*n; j++ )
376 for( i=0; i<j+1; i++, k++ )
377 A[ LDA*i+j ]=ap[ k ];
378 for( i=0, k=0; i<*n; i++ )
379 for( j=i; j<*n; j++, k++ )
380 AP[ k ]=A[ LDA*i+j ];
383 for( j=0, k=0; j<*n; j++ )
384 for( i=j; i<*n; i++, k++ )
385 A[ LDA*i+j ]=ap[ k ];
386 for( i=0, k=0; i<*n; i++ )
387 for( j=0; j<i+1; j++, k++ )
388 AP[ k ]=A[ LDA*i+j ];
void cblas_sspmv(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const int N, const float alpha, const float *Ap, const float *X, const int incX, const float beta, float *Y, const int incY)
void get_uplo_type(char *type, CBLAS_UPLO *uplo)