libflame revision_anchor
Functions
FLA_Herk_external.c File Reference

(r)

Functions

FLA_Error FLA_Herk_external (FLA_Uplo uplo, FLA_Trans trans, FLA_Obj alpha, FLA_Obj A, FLA_Obj beta, FLA_Obj C)
 

Function Documentation

◆ FLA_Herk_external()

FLA_Error FLA_Herk_external ( FLA_Uplo uplo,
FLA_Trans trans,
FLA_Obj alpha,
FLA_Obj A,
FLA_Obj beta,
FLA_Obj C )
14{
15 FLA_Datatype datatype;
16 int k_A;
17 int m_A, n_A;
18 int m_C;
19 int rs_A, cs_A;
20 int rs_C, cs_C;
21 uplo1_t blis_uplo;
22 trans1_t blis_trans;
23
24 if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
25 FLA_Herk_check( uplo, trans, alpha, A, beta, C );
26
27 if ( FLA_Obj_has_zero_dim( C ) ) return FLA_SUCCESS;
28
29 if ( FLA_Obj_has_zero_dim( A ) )
30 {
31 FLA_Scal_external( beta, C );
32 return FLA_SUCCESS;
33 }
34
35 datatype = FLA_Obj_datatype( A );
36
37 m_A = FLA_Obj_length( A );
38 n_A = FLA_Obj_width( A );
39 rs_A = FLA_Obj_row_stride( A );
40 cs_A = FLA_Obj_col_stride( A );
41
42 m_C = FLA_Obj_length( C );
43 rs_C = FLA_Obj_row_stride( C );
44 cs_C = FLA_Obj_col_stride( C );
45
46 if ( trans == FLA_NO_TRANSPOSE )
47 k_A = n_A;
48 else
49 k_A = m_A;
50
51 FLA_Param_map_flame_to_blis_uplo( uplo, &blis_uplo );
52 FLA_Param_map_flame_to_blis_trans( trans, &blis_trans );
53
54
55 switch( datatype ){
56
57 case FLA_FLOAT:
58 {
59 float *buff_A = ( float * ) FLA_FLOAT_PTR( A );
60 float *buff_C = ( float * ) FLA_FLOAT_PTR( C );
61 float *buff_alpha = ( float * ) FLA_FLOAT_PTR( alpha );
62 float *buff_beta = ( float * ) FLA_FLOAT_PTR( beta );
63
64 bl1_ssyrk( blis_uplo,
65 blis_trans,
66 m_C,
67 k_A,
68 buff_alpha,
69 buff_A, rs_A, cs_A,
70 buff_beta,
71 buff_C, rs_C, cs_C );
72
73 break;
74 }
75
76 case FLA_DOUBLE:
77 {
78 double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
79 double *buff_C = ( double * ) FLA_DOUBLE_PTR( C );
80 double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
81 double *buff_beta = ( double * ) FLA_DOUBLE_PTR( beta );
82
83 bl1_dsyrk( blis_uplo,
84 blis_trans,
85 m_C,
86 k_A,
87 buff_alpha,
88 buff_A, rs_A, cs_A,
89 buff_beta,
90 buff_C, rs_C, cs_C );
91
92 break;
93 }
94
95 case FLA_COMPLEX:
96 {
97 scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
98 scomplex *buff_C = ( scomplex * ) FLA_COMPLEX_PTR( C );
99 float *buff_alpha = ( float * ) FLA_FLOAT_PTR( alpha );
100 float *buff_beta = ( float * ) FLA_FLOAT_PTR( beta );
101
102 bl1_cherk( blis_uplo,
103 blis_trans,
104 m_C,
105 k_A,
106 buff_alpha,
107 buff_A, rs_A, cs_A,
108 buff_beta,
109 buff_C, rs_C, cs_C );
110
111 break;
112 }
113
114 case FLA_DOUBLE_COMPLEX:
115 {
116 dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
117 dcomplex *buff_C = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( C );
118 double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
119 double *buff_beta = ( double * ) FLA_DOUBLE_PTR( beta );
120
121 bl1_zherk( blis_uplo,
122 blis_trans,
123 m_C,
124 k_A,
125 buff_alpha,
126 buff_A, rs_A, cs_A,
127 buff_beta,
128 buff_C, rs_C, cs_C );
129
130 break;
131 }
132
133 }
134
135 return FLA_SUCCESS;
136}
FLA_Herk_check
FLA_Error FLA_Herk_check(FLA_Uplo uplo, FLA_Trans trans, FLA_Obj alpha, FLA_Obj A, FLA_Obj beta, FLA_Obj C)
Definition FLA_Herk_check.c:13
FLA_Scal_external
FLA_Error FLA_Scal_external(FLA_Obj alpha, FLA_Obj A)
Definition FLA_Scal_external.c:13
FLA_Obj_width
dim_t FLA_Obj_width(FLA_Obj obj)
Definition FLA_Query.c:123
FLA_Obj_has_zero_dim
FLA_Bool FLA_Obj_has_zero_dim(FLA_Obj A)
Definition FLA_Query.c:400
FLA_Obj_row_stride
dim_t FLA_Obj_row_stride(FLA_Obj obj)
Definition FLA_Query.c:167
FLA_Obj_length
dim_t FLA_Obj_length(FLA_Obj obj)
Definition FLA_Query.c:116
FLA_Param_map_flame_to_blis_trans
void FLA_Param_map_flame_to_blis_trans(FLA_Trans trans, trans1_t *blis_trans)
Definition FLA_Param.c:245
FLA_Obj_col_stride
dim_t FLA_Obj_col_stride(FLA_Obj obj)
Definition FLA_Query.c:174
FLA_Check_error_level
unsigned int FLA_Check_error_level(void)
Definition FLA_Check.c:18
FLA_Param_map_flame_to_blis_uplo
void FLA_Param_map_flame_to_blis_uplo(FLA_Uplo uplo, uplo1_t *blis_uplo)
Definition FLA_Param.c:285
FLA_Obj_datatype
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition FLA_Query.c:13
FLA_Datatype
int FLA_Datatype
Definition FLA_type_defs.h:49
bl1_cherk
void bl1_cherk(uplo1_t uplo, trans1_t trans, int m, int k, float *alpha, scomplex *a, int a_rs, int a_cs, float *beta, scomplex *c, int c_rs, int c_cs)
Definition bl1_herk.c:36
bl1_zherk
void bl1_zherk(uplo1_t uplo, trans1_t trans, int m, int k, double *alpha, dcomplex *a, int a_rs, int a_cs, double *beta, dcomplex *c, int c_rs, int c_cs)
Definition bl1_herk.c:184
bl1_dsyrk
void bl1_dsyrk(uplo1_t uplo, trans1_t trans, int m, int k, double *alpha, double *a, int a_rs, int a_cs, double *beta, double *c, int c_rs, int c_cs)
Definition bl1_syrk.c:109
bl1_ssyrk
void bl1_ssyrk(uplo1_t uplo, trans1_t trans, int m, int k, float *alpha, float *a, int a_rs, int a_cs, float *beta, float *c, int c_rs, int c_cs)
Definition bl1_syrk.c:13
uplo1_t
uplo1_t
Definition blis_type_defs.h:61
trans1_t
trans1_t
Definition blis_type_defs.h:53
dcomplex
Definition blis_type_defs.h:138
scomplex
Definition blis_type_defs.h:133

References bl1_cherk(), bl1_dsyrk(), bl1_ssyrk(), bl1_zherk(), FLA_Check_error_level(), FLA_Herk_check(), FLA_Obj_col_stride(), FLA_Obj_datatype(), FLA_Obj_has_zero_dim(), FLA_Obj_length(), FLA_Obj_row_stride(), FLA_Obj_width(), FLA_Param_map_flame_to_blis_trans(), FLA_Param_map_flame_to_blis_uplo(), and FLA_Scal_external().

Referenced by FLA_Herk(), FLA_Herk_lh_task(), FLA_Herk_ln_task(), FLA_Herk_task(), FLA_Herk_uh_task(), FLA_Herk_un_task(), FLA_Random_spd_matrix(), and FLA_UDdate_UT_unb_var1().

Generated on Thu Sep 23 2021 17:46:17 for libflame by doxygen 1.12.0