31 int i, m, inc_a, inc_b, inc_d, inc_e;
47 float* buff_d = FLA_FLOAT_PTR( d );
48 float* buff_e = FLA_FLOAT_PTR( e );
49 float* buff_1 = FLA_FLOAT_PTR(
FLA_ONE );
63 double* buff_d = FLA_DOUBLE_PTR( d );
64 double* buff_e = FLA_DOUBLE_PTR( e );
65 double* buff_1 = FLA_DOUBLE_PTR(
FLA_ONE );
79 scomplex* buff_a = FLA_COMPLEX_PTR( a );
80 scomplex* buff_b = ( m > 1 ? FLA_COMPLEX_PTR( b ) : NULL );
81 scomplex* buff_d = FLA_COMPLEX_PTR( d );
82 scomplex* buff_e = FLA_COMPLEX_PTR( e );
84 float* buff_0 = FLA_FLOAT_PTR(
FLA_ZERO );
86 for ( i = 0; i < m; ++
i )
100 scomplex* beta1 = buff_b + (i-1)*inc_b;
101 if ( beta1->
imag == 0.0F )
106 bl1_cabsval2( beta1, &absv );
107 bl1_cinvscals( &absv, delta1 );
109 bl1_cscals( delta1, beta1 );
110 beta1->
imag = *buff_0;
112 bl1_cscals( delta1, alpha1 );
116 if ( alpha1->
imag == 0.0F )
121 bl1_cabsval2( alpha1, &absv );
122 bl1_cinvscals( &absv, epsilon1 );
124 bl1_cscals( epsilon1, alpha1 );
125 alpha1->
imag = *buff_0;
131 bl1_cscals( epsilon1, beta2 );
136 case FLA_DOUBLE_COMPLEX:
138 dcomplex* buff_a = FLA_DOUBLE_COMPLEX_PTR( a );
139 dcomplex* buff_b = ( m > 1 ? FLA_DOUBLE_COMPLEX_PTR( b ) : NULL );
140 dcomplex* buff_d = FLA_DOUBLE_COMPLEX_PTR( d );
141 dcomplex* buff_e = FLA_DOUBLE_COMPLEX_PTR( e );
143 double* buff_0 = FLA_DOUBLE_PTR(
FLA_ZERO );
145 for ( i = 0; i < m; ++
i )
149 dcomplex* epsilon1 = buff_e + (
i )*inc_e;
159 dcomplex* beta1 = buff_b + (i-1)*inc_b;
161 bl1_zabsval2( beta1, &absv );
162 bl1_zinvscals( &absv, delta1 );
165 beta1->
imag = *buff_0;
171 bl1_zabsval2( alpha1, &absv );
172 bl1_zinvscals( &absv, epsilon1 );
175 alpha1->
imag = *buff_0;
double imag
Definition: blis_type_defs.h:139
double *restrict alpha1
Definition: bl1_axpyv2bdotaxpy.c:198
Definition: blis_type_defs.h:82
FLA_Obj FLA_ONE
Definition: FLA_Init.c:18
void bl1_dsetv(int m, double *sigma, double *x, int incx)
Definition: bl1_setv.c:39
FLA_Datatype FLA_Obj_datatype(FLA_Obj obj)
Definition: FLA_Query.c:13
void bl1_ssetv(int m, float *sigma, float *x, int incx)
Definition: bl1_setv.c:26
Definition: blis_type_defs.h:132
dim_t FLA_Obj_vector_dim(FLA_Obj obj)
Definition: FLA_Query.c:137
int FLA_Datatype
Definition: FLA_type_defs.h:49
int i
Definition: bl1_axmyv2.c:145
dim_t FLA_Obj_vector_inc(FLA_Obj obj)
Definition: FLA_Query.c:145
float imag
Definition: blis_type_defs.h:134
FLA_Obj FLA_ZERO
Definition: FLA_Init.c:20
Definition: blis_type_defs.h:137