libflame revision_anchor
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().