3 <meta http-equiv="Content-Type" content="text/html; charset=US-ASCII">
4 <title>Derivatives of the Bessel Functions</title>
5 <link rel="stylesheet" href="../../math.css" type="text/css">
6 <meta name="generator" content="DocBook XSL Stylesheets V1.79.1">
7 <link rel="home" href="../../index.html" title="Math Toolkit 2.11.0">
8 <link rel="up" href="../bessel.html" title="Bessel Functions">
9 <link rel="prev" href="sph_bessel.html" title="Spherical Bessel Functions of the First and Second Kinds">
10 <link rel="next" href="../hankel.html" title="Hankel Functions">
12 <body bgcolor="white" text="black" link="#0000FF" vlink="#840084" alink="#0000FF">
13 <table cellpadding="2" width="100%"><tr>
14 <td valign="top"><img alt="Boost C++ Libraries" width="277" height="86" src="../../../../../../boost.png"></td>
15 <td align="center"><a href="../../../../../../index.html">Home</a></td>
16 <td align="center"><a href="../../../../../../libs/libraries.htm">Libraries</a></td>
17 <td align="center"><a href="http://www.boost.org/users/people.html">People</a></td>
18 <td align="center"><a href="http://www.boost.org/users/faq.html">FAQ</a></td>
19 <td align="center"><a href="../../../../../../more/index.htm">More</a></td>
22 <div class="spirit-nav">
23 <a accesskey="p" href="sph_bessel.html"><img src="../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../bessel.html"><img src="../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../../index.html"><img src="../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="../hankel.html"><img src="../../../../../../doc/src/images/next.png" alt="Next"></a>
26 <div class="titlepage"><div><div><h3 class="title">
27 <a name="math_toolkit.bessel.bessel_derivatives"></a><a class="link" href="bessel_derivatives.html" title="Derivatives of the Bessel Functions">Derivatives of
28 the Bessel Functions</a>
29 </h3></div></div></div>
31 <a name="math_toolkit.bessel.bessel_derivatives.h0"></a>
32 <span class="phrase"><a name="math_toolkit.bessel.bessel_derivatives.synopsis"></a></span><a class="link" href="bessel_derivatives.html#math_toolkit.bessel.bessel_derivatives.synopsis">Synopsis</a>
35 <code class="computeroutput"><span class="preprocessor">#include</span> <span class="special"><</span><span class="identifier">boost</span><span class="special">/</span><span class="identifier">math</span><span class="special">/</span><span class="identifier">special_functions</span><span class="special">/</span><span class="identifier">bessel_prime</span><span class="special">.</span><span class="identifier">hpp</span><span class="special">></span></code>
37 <pre class="programlisting"><span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">></span>
38 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cyl_bessel_j_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">);</span>
40 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">,</span> <span class="keyword">class</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">></span>
41 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cyl_bessel_j_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">,</span> <span class="keyword">const</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">&);</span>
43 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">></span>
44 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cyl_neumann_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">);</span>
46 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">,</span> <span class="keyword">class</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">></span>
47 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cyl_neumann_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">,</span> <span class="keyword">const</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">&);</span>
49 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">></span>
50 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cyl_bessel_i_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">);</span>
52 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">,</span> <span class="keyword">class</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">></span>
53 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cyl_bessel_i_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">,</span> <span class="keyword">const</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">&);</span>
55 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">></span>
56 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cyl_bessel_k_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">);</span>
58 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">,</span> <span class="keyword">class</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">></span>
59 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">cyl_bessel_k_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">,</span> <span class="keyword">const</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">&);</span>
61 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">></span>
62 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">sph_bessel_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">);</span>
64 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">,</span> <span class="keyword">class</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">></span>
65 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">sph_bessel_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">,</span> <span class="keyword">const</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">&);</span>
67 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">></span>
68 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">sph_neumann_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">);</span>
70 <span class="keyword">template</span> <span class="special"><</span><span class="keyword">class</span> <span class="identifier">T1</span><span class="special">,</span> <span class="keyword">class</span> <span class="identifier">T2</span><span class="special">,</span> <span class="keyword">class</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">></span>
71 <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>calculated-result-type</em></span></a> <span class="identifier">sph_neumann_prime</span><span class="special">(</span><span class="identifier">T1</span> <span class="identifier">v</span><span class="special">,</span> <span class="identifier">T2</span> <span class="identifier">x</span><span class="special">,</span> <span class="keyword">const</span> <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a><span class="special">&);</span>
74 <a name="math_toolkit.bessel.bessel_derivatives.h1"></a>
75 <span class="phrase"><a name="math_toolkit.bessel.bessel_derivatives.description"></a></span><a class="link" href="bessel_derivatives.html#math_toolkit.bessel.bessel_derivatives.description">Description</a>
78 These functions return the first derivative with respect to <span class="emphasis"><em>x</em></span>
79 of the corresponding Bessel function.
82 The return type of these functions is computed using the <a class="link" href="../result_type.html" title="Calculation of the Type of the Result"><span class="emphasis"><em>result
83 type calculation rules</em></span></a> when T1 and T2 are different types.
84 The functions are also optimised for the relatively common case that T1 is
88 The final <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">Policy</a> argument is optional and can
89 be used to control the behaviour of the function: how it handles errors,
90 what level of precision to use etc. Refer to the <a class="link" href="../../policy.html" title="Chapter 20. Policies: Controlling Precision, Error Handling etc">policy
91 documentation for more details</a>.
94 The functions return the result of <a class="link" href="../error_handling.html#math_toolkit.error_handling.domain_error">domain_error</a>
95 whenever the result is undefined or complex.
98 <a name="math_toolkit.bessel.bessel_derivatives.h2"></a>
99 <span class="phrase"><a name="math_toolkit.bessel.bessel_derivatives.testing"></a></span><a class="link" href="bessel_derivatives.html#math_toolkit.bessel.bessel_derivatives.testing">Testing</a>
102 There are two sets of test values: spot values calculated using <a href="http://www.wolframalpha.com/" target="_top">wolframalpha.com</a>,
103 and a much larger set of tests computed using a relation to the underlying
104 Bessel functions that the implementation does not use.
107 <a name="math_toolkit.bessel.bessel_derivatives.h3"></a>
108 <span class="phrase"><a name="math_toolkit.bessel.bessel_derivatives.accuracy"></a></span><a class="link" href="bessel_derivatives.html#math_toolkit.bessel.bessel_derivatives.accuracy">Accuracy</a>
111 The accuracy of these functions is broadly similar to the underlying Bessel
115 <a name="math_toolkit.bessel.bessel_derivatives.table_cyl_bessel_i_prime_integer_orders_"></a><p class="title"><b>Table 8.50. Error rates for cyl_bessel_i_prime (integer orders)</b></p>
116 <div class="table-contents"><table class="table" summary="Error rates for cyl_bessel_i_prime (integer orders)">
129 GNU C++ version 7.1.0<br> linux<br> double
134 GNU C++ version 7.1.0<br> linux<br> long double
139 Sun compiler version 0x5150<br> Sun Solaris<br> long double
144 Microsoft Visual C++ version 14.1<br> Win32<br> double
152 Bessel I'0: Mathworld Data (Integer Version)
157 <span class="blue">Max = 0ε (Mean = 0ε)</span>
162 <span class="blue">Max = 0.82ε (Mean = 0.259ε)</span>
167 <span class="blue">Max = 0.82ε (Mean = 0.259ε)</span>
172 <span class="blue">Max = 0.82ε (Mean = 0.354ε)</span>
179 Bessel I'1: Mathworld Data (Integer Version)
184 <span class="blue">Max = 0ε (Mean = 0ε)</span>
189 <span class="blue">Max = 1.97ε (Mean = 0.757ε)</span>
194 <span class="blue">Max = 1.97ε (Mean = 0.757ε)</span>
199 <span class="blue">Max = 1.36ε (Mean = 0.782ε)</span>
206 Bessel I'n: Mathworld Data (Integer Version)
211 <span class="blue">Max = 0ε (Mean = 0ε)</span>
216 <span class="blue">Max = 2.31ε (Mean = 1.41ε)</span>
221 <span class="blue">Max = 701ε (Mean = 212ε)</span>
226 <span class="blue">Max = 3.61ε (Mean = 1.22ε)</span>
233 <br class="table-break"><div class="table">
234 <a name="math_toolkit.bessel.bessel_derivatives.table_cyl_bessel_i_prime"></a><p class="title"><b>Table 8.51. Error rates for cyl_bessel_i_prime</b></p>
235 <div class="table-contents"><table class="table" summary="Error rates for cyl_bessel_i_prime">
248 GNU C++ version 7.1.0<br> linux<br> double
253 GNU C++ version 7.1.0<br> linux<br> long double
258 Sun compiler version 0x5150<br> Sun Solaris<br> long double
263 Microsoft Visual C++ version 14.1<br> Win32<br> double
271 Bessel I'0: Mathworld Data
276 <span class="blue">Max = 0ε (Mean = 0ε)</span>
281 <span class="blue">Max = 0.82ε (Mean = 0.259ε)</span>
286 <span class="blue">Max = 0.82ε (Mean = 0.259ε)</span>
291 <span class="blue">Max = 0.82ε (Mean = 0.354ε)</span>
298 Bessel I'1: Mathworld Data
303 <span class="blue">Max = 0ε (Mean = 0ε)</span>
308 <span class="blue">Max = 1.97ε (Mean = 0.757ε)</span>
313 <span class="blue">Max = 1.97ε (Mean = 0.757ε)</span>
318 <span class="blue">Max = 1.36ε (Mean = 0.782ε)</span>
325 Bessel I'n: Mathworld Data
330 <span class="blue">Max = 0ε (Mean = 0ε)</span>
335 <span class="blue">Max = 2.31ε (Mean = 1.41ε)</span>
340 <span class="blue">Max = 701ε (Mean = 212ε)</span>
345 <span class="blue">Max = 3.61ε (Mean = 1.22ε)</span>
352 Bessel I'v: Mathworld Data
357 <span class="blue">Max = 1.62ε (Mean = 0.512ε)</span>
362 <span class="blue">Max = 2.89e+03ε (Mean = 914ε)</span>
367 <span class="blue">Max = 2.89e+03ε (Mean = 914ε)</span>
372 <span class="blue">Max = 3.76e+03ε (Mean = 1.19e+03ε)</span>
379 Bessel I'n: Random Data
384 <span class="blue">Max = 0ε (Mean = 0ε)</span>
389 <span class="blue">Max = 3.95ε (Mean = 1.06ε)</span>
394 <span class="blue">Max = 195ε (Mean = 37.1ε)</span>
399 <span class="blue">Max = 9.85ε (Mean = 1.82ε)</span>
406 Bessel I'v: Random Data
411 <span class="blue">Max = 0ε (Mean = 0ε)</span>
416 <span class="blue">Max = 14.1ε (Mean = 2.93ε)</span>
421 <span class="blue">Max = 336ε (Mean = 68.5ε)</span>
426 <span class="blue">Max = 14ε (Mean = 2.5ε)</span>
433 Bessel I'v: Mathworld Data (large values)
438 <span class="blue">Max = 0ε (Mean = 0ε)</span>
443 <span class="blue">Max = 42.6ε (Mean = 20.2ε)</span>
448 <span class="blue">Max = 42.6ε (Mean = 20.2ε)</span>
453 <span class="blue">Max = 59.5ε (Mean = 26.6ε)</span>
460 <br class="table-break"><div class="table">
461 <a name="math_toolkit.bessel.bessel_derivatives.table_cyl_bessel_j_prime_integer_orders_"></a><p class="title"><b>Table 8.52. Error rates for cyl_bessel_j_prime (integer orders)</b></p>
462 <div class="table-contents"><table class="table" summary="Error rates for cyl_bessel_j_prime (integer orders)">
475 GNU C++ version 7.1.0<br> linux<br> double
480 GNU C++ version 7.1.0<br> linux<br> long double
485 Sun compiler version 0x5150<br> Sun Solaris<br> long double
490 Microsoft Visual C++ version 14.1<br> Win32<br> double
498 Bessel J0': Mathworld Data (Integer Version)
503 <span class="blue">Max = 0ε (Mean = 0ε)</span>
508 <span class="blue">Max = 18.9ε (Mean = 6.82ε)</span>
513 <span class="blue">Max = 18.9ε (Mean = 6.72ε)</span>
518 <span class="blue">Max = 6.62ε (Mean = 2.55ε)</span>
525 Bessel J0': Mathworld Data (Tricky cases) (Integer Version)
530 <span class="blue">Max = 0ε (Mean = 0ε)</span>
535 <span class="blue">Max = 7.44ε (Mean = 3.34ε)</span>
540 <span class="blue">Max = 7.44ε (Mean = 3.31ε)</span>
545 <span class="blue">Max = 3.67ε (Mean = 1.74ε)</span>
552 Bessel J1': Mathworld Data (Integer Version)
557 <span class="blue">Max = 0ε (Mean = 0ε)</span>
562 <span class="blue">Max = 7.9ε (Mean = 3.37ε)</span>
567 <span class="blue">Max = 7.9ε (Mean = 3.37ε)</span>
572 <span class="blue">Max = 0.999ε (Mean = 0.627ε)</span>
579 Bessel J1': Mathworld Data (tricky cases) (Integer Version)
584 <span class="blue">Max = 287ε (Mean = 129ε)</span>
589 <span class="blue">Max = 5.88e+05ε (Mean = 2.63e+05ε)</span>
594 <span class="blue">Max = 5.88e+05ε (Mean = 2.63e+05ε)</span>
599 <span class="blue">Max = 288ε (Mean = 129ε)</span>
606 Bessel JN': Mathworld Data (Integer Version)
611 <span class="blue">Max = 0.527ε (Mean = 0.128ε)</span>
616 <span class="blue">Max = 1.29e+03ε (Mean = 312ε)</span>
621 <span class="blue">Max = 1.29e+03ε (Mean = 355ε)</span>
626 <span class="blue">Max = 14ε (Mean = 6.13ε)</span>
633 <br class="table-break"><div class="table">
634 <a name="math_toolkit.bessel.bessel_derivatives.table_cyl_bessel_j_prime"></a><p class="title"><b>Table 8.53. Error rates for cyl_bessel_j_prime</b></p>
635 <div class="table-contents"><table class="table" summary="Error rates for cyl_bessel_j_prime">
648 GNU C++ version 7.1.0<br> linux<br> double
653 GNU C++ version 7.1.0<br> linux<br> long double
658 Sun compiler version 0x5150<br> Sun Solaris<br> long double
663 Microsoft Visual C++ version 14.1<br> Win32<br> double
671 Bessel J0': Mathworld Data
676 <span class="blue">Max = 0ε (Mean = 0ε)</span>
681 <span class="blue">Max = 18.9ε (Mean = 6.82ε)</span>
686 <span class="blue">Max = 18.9ε (Mean = 6.72ε)</span>
691 <span class="blue">Max = 6.62ε (Mean = 2.55ε)</span>
698 Bessel J0': Mathworld Data (Tricky cases)
703 <span class="blue">Max = 0ε (Mean = 0ε)</span>
708 <span class="blue">Max = 7.44ε (Mean = 3.34ε)</span>
713 <span class="blue">Max = 7.44ε (Mean = 3.31ε)</span>
718 <span class="blue">Max = 3.67ε (Mean = 1.74ε)</span>
725 Bessel J1': Mathworld Data
730 <span class="blue">Max = 0ε (Mean = 0ε)</span>
735 <span class="blue">Max = 7.9ε (Mean = 3.37ε)</span>
740 <span class="blue">Max = 7.9ε (Mean = 3.37ε)</span>
745 <span class="blue">Max = 0.999ε (Mean = 0.627ε)</span>
752 Bessel J1': Mathworld Data (tricky cases)
757 <span class="blue">Max = 287ε (Mean = 129ε)</span>
762 <span class="blue">Max = 5.88e+05ε (Mean = 2.63e+05ε)</span>
767 <span class="blue">Max = 5.88e+05ε (Mean = 2.63e+05ε)</span>
772 <span class="blue">Max = 288ε (Mean = 129ε)</span>
779 Bessel JN': Mathworld Data
784 <span class="blue">Max = 0.527ε (Mean = 0.128ε)</span>
789 <span class="blue">Max = 1.29e+03ε (Mean = 312ε)</span>
794 <span class="blue">Max = 1.29e+03ε (Mean = 355ε)</span>
799 <span class="blue">Max = 14ε (Mean = 6.13ε)</span>
806 Bessel J': Mathworld Data
811 <span class="blue">Max = 21.5ε (Mean = 4.7ε)</span>
816 <span class="blue">Max = 42.5ε (Mean = 9.31ε)</span>
821 <span class="blue">Max = 42.5ε (Mean = 9.32ε)</span>
826 <span class="blue">Max = 23.7ε (Mean = 8ε)</span>
833 Bessel J': Mathworld Data (large values)
838 <span class="blue">Max = 0ε (Mean = 0ε)</span>
843 <span class="blue">Max = 989ε (Mean = 495ε)</span>
848 <span class="blue">Max = 989ε (Mean = 495ε)</span>
853 <span class="blue">Max = 2.9ε (Mean = 1.61ε)</span>
860 Bessel JN': Random Data
865 <span class="blue">Max = 0.593ε (Mean = 0.0396ε)</span>
870 <span class="blue">Max = 11.3ε (Mean = 1.85ε)</span>
875 <span class="blue">Max = 79.4ε (Mean = 16.2ε)</span>
880 <span class="blue">Max = 6.34ε (Mean = 0.999ε)</span>
887 Bessel J': Random Data
892 <span class="blue">Max = 0.885ε (Mean = 0.033ε)</span>
897 <span class="blue">Max = 139ε (Mean = 6.47ε)</span>
902 <span class="blue">Max = 279ε (Mean = 27.2ε)</span>
907 <span class="blue">Max = 176ε (Mean = 9.75ε)</span>
914 Bessel J': Random Data (Tricky large values)
919 <span class="blue">Max = 0ε (Mean = 0ε)</span>
924 <span class="blue">Max = 474ε (Mean = 62.2ε)</span>
929 <span class="blue">Max = 474ε (Mean = 64.5ε)</span>
934 <span class="blue">Max = 379ε (Mean = 45.4ε)</span>
941 <br class="table-break"><div class="table">
942 <a name="math_toolkit.bessel.bessel_derivatives.table_cyl_bessel_k_prime_integer_orders_"></a><p class="title"><b>Table 8.54. Error rates for cyl_bessel_k_prime (integer orders)</b></p>
943 <div class="table-contents"><table class="table" summary="Error rates for cyl_bessel_k_prime (integer orders)">
956 GNU C++ version 7.1.0<br> linux<br> double
961 GNU C++ version 7.1.0<br> linux<br> long double
966 Sun compiler version 0x5150<br> Sun Solaris<br> long double
971 Microsoft Visual C++ version 14.1<br> Win32<br> double
979 Bessel K'0: Mathworld Data (Integer Version)
984 <span class="blue">Max = 0ε (Mean = 0ε)</span>
989 <span class="blue">Max = 0.786ε (Mean = 0.329ε)</span>
994 <span class="blue">Max = 0.786ε (Mean = 0.329ε)</span>
999 <span class="blue">Max = 0.786ε (Mean = 0.39ε)</span>
1006 Bessel K'1: Mathworld Data (Integer Version)
1011 <span class="blue">Max = 0ε (Mean = 0ε)</span>
1016 <span class="blue">Max = 0.736ε (Mean = 0.389ε)</span>
1021 <span class="blue">Max = 0.736ε (Mean = 0.389ε)</span>
1026 <span class="blue">Max = 0.761ε (Mean = 0.444ε)</span>
1033 Bessel K'n: Mathworld Data (Integer Version)
1038 <span class="blue">Max = 0ε (Mean = 0ε)</span>
1043 <span class="blue">Max = 2.16ε (Mean = 1.08ε)</span>
1048 <span class="blue">Max = 2.16ε (Mean = 1.08ε)</span>
1053 <span class="blue">Max = 4.17ε (Mean = 1.75ε)</span>
1060 <br class="table-break"><div class="table">
1061 <a name="math_toolkit.bessel.bessel_derivatives.table_cyl_bessel_k_prime"></a><p class="title"><b>Table 8.55. Error rates for cyl_bessel_k_prime</b></p>
1062 <div class="table-contents"><table class="table" summary="Error rates for cyl_bessel_k_prime">
1075 GNU C++ version 7.1.0<br> linux<br> double
1080 GNU C++ version 7.1.0<br> linux<br> long double
1085 Sun compiler version 0x5150<br> Sun Solaris<br> long double
1090 Microsoft Visual C++ version 14.1<br> Win32<br> double
1098 Bessel K'0: Mathworld Data
1103 <span class="blue">Max = 0ε (Mean = 0ε)</span>
1108 <span class="blue">Max = 0.786ε (Mean = 0.329ε)</span>
1113 <span class="blue">Max = 0.786ε (Mean = 0.329ε)</span>
1118 <span class="blue">Max = 0.786ε (Mean = 0.39ε)</span>
1125 Bessel K'1: Mathworld Data
1130 <span class="blue">Max = 0ε (Mean = 0ε)</span>
1135 <span class="blue">Max = 0.736ε (Mean = 0.389ε)</span>
1140 <span class="blue">Max = 0.736ε (Mean = 0.389ε)</span>
1145 <span class="blue">Max = 0.761ε (Mean = 0.444ε)</span>
1152 Bessel K'n: Mathworld Data
1157 <span class="blue">Max = 0ε (Mean = 0ε)</span>
1162 <span class="blue">Max = 2.16ε (Mean = 1.08ε)</span>
1167 <span class="blue">Max = 2.16ε (Mean = 1.08ε)</span>
1172 <span class="blue">Max = 4.17ε (Mean = 1.75ε)</span>
1179 Bessel K'v: Mathworld Data
1184 <span class="blue">Max = 0ε (Mean = 0ε)</span>
1189 <span class="blue">Max = 3.94ε (Mean = 2.44ε)</span>
1194 <span class="blue">Max = 3.94ε (Mean = 2.34ε)</span>
1199 <span class="blue">Max = 3.94ε (Mean = 1.47ε)</span>
1206 Bessel K'v: Mathworld Data (large values)
1211 <span class="blue">Max = 0ε (Mean = 0ε)</span>
1216 <span class="blue">Max = 59.2ε (Mean = 42.9ε)</span>
1221 <span class="blue">Max = 58.7ε (Mean = 42.6ε)</span>
1226 <span class="blue">Max = 18.6ε (Mean = 11.8ε)</span>
1233 Bessel K'n: Random Data
1238 <span class="blue">Max = 0ε (Mean = 0ε)</span>
1243 <span class="blue">Max = 4.45ε (Mean = 1.19ε)</span>
1248 <span class="blue">Max = 4.45ε (Mean = 1.19ε)</span>
1253 <span class="blue">Max = 9.67ε (Mean = 1.73ε)</span>
1260 Bessel K'v: Random Data
1265 <span class="blue">Max = 0ε (Mean = 0ε)</span>
1270 <span class="blue">Max = 7.95ε (Mean = 1.53ε)</span>
1275 <span class="blue">Max = 7.95ε (Mean = 1.52ε)</span>
1280 <span class="blue">Max = 8.32ε (Mean = 1.65ε)</span>
1287 <br class="table-break"><div class="table">
1288 <a name="math_toolkit.bessel.bessel_derivatives.table_sph_bessel_prime"></a><p class="title"><b>Table 8.56. Error rates for sph_bessel_prime</b></p>
1289 <div class="table-contents"><table class="table" summary="Error rates for sph_bessel_prime">
1302 GNU C++ version 7.1.0<br> linux<br> double
1307 GNU C++ version 7.1.0<br> linux<br> long double
1312 Sun compiler version 0x5150<br> Sun Solaris<br> long double
1317 Microsoft Visual C++ version 14.1<br> Win32<br> double
1324 Bessel j': Random Data
1329 <span class="blue">Max = 0.753ε (Mean = 0.0343ε)</span>
1334 <span class="blue">Max = 167ε (Mean = 12ε)</span>
1339 <span class="blue">Max = 167ε (Mean = 33.2ε)</span>
1344 <span class="blue">Max = 307ε (Mean = 25.2ε)</span>
1350 <br class="table-break"><div class="table">
1351 <a name="math_toolkit.bessel.bessel_derivatives.table_sph_neumann_prime"></a><p class="title"><b>Table 8.57. Error rates for sph_neumann_prime</b></p>
1352 <div class="table-contents"><table class="table" summary="Error rates for sph_neumann_prime">
1365 GNU C++ version 7.1.0<br> linux<br> double
1370 GNU C++ version 7.1.0<br> linux<br> long double
1375 Sun compiler version 0x5150<br> Sun Solaris<br> long double
1380 Microsoft Visual C++ version 14.1<br> Win32<br> double
1392 <span class="blue">Max = 0.988ε (Mean = 0.0869ε)</span>
1397 <span class="blue">Max = 158ε (Mean = 18.8ε)</span>
1402 <span class="blue">Max = 158ε (Mean = 20.2ε)</span>
1407 <span class="blue">Max = 296ε (Mean = 25.6ε)</span>
1413 <br class="table-break"><h5>
1414 <a name="math_toolkit.bessel.bessel_derivatives.h4"></a>
1415 <span class="phrase"><a name="math_toolkit.bessel.bessel_derivatives.implementation"></a></span><a class="link" href="bessel_derivatives.html#math_toolkit.bessel.bessel_derivatives.implementation">Implementation</a>
1418 In the general case, the derivatives are calculated using the relations:
1420 <div class="blockquote"><blockquote class="blockquote"><p>
1421 <span class="inlinemediaobject"><img src="../../../equations/bessel_derivatives1.svg"></span>
1423 </p></blockquote></div>
1425 There are also a number of special cases, for large x we have:
1427 <div class="blockquote"><blockquote class="blockquote"><p>
1428 <span class="inlinemediaobject"><img src="../../../equations/bessel_derivatives4.svg"></span>
1430 </p></blockquote></div>
1434 <div class="blockquote"><blockquote class="blockquote"><p>
1435 <span class="inlinemediaobject"><img src="../../../equations/bessel_derivatives5.svg"></span>
1437 </p></blockquote></div>
1439 <table xmlns:rev="http://www.cs.rpi.edu/~gregod/boost/tools/doc/revision" width="100%"><tr>
1440 <td align="left"></td>
1441 <td align="right"><div class="copyright-footer">Copyright © 2006-2019 Nikhar
1442 Agrawal, Anton Bikineev, Paul A. Bristow, Marco Guazzone, Christopher Kormanyos,
1443 Hubert Holin, Bruno Lalande, John Maddock, Jeremy Murphy, Matthew Pulver, Johan
1444 Råde, Gautam Sewani, Benjamin Sobotta, Nicholas Thompson, Thijs van den Berg,
1445 Daryle Walker and Xiaogang Zhang<p>
1446 Distributed under the Boost Software License, Version 1.0. (See accompanying
1447 file LICENSE_1_0.txt or copy at <a href="http://www.boost.org/LICENSE_1_0.txt" target="_top">http://www.boost.org/LICENSE_1_0.txt</a>)
1452 <div class="spirit-nav">
1453 <a accesskey="p" href="sph_bessel.html"><img src="../../../../../../doc/src/images/prev.png" alt="Prev"></a><a accesskey="u" href="../bessel.html"><img src="../../../../../../doc/src/images/up.png" alt="Up"></a><a accesskey="h" href="../../index.html"><img src="../../../../../../doc/src/images/home.png" alt="Home"></a><a accesskey="n" href="../hankel.html"><img src="../../../../../../doc/src/images/next.png" alt="Next"></a>