libflame revision_anchor
FLA_Her_external.c File Reference

(r)

Functions

FLA_Error FLA_Her_external (FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj x, FLA_Obj A)
 

Function Documentation

◆ FLA_Her_external()

FLA_Error FLA_Her_external ( FLA_Uplo uplo,
FLA_Obj alpha,
FLA_Obj x,
FLA_Obj A )
14{
15 FLA_Datatype datatype;
16 int m_A;
17 int rs_A, cs_A;
18 int inc_x;
19 uplo1_t blis_uplo;
20 conj1_t blis_conj;
21
22 if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
23 FLA_Her_check( uplo, alpha, x, 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
35 FLA_Param_map_flame_to_blis_uplo( uplo, &blis_uplo );
36 FLA_Param_map_flame_to_blis_conj( FLA_NO_CONJUGATE, &blis_conj );
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_alpha = ( float * ) FLA_FLOAT_PTR( alpha );
46
47 bl1_ssyr( blis_uplo,
48 m_A,
49 buff_alpha,
50 buff_x, inc_x,
51 buff_A, rs_A, cs_A );
52
53 break;
54 }
55
56 case FLA_DOUBLE:
57 {
58 double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
59 double *buff_x = ( double * ) FLA_DOUBLE_PTR( x );
60 double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
61
62 bl1_dsyr( blis_uplo,
63 m_A,
64 buff_alpha,
65 buff_x, inc_x,
66 buff_A, rs_A, cs_A );
67
68 break;
69 }
70
71 case FLA_COMPLEX:
72 {
73 scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
74 scomplex *buff_x = ( scomplex * ) FLA_COMPLEX_PTR( x );
75 float *buff_alpha = ( float * ) FLA_FLOAT_PTR( alpha );
76
77 bl1_cher( blis_uplo,
78 blis_conj,
79 m_A,
80 buff_alpha,
81 buff_x, inc_x,
82 buff_A, rs_A, cs_A );
83
84 break;
85 }
86
87 case FLA_DOUBLE_COMPLEX:
88 {
89 dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
90 dcomplex *buff_x = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( x );
91 double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
92
93 bl1_zher( blis_uplo,
94 blis_conj,
95 m_A,
96 buff_alpha,
97 buff_x, inc_x,
98 buff_A, rs_A, cs_A );
99
100 break;
101 }
102
103 }
104
105 return FLA_SUCCESS;
106}
FLA_Error FLA_Her_check(FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj x, FLA_Obj A)
Definition FLA_Her_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_cher(uplo1_t uplo, conj1_t conj, int m, float *alpha, scomplex *x, int incx, scomplex *a, int a_rs, int a_cs)
Definition bl1_her.c:31
void bl1_zher(uplo1_t uplo, conj1_t conj, int m, double *alpha, dcomplex *x, int incx, dcomplex *a, int a_rs, int a_cs)
Definition bl1_her.c:101
void bl1_ssyr(uplo1_t uplo, int m, float *alpha, float *x, int incx, float *a, int a_rs, int a_cs)
Definition bl1_syr.c:13
void bl1_dsyr(uplo1_t uplo, int m, double *alpha, double *x, int incx, double *a, int a_rs, int a_cs)
Definition bl1_syr.c:57
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_cher(), bl1_dsyr(), bl1_ssyr(), bl1_zher(), FLA_Check_error_level(), FLA_Her_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_Chol_l_unb_var3(), FLA_Her(), FLA_Herk_ln_unb_var5(), FLA_Herk_ln_unb_var6(), FLA_Herk_un_unb_var5(), FLA_Herk_un_unb_var6(), and FLA_Ttmm_u_unb_var1().