libflame
revision_anchor
|
Functions | |
void | bli_ssymv (uplo_t uplo, int m, float *alpha, float *a, int a_rs, int a_cs, float *x, int incx, float *beta, float *y, int incy) |
void | bli_dsymv (uplo_t uplo, int m, double *alpha, double *a, int a_rs, int a_cs, double *x, int incx, double *beta, double *y, int incy) |
void | bli_csymv (uplo_t uplo, int m, scomplex *alpha, scomplex *a, int a_rs, int a_cs, scomplex *x, int incx, scomplex *beta, scomplex *y, int incy) |
void | bli_zsymv (uplo_t uplo, int m, dcomplex *alpha, dcomplex *a, int a_rs, int a_cs, dcomplex *x, int incx, dcomplex *beta, dcomplex *y, int incy) |
void | bli_ssymv_blas (uplo_t uplo, int m, float *alpha, float *a, int lda, float *x, int incx, float *beta, float *y, int incy) |
void | bli_dsymv_blas (uplo_t uplo, int m, double *alpha, double *a, int lda, double *x, int incx, double *beta, double *y, int incy) |
void | bli_csymv_blas (uplo_t uplo, int m, scomplex *alpha, scomplex *a, int lda, scomplex *x, int incx, scomplex *beta, scomplex *y, int incy) |
void | bli_zsymv_blas (uplo_t uplo, int m, dcomplex *alpha, dcomplex *a, int lda, dcomplex *x, int incx, dcomplex *beta, dcomplex *y, int incy) |
void bli_csymv | ( | uplo_t | uplo, |
int | m, | ||
scomplex * | alpha, | ||
scomplex * | a, | ||
int | a_rs, | ||
int | a_cs, | ||
scomplex * | x, | ||
int | incx, | ||
scomplex * | beta, | ||
scomplex * | y, | ||
int | incy | ||
) |
References bli_ccreate_contigmr(), bli_cfree_contigm(), bli_csymv_blas(), bli_is_row_storage(), and bli_zero_dim1().
Referenced by FLA_Symv_external().
{ scomplex* a_save = a; int a_rs_save = a_rs; int a_cs_save = a_cs; int lda, inca; // Return early if possible. if ( bli_zero_dim1( m ) ) return; // If necessary, allocate, initialize, and use a temporary contiguous // copy of the matrix rather than the original matrix. bli_ccreate_contigmr( uplo, m, m, a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); // Initialize with values assuming column-major storage. lda = a_cs; inca = a_rs; // If A is a row-major matrix, then we can use the underlying column-major // BLAS implementation by fiddling with the parameters. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( lda, inca ); bli_toggle_uplo( uplo ); } bli_csymv_blas( uplo, m, alpha, a, lda, x, incx, beta, y, incy ); // Free the temporary contiguous matrix. bli_cfree_contigm( a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); }
void bli_csymv_blas | ( | uplo_t | uplo, |
int | m, | ||
scomplex * | alpha, | ||
scomplex * | a, | ||
int | lda, | ||
scomplex * | x, | ||
int | incx, | ||
scomplex * | beta, | ||
scomplex * | y, | ||
int | incy | ||
) |
References bli_callocv(), bli_ccopyv(), bli_cfree(), bli_param_map_to_netlib_side(), bli_param_map_to_netlib_uplo(), BLIS_LEFT, BLIS_NO_CONJUGATE, cblas_csymm(), CblasColMajor, and F77_csymm().
Referenced by bli_csymv().
{ scomplex* x_copy; scomplex* y_copy; int n = 1; int ldx = m; int ldy = m; #ifdef BLIS_ENABLE_CBLAS_INTERFACES enum CBLAS_ORDER cblas_order = CblasColMajor; enum CBLAS_SIDE cblas_side; enum CBLAS_UPLO cblas_uplo; bli_param_map_to_netlib_side( BLIS_LEFT, &cblas_side ); bli_param_map_to_netlib_uplo( uplo, &cblas_uplo ); x_copy = bli_callocv( m ); y_copy = bli_callocv( m ); bli_ccopyv( BLIS_NO_CONJUGATE, m, x, incx, x_copy, 1 ); bli_ccopyv( BLIS_NO_CONJUGATE, m, y, incy, y_copy, 1 ); cblas_csymm( cblas_order, cblas_side, cblas_uplo, m, n, alpha, a, lda, x_copy, ldx, beta, y_copy, ldy ); bli_ccopyv( BLIS_NO_CONJUGATE, m, y_copy, 1, y, incy ); bli_cfree( x_copy ); bli_cfree( y_copy ); #else char blas_side; char blas_uplo; bli_param_map_to_netlib_side( BLIS_LEFT, &blas_side ); bli_param_map_to_netlib_uplo( uplo, &blas_uplo ); x_copy = bli_callocv( m ); y_copy = bli_callocv( m ); bli_ccopyv( BLIS_NO_CONJUGATE, m, x, incx, x_copy, 1 ); bli_ccopyv( BLIS_NO_CONJUGATE, m, y, incy, y_copy, 1 ); F77_csymm ( &blas_side, &blas_uplo, &m, &n, alpha, a, &lda, x_copy, &ldx, beta, y_copy, &ldy ); bli_ccopyv( BLIS_NO_CONJUGATE, m, y_copy, 1, y, incy ); bli_cfree( x_copy ); bli_cfree( y_copy ); #endif }
void bli_dsymv | ( | uplo_t | uplo, |
int | m, | ||
double * | alpha, | ||
double * | a, | ||
int | a_rs, | ||
int | a_cs, | ||
double * | x, | ||
int | incx, | ||
double * | beta, | ||
double * | y, | ||
int | incy | ||
) |
References bli_dcreate_contigmr(), bli_dfree_contigm(), bli_dsymv_blas(), bli_is_row_storage(), and bli_zero_dim1().
Referenced by bli_dhemv(), FLA_Hemv_external(), FLA_Hemvc_external(), FLA_Symv_external(), FLA_Tridiag_UT_l_step_ofd_var2(), FLA_Tridiag_UT_l_step_ofd_var3(), FLA_Tridiag_UT_l_step_opd_var1(), FLA_Tridiag_UT_l_step_opd_var2(), and FLA_Tridiag_UT_l_step_opd_var3().
{ double* a_save = a; int a_rs_save = a_rs; int a_cs_save = a_cs; int lda, inca; // Return early if possible. if ( bli_zero_dim1( m ) ) return; // If necessary, allocate, initialize, and use a temporary contiguous // copy of the matrix rather than the original matrix. bli_dcreate_contigmr( uplo, m, m, a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); // Initialize with values assuming column-major storage. lda = a_cs; inca = a_rs; // If A is a row-major matrix, then we can use the underlying column-major // BLAS implementation by fiddling with the parameters. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( lda, inca ); bli_toggle_uplo( uplo ); } bli_dsymv_blas( uplo, m, alpha, a, lda, x, incx, beta, y, incy ); // Free the temporary contiguous matrix. bli_dfree_contigm( a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); }
void bli_dsymv_blas | ( | uplo_t | uplo, |
int | m, | ||
double * | alpha, | ||
double * | a, | ||
int | lda, | ||
double * | x, | ||
int | incx, | ||
double * | beta, | ||
double * | y, | ||
int | incy | ||
) |
References bli_param_map_to_netlib_uplo(), cblas_dsymv(), CblasColMajor, and F77_dsymv().
Referenced by bli_dsymv().
{ #ifdef BLIS_ENABLE_CBLAS_INTERFACES enum CBLAS_ORDER cblas_order = CblasColMajor; enum CBLAS_UPLO cblas_uplo; bli_param_map_to_netlib_uplo( uplo, &cblas_uplo ); cblas_dsymv( cblas_order, cblas_uplo, m, *alpha, a, lda, x, incx, *beta, y, incy ); #else char blas_uplo; bli_param_map_to_netlib_uplo( uplo, &blas_uplo ); F77_dsymv( &blas_uplo, &m, alpha, a, &lda, x, &incx, beta, y, &incy ); #endif }
void bli_ssymv | ( | uplo_t | uplo, |
int | m, | ||
float * | alpha, | ||
float * | a, | ||
int | a_rs, | ||
int | a_cs, | ||
float * | x, | ||
int | incx, | ||
float * | beta, | ||
float * | y, | ||
int | incy | ||
) |
References bli_is_row_storage(), bli_screate_contigmr(), bli_sfree_contigm(), bli_ssymv_blas(), and bli_zero_dim1().
Referenced by bli_shemv(), FLA_Hemv_external(), FLA_Hemvc_external(), FLA_Symv_external(), FLA_Tridiag_UT_l_step_ofs_var2(), FLA_Tridiag_UT_l_step_ofs_var3(), FLA_Tridiag_UT_l_step_ops_var1(), FLA_Tridiag_UT_l_step_ops_var2(), and FLA_Tridiag_UT_l_step_ops_var3().
{ float* a_save = a; int a_rs_save = a_rs; int a_cs_save = a_cs; int lda, inca; // Return early if possible. if ( bli_zero_dim1( m ) ) return; // If necessary, allocate, initialize, and use a temporary contiguous // copy of the matrix rather than the original matrix. bli_screate_contigmr( uplo, m, m, a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); // Initialize with values assuming column-major storage. lda = a_cs; inca = a_rs; // If A is a row-major matrix, then we can use the underlying column-major // BLAS implementation by fiddling with the parameters. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( lda, inca ); bli_toggle_uplo( uplo ); } bli_ssymv_blas( uplo, m, alpha, a, lda, x, incx, beta, y, incy ); // Free the temporary contiguous matrix. bli_sfree_contigm( a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); }
void bli_ssymv_blas | ( | uplo_t | uplo, |
int | m, | ||
float * | alpha, | ||
float * | a, | ||
int | lda, | ||
float * | x, | ||
int | incx, | ||
float * | beta, | ||
float * | y, | ||
int | incy | ||
) |
References bli_param_map_to_netlib_uplo(), cblas_ssymv(), CblasColMajor, and F77_ssymv().
Referenced by bli_ssymv().
{ #ifdef BLIS_ENABLE_CBLAS_INTERFACES enum CBLAS_ORDER cblas_order = CblasColMajor; enum CBLAS_UPLO cblas_uplo; bli_param_map_to_netlib_uplo( uplo, &cblas_uplo ); cblas_ssymv( cblas_order, cblas_uplo, m, *alpha, a, lda, x, incx, *beta, y, incy ); #else char blas_uplo; bli_param_map_to_netlib_uplo( uplo, &blas_uplo ); F77_ssymv( &blas_uplo, &m, alpha, a, &lda, x, &incx, beta, y, &incy ); #endif }
void bli_zsymv | ( | uplo_t | uplo, |
int | m, | ||
dcomplex * | alpha, | ||
dcomplex * | a, | ||
int | a_rs, | ||
int | a_cs, | ||
dcomplex * | x, | ||
int | incx, | ||
dcomplex * | beta, | ||
dcomplex * | y, | ||
int | incy | ||
) |
References bli_is_row_storage(), bli_zcreate_contigmr(), bli_zero_dim1(), bli_zfree_contigm(), and bli_zsymv_blas().
Referenced by FLA_Symv_external().
{ dcomplex* a_save = a; int a_rs_save = a_rs; int a_cs_save = a_cs; int lda, inca; // Return early if possible. if ( bli_zero_dim1( m ) ) return; // If necessary, allocate, initialize, and use a temporary contiguous // copy of the matrix rather than the original matrix. bli_zcreate_contigmr( uplo, m, m, a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); // Initialize with values assuming column-major storage. lda = a_cs; inca = a_rs; // If A is a row-major matrix, then we can use the underlying column-major // BLAS implementation by fiddling with the parameters. if ( bli_is_row_storage( a_rs, a_cs ) ) { bli_swap_ints( lda, inca ); bli_toggle_uplo( uplo ); } bli_zsymv_blas( uplo, m, alpha, a, lda, x, incx, beta, y, incy ); // Free the temporary contiguous matrix. bli_zfree_contigm( a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); }
void bli_zsymv_blas | ( | uplo_t | uplo, |
int | m, | ||
dcomplex * | alpha, | ||
dcomplex * | a, | ||
int | lda, | ||
dcomplex * | x, | ||
int | incx, | ||
dcomplex * | beta, | ||
dcomplex * | y, | ||
int | incy | ||
) |
References bli_param_map_to_netlib_side(), bli_param_map_to_netlib_uplo(), bli_zallocv(), bli_zcopyv(), bli_zfree(), BLIS_LEFT, BLIS_NO_CONJUGATE, cblas_zsymm(), CblasColMajor, and F77_zsymm().
Referenced by bli_zsymv().
{ dcomplex* x_copy; dcomplex* y_copy; int n = 1; int ldx = m; int ldy = m; #ifdef BLIS_ENABLE_CBLAS_INTERFACES enum CBLAS_ORDER cblas_order = CblasColMajor; enum CBLAS_SIDE cblas_side; enum CBLAS_UPLO cblas_uplo; bli_param_map_to_netlib_side( BLIS_LEFT, &cblas_side ); bli_param_map_to_netlib_uplo( uplo, &cblas_uplo ); x_copy = bli_zallocv( m ); y_copy = bli_zallocv( m ); bli_zcopyv( BLIS_NO_CONJUGATE, m, x, incx, x_copy, 1 ); bli_zcopyv( BLIS_NO_CONJUGATE, m, y, incy, y_copy, 1 ); cblas_zsymm( cblas_order, cblas_side, cblas_uplo, m, n, alpha, a, lda, x_copy, ldx, beta, y_copy, ldy ); bli_zcopyv( BLIS_NO_CONJUGATE, m, y_copy, 1, y, incy ); bli_zfree( x_copy ); bli_zfree( y_copy ); #else char blas_side; char blas_uplo; bli_param_map_to_netlib_side( BLIS_LEFT, &blas_side ); bli_param_map_to_netlib_uplo( uplo, &blas_uplo ); x_copy = bli_zallocv( m ); y_copy = bli_zallocv( m ); bli_zcopyv( BLIS_NO_CONJUGATE, m, x, incx, x_copy, 1 ); bli_zcopyv( BLIS_NO_CONJUGATE, m, y, incy, y_copy, 1 ); F77_zsymm ( &blas_side, &blas_uplo, &m, &n, alpha, a, &lda, x_copy, &ldx, beta, y_copy, &ldy ); bli_zcopyv( BLIS_NO_CONJUGATE, m, y_copy, 1, y, incy ); bli_zfree( x_copy ); bli_zfree( y_copy ); #endif }