libflame revision_anchor
Functions
FLA_Fused_Gerc2_Ahx_Ax_opt_var1.c File Reference

(r)

Functions

FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opt_var1 (FLA_Obj alpha, FLA_Obj u, FLA_Obj y, FLA_Obj z, FLA_Obj A, FLA_Obj x, FLA_Obj v, FLA_Obj w)
 
FLA_Error FLA_Fused_Gerc2_Ahx_Ax_ops_var1 (int m_A, int n_A, float *buff_alpha, float *buff_u, int inc_u, float *buff_y, int inc_y, float *buff_z, int inc_z, float *buff_A, int rs_A, int cs_A, float *buff_x, int inc_x, float *buff_v, int inc_v, float *buff_w, int inc_w)
 
FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opd_var1 (int m_A, int n_A, double *buff_alpha, double *buff_u, int inc_u, double *buff_y, int inc_y, double *buff_z, int inc_z, double *buff_A, int rs_A, int cs_A, double *buff_x, int inc_x, double *buff_v, int inc_v, double *buff_w, int inc_w)
 
FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opc_var1 (int m_A, int n_A, scomplex *buff_alpha, scomplex *buff_u, int inc_u, scomplex *buff_y, int inc_y, scomplex *buff_z, int inc_z, scomplex *buff_A, int rs_A, int cs_A, scomplex *buff_x, int inc_x, scomplex *buff_v, int inc_v, scomplex *buff_w, int inc_w)
 
FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opz_var1 (int m_A, int n_A, dcomplex *buff_alpha, dcomplex *buff_u, int inc_u, dcomplex *buff_y, int inc_y, dcomplex *buff_z, int inc_z, dcomplex *buff_A, int rs_A, int cs_A, dcomplex *buff_x, int inc_x, dcomplex *buff_v, int inc_v, dcomplex *buff_w, int inc_w)
 

Function Documentation

◆ FLA_Fused_Gerc2_Ahx_Ax_opc_var1()

FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opc_var1 ( int  m_A,
int  n_A,
scomplex *  buff_alpha,
scomplex *  buff_u,
int  inc_u,
scomplex *  buff_y,
int  inc_y,
scomplex *  buff_z,
int  inc_z,
scomplex *  buff_A,
int  rs_A,
int  cs_A,
scomplex *  buff_x,
int  inc_x,
scomplex *  buff_v,
int  inc_v,
scomplex *  buff_w,
int  inc_w 
)
337{
338 scomplex* buff_0 = FLA_COMPLEX_PTR( FLA_ZERO );
339 int i;
340
341 bl1_csetv( m_A,
342 buff_0,
343 buff_w, inc_w );
344
345 for ( i = 0; i < n_A; ++i )
346 {
347 scomplex* a1 = buff_A + (i )*cs_A + (0 )*rs_A;
348 scomplex* nu1 = buff_v + (i )*inc_v;
349 scomplex* x = buff_x;
350 scomplex* chi1 = buff_x + (i )*inc_x;
351 scomplex* psi1 = buff_y + (i )*inc_y;
352 scomplex* u = buff_u;
353 scomplex* upsilon1 = buff_u + (i )*inc_u;
354 scomplex* w = buff_w;
355 scomplex* z = buff_z;
356 scomplex* alpha = buff_alpha;
357 scomplex temp1;
358 scomplex temp2;
359 scomplex conj_psi1;
360 scomplex conj_upsilon1;
361
362 /*------------------------------------------------------------*/
363
364 bl1_ccopyconj( psi1, &conj_psi1 );
365 bl1_cmult3( alpha, &conj_psi1, &temp1 );
366
367 bl1_ccopyconj( upsilon1, &conj_upsilon1 );
368 bl1_cmult3( alpha, &conj_upsilon1, &temp2 );
369
370 bl1_caxpyv( BLIS1_NO_CONJUGATE,
371 m_A,
372 &temp1,
373 u, inc_u,
374 a1, rs_A );
375/*
376 F77_caxpy( &m_A,
377 &temp1,
378 u, &inc_u,
379 a1, &rs_A );
380*/
381
382 bl1_caxpyv( BLIS1_NO_CONJUGATE,
383 m_A,
384 &temp2,
385 z, inc_z,
386 a1, rs_A );
387/*
388 F77_caxpy( &m_A,
389 &temp2,
390 z, &inc_z,
391 a1, &rs_A );
392*/
393
394 bl1_cdot( BLIS1_CONJUGATE,
395 m_A,
396 a1, rs_A,
397 x, inc_x,
398 nu1 );
399
400 bl1_caxpyv( BLIS1_NO_CONJUGATE,
401 m_A,
402 chi1,
403 a1, rs_A,
404 w, inc_w );
405/*
406 F77_caxpy( &m_A,
407 chi1,
408 a1, &rs_A,
409 w, &inc_w );
410*/
411
412 /*------------------------------------------------------------*/
413
414 }
415
416 return FLA_SUCCESS;
417}
FLA_ZERO
FLA_Obj FLA_ZERO
Definition FLA_Init.c:20
psi1
double *restrict psi1
Definition bl1_axmyv2.c:141
i
int i
Definition bl1_axmyv2.c:145
chi1
chi1
Definition bl1_axmyv2.c:366
temp2
double temp2
Definition bl1_axpyv2b.c:147
temp1
double temp1
Definition bl1_axpyv2b.c:146
upsilon1
upsilon1
Definition bl1_axpyv2bdotaxpy.c:225
bl1_caxpyv
void bl1_caxpyv(conj1_t conj, int n, scomplex *alpha, scomplex *x, int incx, scomplex *y, int incy)
Definition bl1_axpyv.c:29
bl1_cdot
void bl1_cdot(conj1_t conj, int n, scomplex *x, int incx, scomplex *y, int incy, scomplex *rho)
Definition bl1_dot.c:39
bl1_csetv
void bl1_csetv(int m, scomplex *sigma, scomplex *x, int incx)
Definition bl1_setv.c:52
BLIS1_CONJUGATE
@ BLIS1_CONJUGATE
Definition blis_type_defs.h:82
BLIS1_NO_CONJUGATE
@ BLIS1_NO_CONJUGATE
Definition blis_type_defs.h:81
scomplex
Definition blis_type_defs.h:133

References bl1_caxpyv(), bl1_cdot(), bl1_csetv(), BLIS1_CONJUGATE, BLIS1_NO_CONJUGATE, chi1, FLA_ZERO, i, psi1, temp1, temp2, and upsilon1.

Referenced by FLA_Fused_Gerc2_Ahx_Ax_opt_var1(), and FLA_Hess_UT_step_ofc_var3().

◆ FLA_Fused_Gerc2_Ahx_Ax_opd_var1()

FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opd_var1 ( int  m_A,
int  n_A,
double *  buff_alpha,
double *  buff_u,
int  inc_u,
double *  buff_y,
int  inc_y,
double *  buff_z,
int  inc_z,
double *  buff_A,
int  rs_A,
int  cs_A,
double *  buff_x,
int  inc_x,
double *  buff_v,
int  inc_v,
double *  buff_w,
int  inc_w 
)
258{
259 double zero = bl1_d0();
260
261 double* restrict alpha = buff_alpha;
262 double* restrict u = buff_u;
263 double* restrict z = buff_z;
264 double* restrict x = buff_x;
265 double* restrict w = buff_w;
266
267 double* restrict a1;
268 double* restrict nu1;
269 double* restrict chi1;
270 double* restrict psi1;
271 double* restrict upsilon1;
272
273 double alpha_psi1;
274 double alpha_upsilon1;
275
276 int n_run = n_A / 1;
277 //int n_left = n_A % 1;
278 int step_a1 = 1*cs_A;
279 int step_nu1 = 1*inc_v;
280 int step_chi1 = 1*inc_x;
281 int step_psi1 = 1*inc_y;
282 int step_upsilon1 = 1*inc_u;
283 int i;
284
285 bl1_dsetv( m_A,
286 &zero,
287 buff_w, inc_w );
288
289 a1 = buff_A;
290 nu1 = buff_v;
291 chi1 = buff_x;
292 psi1 = buff_y;
293 upsilon1 = buff_u;
294
295 for ( i = 0; i < n_run; ++i )
296 {
297 /*------------------------------------------------------------*/
298
299 bl1_dmult3( alpha, psi1, &alpha_psi1 );
300 bl1_dmult3( alpha, upsilon1, &alpha_upsilon1 );
301
302 bl1_daxpyv2bdotaxpy( m_A,
303 &alpha_psi1,
304 u, inc_u,
305 &alpha_upsilon1,
306 z, inc_z,
307 a1, rs_A,
308 x, inc_x,
309 chi1,
310 nu1,
311 w, inc_w );
312
313 /*------------------------------------------------------------*/
314
315 a1 += step_a1;
316 nu1 += step_nu1;
317 chi1 += step_chi1;
318 psi1 += step_psi1;
319 upsilon1 += step_upsilon1;
320 }
321
322 return FLA_SUCCESS;
323}
n_run
int n_run
Definition bl1_axmyv2.c:148
bl1_daxpyv2bdotaxpy
void bl1_daxpyv2bdotaxpy(int n, double *beta, double *u, int inc_u, double *gamma, double *z, int inc_z, double *a, int inc_a, double *x, int inc_x, double *kappa, double *rho, double *w, int inc_w)
Definition bl1_axpyv2bdotaxpy.c:36
bl1_dsetv
void bl1_dsetv(int m, double *sigma, double *x, int incx)
Definition bl1_setv.c:39
bl1_d0
double bl1_d0(void)
Definition bl1_constants.c:118

References bl1_d0(), bl1_daxpyv2bdotaxpy(), bl1_dsetv(), chi1, i, n_run, psi1, and upsilon1.

Referenced by FLA_Fused_Gerc2_Ahx_Ax_opt_var1(), and FLA_Hess_UT_step_ofd_var3().

◆ FLA_Fused_Gerc2_Ahx_Ax_ops_var1()

FLA_Error FLA_Fused_Gerc2_Ahx_Ax_ops_var1 ( int  m_A,
int  n_A,
float *  buff_alpha,
float *  buff_u,
int  inc_u,
float *  buff_y,
int  inc_y,
float *  buff_z,
int  inc_z,
float *  buff_A,
int  rs_A,
int  cs_A,
float *  buff_x,
int  inc_x,
float *  buff_v,
int  inc_v,
float *  buff_w,
int  inc_w 
)
160{
161 float* buff_0 = FLA_FLOAT_PTR( FLA_ZERO );
162 int i;
163
164 bl1_ssetv( m_A,
165 buff_0,
166 buff_w, inc_w );
167
168 for ( i = 0; i < n_A; ++i )
169 {
170 float* a1 = buff_A + (i )*cs_A + (0 )*rs_A;
171 float* nu1 = buff_v + (i )*inc_v;
172 float* x = buff_x;
173 float* chi1 = buff_x + (i )*inc_x;
174 float* psi1 = buff_y + (i )*inc_y;
175 float* u = buff_u;
176 float* upsilon1 = buff_u + (i )*inc_u;
177 float* w = buff_w;
178 float* z = buff_z;
179 float* alpha = buff_alpha;
180 float temp1;
181 float temp2;
182
183 /*------------------------------------------------------------*/
184
185 // bl1_scopyconj( psi1, &conj_psi1 );
186 // bl1_smult3( alpha, &conj_psi1, &temp1 );
187 temp1 = *alpha * *psi1;
188
189 // bl1_scopyconj( upsilon1, &conj_upsilon1 );
190 // bl1_smult3( alpha, &conj_upsilon1, &temp2 );
191 temp2 = *alpha * *upsilon1;
192
193 bl1_saxpyv( BLIS1_NO_CONJUGATE,
194 m_A,
195 &temp1,
196 u, inc_u,
197 a1, rs_A );
198/*
199 F77_saxpy( &m_A,
200 &temp1,
201 u, &inc_u,
202 a1, &rs_A );
203*/
204
205 bl1_saxpyv( BLIS1_NO_CONJUGATE,
206 m_A,
207 &temp2,
208 z, inc_z,
209 a1, rs_A );
210/*
211 F77_saxpy( &m_A,
212 &temp2,
213 z, &inc_z,
214 a1, &rs_A );
215*/
216
217 bl1_sdot( BLIS1_CONJUGATE,
218 m_A,
219 a1, rs_A,
220 x, inc_x,
221 nu1 );
222/*
223 *nu1 = F77_sdot( &m_A,
224 a1, &rs_A,
225 x, &inc_x );
226*/
227
228 bl1_saxpyv( BLIS1_NO_CONJUGATE,
229 m_A,
230 chi1,
231 a1, rs_A,
232 w, inc_w );
233/*
234 F77_saxpy( &m_A,
235 chi1,
236 a1, &rs_A,
237 w, &inc_w );
238*/
239 /*------------------------------------------------------------*/
240
241 }
242
243 return FLA_SUCCESS;
244}
bl1_saxpyv
void bl1_saxpyv(conj1_t conj, int n, float *alpha, float *x, int incx, float *y, int incy)
Definition bl1_axpyv.c:13
bl1_sdot
void bl1_sdot(conj1_t conj, int n, float *x, int incx, float *y, int incy, float *rho)
Definition bl1_dot.c:13
bl1_ssetv
void bl1_ssetv(int m, float *sigma, float *x, int incx)
Definition bl1_setv.c:26

References bl1_saxpyv(), bl1_sdot(), bl1_ssetv(), BLIS1_CONJUGATE, BLIS1_NO_CONJUGATE, chi1, FLA_ZERO, i, psi1, temp1, temp2, and upsilon1.

Referenced by FLA_Fused_Gerc2_Ahx_Ax_opt_var1(), and FLA_Hess_UT_step_ofs_var3().

◆ FLA_Fused_Gerc2_Ahx_Ax_opt_var1()

FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opt_var1 ( FLA_Obj  alpha,
FLA_Obj  u,
FLA_Obj  y,
FLA_Obj  z,
FLA_Obj  A,
FLA_Obj  x,
FLA_Obj  v,
FLA_Obj  w 
)
14{
15/*
16 Effective computation:
17 A = A + alpha * ( u * y' + z * u' );
18 v = A' * x;
19 w = A * x;
20*/
21 FLA_Datatype datatype;
22 int m_A, n_A;
23 int rs_A, cs_A;
24 int inc_u, inc_y, inc_z, inc_x, inc_v, inc_w;
25
26 datatype = FLA_Obj_datatype( A );
27
28 m_A = FLA_Obj_length( A );
29 n_A = FLA_Obj_width( A );
30
31 rs_A = FLA_Obj_row_stride( A );
32 cs_A = FLA_Obj_col_stride( A );
33
34 inc_u = FLA_Obj_vector_inc( u );
35 inc_y = FLA_Obj_vector_inc( y );
36 inc_z = FLA_Obj_vector_inc( z );
37 inc_x = FLA_Obj_vector_inc( x );
38 inc_v = FLA_Obj_vector_inc( v );
39 inc_w = FLA_Obj_vector_inc( w );
40
41
42 switch ( datatype )
43 {
44 case FLA_FLOAT:
45 {
46 float* buff_A = FLA_FLOAT_PTR( A );
47 float* buff_u = FLA_FLOAT_PTR( u );
48 float* buff_y = FLA_FLOAT_PTR( y );
49 float* buff_z = FLA_FLOAT_PTR( z );
50 float* buff_x = FLA_FLOAT_PTR( x );
51 float* buff_v = FLA_FLOAT_PTR( v );
52 float* buff_w = FLA_FLOAT_PTR( w );
53 float* buff_alpha = FLA_FLOAT_PTR( alpha );
54
55 FLA_Fused_Gerc2_Ahx_Ax_ops_var1( m_A,
56 n_A,
57 buff_alpha,
58 buff_u, inc_u,
59 buff_y, inc_y,
60 buff_z, inc_z,
61 buff_A, rs_A, cs_A,
62 buff_x, inc_x,
63 buff_v, inc_v,
64 buff_w, inc_w );
65
66 break;
67 }
68
69 case FLA_DOUBLE:
70 {
71 double* buff_A = FLA_DOUBLE_PTR( A );
72 double* buff_u = FLA_DOUBLE_PTR( u );
73 double* buff_y = FLA_DOUBLE_PTR( y );
74 double* buff_z = FLA_DOUBLE_PTR( z );
75 double* buff_x = FLA_DOUBLE_PTR( x );
76 double* buff_v = FLA_DOUBLE_PTR( v );
77 double* buff_w = FLA_DOUBLE_PTR( w );
78 double* buff_alpha = FLA_DOUBLE_PTR( alpha );
79
80 FLA_Fused_Gerc2_Ahx_Ax_opd_var1( m_A,
81 n_A,
82 buff_alpha,
83 buff_u, inc_u,
84 buff_y, inc_y,
85 buff_z, inc_z,
86 buff_A, rs_A, cs_A,
87 buff_x, inc_x,
88 buff_v, inc_v,
89 buff_w, inc_w );
90
91 break;
92 }
93
94 case FLA_COMPLEX:
95 {
96 scomplex* buff_A = FLA_COMPLEX_PTR( A );
97 scomplex* buff_u = FLA_COMPLEX_PTR( u );
98 scomplex* buff_y = FLA_COMPLEX_PTR( y );
99 scomplex* buff_z = FLA_COMPLEX_PTR( z );
100 scomplex* buff_x = FLA_COMPLEX_PTR( x );
101 scomplex* buff_v = FLA_COMPLEX_PTR( v );
102 scomplex* buff_w = FLA_COMPLEX_PTR( w );
103 scomplex* buff_alpha = FLA_COMPLEX_PTR( alpha );
104
105 FLA_Fused_Gerc2_Ahx_Ax_opc_var1( m_A,
106 n_A,
107 buff_alpha,
108 buff_u, inc_u,
109 buff_y, inc_y,
110 buff_z, inc_z,
111 buff_A, rs_A, cs_A,
112 buff_x, inc_x,
113 buff_v, inc_v,
114 buff_w, inc_w );
115
116 break;
117 }
118
119 case FLA_DOUBLE_COMPLEX:
120 {
121 dcomplex* buff_A = FLA_DOUBLE_COMPLEX_PTR( A );
122 dcomplex* buff_u = FLA_DOUBLE_COMPLEX_PTR( u );
123 dcomplex* buff_y = FLA_DOUBLE_COMPLEX_PTR( y );
124 dcomplex* buff_z = FLA_DOUBLE_COMPLEX_PTR( z );
125 dcomplex* buff_x = FLA_DOUBLE_COMPLEX_PTR( x );
126 dcomplex* buff_v = FLA_DOUBLE_COMPLEX_PTR( v );
127 dcomplex* buff_w = FLA_DOUBLE_COMPLEX_PTR( w );
128 dcomplex* buff_alpha = FLA_DOUBLE_COMPLEX_PTR( alpha );
129
130 FLA_Fused_Gerc2_Ahx_Ax_opz_var1( m_A,
131 n_A,
132 buff_alpha,
133 buff_u, inc_u,
134 buff_y, inc_y,
135 buff_z, inc_z,
136 buff_A, rs_A, cs_A,
137 buff_x, inc_x,
138 buff_v, inc_v,
139 buff_w, inc_w );
140
141 break;
142 }
143 }
144
145 return FLA_SUCCESS;
146}
FLA_Fused_Gerc2_Ahx_Ax_opz_var1
FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opz_var1(int m_A, int n_A, dcomplex *buff_alpha, dcomplex *buff_u, int inc_u, dcomplex *buff_y, int inc_y, dcomplex *buff_z, int inc_z, dcomplex *buff_A, int rs_A, int cs_A, dcomplex *buff_x, int inc_x, dcomplex *buff_v, int inc_v, dcomplex *buff_w, int inc_w)
Definition FLA_Fused_Gerc2_Ahx_Ax_opt_var1.c:421
FLA_Fused_Gerc2_Ahx_Ax_opd_var1
FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opd_var1(int m_A, int n_A, double *buff_alpha, double *buff_u, int inc_u, double *buff_y, int inc_y, double *buff_z, int inc_z, double *buff_A, int rs_A, int cs_A, double *buff_x, int inc_x, double *buff_v, int inc_v, double *buff_w, int inc_w)
Definition FLA_Fused_Gerc2_Ahx_Ax_opt_var1.c:248
FLA_Fused_Gerc2_Ahx_Ax_opc_var1
FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opc_var1(int m_A, int n_A, scomplex *buff_alpha, scomplex *buff_u, int inc_u, scomplex *buff_y, int inc_y, scomplex *buff_z, int inc_z, scomplex *buff_A, int rs_A, int cs_A, scomplex *buff_x, int inc_x, scomplex *buff_v, int inc_v, scomplex *buff_w, int inc_w)
Definition FLA_Fused_Gerc2_Ahx_Ax_opt_var1.c:327
FLA_Fused_Gerc2_Ahx_Ax_ops_var1
FLA_Error FLA_Fused_Gerc2_Ahx_Ax_ops_var1(int m_A, int n_A, float *buff_alpha, float *buff_u, int inc_u, float *buff_y, int inc_y, float *buff_z, int inc_z, float *buff_A, int rs_A, int cs_A, float *buff_x, int inc_x, float *buff_v, int inc_v, float *buff_w, int inc_w)
Definition FLA_Fused_Gerc2_Ahx_Ax_opt_var1.c:150
FLA_Obj_width
dim_t FLA_Obj_width(FLA_Obj obj)
Definition FLA_Query.c:123
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_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
dcomplex
Definition blis_type_defs.h:138

References FLA_Fused_Gerc2_Ahx_Ax_opc_var1(), FLA_Fused_Gerc2_Ahx_Ax_opd_var1(), FLA_Fused_Gerc2_Ahx_Ax_ops_var1(), FLA_Fused_Gerc2_Ahx_Ax_opz_var1(), FLA_Obj_col_stride(), FLA_Obj_datatype(), FLA_Obj_length(), FLA_Obj_row_stride(), FLA_Obj_vector_inc(), FLA_Obj_width(), and i.

◆ FLA_Fused_Gerc2_Ahx_Ax_opz_var1()

FLA_Error FLA_Fused_Gerc2_Ahx_Ax_opz_var1 ( int  m_A,
int  n_A,
dcomplex *  buff_alpha,
dcomplex *  buff_u,
int  inc_u,
dcomplex *  buff_y,
int  inc_y,
dcomplex *  buff_z,
int  inc_z,
dcomplex *  buff_A,
int  rs_A,
int  cs_A,
dcomplex *  buff_x,
int  inc_x,
dcomplex *  buff_v,
int  inc_v,
dcomplex *  buff_w,
int  inc_w 
)
431{
432 dcomplex zero = bl1_z0();
433
434 dcomplex* restrict alpha = buff_alpha;
435 dcomplex* restrict u = buff_u;
436 dcomplex* restrict z = buff_z;
437 dcomplex* restrict x = buff_x;
438 dcomplex* restrict w = buff_w;
439
440 dcomplex* restrict a1;
441 dcomplex* restrict nu1;
442 dcomplex* restrict chi1;
443 dcomplex* restrict psi1;
444 dcomplex* restrict upsilon1;
445
446 dcomplex temp1;
447 dcomplex temp2;
448 dcomplex conj_psi1;
449 dcomplex conj_upsilon1;
450
451 int n_run = n_A / 1;
452 //int n_left = n_A % 1;
453 int step_a1 = 1*cs_A;
454 int step_nu1 = 1*inc_v;
455 int step_chi1 = 1*inc_x;
456 int step_psi1 = 1*inc_y;
457 int step_upsilon1 = 1*inc_u;
458 int i;
459
460 bl1_zsetv( m_A,
461 &zero,
462 buff_w, inc_w );
463
464 a1 = buff_A;
465 nu1 = buff_v;
466 chi1 = buff_x;
467 psi1 = buff_y;
468 upsilon1 = buff_u;
469
470 for ( i = 0; i < n_run; ++i )
471 {
472 /*------------------------------------------------------------*/
473
474 bl1_zcopyconj( psi1, &conj_psi1 );
475 bl1_zmult3( alpha, &conj_psi1, &temp1 );
476
477 bl1_zcopyconj( upsilon1, &conj_upsilon1 );
478 bl1_zmult3( alpha, &conj_upsilon1, &temp2 );
479
480/*
481 bl1_zaxpyv2bdotaxpy( m_A,
482 &temp1,
483 u, inc_u,
484 &temp2,
485 z, inc_z,
486 a1, rs_A,
487 x, inc_x,
488 chi1,
489 nu1,
490 w, inc_w );
491*/
492
493 bl1_zaxpyv2b( m_A,
494 &temp1,
495 &temp2,
496 u, inc_u,
497 z, inc_z,
498 a1, rs_A );
499 bl1_zdotaxpy( m_A,
500 a1, rs_A,
501 x, inc_x,
502 chi1,
503 nu1,
504 w, inc_w );
505
506 /*------------------------------------------------------------*/
507
508 a1 += step_a1;
509 nu1 += step_nu1;
510 chi1 += step_chi1;
511 psi1 += step_psi1;
512 upsilon1 += step_upsilon1;
513 }
514
515 return FLA_SUCCESS;
516}
bl1_zaxpyv2b
void bl1_zaxpyv2b(int n, dcomplex *alpha1, dcomplex *alpha2, dcomplex *x1, int inc_x1, dcomplex *x2, int inc_x2, dcomplex *y, int inc_y)
Definition bl1_axpyv2b.c:210
bl1_zdotaxpy
void bl1_zdotaxpy(int n, dcomplex *a, int inc_a, dcomplex *x, int inc_x, dcomplex *kappa, dcomplex *rho, dcomplex *w, int inc_w)
Definition bl1_dotaxpy.c:258
bl1_z0
dcomplex bl1_z0(void)
Definition bl1_constants.c:133
bl1_zsetv
void bl1_zsetv(int m, dcomplex *sigma, dcomplex *x, int incx)
Definition bl1_setv.c:66

References bl1_z0(), bl1_zaxpyv2b(), bl1_zdotaxpy(), bl1_zsetv(), chi1, i, n_run, psi1, temp1, temp2, and upsilon1.

Referenced by FLA_Fused_Gerc2_Ahx_Ax_opt_var1(), and FLA_Hess_UT_step_ofz_var3().

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