libflame  revision_anchor
Functions
FLA_Syr2_external.c File Reference

(r)

Functions

FLA_Error FLA_Syr2_external (FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj x, FLA_Obj y, FLA_Obj A)
 

Function Documentation

◆ FLA_Syr2_external()

FLA_Error FLA_Syr2_external ( FLA_Uplo  uplo,
FLA_Obj  alpha,
FLA_Obj  x,
FLA_Obj  y,
FLA_Obj  A 
)

References bl1_csyr2(), bl1_dsyr2(), bl1_ssyr2(), bl1_zsyr2(), 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_Syr2_check().

Referenced by FLA_Syr2(), FLA_Syr2k_ln_unb_var10(), FLA_Syr2k_ln_unb_var9(), FLA_Syr2k_lt_unb_var10(), FLA_Syr2k_lt_unb_var9(), FLA_Syr2k_un_unb_var10(), FLA_Syr2k_un_unb_var9(), FLA_Syr2k_ut_unb_var10(), and FLA_Syr2k_ut_unb_var9().

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_Syr2_check( uplo, alpha, x, y, A );
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 
48  bl1_ssyr2( blis_uplo,
49  m_A,
50  buff_alpha,
51  buff_x, inc_x,
52  buff_y, inc_y,
53  buff_A, rs_A, cs_A );
54 
55  break;
56  }
57 
58  case FLA_DOUBLE:
59  {
60  double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
61  double *buff_x = ( double * ) FLA_DOUBLE_PTR( x );
62  double *buff_y = ( double * ) FLA_DOUBLE_PTR( y );
63  double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
64 
65  bl1_dsyr2( blis_uplo,
66  m_A,
67  buff_alpha,
68  buff_x, inc_x,
69  buff_y, inc_y,
70  buff_A, rs_A, cs_A );
71 
72  break;
73  }
74 
75  case FLA_COMPLEX:
76  {
77  scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
78  scomplex *buff_x = ( scomplex * ) FLA_COMPLEX_PTR( x );
79  scomplex *buff_y = ( scomplex * ) FLA_COMPLEX_PTR( y );
80  scomplex *buff_alpha = ( scomplex * ) FLA_COMPLEX_PTR( alpha );
81 
82  bl1_csyr2( blis_uplo,
83  m_A,
84  buff_alpha,
85  buff_x, inc_x,
86  buff_y, inc_y,
87  buff_A, rs_A, cs_A );
88 
89  break;
90  }
91 
92  case FLA_DOUBLE_COMPLEX:
93  {
94  dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
95  dcomplex *buff_x = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( x );
96  dcomplex *buff_y = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( y );
97  dcomplex *buff_alpha = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( alpha );
98 
99  bl1_zsyr2( blis_uplo,
100  m_A,
101  buff_alpha,
102  buff_x, inc_x,
103  buff_y, inc_y,
104  buff_A, rs_A, cs_A );
105 
106  break;
107  }
108 
109  }
110 
111  return FLA_SUCCESS;
112 }
uplo1_t
Definition: blis_type_defs.h:60
dim_t FLA_Obj_row_stride(FLA_Obj obj)
Definition: FLA_Query.c:167
void bl1_zsyr2(uplo1_t uplo, int m, dcomplex *alpha, dcomplex *x, int incx, dcomplex *y, int incy, dcomplex *a, int a_rs, int a_cs)
Definition: bl1_syr2.c:148
void bl1_ssyr2(uplo1_t uplo, int m, float *alpha, float *x, int incx, float *y, int incy, float *a, int a_rs, int a_cs)
Definition: bl1_syr2.c:13
void bl1_dsyr2(uplo1_t uplo, int m, double *alpha, double *x, int incx, double *y, int incy, double *a, int a_rs, int a_cs)
Definition: bl1_syr2.c:58
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_csyr2(uplo1_t uplo, int m, scomplex *alpha, scomplex *x, int incx, scomplex *y, int incy, scomplex *a, int a_rs, int a_cs)
Definition: bl1_syr2.c:103
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
FLA_Error FLA_Syr2_check(FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj x, FLA_Obj y, FLA_Obj A)
Definition: FLA_Syr2_check.c:13
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
Definition: blis_type_defs.h:137