libflame revision_anchor
FLA_Symm_external.c File Reference

(r)

Functions

FLA_Error FLA_Symm_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_Symm_external()

FLA_Error FLA_Symm_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_Symm_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_csymm( 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_zsymm( 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_Symm_check(FLA_Side side, FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj A, FLA_Obj B, FLA_Obj beta, FLA_Obj C)
Definition FLA_Symm_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_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_zsymm(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_symm.c:796
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
void bl1_csymm(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_symm.c:535
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_csymm(), bl1_dsymm(), bl1_ssymm(), bl1_zsymm(), 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_side(), FLA_Param_map_flame_to_blis_uplo(), and FLA_Symm_check().

Referenced by FLA_Symm(), FLA_Symm_ll_task(), FLA_Symm_lu_task(), FLA_Symm_rl_task(), FLA_Symm_ru_task(), and FLA_Symm_task().