Lines Matching refs:delta_u
109 double delta_u = 0.0; // dummy initialize, first used with * 0 in BesselJ() local
120 delta_u = g_bar_delta_u / g_bar; in BesselJ()
121 u = u + delta_u ; // u_k = u_{k-1} + delta_u_k in BesselJ()
133 g_bar_delta_u = - g * delta_u - m_bar * u; // alpha_k = 0.0 in BesselJ()
135 delta_u = g_bar_delta_u / g_bar; in BesselJ()
136 u = u + delta_u; in BesselJ()
142 g_bar_delta_u = f_bar - g * delta_u - m_bar * u; // alpha_k = 1.0 in BesselJ()
144 delta_u = g_bar_delta_u / g_bar; in BesselJ()
145 u = u + delta_u; in BesselJ()
154 g_bar_delta_u = - g * delta_u - m_bar * u; in BesselJ()
156 delta_u = g_bar_delta_u / g_bar; in BesselJ()
157 u = u + delta_u; in BesselJ()
160 bHasfound = (fabs(delta_u)<=fabs(u)*epsilon); in BesselJ()
358 double delta_u = g_bar_delta_u / g_bar; in Bessely0() local
377 g_bar_delta_u = f_bar * alpha - g * delta_u - m_bar * u; in Bessely0()
379 delta_u = g_bar_delta_u / g_bar; in Bessely0()
380 u = u+delta_u; in Bessely0()
383 bHasFound = (fabs(delta_u)<=fabs(u)*epsilon); in Bessely0()
414 double delta_u = g_bar_delta_u / g_bar; in Bessely1() local
415 u = u + delta_u; in Bessely1()
435 g_bar_delta_u = f_bar * alpha - g * delta_u - m_bar * u; in Bessely1()
437 delta_u = g_bar_delta_u / g_bar; in Bessely1()
438 u = u+delta_u; in Bessely1()
441 bHasFound = (fabs(delta_u)<=fabs(u)*epsilon); in Bessely1()