libflame  revision_anchor
Functions
FLA_Syrk_external.c File Reference

(r)

Functions

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

Function Documentation

◆ FLA_Syrk_external()

FLA_Error FLA_Syrk_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_Syrk_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  scomplex *buff_alpha = ( scomplex * ) FLA_COMPLEX_PTR( alpha );
100  scomplex *buff_beta = ( scomplex * ) FLA_COMPLEX_PTR( beta );
101 
102  bl1_csyrk( 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  dcomplex *buff_alpha = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( alpha );
119  dcomplex *buff_beta = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( beta );
120 
121  bl1_zsyrk( 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_Error FLA_Syrk_check(FLA_Uplo uplo, FLA_Trans trans, FLA_Obj alpha, FLA_Obj A, FLA_Obj beta, FLA_Obj C)
Definition: FLA_Syrk_check.c:13
FLA_Error FLA_Scal_external(FLA_Obj alpha, FLA_Obj A)
Definition: FLA_Scal_external.c:13
dim_t FLA_Obj_width(FLA_Obj obj)
Definition: FLA_Query.c:123
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
void FLA_Param_map_flame_to_blis_trans(FLA_Trans trans, trans1_t *blis_trans)
Definition: FLA_Param.c:245
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_uplo(FLA_Uplo uplo, uplo1_t *blis_uplo)
Definition: FLA_Param.c:285
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
int FLA_Datatype
Definition: FLA_type_defs.h:49
void bl1_csyrk(uplo1_t uplo, trans1_t trans, int m, int k, scomplex *alpha, scomplex *a, int a_rs, int a_cs, scomplex *beta, scomplex *c, int c_rs, int c_cs)
Definition: bl1_syrk.c:205
void bl1_zsyrk(uplo1_t uplo, trans1_t trans, int m, int k, dcomplex *alpha, dcomplex *a, int a_rs, int a_cs, dcomplex *beta, dcomplex *c, int c_rs, int c_cs)
Definition: bl1_syrk.c:301
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
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
Definition: blis_type_defs.h:61
trans1_t
Definition: blis_type_defs.h:53
Definition: blis_type_defs.h:138
Definition: blis_type_defs.h:133

References bl1_csyrk(), bl1_dsyrk(), bl1_ssyrk(), bl1_zsyrk(), 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_width(), FLA_Param_map_flame_to_blis_trans(), FLA_Param_map_flame_to_blis_uplo(), FLA_Scal_external(), and FLA_Syrk_check().

Referenced by FLA_Syrk(), FLA_Syrk_ln_task(), FLA_Syrk_lt_task(), FLA_Syrk_task(), FLA_Syrk_un_task(), and FLA_Syrk_ut_task().