libflame
revision_anchor
|
Functions | |
void | bli_ssyr (uplo_t uplo, int m, float *alpha, float *x, int incx, float *a, int a_rs, int a_cs) |
void | bli_dsyr (uplo_t uplo, int m, double *alpha, double *x, int incx, double *a, int a_rs, int a_cs) |
void | bli_csyr (uplo_t uplo, int m, scomplex *alpha, scomplex *x, int incx, scomplex *a, int a_rs, int a_cs) |
void | bli_zsyr (uplo_t uplo, int m, dcomplex *alpha, dcomplex *x, int incx, dcomplex *a, int a_rs, int a_cs) |
void | bli_ssyr_blas (uplo_t uplo, int m, float *alpha, float *x, int incx, float *a, int lda) |
void | bli_dsyr_blas (uplo_t uplo, int m, double *alpha, double *x, int incx, double *a, int lda) |
void | bli_csyr_blas (uplo_t uplo, int m, scomplex *alpha, scomplex *x, int incx, scomplex *a, int lda) |
void | bli_zsyr_blas (uplo_t uplo, int m, dcomplex *alpha, dcomplex *x, int incx, dcomplex *a, int lda) |
void bli_csyr | ( | uplo_t | uplo, |
int | m, | ||
scomplex * | alpha, | ||
scomplex * | x, | ||
int | incx, | ||
scomplex * | a, | ||
int | a_rs, | ||
int | a_cs | ||
) |
References bli_ccreate_contigmr(), bli_cfree_saved_contigm(), bli_csyr_blas(), bli_is_row_storage(), and bli_zero_dim1().
Referenced by FLA_Syr_external().
{ int m_save = m; 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_csyr_blas( uplo, m, alpha, x, incx, a, lda ); // Free the temporary contiguous matrix. bli_cfree_saved_contigm( m_save, m_save, a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); }
void bli_csyr_blas | ( | uplo_t | uplo, |
int | m, | ||
scomplex * | alpha, | ||
scomplex * | x, | ||
int | incx, | ||
scomplex * | a, | ||
int | lda | ||
) |
References bli_callocv(), bli_ccopyv(), bli_cfree(), bli_param_map_to_netlib_trans(), bli_param_map_to_netlib_uplo(), BLIS_NO_CONJUGATE, BLIS_NO_TRANSPOSE, cblas_csyrk(), CblasColMajor, F77_csyrk(), scomplex::imag, and scomplex::real.
Referenced by bli_csyr().
{ scomplex* x_copy; scomplex beta; int k = 1; int ldx = m; #ifdef BLIS_ENABLE_CBLAS_INTERFACES enum CBLAS_ORDER cblas_order = CblasColMajor; enum CBLAS_UPLO cblas_uplo; enum CBLAS_TRANSPOSE cblas_trans; bli_param_map_to_netlib_uplo( uplo, &cblas_uplo ); bli_param_map_to_netlib_trans( BLIS_NO_TRANSPOSE, &cblas_trans ); x_copy = bli_callocv( m ); bli_ccopyv( BLIS_NO_CONJUGATE, m, x, incx, x_copy, 1 ); beta.real = 1.0; beta.imag = 0.0; cblas_csyrk( cblas_order, cblas_uplo, cblas_trans, m, k, alpha, x_copy, ldx, &beta, a, lda ); bli_cfree( x_copy ); #else char blas_uplo; char blas_trans; bli_param_map_to_netlib_uplo( uplo, &blas_uplo ); bli_param_map_to_netlib_trans( BLIS_NO_TRANSPOSE, &blas_trans ); x_copy = bli_callocv( m ); bli_ccopyv( BLIS_NO_CONJUGATE, m, x, incx, x_copy, 1 ); beta.real = 1.0; beta.imag = 0.0; F77_csyrk ( &blas_uplo, &blas_trans, &m, &k, alpha, x_copy, &ldx, &beta, a, &lda ); bli_cfree( x_copy ); #endif }
void bli_dsyr | ( | uplo_t | uplo, |
int | m, | ||
double * | alpha, | ||
double * | x, | ||
int | incx, | ||
double * | a, | ||
int | a_rs, | ||
int | a_cs | ||
) |
References bli_dcreate_contigmr(), bli_dfree_saved_contigm(), bli_dsyr_blas(), bli_is_row_storage(), and bli_zero_dim1().
Referenced by bli_dher(), FLA_Chol_l_opd_var3(), FLA_Chol_u_opd_var3(), FLA_Her_external(), FLA_Herc_external(), FLA_Syr_external(), FLA_Ttmm_l_opd_var1(), and FLA_Ttmm_u_opd_var1().
{ int m_save = m; 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_dsyr_blas( uplo, m, alpha, x, incx, a, lda ); // Free the temporary contiguous matrix. bli_dfree_saved_contigm( m_save, m_save, a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); }
void bli_dsyr_blas | ( | uplo_t | uplo, |
int | m, | ||
double * | alpha, | ||
double * | x, | ||
int | incx, | ||
double * | a, | ||
int | lda | ||
) |
References bli_param_map_to_netlib_uplo(), cblas_dsyr(), CblasColMajor, and F77_dsyr().
Referenced by bli_dsyr().
{ #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_dsyr( cblas_order, cblas_uplo, m, *alpha, x, incx, a, lda ); #else char blas_uplo; bli_param_map_to_netlib_uplo( uplo, &blas_uplo ); F77_dsyr( &blas_uplo, &m, alpha, x, &incx, a, &lda ); #endif }
void bli_ssyr | ( | uplo_t | uplo, |
int | m, | ||
float * | alpha, | ||
float * | x, | ||
int | incx, | ||
float * | a, | ||
int | a_rs, | ||
int | a_cs | ||
) |
References bli_is_row_storage(), bli_screate_contigmr(), bli_sfree_saved_contigm(), bli_ssyr_blas(), and bli_zero_dim1().
Referenced by bli_sher(), FLA_Chol_l_ops_var3(), FLA_Chol_u_ops_var3(), FLA_Her_external(), FLA_Herc_external(), FLA_Syr_external(), FLA_Ttmm_l_ops_var1(), and FLA_Ttmm_u_ops_var1().
{ int m_save = m; 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_ssyr_blas( uplo, m, alpha, x, incx, a, lda ); // Free the temporary contiguous matrix. bli_sfree_saved_contigm( m_save, m_save, a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); }
void bli_ssyr_blas | ( | uplo_t | uplo, |
int | m, | ||
float * | alpha, | ||
float * | x, | ||
int | incx, | ||
float * | a, | ||
int | lda | ||
) |
References bli_param_map_to_netlib_uplo(), cblas_ssyr(), CblasColMajor, and F77_ssyr().
Referenced by bli_ssyr().
{ #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_ssyr( cblas_order, cblas_uplo, m, *alpha, x, incx, a, lda ); #else char blas_uplo; bli_param_map_to_netlib_uplo( uplo, &blas_uplo ); F77_ssyr( &blas_uplo, &m, alpha, x, &incx, a, &lda ); #endif }
void bli_zsyr | ( | uplo_t | uplo, |
int | m, | ||
dcomplex * | alpha, | ||
dcomplex * | x, | ||
int | incx, | ||
dcomplex * | a, | ||
int | a_rs, | ||
int | a_cs | ||
) |
References bli_is_row_storage(), bli_zcreate_contigmr(), bli_zero_dim1(), bli_zfree_saved_contigm(), and bli_zsyr_blas().
Referenced by FLA_Syr_external().
{ int m_save = m; 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_zsyr_blas( uplo, m, alpha, x, incx, a, lda ); // Free the temporary contiguous matrix. bli_zfree_saved_contigm( m_save, m_save, a_save, a_rs_save, a_cs_save, &a, &a_rs, &a_cs ); }
void bli_zsyr_blas | ( | uplo_t | uplo, |
int | m, | ||
dcomplex * | alpha, | ||
dcomplex * | x, | ||
int | incx, | ||
dcomplex * | a, | ||
int | lda | ||
) |
References bli_param_map_to_netlib_trans(), bli_param_map_to_netlib_uplo(), bli_zallocv(), bli_zcopyv(), bli_zfree(), BLIS_NO_CONJUGATE, BLIS_NO_TRANSPOSE, cblas_zsyrk(), CblasColMajor, F77_zsyrk(), dcomplex::imag, and dcomplex::real.
Referenced by bli_zsyr().
{ dcomplex* x_copy; dcomplex beta; int k = 1; int ldx = m; #ifdef BLIS_ENABLE_CBLAS_INTERFACES enum CBLAS_ORDER cblas_order = CblasColMajor; enum CBLAS_UPLO cblas_uplo; enum CBLAS_TRANSPOSE cblas_trans; bli_param_map_to_netlib_uplo( uplo, &cblas_uplo ); bli_param_map_to_netlib_trans( BLIS_NO_TRANSPOSE, &cblas_trans ); x_copy = bli_zallocv( m ); bli_zcopyv( BLIS_NO_CONJUGATE, m, x, incx, x_copy, 1 ); beta.real = 1.0; beta.imag = 0.0; cblas_zsyrk( cblas_order, cblas_uplo, cblas_trans, m, k, alpha, x_copy, ldx, &beta, a, lda ); bli_zfree( x_copy ); #else char blas_uplo; char blas_trans; bli_param_map_to_netlib_uplo( uplo, &blas_uplo ); bli_param_map_to_netlib_trans( BLIS_NO_TRANSPOSE, &blas_trans ); x_copy = bli_zallocv( m ); bli_zcopyv( BLIS_NO_CONJUGATE, m, x, incx, x_copy, 1 ); beta.real = 1.0; beta.imag = 0.0; F77_zsyrk ( &blas_uplo, &blas_trans, &m, &k, alpha, x_copy, &ldx, &beta, a, &lda ); bli_zfree( x_copy ); #endif }