libflame revision_anchor
Functions
FLA_Hemv_external.c File Reference

(r)

Functions

FLA_Error FLA_Hemv_external (FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj A, FLA_Obj x, FLA_Obj beta, FLA_Obj y)
 

Function Documentation

◆ FLA_Hemv_external()

FLA_Error FLA_Hemv_external ( FLA_Uplo uplo,
FLA_Obj alpha,
FLA_Obj A,
FLA_Obj x,
FLA_Obj beta,
FLA_Obj y )
14{
15 FLA_Datatype datatype;
16 int m_A;
17 int rs_A, cs_A;
18 int inc_x;
19 int inc_y;
20 uplo1_t blis_uplo;
21 conj1_t blis_conj;
22
23 if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
24 FLA_Hemv_check( uplo, alpha, A, x, beta, y );
25
26 if ( FLA_Obj_has_zero_dim( A ) ) return FLA_SUCCESS;
27
28 datatype = FLA_Obj_datatype( A );
29
30 m_A = FLA_Obj_length( A );
31 rs_A = FLA_Obj_row_stride( A );
32 cs_A = FLA_Obj_col_stride( A );
33
34 inc_x = FLA_Obj_vector_inc( x );
35 inc_y = FLA_Obj_vector_inc( y );
36
37 FLA_Param_map_flame_to_blis_uplo( uplo, &blis_uplo );
38 FLA_Param_map_flame_to_blis_conj( FLA_NO_CONJUGATE, &blis_conj );
39
40
41 switch( datatype ){
42
43 case FLA_FLOAT:
44 {
45 float *buff_A = ( float * ) FLA_FLOAT_PTR( A );
46 float *buff_x = ( float * ) FLA_FLOAT_PTR( x );
47 float *buff_y = ( float * ) FLA_FLOAT_PTR( y );
48 float *buff_alpha = ( float * ) FLA_FLOAT_PTR( alpha );
49 float *buff_beta = ( float * ) FLA_FLOAT_PTR( beta );
50
51 bl1_ssymv( blis_uplo,
52 m_A,
53 buff_alpha,
54 buff_A, rs_A, cs_A,
55 buff_x, inc_x,
56 buff_beta,
57 buff_y, inc_y );
58
59 break;
60 }
61
62 case FLA_DOUBLE:
63 {
64 double *buff_A = ( double * ) FLA_DOUBLE_PTR( A );
65 double *buff_x = ( double * ) FLA_DOUBLE_PTR( x );
66 double *buff_y = ( double * ) FLA_DOUBLE_PTR( y );
67 double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
68 double *buff_beta = ( double * ) FLA_DOUBLE_PTR( beta );
69
70 bl1_dsymv( blis_uplo,
71 m_A,
72 buff_alpha,
73 buff_A, rs_A, cs_A,
74 buff_x, inc_x,
75 buff_beta,
76 buff_y, inc_y );
77
78 break;
79 }
80
81 case FLA_COMPLEX:
82 {
83 scomplex *buff_A = ( scomplex * ) FLA_COMPLEX_PTR( A );
84 scomplex *buff_x = ( scomplex * ) FLA_COMPLEX_PTR( x );
85 scomplex *buff_y = ( scomplex * ) FLA_COMPLEX_PTR( y );
86 scomplex *buff_alpha = ( scomplex * ) FLA_COMPLEX_PTR( alpha );
87 scomplex *buff_beta = ( scomplex * ) FLA_COMPLEX_PTR( beta );
88
89 bl1_chemv( blis_uplo,
90 blis_conj,
91 m_A,
92 buff_alpha,
93 buff_A, rs_A, cs_A,
94 buff_x, inc_x,
95 buff_beta,
96 buff_y, inc_y );
97
98 break;
99 }
100
101 case FLA_DOUBLE_COMPLEX:
102 {
103 dcomplex *buff_A = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( A );
104 dcomplex *buff_x = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( x );
105 dcomplex *buff_y = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( y );
106 dcomplex *buff_alpha = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( alpha );
107 dcomplex *buff_beta = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( beta );
108
109 bl1_zhemv( blis_uplo,
110 blis_conj,
111 m_A,
112 buff_alpha,
113 buff_A, rs_A, cs_A,
114 buff_x, inc_x,
115 buff_beta,
116 buff_y, inc_y );
117
118 break;
119 }
120
121 }
122
123 return FLA_SUCCESS;
124}
FLA_Hemv_check
FLA_Error FLA_Hemv_check(FLA_Uplo uplo, FLA_Obj alpha, FLA_Obj A, FLA_Obj x, FLA_Obj beta, FLA_Obj y)
Definition FLA_Hemv_check.c:13
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_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_conj
void FLA_Param_map_flame_to_blis_conj(FLA_Conj conj, conj1_t *blis_conj)
Definition FLA_Param.c:269
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_vector_inc
dim_t FLA_Obj_vector_inc(FLA_Obj obj)
Definition FLA_Query.c:145
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_zhemv
void bl1_zhemv(uplo1_t uplo, conj1_t conj, int m, dcomplex *alpha, dcomplex *a, int a_rs, int a_cs, dcomplex *x, int incx, dcomplex *beta, dcomplex *y, int incy)
Definition bl1_hemv.c:134
bl1_chemv
void bl1_chemv(uplo1_t uplo, conj1_t conj, int m, scomplex *alpha, scomplex *a, int a_rs, int a_cs, scomplex *x, int incx, scomplex *beta, scomplex *y, int incy)
Definition bl1_hemv.c:35
bl1_dsymv
void bl1_dsymv(uplo1_t uplo, int m, double *alpha, double *a, int a_rs, int a_cs, double *x, int incx, double *beta, double *y, int incy)
Definition bl1_symv.c:56
bl1_ssymv
void bl1_ssymv(uplo1_t uplo, int m, float *alpha, float *a, int a_rs, int a_cs, float *x, int incx, float *beta, float *y, int incy)
Definition bl1_symv.c:13
uplo1_t
uplo1_t
Definition blis_type_defs.h:61
conj1_t
conj1_t
Definition blis_type_defs.h:80
dcomplex
Definition blis_type_defs.h:138
scomplex
Definition blis_type_defs.h:133

References bl1_chemv(), bl1_dsymv(), bl1_ssymv(), bl1_zhemv(), FLA_Check_error_level(), FLA_Hemv_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_Eig_gest_nl_unb_var1(), FLA_Eig_gest_nl_unb_var2(), FLA_Hemm_ll_unb_var10(), FLA_Hemm_ll_unb_var9(), FLA_Hemm_lu_unb_var10(), FLA_Hemm_lu_unb_var9(), and FLA_Hemv().

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