libflame  revision_anchor
Functions
FLA_Dotcs_external.c File Reference

(r)

Functions

FLA_Error FLA_Dotcs_external (FLA_Conj conj, FLA_Obj alpha, FLA_Obj x, FLA_Obj y, FLA_Obj beta, FLA_Obj rho)
 

Function Documentation

◆ FLA_Dotcs_external()

FLA_Error FLA_Dotcs_external ( FLA_Conj  conj,
FLA_Obj  alpha,
FLA_Obj  x,
FLA_Obj  y,
FLA_Obj  beta,
FLA_Obj  rho 
)

References bl1_cdots(), bl1_ddots(), bl1_sdots(), bl1_zdots(), FLA_Check_error_level(), FLA_Dotcs_check(), FLA_Obj_datatype(), FLA_Obj_has_zero_dim(), FLA_Obj_vector_dim(), FLA_Obj_vector_inc(), FLA_Param_map_flame_to_blis_conj(), and FLA_Scal_external().

Referenced by FLA_Chol_l_unb_var1(), FLA_Chol_l_unb_var2(), FLA_Chol_u_unb_var1(), FLA_Chol_u_unb_var2(), FLA_Dotcs(), FLA_Herk_lh_unb_var1(), FLA_Herk_lh_unb_var2(), FLA_Herk_lh_unb_var3(), FLA_Herk_lh_unb_var4(), FLA_Herk_ln_unb_var1(), FLA_Herk_ln_unb_var2(), FLA_Herk_ln_unb_var3(), FLA_Herk_ln_unb_var4(), FLA_Herk_uh_unb_var1(), FLA_Herk_uh_unb_var2(), FLA_Herk_uh_unb_var3(), FLA_Herk_uh_unb_var4(), FLA_Herk_un_unb_var1(), FLA_Herk_un_unb_var2(), FLA_Herk_un_unb_var3(), FLA_Herk_un_unb_var4(), FLA_Ttmm_l_unb_var2(), FLA_Ttmm_l_unb_var3(), FLA_Ttmm_u_unb_var2(), and FLA_Ttmm_u_unb_var3().

14 {
15  FLA_Datatype datatype;
16  int num_elem;
17  int inc_x;
18  int inc_y;
19  conj1_t blis_conj;
20 
21  if ( FLA_Check_error_level() == FLA_FULL_ERROR_CHECKING )
22  FLA_Dotcs_check( conj, alpha, x, y, beta, rho );
23 
24  if ( FLA_Obj_has_zero_dim( x ) )
25  {
26  FLA_Scal_external( beta, rho );
27  return FLA_SUCCESS;
28  }
29 
30  datatype = FLA_Obj_datatype( x );
31 
32  inc_x = FLA_Obj_vector_inc( x );
33  inc_y = FLA_Obj_vector_inc( y );
34  num_elem = FLA_Obj_vector_dim( x );
35 
36  FLA_Param_map_flame_to_blis_conj( conj, &blis_conj );
37 
38  switch ( datatype ){
39 
40  case FLA_FLOAT:
41  {
42  float *buff_x = ( float * ) FLA_FLOAT_PTR( x );
43  float *buff_y = ( float * ) FLA_FLOAT_PTR( y );
44  float *buff_rho = ( float * ) FLA_FLOAT_PTR( rho );
45  float *buff_alpha = ( float * ) FLA_FLOAT_PTR( alpha );
46  float *buff_beta = ( float * ) FLA_FLOAT_PTR( beta );
47 
48  bl1_sdots( blis_conj,
49  num_elem,
50  buff_alpha,
51  buff_x, inc_x,
52  buff_y, inc_y,
53  buff_beta,
54  buff_rho );
55 
56  break;
57  }
58 
59  case FLA_DOUBLE:
60  {
61  double *buff_x = ( double * ) FLA_DOUBLE_PTR( x );
62  double *buff_y = ( double * ) FLA_DOUBLE_PTR( y );
63  double *buff_rho = ( double * ) FLA_DOUBLE_PTR( rho );
64  double *buff_alpha = ( double * ) FLA_DOUBLE_PTR( alpha );
65  double *buff_beta = ( double * ) FLA_DOUBLE_PTR( beta );
66 
67  bl1_ddots( blis_conj,
68  num_elem,
69  buff_alpha,
70  buff_x, inc_x,
71  buff_y, inc_y,
72  buff_beta,
73  buff_rho );
74 
75  break;
76  }
77 
78  case FLA_COMPLEX:
79  {
80  scomplex *buff_x = ( scomplex * ) FLA_COMPLEX_PTR( x );
81  scomplex *buff_y = ( scomplex * ) FLA_COMPLEX_PTR( y );
82  scomplex *buff_rho = ( scomplex * ) FLA_COMPLEX_PTR( rho );
83  scomplex *buff_alpha = ( scomplex * ) FLA_COMPLEX_PTR( alpha );
84  scomplex *buff_beta = ( scomplex * ) FLA_COMPLEX_PTR( beta );
85 
86  bl1_cdots( blis_conj,
87  num_elem,
88  buff_alpha,
89  buff_x, inc_x,
90  buff_y, inc_y,
91  buff_beta,
92  buff_rho );
93 
94  break;
95  }
96 
97  case FLA_DOUBLE_COMPLEX:
98  {
99  dcomplex *buff_x = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( x );
100  dcomplex *buff_y = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( y );
101  dcomplex *buff_rho = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( rho );
102  dcomplex *buff_alpha = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( alpha );
103  dcomplex *buff_beta = ( dcomplex * ) FLA_DOUBLE_COMPLEX_PTR( beta );
104 
105  bl1_zdots( blis_conj,
106  num_elem,
107  buff_alpha,
108  buff_x, inc_x,
109  buff_y, inc_y,
110  buff_beta,
111  buff_rho );
112 
113  break;
114  }
115 
116  }
117 
118  return FLA_SUCCESS;
119 }
FLA_Error FLA_Scal_external(FLA_Obj alpha, FLA_Obj A)
Definition: FLA_Scal_external.c:13
conj1_t
Definition: blis_type_defs.h:79
void bl1_cdots(conj1_t conj, int n, scomplex *alpha, scomplex *x, int incx, scomplex *y, int incy, scomplex *beta, scomplex *rho)
Definition: bl1_dots.c:39
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
FLA_Bool FLA_Obj_has_zero_dim(FLA_Obj A)
Definition: FLA_Query.c:400
void bl1_sdots(conj1_t conj, int n, float *alpha, float *x, int incx, float *y, int incy, float *beta, float *rho)
Definition: bl1_dots.c:13
Definition: blis_type_defs.h:132
dim_t FLA_Obj_vector_dim(FLA_Obj obj)
Definition: FLA_Query.c:137
unsigned int FLA_Check_error_level(void)
Definition: FLA_Check.c:18
void FLA_Param_map_flame_to_blis_conj(FLA_Conj conj, conj1_t *blis_conj)
Definition: FLA_Param.c:269
int FLA_Datatype
Definition: FLA_type_defs.h:49
void bl1_ddots(conj1_t conj, int n, double *alpha, double *x, int incx, double *y, int incy, double *beta, double *rho)
Definition: bl1_dots.c:26
FLA_Error FLA_Dotcs_check(FLA_Conj conj, FLA_Obj alpha, FLA_Obj x, FLA_Obj y, FLA_Obj beta, FLA_Obj rho)
Definition: FLA_Dotcs_check.c:13
dim_t FLA_Obj_vector_inc(FLA_Obj obj)
Definition: FLA_Query.c:145
void bl1_zdots(conj1_t conj, int n, dcomplex *alpha, dcomplex *x, int incx, dcomplex *y, int incy, dcomplex *beta, dcomplex *rho)
Definition: bl1_dots.c:56
Definition: blis_type_defs.h:137