48 d1p = buff_d + (m_d-2)*inc_d;
49 e1p = buff_e + (m_d-2)*inc_e;
50 d2p = buff_d + (m_d-1)*inc_d;
77 for (
i = 0;
i < m_l; ++
i )
79 int* status = buff_ls + (
i )*inc_ls;
83 double* lambda1 = buff_l + (
i )*inc_l;
85 pshift_cand = *lambda1;
86 dist_cand = fabs( wilkshift - pshift_cand );
99 for (
i = 0;
i < m_l; ++
i )
101 double* lambda1 = buff_l + (
i )*inc_l;
102 int* status = buff_ls + (
i )*inc_ls;
103 double dist = fabs( wilkshift - *lambda1 );
105 if ( *status == 1 )
continue;
107 if ( dist < dist_cand )
110 pshift_cand = *lambda1;
FLA_Error FLA_Wilkshift_tridiag_opd(double delta1, double epsilon, double delta2, double *kappa)
Definition: FLA_Wilkshift_tridiag.c:155
int i
Definition: bl1_axmyv2.c:145