libflame  revision_anchor
Functions
FLA_Hemv_external.c File Reference

(r)

Functions

FLA_Error FLA_Hemv_external (FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj A, FLA_Obj x, FLA_Obj beta, FLA_Obj y)
 

Function Documentation

◆ FLA_Hemv_external()

FLA_Error FLA_Hemv_external ( FLA_Uplo  uplo,
FLA_Obj  alpha,
FLA_Obj  A,
FLA_Obj  x,
FLA_Obj  beta,
FLA_Obj  y 
)

References bl1_chemv(), bl1_dsymv(), bl1_ssymv(), bl1_zhemv(), FLA_Check_error_level(), FLA_Hemv_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_Param_map_flame_to_blis_conj(), and FLA_Param_map_flame_to_blis_uplo().

Referenced by FLA_Eig_gest_nl_unb_var1(), FLA_Eig_gest_nl_unb_var2(), FLA_Hemm_ll_unb_var10(), FLA_Hemm_ll_unb_var9(), FLA_Hemm_lu_unb_var10(), FLA_Hemm_lu_unb_var9(), and FLA_Hemv().

14 {
15  FLA_Datatype datatype;
16  int m_A;
17  int rs_A, cs_A;
18  int inc_x;
19  int inc_y;
20  uplo1_t blis_uplo;
21  conj1_t blis_conj;
22 
23  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
24  FLA_Hemv_check( uplo, alpha, A, x, beta, y );
25 
26  if ( FLA_Obj_has_zero_dim( A ) ) return FLA_SUCCESS;
27 
28  datatype = FLA_Obj_datatype( A );
29 
30  m_A = FLA_Obj_length( A );
31  rs_A = FLA_Obj_row_stride( A );
32  cs_A = FLA_Obj_col_stride( A );
33 
34  inc_x = FLA_Obj_vector_inc( x );
35  inc_y = FLA_Obj_vector_inc( y );
36 
37  FLA_Param_map_flame_to_blis_uplo( uplo, &blis_uplo );
38  FLA_Param_map_flame_to_blis_conj( FLA_NO_CONJUGATE, &blis_conj );
39 
40 
41  switch( datatype ){
42 
43  case FLA_FLOAT:
44  {
45  float *buff_A = ( float * ) FLA_FLOAT_PTR( A );
46  float *buff_x = ( float * ) FLA_FLOAT_PTR( x );
47  float *buff_y = ( float * ) FLA_FLOAT_PTR( y );
48  float *buff_alpha = ( float * ) FLA_FLOAT_PTR( alpha );
49  float *buff_beta = ( float * ) FLA_FLOAT_PTR( beta );
50 
51  bl1_ssymv( blis_uplo,
52  m_A,
53  buff_alpha,
54  buff_A, rs_A, cs_A,
55  buff_x, inc_x,
56  buff_beta,
57  buff_y, inc_y );
58 
59  break;
60  }
61 
62  case FLA_DOUBLE:
63  {
64  double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
65  double *buff_x = ( double * ) FLA_DOUBLE_PTR( x );
66  double *buff_y = ( double * ) FLA_DOUBLE_PTR( y );
67  double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
68  double *buff_beta = ( double * ) FLA_DOUBLE_PTR( beta );
69 
70  bl1_dsymv( blis_uplo,
71  m_A,
72  buff_alpha,
73  buff_A, rs_A, cs_A,
74  buff_x, inc_x,
75  buff_beta,
76  buff_y, inc_y );
77 
78  break;
79  }
80 
81  case FLA_COMPLEX:
82  {
83  scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
84  scomplex *buff_x = ( scomplex * ) FLA_COMPLEX_PTR( x );
85  scomplex *buff_y = ( scomplex * ) FLA_COMPLEX_PTR( y );
86  scomplex *buff_alpha = ( scomplex * ) FLA_COMPLEX_PTR( alpha );
87  scomplex *buff_beta = ( scomplex * ) FLA_COMPLEX_PTR( beta );
88 
89  bl1_chemv( blis_uplo,
90  blis_conj,
91  m_A,
92  buff_alpha,
93  buff_A, rs_A, cs_A,
94  buff_x, inc_x,
95  buff_beta,
96  buff_y, inc_y );
97 
98  break;
99  }
100 
101  case FLA_DOUBLE_COMPLEX:
102  {
103  dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
104  dcomplex *buff_x = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( x );
105  dcomplex *buff_y = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( y );
106  dcomplex *buff_alpha = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( alpha );
107  dcomplex *buff_beta = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( beta );
108 
109  bl1_zhemv( blis_uplo,
110  blis_conj,
111  m_A,
112  buff_alpha,
113  buff_A, rs_A, cs_A,
114  buff_x, inc_x,
115  buff_beta,
116  buff_y, inc_y );
117 
118  break;
119  }
120 
121  }
122 
123  return FLA_SUCCESS;
124 }
uplo1_t
Definition: blis_type_defs.h:60
conj1_t
Definition: blis_type_defs.h:79
dim_t FLA_Obj_row_stride(FLA_Obj obj)
Definition: FLA_Query.c:167
void bl1_ssymv(uplo1_t uplo, int m, float *alpha, float *a, int a_rs, int a_cs, float *x, int incx, float *beta, float *y, int incy)
Definition: bl1_symv.c:13
void bl1_chemv(uplo1_t uplo, conj1_t conj, int m, scomplex *alpha, scomplex *a, int a_rs, int a_cs, scomplex *x, int incx, scomplex *beta, scomplex *y, int incy)
Definition: bl1_hemv.c:35
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
void FLA_Param_map_flame_to_blis_uplo(FLA_Uplo uplo, uplo1_t *blis_uplo)
Definition: FLA_Param.c:285
Definition: blis_type_defs.h:132
void bl1_dsymv(uplo1_t uplo, int m, double *alpha, double *a, int a_rs, int a_cs, double *x, int incx, double *beta, double *y, int incy)
Definition: bl1_symv.c:56
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
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition: FLA_Query.c:174
dim_t FLA_Obj_vector_inc(FLA_Obj obj)
Definition: FLA_Query.c:145
void bl1_zhemv(uplo1_t uplo, conj1_t conj, int m, dcomplex *alpha, dcomplex *a, int a_rs, int a_cs, dcomplex *x, int incx, dcomplex *beta, dcomplex *y, int incy)
Definition: bl1_hemv.c:134
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116
FLA_Error FLA_Hemv_check(FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj A, FLA_Obj x, FLA_Obj beta, FLA_Obj y)
Definition: FLA_Hemv_check.c:13
Definition: blis_type_defs.h:137