libflame  revision_anchor
Functions
FLA_Her2c_external.c File Reference

(r)

Functions

FLA_Error FLA_Her2c_external (FLA_Uplo uplo, FLA_Conj conj, FLA_Obj alpha, FLA_Obj x, FLA_Obj y, FLA_Obj A)
 

Function Documentation

◆ FLA_Her2c_external()

FLA_Error FLA_Her2c_external ( FLA_Uplo  uplo,
FLA_Conj  conj,
FLA_Obj  alpha,
FLA_Obj  x,
FLA_Obj  y,
FLA_Obj  A 
)
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_Her2c_check( uplo, conj, alpha, x, y, A );
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( conj, &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 
50  bl1_ssyr2( blis_uplo,
51  m_A,
52  buff_alpha,
53  buff_x, inc_x,
54  buff_y, inc_y,
55  buff_A, rs_A, cs_A );
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 
67  bl1_dsyr2( blis_uplo,
68  m_A,
69  buff_alpha,
70  buff_x, inc_x,
71  buff_y, inc_y,
72  buff_A, rs_A, cs_A );
73  break;
74  }
75 
76  case FLA_COMPLEX:
77  {
78  scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
79  scomplex *buff_x = ( scomplex * ) FLA_COMPLEX_PTR( x );
80  scomplex *buff_y = ( scomplex * ) FLA_COMPLEX_PTR( y );
81  scomplex *buff_alpha = ( scomplex * ) FLA_COMPLEX_PTR( alpha );
82 
83  bl1_cher2( blis_uplo,
84  blis_conj,
85  m_A,
86  buff_alpha,
87  buff_x, inc_x,
88  buff_y, inc_y,
89  buff_A, rs_A, cs_A );
90 
91  break;
92  }
93 
94  case FLA_DOUBLE_COMPLEX:
95  {
96  dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
97  dcomplex *buff_x = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( x );
98  dcomplex *buff_y = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( y );
99  dcomplex *buff_alpha = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( alpha );
100 
101  bl1_zher2( blis_uplo,
102  blis_conj,
103  m_A,
104  buff_alpha,
105  buff_x, inc_x,
106  buff_y, inc_y,
107  buff_A, rs_A, cs_A );
108 
109  break;
110  }
111 
112  }
113 
114  return FLA_SUCCESS;
115 }
FLA_Error FLA_Her2c_check(FLA_Uplo uplo, FLA_Conj conj, FLA_Obj alpha, FLA_Obj x, FLA_Obj y, FLA_Obj A)
Definition: FLA_Her2c_check.c:13
FLA_Bool FLA_Obj_has_zero_dim(FLA_Obj A)
Definition: FLA_Query.c:400
dim_t FLA_Obj_row_stride(FLA_Obj obj)
Definition: FLA_Query.c:167
dim_t FLA_Obj_length(FLA_Obj obj)
Definition: FLA_Query.c:116
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition: FLA_Query.c:174
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
void FLA_Param_map_flame_to_blis_uplo(FLA_Uplo uplo, uplo1_t *blis_uplo)
Definition: FLA_Param.c:285
dim_t FLA_Obj_vector_inc(FLA_Obj obj)
Definition: FLA_Query.c:145
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
int FLA_Datatype
Definition: FLA_type_defs.h:49
void bl1_zher2(uplo1_t uplo, conj1_t conj, int m, dcomplex *alpha, dcomplex *x, int incx, dcomplex *y, int incy, dcomplex *a, int a_rs, int a_cs)
Definition: bl1_her2.c:121
void bl1_cher2(uplo1_t uplo, conj1_t conj, int m, scomplex *alpha, scomplex *x, int incx, scomplex *y, int incy, scomplex *a, int a_rs, int a_cs)
Definition: bl1_her2.c:33
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
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
uplo1_t
Definition: blis_type_defs.h:61
conj1_t
Definition: blis_type_defs.h:80
Definition: blis_type_defs.h:138
Definition: blis_type_defs.h:133

References bl1_cher2(), bl1_dsyr2(), bl1_ssyr2(), bl1_zher2(), FLA_Check_error_level(), FLA_Her2c_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_il_unb_var4(), FLA_Eig_gest_il_unb_var5(), FLA_Eig_gest_iu_unb_var4(), FLA_Eig_gest_iu_unb_var5(), FLA_Eig_gest_nl_unb_var4(), FLA_Eig_gest_nl_unb_var5(), FLA_Eig_gest_nu_unb_var4(), FLA_Eig_gest_nu_unb_var5(), FLA_Her2c(), FLA_Her2k_lh_unb_var10(), FLA_Her2k_lh_unb_var9(), FLA_Her2k_uh_unb_var10(), and FLA_Her2k_uh_unb_var9().