libflame  revision_anchor
Functions
FLA_Gemvc_external.c File Reference

(r)

Functions

FLA_Error FLA_Gemvc_external (FLA_Trans transa, FLA_Conj conjx, FLA_Obj alpha, FLA_Obj A, FLA_Obj x, FLA_Obj beta, FLA_Obj y)
 

Function Documentation

◆ FLA_Gemvc_external()

FLA_Error FLA_Gemvc_external ( FLA_Trans  transa,
FLA_Conj  conjx,
FLA_Obj  alpha,
FLA_Obj  A,
FLA_Obj  x,
FLA_Obj  beta,
FLA_Obj  y 
)

References bl1_cgemv(), bl1_dgemv(), bl1_sgemv(), bl1_zgemv(), FLA_Check_error_level(), FLA_Gemvc_check(), FLA_Obj_col_stride(), FLA_Obj_datatype(), FLA_Obj_has_zero_dim(), FLA_Obj_length(), FLA_Obj_row_stride(), FLA_Obj_vector_inc(), FLA_Obj_width(), FLA_Param_map_flame_to_blis_conj(), FLA_Param_map_flame_to_blis_trans(), and FLA_Scal_external().

Referenced by FLA_Accum_T_UT_fr_unb_var1(), FLA_Apply_H2_UT_l_unb_var1(), FLA_Apply_H2_UT_piv_row(), FLA_Apply_HUD_UT_lh_unb_var1(), FLA_Chol_l_unb_var2(), FLA_Chol_u_unb_var2(), FLA_Eig_gest_il_unb_var2(), FLA_Eig_gest_il_unb_var3(), FLA_Eig_gest_iu_unb_var2(), FLA_Eig_gest_iu_unb_var3(), FLA_Eig_gest_nl_unb_var2(), FLA_Eig_gest_nu_unb_var2(), FLA_Gemm_cc_unb_var1(), FLA_Gemm_cc_unb_var2(), FLA_Gemm_cc_unb_var3(), FLA_Gemm_cc_unb_var4(), FLA_Gemm_ch_unb_var1(), FLA_Gemm_ch_unb_var2(), FLA_Gemm_ch_unb_var3(), FLA_Gemm_ch_unb_var4(), FLA_Gemm_cn_unb_var1(), FLA_Gemm_cn_unb_var2(), FLA_Gemm_ct_unb_var1(), FLA_Gemm_ct_unb_var2(), FLA_Gemm_hc_unb_var1(), FLA_Gemm_hc_unb_var2(), FLA_Gemm_hc_unb_var3(), FLA_Gemm_hc_unb_var4(), FLA_Gemm_hh_unb_var1(), FLA_Gemm_hh_unb_var2(), FLA_Gemm_hh_unb_var3(), FLA_Gemm_hh_unb_var4(), FLA_Gemm_hn_unb_var1(), FLA_Gemm_hn_unb_var2(), FLA_Gemm_ht_unb_var1(), FLA_Gemm_ht_unb_var2(), FLA_Gemm_nc_unb_var3(), FLA_Gemm_nc_unb_var4(), FLA_Gemm_nh_unb_var3(), FLA_Gemm_nh_unb_var4(), FLA_Gemm_tc_unb_var3(), FLA_Gemm_tc_unb_var4(), FLA_Gemm_th_unb_var3(), FLA_Gemm_th_unb_var4(), FLA_Gemvc(), FLA_Hemm_ll_unb_var2(), FLA_Hemm_ll_unb_var3(), FLA_Hemm_ll_unb_var5(), FLA_Hemm_ll_unb_var6(), FLA_Hemm_lu_unb_var1(), FLA_Hemm_lu_unb_var2(), FLA_Hemm_lu_unb_var6(), FLA_Hemm_lu_unb_var7(), FLA_Hemm_rl_unb_var1(), FLA_Hemm_rl_unb_var2(), FLA_Hemm_rl_unb_var6(), FLA_Hemm_rl_unb_var7(), FLA_Hemm_ru_unb_var2(), FLA_Hemm_ru_unb_var3(), FLA_Hemm_ru_unb_var5(), FLA_Hemm_ru_unb_var6(), FLA_Her2k_lh_unb_var1(), FLA_Her2k_lh_unb_var2(), FLA_Her2k_lh_unb_var3(), FLA_Her2k_lh_unb_var6(), FLA_Her2k_lh_unb_var7(), FLA_Her2k_lh_unb_var8(), FLA_Her2k_ln_unb_var2(), FLA_Her2k_ln_unb_var3(), FLA_Her2k_ln_unb_var4(), FLA_Her2k_ln_unb_var5(), FLA_Her2k_ln_unb_var6(), FLA_Her2k_ln_unb_var7(), FLA_Her2k_uh_unb_var1(), FLA_Her2k_uh_unb_var2(), FLA_Her2k_uh_unb_var3(), FLA_Her2k_uh_unb_var6(), FLA_Her2k_uh_unb_var7(), FLA_Her2k_uh_unb_var8(), FLA_Her2k_un_unb_var2(), FLA_Her2k_un_unb_var3(), FLA_Her2k_un_unb_var4(), FLA_Her2k_un_unb_var5(), FLA_Her2k_un_unb_var6(), FLA_Her2k_un_unb_var7(), FLA_Herk_lh_unb_var1(), FLA_Herk_lh_unb_var4(), FLA_Herk_ln_unb_var1(), FLA_Herk_ln_unb_var2(), FLA_Herk_ln_unb_var3(), FLA_Herk_uh_unb_var1(), FLA_Herk_uh_unb_var3(), FLA_Herk_un_unb_var2(), FLA_Herk_un_unb_var4(), FLA_LQ_UT_unb_var2(), FLA_Trmm_llc_unb_var1(), FLA_Trmm_llh_unb_var1(), FLA_Trmm_luc_unb_var1(), FLA_Trmm_luh_unb_var1(), FLA_Trmm_rlc_unb_var1(), FLA_Trmm_rlh_unb_var1(), FLA_Trmm_ruc_unb_var1(), FLA_Trmm_ruh_unb_var1(), FLA_Trsm_llc_unb_var1(), FLA_Trsm_llh_unb_var1(), FLA_Trsm_luc_unb_var1(), FLA_Trsm_luh_unb_var1(), FLA_Trsm_rlc_unb_var1(), FLA_Trsm_rlh_unb_var1(), FLA_Trsm_ruc_unb_var1(), FLA_Trsm_ruh_unb_var1(), FLA_Ttmm_l_unb_var2(), and FLA_Ttmm_u_unb_var2().

14 {
15  FLA_Datatype datatype;
16  int m_A, n_A;
17  int rs_A, cs_A;
18  int inc_x;
19  int inc_y;
20  trans1_t blis_transa;
21  conj1_t blis_conjx;
22 
23  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
24  FLA_Gemvc_check( transa, conjx, alpha, A, x, beta, y );
25 
26  if ( FLA_Obj_has_zero_dim( A ) )
27  {
28  FLA_Scal_external( beta, y );
29  return FLA_SUCCESS;
30  }
31 
32  datatype = FLA_Obj_datatype( A );
33 
34  m_A = FLA_Obj_length( A );
35  n_A = FLA_Obj_width( A );
36  rs_A = FLA_Obj_row_stride( A );
37  cs_A = FLA_Obj_col_stride( A );
38 
39  inc_x = FLA_Obj_vector_inc( x );
40  inc_y = FLA_Obj_vector_inc( y );
41 
42  FLA_Param_map_flame_to_blis_trans( transa, &blis_transa );
43  FLA_Param_map_flame_to_blis_conj( conjx, &blis_conjx );
44 
45 
46  switch( datatype ){
47 
48  case FLA_FLOAT:
49  {
50  float *buff_A = ( float * ) FLA_FLOAT_PTR( A );
51  float *buff_x = ( float * ) FLA_FLOAT_PTR( x );
52  float *buff_y = ( float * ) FLA_FLOAT_PTR( y );
53  float *buff_alpha = ( float * ) FLA_FLOAT_PTR( alpha );
54  float *buff_beta = ( float * ) FLA_FLOAT_PTR( beta );
55 
56  bl1_sgemv( blis_transa,
57  blis_conjx,
58  m_A,
59  n_A,
60  buff_alpha,
61  buff_A, rs_A, cs_A,
62  buff_x, inc_x,
63  buff_beta,
64  buff_y, inc_y );
65 
66  break;
67  }
68 
69  case FLA_DOUBLE:
70  {
71  double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
72  double *buff_x = ( double * ) FLA_DOUBLE_PTR( x );
73  double *buff_y = ( double * ) FLA_DOUBLE_PTR( y );
74  double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
75  double *buff_beta = ( double * ) FLA_DOUBLE_PTR( beta );
76 
77  bl1_dgemv( blis_transa,
78  blis_conjx,
79  m_A,
80  n_A,
81  buff_alpha,
82  buff_A, rs_A, cs_A,
83  buff_x, inc_x,
84  buff_beta,
85  buff_y, inc_y );
86 
87  break;
88  }
89 
90  case FLA_COMPLEX:
91  {
92  scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
93  scomplex *buff_x = ( scomplex * ) FLA_COMPLEX_PTR( x );
94  scomplex *buff_y = ( scomplex * ) FLA_COMPLEX_PTR( y );
95  scomplex *buff_alpha = ( scomplex * ) FLA_COMPLEX_PTR( alpha );
96  scomplex *buff_beta = ( scomplex * ) FLA_COMPLEX_PTR( beta );
97 
98  bl1_cgemv( blis_transa,
99  blis_conjx,
100  m_A,
101  n_A,
102  buff_alpha,
103  buff_A, rs_A, cs_A,
104  buff_x, inc_x,
105  buff_beta,
106  buff_y, inc_y );
107 
108  break;
109  }
110 
111  case FLA_DOUBLE_COMPLEX:
112  {
113  dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
114  dcomplex *buff_x = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( x );
115  dcomplex *buff_y = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( y );
116  dcomplex *buff_alpha = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( alpha );
117  dcomplex *buff_beta = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( beta );
118 
119  bl1_zgemv( blis_transa,
120  blis_conjx,
121  m_A,
122  n_A,
123  buff_alpha,
124  buff_A, rs_A, cs_A,
125  buff_x, inc_x,
126  buff_beta,
127  buff_y, inc_y );
128 
129  break;
130  }
131 
132  }
133 
134  return FLA_SUCCESS;
135 }
void bl1_dgemv(trans1_t transa, conj1_t conjx, int m, int n, double *alpha, double *a, int a_rs, int a_cs, double *x, int incx, double *beta, double *y, int incy)
Definition: bl1_gemv.c:69
void FLA_Param_map_flame_to_blis_trans(FLA_Trans trans, trans1_t *blis_trans)
Definition: FLA_Param.c:245
FLA_Error FLA_Scal_external(FLA_Obj alpha, FLA_Obj A)
Definition: FLA_Scal_external.c:13
conj1_t
Definition: blis_type_defs.h:79
dim_t FLA_Obj_row_stride(FLA_Obj obj)
Definition: FLA_Query.c:167
trans1_t
Definition: blis_type_defs.h:52
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
FLA_Bool FLA_Obj_has_zero_dim(FLA_Obj A)
Definition: FLA_Query.c:400
dim_t FLA_Obj_width(FLA_Obj obj)
Definition: FLA_Query.c:123
void bl1_sgemv(trans1_t transa, conj1_t conjx, int m, int n, float *alpha, float *a, int a_rs, int a_cs, float *x, int incx, float *beta, float *y, int incy)
Definition: bl1_gemv.c:13
Definition: blis_type_defs.h:132
void bl1_zgemv(trans1_t transa, conj1_t conjx, int m, int n, dcomplex *alpha, dcomplex *a, int a_rs, int a_cs, dcomplex *x, int incx, dcomplex *beta, dcomplex *y, int incy)
Definition: bl1_gemv.c:255
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
void FLA_Param_map_flame_to_blis_conj(FLA_Conj conj, conj1_t *blis_conj)
Definition: FLA_Param.c:269
int FLA_Datatype
Definition: FLA_type_defs.h:49
void bl1_cgemv(trans1_t transa, conj1_t conjx, int m, int n, scomplex *alpha, scomplex *a, int a_rs, int a_cs, scomplex *x, int incx, scomplex *beta, scomplex *y, int incy)
Definition: bl1_gemv.c:125
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition: FLA_Query.c:174
FLA_Error FLA_Gemvc_check(FLA_Trans transa, FLA_Conj conjx, FLA_Obj alpha, FLA_Obj A, FLA_Obj x, FLA_Obj beta, FLA_Obj y)
Definition: FLA_Gemvc_check.c:13
dim_t FLA_Obj_vector_inc(FLA_Obj obj)
Definition: FLA_Query.c:145
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116
Definition: blis_type_defs.h:137