libflame  revision_anchor
Functions
FLA_Hemm_external.c File Reference

(r)

Functions

FLA_Error FLA_Hemm_external (FLA_Side side, FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj A, FLA_Obj B, FLA_Obj beta, FLA_Obj C)
 

Function Documentation

◆ FLA_Hemm_external()

FLA_Error FLA_Hemm_external ( FLA_Side  side,
FLA_Uplo  uplo,
FLA_Obj  alpha,
FLA_Obj  A,
FLA_Obj  B,
FLA_Obj  beta,
FLA_Obj  C 
)
14 {
15  FLA_Datatype datatype;
16  int m_C, n_C;
17  int rs_A, cs_A;
18  int rs_B, cs_B;
19  int rs_C, cs_C;
20  side1_t blis_side;
21  uplo1_t blis_uplo;
22 
23  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
24  FLA_Hemm_check( side, uplo, alpha, A, B, beta, C );
25 
26  if ( FLA_Obj_has_zero_dim( C ) ) return FLA_SUCCESS;
27 
28  datatype = FLA_Obj_datatype( A );
29 
30  rs_A = FLA_Obj_row_stride( A );
31  cs_A = FLA_Obj_col_stride( A );
32 
33  rs_B = FLA_Obj_row_stride( B );
34  cs_B = FLA_Obj_col_stride( B );
35 
36  m_C = FLA_Obj_length( C );
37  n_C = FLA_Obj_width( C );
38  rs_C = FLA_Obj_row_stride( C );
39  cs_C = FLA_Obj_col_stride( C );
40 
41  FLA_Param_map_flame_to_blis_side( side, &blis_side );
42  FLA_Param_map_flame_to_blis_uplo( uplo, &blis_uplo );
43 
44 
45  switch( datatype ){
46 
47  case FLA_FLOAT:
48  {
49  float *buff_A = ( float * ) FLA_FLOAT_PTR( A );
50  float *buff_B = ( float * ) FLA_FLOAT_PTR( B );
51  float *buff_C = ( float * ) FLA_FLOAT_PTR( C );
52  float *buff_alpha = ( float * ) FLA_FLOAT_PTR( alpha );
53  float *buff_beta = ( float * ) FLA_FLOAT_PTR( beta );
54 
55  bl1_ssymm( blis_side,
56  blis_uplo,
57  m_C,
58  n_C,
59  buff_alpha,
60  buff_A, rs_A, cs_A,
61  buff_B, rs_B, cs_B,
62  buff_beta,
63  buff_C, rs_C, cs_C );
64 
65  break;
66  }
67 
68  case FLA_DOUBLE:
69  {
70  double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
71  double *buff_B = ( double * ) FLA_DOUBLE_PTR( B );
72  double *buff_C = ( double * ) FLA_DOUBLE_PTR( C );
73  double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
74  double *buff_beta = ( double * ) FLA_DOUBLE_PTR( beta );
75 
76  bl1_dsymm( blis_side,
77  blis_uplo,
78  m_C,
79  n_C,
80  buff_alpha,
81  buff_A, rs_A, cs_A,
82  buff_B, rs_B, cs_B,
83  buff_beta,
84  buff_C, rs_C, cs_C );
85 
86  break;
87  }
88 
89  case FLA_COMPLEX:
90  {
91  scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
92  scomplex *buff_B = ( scomplex * ) FLA_COMPLEX_PTR( B );
93  scomplex *buff_C = ( scomplex * ) FLA_COMPLEX_PTR( C );
94  scomplex *buff_alpha = ( scomplex * ) FLA_COMPLEX_PTR( alpha );
95  scomplex *buff_beta = ( scomplex * ) FLA_COMPLEX_PTR( beta );
96 
97  bl1_chemm( blis_side,
98  blis_uplo,
99  m_C,
100  n_C,
101  buff_alpha,
102  buff_A, rs_A, cs_A,
103  buff_B, rs_B, cs_B,
104  buff_beta,
105  buff_C, rs_C, cs_C );
106 
107  break;
108  }
109 
110  case FLA_DOUBLE_COMPLEX:
111  {
112  dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
113  dcomplex *buff_B = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( B );
114  dcomplex *buff_C = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( C );
115  dcomplex *buff_alpha = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( alpha );
116  dcomplex *buff_beta = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( beta );
117 
118  bl1_zhemm( blis_side,
119  blis_uplo,
120  m_C,
121  n_C,
122  buff_alpha,
123  buff_A, rs_A, cs_A,
124  buff_B, rs_B, cs_B,
125  buff_beta,
126  buff_C, rs_C, cs_C );
127 
128  break;
129  }
130 
131  }
132 
133  return FLA_SUCCESS;
134 }
FLA_Error FLA_Hemm_check(FLA_Side side, FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj A, FLA_Obj B, FLA_Obj beta, FLA_Obj C)
Definition: FLA_Hemm_check.c:13
void FLA_Param_map_flame_to_blis_side(FLA_Uplo side, side1_t *blis_side)
Definition: FLA_Param.c:301
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
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_chemm(side1_t side, uplo1_t uplo, int m, int n, scomplex *alpha, scomplex *a, int a_rs, int a_cs, scomplex *b, int b_rs, int b_cs, scomplex *beta, scomplex *c, int c_rs, int c_cs)
Definition: bl1_hemm.c:38
void bl1_zhemm(side1_t side, uplo1_t uplo, int m, int n, dcomplex *alpha, dcomplex *a, int a_rs, int a_cs, dcomplex *b, int b_rs, int b_cs, dcomplex *beta, dcomplex *c, int c_rs, int c_cs)
Definition: bl1_hemm.c:348
void bl1_ssymm(side1_t side, uplo1_t uplo, int m, int n, float *alpha, float *a, int a_rs, int a_cs, float *b, int b_rs, int b_cs, float *beta, float *c, int c_rs, int c_cs)
Definition: bl1_symm.c:13
void bl1_dsymm(side1_t side, uplo1_t uplo, int m, int n, double *alpha, double *a, int a_rs, int a_cs, double *b, int b_rs, int b_cs, double *beta, double *c, int c_rs, int c_cs)
Definition: bl1_symm.c:274
uplo1_t
Definition: blis_type_defs.h:61
side1_t
Definition: blis_type_defs.h:67
Definition: blis_type_defs.h:138
Definition: blis_type_defs.h:133

References bl1_chemm(), bl1_dsymm(), bl1_ssymm(), bl1_zhemm(), FLA_Check_error_level(), FLA_Hemm_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_side(), and FLA_Param_map_flame_to_blis_uplo().

Referenced by FLA_Hemm(), FLA_Hemm_ll_task(), FLA_Hemm_lu_task(), FLA_Hemm_rl_task(), FLA_Hemm_ru_task(), and FLA_Hemm_task().