1 /*
2
3 Copyright (C) 2014, The University of Texas at Austin
4
5 This file is part of libflame and is available under the 3-Clause
6 BSD license, which can be found in the LICENSE file at the top-level
7 directory, or at http://opensource.org/licenses/BSD-3-Clause
8
9 */
10
11 #include "FLAME.h"
12
FLA_Dot2cs_external(FLA_Conj conj,FLA_Obj alpha,FLA_Obj x,FLA_Obj y,FLA_Obj beta,FLA_Obj rho)13 FLA_Error FLA_Dot2cs_external( FLA_Conj conj, FLA_Obj alpha, FLA_Obj x, FLA_Obj y, FLA_Obj beta, FLA_Obj rho )
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_Dot2cs_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_sdot2s( 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_ddot2s( 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_cdot2s( 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_zdot2s( 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 }
120
121