Lines Matching refs:pol
77 inline T sph_bessel_j_small_z_series(unsigned v, T x, const Policy& pol) in sph_bessel_j_small_z_series() argument
88 …heck_series_iterations<T>("boost::math::sph_bessel_j_small_z_series<%1%>(%1%,%1%)", max_iter, pol); in sph_bessel_j_small_z_series()
93 T cyl_bessel_j_imp(T v, T x, const bessel_no_int_tag& t, const Policy& pol)
102 T r = cyl_bessel_j_imp(v, T(-x), t, pol);
103 if(iround(v, pol) & 1)
110 "Got x = %1%, but we need x >= 0", x, pol);
114 bessel_jy(v, x, &j, &y, need_j, pol);
119 inline T cyl_bessel_j_imp(T v, T x, const bessel_maybe_int_tag&, const Policy& pol) in cyl_bessel_j_imp() argument
122 int ival = detail::iconv(v, pol); in cyl_bessel_j_imp()
127 return bessel_jn(ival, x, pol); in cyl_bessel_j_imp()
129 return cyl_bessel_j_imp(v, x, bessel_no_int_tag(), pol); in cyl_bessel_j_imp()
133 inline T cyl_bessel_j_imp(int v, T x, const bessel_int_tag&, const Policy& pol) in cyl_bessel_j_imp() argument
136 return bessel_jn(v, x, pol); in cyl_bessel_j_imp()
140 inline T sph_bessel_j_imp(unsigned n, T x, const Policy& pol) in sph_bessel_j_imp() argument
146 "Got x = %1%, but function requires x > 0.", x, pol); in sph_bessel_j_imp()
151 return boost::math::sinc_pi(x, pol); in sph_bessel_j_imp()
162 return sph_bessel_j_small_z_series(n, x, pol); in sph_bessel_j_imp()
167 * cyl_bessel_j_imp(T(T(n)+T(0.5f)), x, bessel_no_int_tag(), pol); in sph_bessel_j_imp()
171 T cyl_bessel_i_imp(T v, T x, const Policy& pol)
185 T r = cyl_bessel_i_imp(v, T(-x), pol);
186 if(iround(v, pol) & 1)
193 "Got x = %1%, but we need x >= 0", x, pol);
221 return bessel_i_small_z_series(v, x, pol);
223 bessel_ik(v, x, &I, &K, need_i, pol);
228 inline T cyl_bessel_k_imp(T v, T x, const bessel_no_int_tag& /* t */, const Policy& pol) in cyl_bessel_k_imp() argument
236 "Got x = %1%, but we need x > 0", x, pol); in cyl_bessel_k_imp()
240 return (v == 0) ? policies::raise_overflow_error<T>(function, 0, pol) in cyl_bessel_k_imp()
243 "Got x = %1%, but we need x > 0", x, pol); in cyl_bessel_k_imp()
246 bessel_ik(v, x, &I, &K, need_k, pol); in cyl_bessel_k_imp()
251 inline T cyl_bessel_k_imp(T v, T x, const bessel_maybe_int_tag&, const Policy& pol) in cyl_bessel_k_imp() argument
256 return bessel_kn(itrunc(v), x, pol); in cyl_bessel_k_imp()
258 return cyl_bessel_k_imp(v, x, bessel_no_int_tag(), pol); in cyl_bessel_k_imp()
262 inline T cyl_bessel_k_imp(int v, T x, const bessel_int_tag&, const Policy& pol) in cyl_bessel_k_imp() argument
264 return bessel_kn(v, x, pol); in cyl_bessel_k_imp()
268 inline T cyl_neumann_imp(T v, T x, const bessel_no_int_tag&, const Policy& pol) in cyl_neumann_imp() argument
278 policies::raise_overflow_error<T>(function, 0, pol) in cyl_neumann_imp()
281 "Got x = %1%, but result is complex for x <= 0", x, pol); in cyl_neumann_imp()
284 bessel_jy(v, x, &j, &y, need_y, pol); in cyl_neumann_imp()
291 return -policies::raise_overflow_error<T>(function, 0, pol); in cyl_neumann_imp()
296 inline T cyl_neumann_imp(T v, T x, const bessel_maybe_int_tag&, const Policy& pol) in cyl_neumann_imp() argument
305 T r = bessel_yn(itrunc(v, pol), x, pol); in cyl_neumann_imp()
309 T r = cyl_neumann_imp<T>(v, x, bessel_no_int_tag(), pol); in cyl_neumann_imp()
315 inline T cyl_neumann_imp(int v, T x, const bessel_int_tag&, const Policy& pol) in cyl_neumann_imp() argument
317 return bessel_yn(v, x, pol); in cyl_neumann_imp()
321 inline T sph_neumann_imp(unsigned v, T x, const Policy& pol) in sph_neumann_imp() argument
332 "Got x = %1%, but function requires x > 0.", x, pol); in sph_neumann_imp()
335 return -policies::raise_overflow_error<T>(function, 0, pol); in sph_neumann_imp()
337 T result = cyl_neumann_imp(T(T(v)+0.5f), x, bessel_no_int_tag(), pol); in sph_neumann_imp()
341 return -policies::raise_overflow_error<T>(function, 0, pol); in sph_neumann_imp()
347 inline T cyl_bessel_j_zero_imp(T v, int m, const Policy& pol) in cyl_bessel_j_zero_imp() argument
358 …icies::raise_domain_error<T>(function, "Order argument is %1%, but must be finite >= 0 !", v, pol); in cyl_bessel_j_zero_imp()
365 …T>(function, "Requested the %1%'th zero, but the rank must be positive !", static_cast<T>(m), pol); in cyl_bessel_j_zero_imp()
381 …>(function, "Requested the %1%'th zero of J0, but the rank must be > 0 !", static_cast<T>(m), pol); in cyl_bessel_j_zero_imp()
389 …h zero of Jv for negative, non-integer order, but the rank must be > 0 !", static_cast<T>(m), pol); in cyl_bessel_j_zero_imp()
399 …ssel_zero::cyl_bessel_j_zero_detail::initial_guess<T, Policy>((order_is_integer ? vv : v), m, pol); in cyl_bessel_j_zero_imp()
409 …etail::function_object_jv_and_jv_prime<T, Policy>((order_is_integer ? vv : v), order_is_zero, pol), in cyl_bessel_j_zero_imp()
426 inline T cyl_neumann_zero_imp(T v, int m, const Policy& pol) in cyl_neumann_zero_imp() argument
435 …icies::raise_domain_error<T>(function, "Order argument is %1%, but must be finite >= 0 !", v, pol); in cyl_neumann_zero_imp()
441 …T>(function, "Requested the %1%'th zero, but the rank must be positive !", static_cast<T>(m), pol); in cyl_neumann_zero_imp()
454 return boost::math::detail::cyl_neumann_zero_imp(vv, m, pol); in cyl_neumann_zero_imp()
467 …o of Yv for negative, non-half-integer order, but the rank must be > 0 !", static_cast<T>(m), pol); in cyl_neumann_zero_imp()
473 return boost::math::detail::cyl_bessel_j_zero_imp(vv, m, pol); in cyl_neumann_zero_imp()
478 …t = boost::math::detail::bessel_zero::cyl_neumann_zero_detail::initial_guess<T, Policy>(v, m, pol); in cyl_neumann_zero_imp()
488 …::detail::bessel_zero::cyl_neumann_zero_detail::function_object_yv_and_yv_prime<T, Policy>(v, pol), in cyl_neumann_zero_imp()
672 const Policy& pol) in cyl_bessel_j_zero() argument
681 *out_it = boost::math::cyl_bessel_j_zero(v, start_index + i, pol); in cyl_bessel_j_zero()
733 const Policy& pol) in cyl_neumann_zero() argument
742 *out_it = boost::math::cyl_neumann_zero(v, start_index + i, pol); in cyl_neumann_zero()