libflame  revision_anchor
Functions
FLA_Symv_external.c File Reference

(r)

Functions

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

Function Documentation

◆ FLA_Symv_external()

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

References bl1_csymv(), bl1_dsymv(), bl1_ssymv(), bl1_zsymv(), FLA_Check_error_level(), 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_uplo(), and FLA_Symv_check().

Referenced by FLA_Symm_ll_unb_var10(), FLA_Symm_ll_unb_var9(), FLA_Symm_lu_unb_var10(), FLA_Symm_lu_unb_var9(), FLA_Symm_rl_unb_var10(), FLA_Symm_rl_unb_var9(), FLA_Symm_ru_unb_var10(), FLA_Symm_ru_unb_var9(), and FLA_Symv().

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