45 for ( ij_tl = ij_begin; ij_tl < m_A - 1; ++ij_tl )
47 double* d1 = buff_d + (ij_tl )*inc_d;
48 double* d2 = buff_d + (ij_tl+1)*inc_d;
49 double* e1 = buff_e + (ij_tl )*inc_e;
50 double abs_e1 = fabs( *e1 );
54 if ( abs_e1 != rzero )
56 if ( abs_e1 <= eps * sqrt( fabs( *d1 ) ) *
60 printf(
"FLA_Tevd_find_submatrix_opd: nudging non-zero subdiagonal element (e1) to zero.\n" );
61 printf(
" d[%3d] = %22.19e\n", ij_tl, *d1 );
62 printf(
" e[%3d] d[%3d] = %22.19e %22.19e\n", ij_tl, ij_tl+1, *e1, *d2 );
73 printf(
"FLA_Tevd_find_submatrix_opd: found non-zero subdiagonal element\n" );
74 printf(
" e[%3d] = %22.19e\n", ij_tl, *e1 );
83 if ( ij_tl == m_A - 1 )
86 printf(
"FLA_Tevd_find_submatrix_opd: no submatrices found.\n" );
95 for ( ij_br = ij_tl; ij_br < m_A - 1; ++ij_br )
97 double* d1 = buff_d + (ij_br )*inc_d;
98 double* d2 = buff_d + (ij_br+1)*inc_d;
99 double* e1 = buff_e + (ij_br )*inc_e;
100 double abs_e1 = fabs( *e1 );
104 if ( abs_e1 != rzero )
106 if ( abs_e1 <= eps * sqrt( fabs( *d1 ) ) *
107 sqrt( fabs( *d2 ) ) )
110 printf(
"FLA_Tevd_find_submatrix_opd: nudging non-zero subdiagonal element (e1) to zero.\n" );
111 printf(
" d[%3d] = %22.19e\n", ij_br, *d1 );
112 printf(
" e[%3d] d[%3d] = %22.19e %22.19e\n", ij_br, ij_br+1, *e1, *d2 );
123 printf(
"FLA_Tevd_find_submatrix_opd: found zero subdiagonal element\n" );
124 printf(
" e[%3d] = %22.19e\n", ij_br, *e1 );
double FLA_Mach_params_opd(FLA_Machval machval)
Definition: FLA_Mach_params.c:74
double bl1_d0(void)
Definition: bl1_constants.c:118