51 for ( i = m_A - 1; i > 1; --
i )
53 dcomplex* G1 = buff_G + (k_iter)*cs_G;
55 int k_left = n_G - k_iter;
133 if ( r_val == FLA_FAILURE && k_iter == n_G )
136 printf(
"FLA_Tevd_iteracc_v_opd_var1: failed to converge (m_A11 = %d) after %2d iters k_total=%d/%d\n", i, k, k_iter, n_G );
138 *n_iter_perf = k_iter;
150 printf(
"FLA_Tevd_iteracc_v_opd_var3: found eig %22.15e in col %3d (n=%d) after %2d iters k_total=%d/%d\n", buff_d[ r_val*inc_d ], ijTL+r_val, m_ATL, k, k_iter, n_G );
152 printf(
"FLA_Tevd_iteracc_v_opd_var3: split occurred in col %3d (n=%d) after %2d iters k_total=%d/%d\n", ijTL+r_val, m_ATL, k, k_iter, n_G );
159 *n_iter_perf = k_iter;
168 int m_TLr = r_val + 1;
169 int m_BRr = m_ATL - m_TLr;
172 int n_Gr = n_G - k_iter;
173 double* dTL = buff_d + (0 )*inc_d;
174 double* eTL = buff_e + (0 )*inc_e;
175 double* puTL = buff_pu+ (0 )*inc_pu;
176 dcomplex* GT = buff_G + (0 )*rs_G + (k_iter)*cs_G;
177 double* dBR = buff_d + (ijBRr)*inc_d;
178 double* eBR = buff_e + (ijBRr)*inc_e;
179 double* puBR = buff_pu+ (ijBRr)*inc_pu;
180 dcomplex* GB = buff_G + (ijBRr)*rs_G + (k_iter)*cs_G;
188 printf(
"FLA_Tevd_iteracc_v_opd_var3: Internal deflation in col %d\n", ijTL+r_val );
189 printf(
"FLA_Tevd_iteracc_v_opd_var3: alpha11 = %23.19e\n", buff_d[r_val*inc_d] );
190 printf(
"FLA_Tevd_iteracc_v_opd_var3: alpha21 alpha22 = %23.19e %23.19e\n", buff_e[r_val*inc_e], buff_d[(r_val+1)*inc_d] );
193 printf(
"FLA_Tevd_iteracc_v_opd_var3: recursing: m_TLr m_BRr: %d %d\n", m_TLr, m_BRr );
194 printf(
"FLA_Tevd_iteracc_v_opd_var3: ijTLr ijBRr: %d %d\n", ijTLr, ijBRr );
195 printf(
"FLA_Tevd_iteracc_v_opd_var3: GB(0,0) i,j: %d %d\n", ijTL + m_TLr+1, k_iter );
220 *n_iter_perf = k_iter + max( n_iter_perfTL, n_iter_perfBR );
223 printf(
"FLA_Tevd_iteracc_v_opd_var3: num deflations: %d = (prev:%d, TL:%d, BR:%d)\n", n_deflations + n_deflationsTL + n_deflationsBR, n_deflations, n_deflationsTL, n_deflationsBR );
224 printf(
"FLA_Tevd_iteracc_v_opd_var3: num iterations: %d = (prev:%d, TL:%d, BR:%d)\n", *n_iter_perf, k_iter, n_iter_perfTL, n_iter_perfBR );
226 return n_deflations + n_deflationsTL + n_deflationsBR;
235 dcomplex* g1 = buff_G + (k_iter)*cs_G;
237 double* alpha11 = buff_d + (0 )*inc_d;
238 double* alpha21 = buff_e + (0 )*inc_e;
239 double* alpha22 = buff_d + (1 )*inc_d;
273 printf(
"FLA_Tevd_iteracc_v_opd_var3: Hevv eig %22.15e in col %3d (n=%d) after %2d iters k_total=%d/%d\n", buff_d[ 1*inc_d ], ijTL+1, 2, 1, k_iter, n_G );
274 printf(
"FLA_Tevd_iteracc_v_opd_var3: Hevv eig %22.15e in col %3d (n=%d) after %2d iters k_total=%d/%d\n", buff_d[ 0*inc_d ], ijTL+0, 2, 0, k_iter, n_G );
279 *n_iter_perf = k_iter;
FLA_Error FLA_Hevv_2x2_opd(double *alpha11, double *alpha21, double *alpha22, double *lambda1, double *lambda2, double *gamma1, double *sigma1)
Definition: FLA_Hevv_2x2.c:249
FLA_Error FLA_Tevd_iteracc_v_opd_var3(int m_A, int m_U, int n_G, int ijTL, double *buff_d, int inc_d, double *buff_e, int inc_e, double *buff_l, int inc_l, int *buff_ls, int inc_ls, double *buff_pu, int inc_pu, dcomplex *buff_G, int rs_G, int cs_G, int *n_iter_perf)
Definition: FLA_Tevd_iteracc_v_opt_var3.c:30
double imag
Definition: blis_type_defs.h:139
int FLA_Error
Definition: FLA_type_defs.h:47
double real
Definition: blis_type_defs.h:139
FLA_Error FLA_Tevd_eigval_v_opd_var3(int m_A, int m_U, int n_G, dcomplex *buff_G, int rs_G, int cs_G, double *buff_d, int inc_d, double *buff_e, int inc_e, double *buff_l, int inc_l, int *buff_ls, int inc_ls, double *buff_pu, int inc_pu, int *n_iter)
Definition: FLA_Tevd_eigval_v_opt_var3.c:30
int i
Definition: bl1_axmyv2.c:145
Definition: blis_type_defs.h:137