annotate libinterp/corefcn/besselj.cc @ 30944:4727765dbc89

maint: Use Octave conventions for coding BIST tests in cset 84944164799e. * besselj.cc (Fairy): Use ';' at end of assert statements within a %!test block. Don't use parentheses around %!error block.
author Rik <rik@octave.org>
date Mon, 18 Apr 2022 13:02:35 -0700
parents 272614f05636
children e0aabe59c1aa
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
27923
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
1 ////////////////////////////////////////////////////////////////////////
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
2 //
30564
796f54d4ddbf update Octave Project Developers copyright for the new year
John W. Eaton <jwe@octave.org>
parents: 29961
diff changeset
3 // Copyright (C) 1997-2022 The Octave Project Developers
27923
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
4 //
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
5 // See the file COPYRIGHT.md in the top-level directory of this
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
6 // distribution or <https://octave.org/copyright/>.
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
7 //
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
8 // This file is part of Octave.
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
9 //
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
10 // Octave is free software: you can redistribute it and/or modify it
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
11 // under the terms of the GNU General Public License as published by
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
12 // the Free Software Foundation, either version 3 of the License, or
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
13 // (at your option) any later version.
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
14 //
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
15 // Octave is distributed in the hope that it will be useful, but
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
16 // WITHOUT ANY WARRANTY; without even the implied warranty of
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
18 // GNU General Public License for more details.
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
19 //
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
20 // You should have received a copy of the GNU General Public License
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
21 // along with Octave; see the file COPYING. If not, see
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
22 // <https://www.gnu.org/licenses/>.
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
23 //
bd51beb6205e update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents: 27919
diff changeset
24 ////////////////////////////////////////////////////////////////////////
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
25
21724
aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef"
John W. Eaton <jwe@octave.org>
parents: 21547
diff changeset
26 #if defined (HAVE_CONFIG_H)
21301
40de9f8f23a6 Use '#include "config.h"' rather than <config.h>.
Rik <rik@octave.org>
parents: 21200
diff changeset
27 # include "config.h"
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
28 #endif
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
29
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
30 #include "lo-specfun.h"
4153
6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe]
jwe
parents: 3810
diff changeset
31 #include "quit.h"
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
32
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14854
diff changeset
33 #include "defun.h"
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
34 #include "error.h"
20940
48b2ad5ee801 maint: Rename oct-obj.[cc|h] to ovl.[cc|h] for clarity.
Rik <rik@octave.org>
parents: 20892
diff changeset
35 #include "ovl.h"
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
36 #include "utils.h"
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
37
29958
32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 29359
diff changeset
38 OCTAVE_NAMESPACE_BEGIN
32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 29359
diff changeset
39
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
40 enum bessel_type
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
41 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
42 BESSEL_J,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
43 BESSEL_Y,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
44 BESSEL_I,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
45 BESSEL_K,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
46 BESSEL_H1,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
47 BESSEL_H2
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
48 };
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
49
22197
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
50 #define DO_BESSEL(type, alpha, x, scaled, ierr, result) \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
51 do \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
52 { \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
53 switch (type) \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
54 { \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
55 case BESSEL_J: \
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
56 result = math::besselj (alpha, x, scaled, ierr); \
22197
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
57 break; \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
58 \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
59 case BESSEL_Y: \
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
60 result = math::bessely (alpha, x, scaled, ierr); \
22197
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
61 break; \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
62 \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
63 case BESSEL_I: \
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
64 result = math::besseli (alpha, x, scaled, ierr); \
22197
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
65 break; \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
66 \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
67 case BESSEL_K: \
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
68 result = math::besselk (alpha, x, scaled, ierr); \
22197
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
69 break; \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
70 \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
71 case BESSEL_H1: \
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
72 result = math::besselh1 (alpha, x, scaled, ierr); \
22197
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
73 break; \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
74 \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
75 case BESSEL_H2: \
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
76 result = math::besselh2 (alpha, x, scaled, ierr); \
22197
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
77 break; \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
78 \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
79 default: \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
80 break; \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
81 } \
e43d83253e28 refill multi-line macro definitions
John W. Eaton <jwe@octave.org>
parents: 22022
diff changeset
82 } \
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
83 while (0)
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
84
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
85 octave_value_list
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
86 do_bessel (enum bessel_type type, const char *fcn,
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
87 const octave_value_list& args, int nargout)
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
88 {
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
89 int nargin = args.length ();
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
90
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
91 if (nargin < 2 || nargin > 3)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
92 print_usage ();
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
93
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
94 octave_value_list retval (nargout > 1 ? 2 : 1);
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
95
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
96 bool scaled = false;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
97 if (nargin == 3)
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
98 {
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
99 octave_value opt_arg = args(2);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
100 bool rpt_error = false;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
101
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
102 if (! opt_arg.is_scalar_type ())
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
103 rpt_error = true;
23586
f6c5db0a02e7 maint: Deprecate is_numeric_type and replace with isnumeric.
Rik <rik@octave.org>
parents: 23579
diff changeset
104 else if (opt_arg.isnumeric ())
20089
93b3757a9623 Use stricter input validation for OPT arg of bessel functions (bug #44867).
Rik <rik@octave.org>
parents: 19697
diff changeset
105 {
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
106 double opt_val = opt_arg.double_value ();
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
107 if (opt_val != 0.0 && opt_val != 1.0)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
108 rpt_error = true;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
109 scaled = (opt_val == 1.0);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
110 }
23579
c20a0fa91c0c maint: Deprecate is_bool_type and replace with islogical.
Rik <rik@octave.org>
parents: 23220
diff changeset
111 else if (opt_arg.islogical ())
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
112 scaled = opt_arg.bool_value ();
20089
93b3757a9623 Use stricter input validation for OPT arg of bessel functions (bug #44867).
Rik <rik@octave.org>
parents: 19697
diff changeset
113
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
114 if (rpt_error)
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
115 error ("%s: OPT must be 0 (or false) or 1 (or true)", fcn);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
116 }
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
117
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
118 octave_value alpha_arg = args(0);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
119 octave_value x_arg = args(1);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
120
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
121 if (alpha_arg.is_single_type () || x_arg.is_single_type ())
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
122 {
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
123 if (alpha_arg.is_scalar_type ())
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
124 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
125 float alpha = args(0).xfloat_value ("%s: ALPHA must be a scalar or matrix", fcn);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
126
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
127 if (x_arg.is_scalar_type ())
20089
93b3757a9623 Use stricter input validation for OPT arg of bessel functions (bug #44867).
Rik <rik@octave.org>
parents: 19697
diff changeset
128 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
129 FloatComplex x = x_arg.xfloat_complex_value ("%s: X must be a scalar or matrix", fcn);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
130
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
131 octave_idx_type ierr;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
132 octave_value result;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
133
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
134 DO_BESSEL (type, alpha, x, scaled, ierr, result);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
135
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
136 retval(0) = result;
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
137 if (nargout > 1)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
138 retval(1) = static_cast<float> (ierr);
20089
93b3757a9623 Use stricter input validation for OPT arg of bessel functions (bug #44867).
Rik <rik@octave.org>
parents: 19697
diff changeset
139 }
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
140 else
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
141 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
142 FloatComplexNDArray x = x_arg.xfloat_complex_array_value ("%s: X must be a scalar or matrix", fcn);
20089
93b3757a9623 Use stricter input validation for OPT arg of bessel functions (bug #44867).
Rik <rik@octave.org>
parents: 19697
diff changeset
143
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
144 Array<octave_idx_type> ierr;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
145 octave_value result;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
146
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
147 DO_BESSEL (type, alpha, x, scaled, ierr, result);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
148
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
149 retval(0) = result;
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
150 if (nargout > 1)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
151 retval(1) = NDArray (ierr);
20089
93b3757a9623 Use stricter input validation for OPT arg of bessel functions (bug #44867).
Rik <rik@octave.org>
parents: 19697
diff changeset
152 }
93b3757a9623 Use stricter input validation for OPT arg of bessel functions (bug #44867).
Rik <rik@octave.org>
parents: 19697
diff changeset
153 }
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
154 else
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
155 {
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
156 dim_vector dv0 = args(0).dims ();
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
157 dim_vector dv1 = args(1).dims ();
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
158
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
159 bool args0_is_row_vector = (dv0(1) == dv0.numel ());
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
160 bool args1_is_col_vector = (dv1(0) == dv1.numel ());
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
161
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
162 if (args0_is_row_vector && args1_is_col_vector)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
163 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
164 FloatRowVector ralpha = args(0).xfloat_row_vector_value ("%s: ALPHA must be a scalar or matrix", fcn);
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
165
27277
db687716fed6 style fixes: generally aim to break long lines before operators, not after
John W. Eaton <jwe@octave.org>
parents: 26376
diff changeset
166 FloatComplexColumnVector cx
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
167 = x_arg.xfloat_complex_column_vector_value ("%s: X must be a scalar or matrix", fcn);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
168
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
169 Array<octave_idx_type> ierr;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
170 octave_value result;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
171
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
172 DO_BESSEL (type, ralpha, cx, scaled, ierr, result);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
173
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
174 retval(0) = result;
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
175 if (nargout > 1)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
176 retval(1) = NDArray (ierr);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
177 }
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
178 else
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
179 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
180 FloatNDArray alpha = args(0).xfloat_array_value ("%s: ALPHA must be a scalar or matrix", fcn);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
181
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
182 if (x_arg.is_scalar_type ())
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
183 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
184 FloatComplex x = x_arg.xfloat_complex_value ("%s: X must be a scalar or matrix", fcn);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
185
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
186 Array<octave_idx_type> ierr;
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
187 octave_value result;
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
188
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
189 DO_BESSEL (type, alpha, x, scaled, ierr, result);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
190
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
191 retval(0) = result;
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
192 if (nargout > 1)
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
193 retval(1) = NDArray (ierr);
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
194 }
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
195 else
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
196 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
197 FloatComplexNDArray x = x_arg.xfloat_complex_array_value ("%s: X must be a scalar or matrix", fcn);
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
198
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
199 Array<octave_idx_type> ierr;
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
200 octave_value result;
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
201
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
202 DO_BESSEL (type, alpha, x, scaled, ierr, result);
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
203
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
204 retval(0) = result;
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
205 if (nargout > 1)
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
206 retval(1) = NDArray (ierr);
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
207 }
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
208 }
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
209 }
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
210 }
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
211 else
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
212 {
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
213 if (alpha_arg.is_scalar_type ())
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
214 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
215 double alpha = args(0).xdouble_value ("%s: ALPHA must be a scalar or matrix", fcn);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
216
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
217 if (x_arg.is_scalar_type ())
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
218 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
219 Complex x = x_arg.xcomplex_value ("%s: X must be a scalar or matrix", fcn);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
220
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
221 octave_idx_type ierr;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
222 octave_value result;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
223
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
224 DO_BESSEL (type, alpha, x, scaled, ierr, result);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
225
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
226 retval(0) = result;
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
227 if (nargout > 1)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
228 retval(1) = static_cast<double> (ierr);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
229 }
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
230 else
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
231 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
232 ComplexNDArray x = x_arg.xcomplex_array_value ("%s: X must be a scalar or matrix", fcn);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
233
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
234 Array<octave_idx_type> ierr;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
235 octave_value result;
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
236
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
237 DO_BESSEL (type, alpha, x, scaled, ierr, result);
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11553
diff changeset
238
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
239 retval(0) = result;
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
240 if (nargout > 1)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
241 retval(1) = NDArray (ierr);
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
242 }
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
243 }
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
244 else
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
245 {
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
246 dim_vector dv0 = args(0).dims ();
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
247 dim_vector dv1 = args(1).dims ();
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
248
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
249 bool args0_is_row_vector = (dv0(1) == dv0.numel ());
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
250 bool args1_is_col_vector = (dv1(0) == dv1.numel ());
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
251
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
252 if (args0_is_row_vector && args1_is_col_vector)
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
253 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
254 RowVector ralpha = args(0).xrow_vector_value ("%s: ALPHA must be a scalar or matrix", fcn);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
255
27277
db687716fed6 style fixes: generally aim to break long lines before operators, not after
John W. Eaton <jwe@octave.org>
parents: 26376
diff changeset
256 ComplexColumnVector cx
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
257 = x_arg.xcomplex_column_vector_value ("%s: X must be a scalar or matrix", fcn);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
258
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
259 Array<octave_idx_type> ierr;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
260 octave_value result;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
261
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
262 DO_BESSEL (type, ralpha, cx, scaled, ierr, result);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
263
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
264 retval(0) = result;
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
265 if (nargout > 1)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
266 retval(1) = NDArray (ierr);
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
267 }
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
268 else
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
269 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
270 NDArray alpha = args(0).xarray_value ("%s: ALPHA must be a scalar or matrix", fcn);
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
271
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
272 if (x_arg.is_scalar_type ())
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
273 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
274 Complex x = x_arg.xcomplex_value ("%s: X must be a scalar or matrix", fcn);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
275
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
276 Array<octave_idx_type> ierr;
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
277 octave_value result;
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
278
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
279 DO_BESSEL (type, alpha, x, scaled, ierr, result);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
280
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
281 retval(0) = result;
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
282 if (nargout > 1)
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
283 retval(1) = NDArray (ierr);
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
284 }
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
285 else
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
286 {
30896
c9788d7f6e65 maint: Use "fcn" as preferred abbreviation for "function" in libinterp/.
Rik <rik@octave.org>
parents: 30564
diff changeset
287 ComplexNDArray x = x_arg.xcomplex_array_value ("%s: X must be a scalar or matrix", fcn);
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
288
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
289 Array<octave_idx_type> ierr;
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
290 octave_value result;
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
291
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
292 DO_BESSEL (type, alpha, x, scaled, ierr, result);
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
293
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
294 retval(0) = result;
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
295 if (nargout > 1)
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
296 retval(1) = NDArray (ierr);
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
297 }
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
298 }
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9732
diff changeset
299 }
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
300 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
301
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
302 return retval;
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
303 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
304
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14854
diff changeset
305 DEFUN (besselj, args, nargout,
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
306 doc: /* -*- texinfo -*-
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
307 @deftypefn {} {@var{J} =} besselj (@var{alpha}, @var{x})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
308 @deftypefnx {} {@var{J} =} besselj (@var{alpha}, @var{x}, @var{opt})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
309 @deftypefnx {} {[@var{J}, @var{ierr}] =} besselj (@dots{})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
310 Compute Bessel functions of the first kind.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
311
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
312 The order of the Bessel function @var{alpha} must be real. The points for
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
313 evaluation @var{x} may be complex.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
314
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
315 If the optional argument @var{opt} is 1 or true, the result @var{J} is
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
316 multiplied by @w{@code{exp (-abs (imag (@var{x})))}}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
317
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
318 If @var{alpha} is a scalar, the result is the same size as @var{x}. If @var{x}
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
319 is a scalar, the result is the same size as @var{alpha}. If @var{alpha} is a
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
320 row vector and @var{x} is a column vector, the result is a matrix with
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
321 @code{length (@var{x})} rows and @code{length (@var{alpha})} columns.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
322 Otherwise, @var{alpha} and @var{x} must conform and the result will be the same
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
323 size.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
324
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
325 If requested, @var{ierr} contains the following status information and is the
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
326 same size as the result.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
327
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
328 @enumerate 0
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
329 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
330 Normal return.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
331
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
332 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
333 Input error, return @code{NaN}.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
334
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
335 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
336 Overflow, return @code{Inf}.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
337
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
338 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
339 Loss of significance by argument reduction results in less than half of machine
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
340 accuracy.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
341
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
342 @item
24906
451f4f291f46 Modified Bessel functions to compute the output with any input
Michele Ginesi <michele.ginesi@gmail.com>
parents: 24594
diff changeset
343 Loss of significance by argument reduction, output may be inaccurate.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
344
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
345 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
346 Error---no computation, algorithm termination condition not met, return
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
347 @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
348 @end enumerate
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
349
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
350 @seealso{bessely, besseli, besselk, besselh}
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
351 @end deftypefn */)
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
352 {
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
353 return do_bessel (BESSEL_J, "besselj", args, nargout);
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
354 }
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
355
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
356 /*
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
357 %!# Function besselj is tested along with other bessels at the end of this file
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
358 */
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
359
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
360 DEFUN (bessely, args, nargout,
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
361 doc: /* -*- texinfo -*-
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
362 @deftypefn {} {@var{Y} =} bessely (@var{alpha}, @var{x})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
363 @deftypefnx {} {@var{Y} =} bessely (@var{alpha}, @var{x}, @var{opt})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
364 @deftypefnx {} {[@var{Y}, @var{ierr}] =} bessely (@dots{})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
365 Compute Bessel functions of the second kind.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
366
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
367 The order of the Bessel function @var{alpha} must be real. The points for
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
368 evaluation @var{x} may be complex.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
369
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
370 If the optional argument @var{opt} is 1 or true, the result @var{Y} is
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
371 multiplied by @w{@code{exp (-abs (imag (@var{x})))}}.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
372
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
373 If @var{alpha} is a scalar, the result is the same size as @var{x}. If @var{x}
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
374 is a scalar, the result is the same size as @var{alpha}. If @var{alpha} is a
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
375 row vector and @var{x} is a column vector, the result is a matrix with
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
376 @code{length (@var{x})} rows and @code{length (@var{alpha})} columns.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
377 Otherwise, @var{alpha} and @var{x} must conform and the result will be the same
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
378 size.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
379
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
380 If requested, @var{ierr} contains the following status information and is the
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
381 same size as the result.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
382
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
383 @enumerate 0
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
384 @item
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
385 Normal return.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
386
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
387 @item
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
388 Input error, return @code{NaN}.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
389
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
390 @item
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
391 Overflow, return @code{Inf}.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
392
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
393 @item
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
394 Loss of significance by argument reduction results in less than half of machine
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
395 accuracy.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
396
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
397 @item
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
398 Complete loss of significance by argument reduction, return @code{NaN}.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
399
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
400 @item
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
401 Error---no computation, algorithm termination condition not met, return
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
402 @code{NaN}.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
403 @end enumerate
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
404
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
405 @seealso{besselj, besseli, besselk, besselh}
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
406 @end deftypefn */)
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
407 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
408 return do_bessel (BESSEL_Y, "bessely", args, nargout);
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
409 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
410
20872
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
411 /*
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
412 %!# Function bessely is tested along with other bessels at the end of this file
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
413 */
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
414
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14854
diff changeset
415 DEFUN (besseli, args, nargout,
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
416 doc: /* -*- texinfo -*-
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
417 @deftypefn {} {@var{I} =} besseli (@var{alpha}, @var{x})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
418 @deftypefnx {} {@var{I} =} besseli (@var{alpha}, @var{x}, @var{opt})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
419 @deftypefnx {} {[@var{I}, @var{ierr}] =} besseli (@dots{})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
420 Compute modified Bessel functions of the first kind.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
421
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
422 The order of the Bessel function @var{alpha} must be real. The points for
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
423 evaluation @var{x} may be complex.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
424
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
425 If the optional argument @var{opt} is 1 or true, the result @var{I} is
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
426 multiplied by @w{@code{exp (-abs (real (@var{x})))}}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
427
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
428 If @var{alpha} is a scalar, the result is the same size as @var{x}. If @var{x}
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
429 is a scalar, the result is the same size as @var{alpha}. If @var{alpha} is a
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
430 row vector and @var{x} is a column vector, the result is a matrix with
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
431 @code{length (@var{x})} rows and @code{length (@var{alpha})} columns.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
432 Otherwise, @var{alpha} and @var{x} must conform and the result will be the same
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
433 size.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
434
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
435 If requested, @var{ierr} contains the following status information and is the
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
436 same size as the result.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
437
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
438 @enumerate 0
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
439 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
440 Normal return.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
441
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
442 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
443 Input error, return @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
444
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
445 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
446 Overflow, return @code{Inf}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
447
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
448 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
449 Loss of significance by argument reduction results in less than half of machine
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
450 accuracy.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
451
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
452 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
453 Complete loss of significance by argument reduction, return @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
454
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
455 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
456 Error---no computation, algorithm termination condition not met, return
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
457 @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
458 @end enumerate
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
459
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
460 @seealso{besselk, besselj, bessely, besselh}
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
461 @end deftypefn */)
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
462
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
463 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
464 return do_bessel (BESSEL_I, "besseli", args, nargout);
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
465 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
466
20872
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
467 /*
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
468 %!# Function besseli is tested along with other bessels at the end of this file
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
469 */
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
470
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14854
diff changeset
471 DEFUN (besselk, args, nargout,
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
472 doc: /* -*- texinfo -*-
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
473 @deftypefn {} {@var{K} =} besselk (@var{alpha}, @var{x})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
474 @deftypefnx {} {@var{K} =} besselk (@var{alpha}, @var{x}, @var{opt})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
475 @deftypefnx {} {[@var{K}, @var{ierr}] =} besselk (@dots{})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
476
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
477 Compute modified Bessel functions of the second kind.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
478
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
479 The order of the Bessel function @var{alpha} must be real. The points for
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
480 evaluation @var{x} may be complex.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
481
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
482 If the optional argument @var{opt} is 1 or true, the result @var{K} is
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
483 multiplied by @w{@code{exp (@var{x})}}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
484
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
485 If @var{alpha} is a scalar, the result is the same size as @var{x}. If @var{x}
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
486 is a scalar, the result is the same size as @var{alpha}. If @var{alpha} is a
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
487 row vector and @var{x} is a column vector, the result is a matrix with
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
488 @code{length (@var{x})} rows and @code{length (@var{alpha})} columns.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
489 Otherwise, @var{alpha} and @var{x} must conform and the result will be the same
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
490 size.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
491
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
492 If requested, @var{ierr} contains the following status information and is the
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
493 same size as the result.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
494
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
495 @enumerate 0
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
496 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
497 Normal return.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
498
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
499 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
500 Input error, return @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
501
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
502 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
503 Overflow, return @code{Inf}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
504
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
505 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
506 Loss of significance by argument reduction results in less than half of machine
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
507 accuracy.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
508
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
509 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
510 Complete loss of significance by argument reduction, return @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
511
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
512 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
513 Error---no computation, algorithm termination condition not met, return
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
514 @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
515 @end enumerate
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
516
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
517 @seealso{besseli, besselj, bessely, besselh}
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
518 @end deftypefn */)
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
519 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
520 return do_bessel (BESSEL_K, "besselk", args, nargout);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
521 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
522
20872
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
523 /*
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
524 %!# Function besselk is tested along with other bessels at the end of this file
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
525 */
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
526
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14854
diff changeset
527 DEFUN (besselh, args, nargout,
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
528 doc: /* -*- texinfo -*-
24594
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
529 @deftypefn {} {@var{H} =} besselh (@var{alpha}, @var{x})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
530 @deftypefnx {} {@var{H} =} besselh (@var{alpha}, @var{k}, @var{x})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
531 @deftypefnx {} {@var{H} =} besselh (@var{alpha}, @var{k}, @var{x}, @var{opt})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
532 @deftypefnx {} {[@var{H}, @var{ierr}] =} besselh (@dots{})
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
533 Compute Bessel functions of the third kind (Hankel functions).
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
534
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
535 The order of the Bessel function @var{alpha} must be real. The kind of Hankel
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
536 function is specified by @var{k} and may be either first (@var{k} = 1) or
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
537 second (@var{k} = 2). The default is Hankel functions of the first kind. The
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
538 points for evaluation @var{x} may be complex.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
539
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
540 If the optional argument @var{opt} is 1 or true, the result is multiplied
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
541 by @code{exp (-I*@var{x})} for @var{k} = 1 or @code{exp (I*@var{x})} for
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
542 @var{k} = 2.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
543
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
544 If @var{alpha} is a scalar, the result is the same size as @var{x}. If @var{x}
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
545 is a scalar, the result is the same size as @var{alpha}. If @var{alpha} is a
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
546 row vector and @var{x} is a column vector, the result is a matrix with
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
547 @code{length (@var{x})} rows and @code{length (@var{alpha})} columns.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
548 Otherwise, @var{alpha} and @var{x} must conform and the result will be the same
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
549 size.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
550
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
551 If requested, @var{ierr} contains the following status information and is the
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
552 same size as the result.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
553
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
554 @enumerate 0
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
555 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
556 Normal return.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
557
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
558 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
559 Input error, return @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
560
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
561 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
562 Overflow, return @code{Inf}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
563
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
564 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
565 Loss of significance by argument reduction results in less than half of machine
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
566 accuracy.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
567
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
568 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
569 Complete loss of significance by argument reduction, return @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
570
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
571 @item
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
572 Error---no computation, algorithm termination condition not met, return
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
573 @code{NaN}.
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
574 @end enumerate
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
575
a81af57f6144 doc: Redo documentation of Bessel functions.
Rik <rik@octave.org>
parents: 24534
diff changeset
576 @seealso{besselj, bessely, besseli, besselk}
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
577 @end deftypefn */)
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
578 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
579 int nargin = args.length ();
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
580
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
581 if (nargin < 2 || nargin > 4)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
582 print_usage ();
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
583
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
584 octave_value_list retval;
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
585
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
586 if (nargin == 2)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
587 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
588 retval = do_bessel (BESSEL_H1, "besselh", args, nargout);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
589 }
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
590 else
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
591 {
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
592 octave_idx_type kind = args(1).xint_value ("besselh: invalid value of K");
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
593
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
594 octave_value_list tmp_args;
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
595
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
596 if (nargin == 4)
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
597 tmp_args(2) = args(3);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
598
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
599 tmp_args(1) = args(2);
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
600 tmp_args(0) = args(0);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
601
20703
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
602 if (kind == 1)
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
603 retval = do_bessel (BESSEL_H1, "besselh", tmp_args, nargout);
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
604 else if (kind == 2)
85e5efae848a eliminate more uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20559
diff changeset
605 retval = do_bessel (BESSEL_H2, "besselh", tmp_args, nargout);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
606 else
20711
7b608fadc663 Make error messages more specific about the variable and problem encountered.
Rik <rik@octave.org>
parents: 20703
diff changeset
607 error ("besselh: K must be 1 or 2");
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
608 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
609
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
610 return retval;
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
611 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
612
20872
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
613 /*
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
614 %!# Function besselh is tested along with other bessels at the end of this file
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
615 */
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
616
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14854
diff changeset
617 DEFUN (airy, args, nargout,
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
618 doc: /* -*- texinfo -*-
30937
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
619 @deftypefn {} {@var{a} =} airy (@var{z})
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
620 @deftypefnx {} {@var{a} =} airy (@var{k}, @var{z})
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
621 @deftypefnx {} {@var{a} =} airy (@var{k}, @var{z}, @var{scale})
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
622 @deftypefnx {} {[@var{a}, @var{ierr}] =} airy (@dots{})
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
623
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
624 Compute Airy functions of the first and second kind, and their derivatives.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
625
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
626 @example
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
627 @group
30937
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
628 K Function Scale factor (if @var{scale} is true)
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
629 --- -------- ---------------------------------------
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
630 0 Ai (Z) exp ((2/3) * Z * sqrt (Z))
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
631 1 dAi(Z)/dZ exp ((2/3) * Z * sqrt (Z))
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
632 2 Bi (Z) exp (-abs (real ((2/3) * Z * sqrt (Z))))
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
633 3 dBi(Z)/dZ exp (-abs (real ((2/3) * Z * sqrt (Z))))
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
634 @end group
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
635 @end example
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
636
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
637 The function call @code{airy (@var{z})} is equivalent to
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
638 @code{airy (0, @var{z})}.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
639
30937
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
640 The optional third input @var{scale} determines whether to
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
641 apply scaling as described above. It is false by default.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
642
30937
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
643 The result @var{a} is the same size as @var{z}.
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
644
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
645 The optional output @var{ierr} contains the following status information and
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
646 is the same size as the result.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
647
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
648 @enumerate 0
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
649 @item
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
650 Normal return.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
651
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
652 @item
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
653 Input error, return @code{NaN}.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
654
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
655 @item
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
656 Overflow, return @code{Inf}.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
657
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
658 @item
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
659 Loss of significance by argument reduction results in less than half
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
660 of machine accuracy.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
661
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
662 @item
24906
451f4f291f46 Modified Bessel functions to compute the output with any input
Michele Ginesi <michele.ginesi@gmail.com>
parents: 24594
diff changeset
663 Loss of significance by argument reduction, output may be inaccurate.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
664
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
665 @item
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
666 Error---no computation, algorithm termination condition not met,
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
667 return @code{NaN}.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
668 @end enumerate
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21786
diff changeset
669 @end deftypefn */)
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
670 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
671 int nargin = args.length ();
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
672
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
673 if (nargin < 1 || nargin > 3)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
674 print_usage ();
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
675
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
676 octave_value_list retval (nargout > 1 ? 2 : 1);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
677
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
678 int kind = 0;
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
679 if (nargin > 1)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
680 {
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
681 kind = args(0).xint_value ("airy: K must be an integer value");
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
682
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
683 if (kind < 0 || kind > 3)
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
684 error ("airy: K must be 0, 1, 2, or 3");
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
685 }
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
686
30937
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
687 bool scale = (nargin == 3) && args(2).xbool_value ("airy: scale option must be a logical value");
20892
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
688
c07bee629973 2015 Code Sprint: use ovl ().
Rik <rik@octave.org>
parents: 20872
diff changeset
689 int idx = (nargin == 1 ? 0 : 1);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
690
30937
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
691 Array<octave_idx_type> ierr;
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
692 octave_value result;
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
693
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
694 if (args(idx).is_single_type ())
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
695 {
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
696 FloatComplexNDArray z = args(idx).xfloat_complex_array_value ("airy: Z must be a complex matrix");
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
697
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
698 if (kind > 1)
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
699 result = math::biry (z, kind == 3, scale, ierr);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
700 else
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
701 result = math::airy (z, kind == 1, scale, ierr);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
702 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
703 else
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
704 {
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
705 ComplexNDArray z = args(idx).xcomplex_array_value ("airy: Z must be a complex matrix");
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
706
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
707 if (kind > 1)
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
708 result = math::biry (z, kind == 3, scale, ierr);
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
709 else
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
710 result = math::airy (z, kind == 1, scale, ierr);
30937
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
711 }
20803
c22206c1a88f eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20741
diff changeset
712
30937
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
713 retval(0) = result;
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
714 if (nargout > 1)
85a67c1a5712 besselj.cc: Improve input validation for function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30896
diff changeset
715 retval(1) = NDArray (ierr);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
716
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
717 return retval;
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
718 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
719
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
720 /*
30940
84944164799e besselj.cc: Add BISTs for the function airy, remove FIXME (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30937
diff changeset
721 %!test <*62321>
30944
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
722 %! assert (airy (0, +1, false), 0.1352924163128814, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
723 %! assert (airy (0, -1, false), 0.5355608832923521, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
724 %! assert (airy (0, +i, false), 0.3314933054321411 - 0.3174498589684437i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
725 %! assert (airy (0, -i, false), 0.3314933054321411 + 0.3174498589684437i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
726 %! assert (airy (0, +1, true), 0.2635136447491401, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
727 %! assert (airy (0, -1, true), 0.4208904755499093 - 0.3311746779333462i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
728 %! assert (airy (0, +i, true), 0.2743053542644657 - 0.0825606941400591i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
729 %! assert (airy (0, -i, true), 0.2743053542644657 + 0.0825606941400591i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
730 %! assert (airy (1, +1, false), -0.1591474412967932, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
731 %! assert (airy (1, -1, false), -0.0101605671166451, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
732 %! assert (airy (1, +i, false), -0.4324926598418070 + 0.0980478562292432i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
733 %! assert (airy (1, -i, false), -0.4324926598418070 - 0.0980478562292432i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
734 %! assert (airy (1, +1, true), -0.3099768889605148, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
735 %! assert (airy (1, -1, true), -0.0079850602592406 + 0.0062829878869967i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
736 %! assert (airy (1, +i, true), -0.2682792317254506 - 0.0680655016581369i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
737 %! assert (airy (1, -i, true), -0.2682792317254506 + 0.0680655016581369i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
738 %! assert (airy (2, +1, false), 1.2074235949528715, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
739 %! assert (airy (2, -1, false), 0.1039973894969446, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
740 %! assert (airy (2, +i, false), 0.6488582083303949 + 0.3449586347680484i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
741 %! assert (airy (2, -i, false), 0.6488582083303949 - 0.3449586347680484i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
742 %! assert (airy (2, +1, true), 0.6199119435726786, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
743 %! assert (airy (2, -1, true), 0.1039973894969446, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
744 %! assert (airy (2, +i, true), 0.4049686654663902 + 0.2152973271658010i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
745 %! assert (airy (2, -i, true), 0.4049686654663902 - 0.2152973271658010i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
746 %! assert (airy (3, +1, false), 0.9324359333927754, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
747 %! assert (airy (3, -1, false), 0.5923756264227923, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
748 %! assert (airy (3, +i, false), 0.1350266467108190 - 0.1288373867812549i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
749 %! assert (airy (3, -i, false), 0.1350266467108190 + 0.1288373867812549i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
750 %! assert (airy (3, +1, true), 0.4787285706049846, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
751 %! assert (airy (3, -1, true), 0.5923756264227923, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
752 %! assert (airy (3, +i, true), 0.0842735134099415 - 0.0804106412111761i, 1e-15);
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
753 %! assert (airy (3, -i, true), 0.0842735134099415 + 0.0804106412111761i, 1e-15);
30940
84944164799e besselj.cc: Add BISTs for the function airy, remove FIXME (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30937
diff changeset
754
30941
272614f05636 besselj.cc: Add more BISTs for the function airy (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30940
diff changeset
755 Input validation tests
30944
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
756 %!error airy ()
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
757 %!error airy (0, 1, 2, 3)
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
758 %!error <K must be an integer value> airy ("foo", 2, false)
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
759 %!error <K must be 0, 1, 2, or 3> airy (3743, 2, false)
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
760 %!error <scale option must be a logical value> airy (0, 2, "foo")
4727765dbc89 maint: Use Octave conventions for coding BIST tests in cset 84944164799e.
Rik <rik@octave.org>
parents: 30941
diff changeset
761 %!error <Z must be a complex matrix> airy (0, "foo", false)
20872
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
762 */
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
763
30940
84944164799e besselj.cc: Add BISTs for the function airy, remove FIXME (bug #62321)
Arun Giridhar <arungiridhar@gmail.com>
parents: 30937
diff changeset
764
20872
ead410040418 2015 Code Sprint: mark bessel functions as tested
Mike Miller <mtmiller@octave.org>
parents: 20853
diff changeset
765 /*
17336
b81b9d079515 Use '##' for comments which stand alone on a line.
Rik <rik@octave.org>
parents: 15195
diff changeset
766 ## Test values computed with GP/PARI version 2.3.3
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
767 %!shared alpha, x, jx, yx, ix, kx, nix
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
768 %!
17336
b81b9d079515 Use '##' for comments which stand alone on a line.
Rik <rik@octave.org>
parents: 15195
diff changeset
769 %! ## Bessel functions, even order, positive and negative x
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
770 %! alpha = 2; x = 1.25;
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
771 %! jx = 0.1710911312405234823613091417;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
772 %! yx = -1.193199310178553861283790424;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
773 %! ix = 0.2220184483766341752692212604;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
774 %! kx = 0.9410016167388185767085460540;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
775 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
776 %!assert (besselj (alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
777 %!assert (bessely (alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
778 %!assert (besseli (alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
779 %!assert (besselk (alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
780 %!assert (besselh (alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
781 %!assert (besselh (alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
782 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
783 %!assert (besselj (alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
784 %!assert (bessely (alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
785 %!assert (besseli (alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
786 %!assert (besselk (alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
787 %!assert (besselh (alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
788 %!assert (besselh (alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
789 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
790 %!assert (besselj (-alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
791 %!assert (bessely (-alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
792 %!assert (besseli (-alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
793 %!assert (besselk (-alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
794 %!assert (besselh (-alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
795 %!assert (besselh (-alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
796 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
797 %!assert (besselj (-alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
798 %!assert (bessely (-alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
799 %!assert (besseli (-alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
800 %!assert (besselk (-alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
801 %!assert (besselh (-alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
802 %!assert (besselh (-alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
803 %!
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
804 %! x *= -1;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
805 %! yx = -1.193199310178553861283790424 + 0.3421822624810469647226182835*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
806 %! kx = 0.9410016167388185767085460540 - 0.6974915263814386815610060884*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
807 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
808 %!assert (besselj (alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
809 %!assert (bessely (alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
810 %!assert (besseli (alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
811 %!assert (besselk (alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
812 %!assert (besselh (alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
813 %!assert (besselh (alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
814 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
815 %!assert (besselj (alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
816 %!assert (bessely (alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
817 %!assert (besseli (alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
818 %!assert (besselk (alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
819 %!assert (besselh (alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
820 %!assert (besselh (alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
821 %!
17336
b81b9d079515 Use '##' for comments which stand alone on a line.
Rik <rik@octave.org>
parents: 15195
diff changeset
822 %! ## Bessel functions, odd order, positive and negative x
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
823 %! alpha = 3; x = 2.5;
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
824 %! jx = 0.2166003910391135247666890035;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
825 %! yx = -0.7560554967536709968379029772;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
826 %! ix = 0.4743704087780355895548240179;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
827 %! kx = 0.2682271463934492027663765197;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
828 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
829 %!assert (besselj (alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
830 %!assert (bessely (alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
831 %!assert (besseli (alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
832 %!assert (besselk (alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
833 %!assert (besselh (alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
834 %!assert (besselh (alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
835 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
836 %!assert (besselj (alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
837 %!assert (bessely (alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
838 %!assert (besseli (alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
839 %!assert (besselk (alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
840 %!assert (besselh (alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
841 %!assert (besselh (alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
842 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
843 %!assert (besselj (-alpha,x), -jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
844 %!assert (bessely (-alpha,x), -yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
845 %!assert (besseli (-alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
846 %!assert (besselk (-alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
847 %!assert (besselh (-alpha,1,x), -(jx + I*yx), 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
848 %!assert (besselh (-alpha,2,x), -(jx - I*yx), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
849 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
850 %!assert (besselj (-alpha,x,1), -jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
851 %!assert (bessely (-alpha,x,1), -yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
852 %!assert (besseli (-alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
853 %!assert (besselk (-alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
854 %!assert (besselh (-alpha,1,x,1), -(jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
855 %!assert (besselh (-alpha,2,x,1), -(jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
856 %!
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
857 %! x *= -1;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
858 %! jx = -jx;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
859 %! yx = 0.7560554967536709968379029772 - 0.4332007820782270495333780070*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
860 %! ix = -ix;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
861 %! kx = -0.2682271463934492027663765197 - 1.490278591297463775542004240*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
862 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
863 %!assert (besselj (alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
864 %!assert (bessely (alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
865 %!assert (besseli (alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
866 %!assert (besselk (alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
867 %!assert (besselh (alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
868 %!assert (besselh (alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
869 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
870 %!assert (besselj (alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
871 %!assert (bessely (alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
872 %!assert (besseli (alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
873 %!assert (besselk (alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
874 %!assert (besselh (alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
875 %!assert (besselh (alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
876 %!
17336
b81b9d079515 Use '##' for comments which stand alone on a line.
Rik <rik@octave.org>
parents: 15195
diff changeset
877 %! ## Bessel functions, fractional order, positive and negative x
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
878 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
879 %! alpha = 3.5; x = 2.75;
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
880 %! jx = 0.1691636439842384154644784389;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
881 %! yx = -0.8301381935499356070267953387;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
882 %! ix = 0.3930540878794826310979363668;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
883 %! kx = 0.2844099013460621170288192503;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
884 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
885 %!assert (besselj (alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
886 %!assert (bessely (alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
887 %!assert (besseli (alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
888 %!assert (besselk (alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
889 %!assert (besselh (alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
890 %!assert (besselh (alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
891 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
892 %!assert (besselj (alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
893 %!assert (bessely (alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
894 %!assert (besseli (alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
895 %!assert (besselk (alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
896 %!assert (besselh (alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
897 %!assert (besselh (alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
898 %!
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
899 %! nix = 0.2119931212254662995364461998;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
900 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
901 %!assert (besselj (-alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
902 %!assert (bessely (-alpha,x), -jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
903 %!assert (besseli (-alpha,x), nix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
904 %!assert (besselk (-alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
905 %!assert (besselh (-alpha,1,x), -I*(jx + I*yx), 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
906 %!assert (besselh (-alpha,2,x), I*(jx - I*yx), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
907 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
908 %!assert (besselj (-alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
909 %!assert (bessely (-alpha,x,1), -jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
910 %!assert (besseli (-alpha,x,1), nix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
911 %!assert (besselk (-alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
912 %!assert (besselh (-alpha,1,x,1), -I*(jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
913 %!assert (besselh (-alpha,2,x,1), I*(jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
914 %!
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
915 %! x *= -1;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
916 %! jx *= -I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
917 %! yx = -0.8301381935499356070267953387*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
918 %! ix *= -I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
919 %! kx = -0.9504059335995575096509874508*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
920 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
921 %!assert (besselj (alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
922 %!assert (bessely (alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
923 %!assert (besseli (alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
924 %!assert (besselk (alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
925 %!assert (besselh (alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
926 %!assert (besselh (alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
927 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
928 %!assert (besselj (alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
929 %!assert (bessely (alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
930 %!assert (besseli (alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
931 %!assert (besselk (alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
932 %!assert (besselh (alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
933 %!assert (besselh (alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
934 %!
17336
b81b9d079515 Use '##' for comments which stand alone on a line.
Rik <rik@octave.org>
parents: 15195
diff changeset
935 %! ## Bessel functions, even order, complex x
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
936 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
937 %! alpha = 2; x = 1.25 + 3.625 * I;
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
938 %! jx = -1.299533366810794494030065917 + 4.370833116012278943267479589*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
939 %! yx = -4.370357232383223896393056727 - 1.283083391453582032688834041*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
940 %! ix = -0.6717801680341515541002273932 - 0.2314623443930774099910228553*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
941 %! kx = -0.01108009888623253515463783379 + 0.2245218229358191588208084197*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
942 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
943 %!assert (besselj (alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
944 %!assert (bessely (alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
945 %!assert (besseli (alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
946 %!assert (besselk (alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
947 %!assert (besselh (alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
948 %!assert (besselh (alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
949 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
950 %!assert (besselj (alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
951 %!assert (bessely (alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
952 %!assert (besseli (alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
953 %!assert (besselk (alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
954 %!assert (besselh (alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
955 %!assert (besselh (alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
956 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
957 %!assert (besselj (-alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
958 %!assert (bessely (-alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
959 %!assert (besseli (-alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
960 %!assert (besselk (-alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
961 %!assert (besselh (-alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
962 %!assert (besselh (-alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
963 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
964 %!assert (besselj (-alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
965 %!assert (bessely (-alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
966 %!assert (besseli (-alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
967 %!assert (besselk (-alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
968 %!assert (besselh (-alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
969 %!assert (besselh (-alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
970 %!
17336
b81b9d079515 Use '##' for comments which stand alone on a line.
Rik <rik@octave.org>
parents: 15195
diff changeset
971 %! ## Bessel functions, odd order, complex x
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
972 %!
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
973 %! alpha = 3; x = 2.5 + 1.875 * I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
974 %! jx = 0.1330721523048277493333458596 + 0.5386295217249660078754395597*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
975 %! yx = -0.6485072392105829901122401551 + 0.2608129289785456797046996987*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
976 %! ix = -0.6182064685486998097516365709 + 0.4677561094683470065767989920*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
977 %! kx = -0.1568585587733540007867882337 - 0.05185853709490846050505141321*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
978 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
979 %!assert (besselj (alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
980 %!assert (bessely (alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
981 %!assert (besseli (alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
982 %!assert (besselk (alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
983 %!assert (besselh (alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
984 %!assert (besselh (alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
985 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
986 %!assert (besselj (alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
987 %!assert (bessely (alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
988 %!assert (besseli (alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
989 %!assert (besselk (alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
990 %!assert (besselh (alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
991 %!assert (besselh (alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
992 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
993 %!assert (besselj (-alpha,x), -jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
994 %!assert (bessely (-alpha,x), -yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
995 %!assert (besseli (-alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
996 %!assert (besselk (-alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
997 %!assert (besselh (-alpha,1,x), -(jx + I*yx), 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
998 %!assert (besselh (-alpha,2,x), -(jx - I*yx), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
999 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1000 %!assert (besselj (-alpha,x,1), -jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1001 %!assert (bessely (-alpha,x,1), -yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1002 %!assert (besseli (-alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1003 %!assert (besselk (-alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1004 %!assert (besselh (-alpha,1,x,1), -(jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1005 %!assert (besselh (-alpha,2,x,1), -(jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1006 %!
17336
b81b9d079515 Use '##' for comments which stand alone on a line.
Rik <rik@octave.org>
parents: 15195
diff changeset
1007 %! ## Bessel functions, fractional order, complex x
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1008 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1009 %! alpha = 3.5; x = 1.75 + 4.125 * I;
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1010 %! jx = -3.018566131370455929707009100 - 0.7585648436793900607704057611*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1011 %! yx = 0.7772278839106298215614791107 - 3.018518722313849782683792010*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1012 %! ix = 0.2100873577220057189038160913 - 0.6551765604618246531254970926*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1013 %! kx = 0.1757147290513239935341488069 + 0.08772348296883849205562558311*I;
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1014 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1015 %!assert (besselj (alpha,x), jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1016 %!assert (bessely (alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1017 %!assert (besseli (alpha,x), ix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1018 %!assert (besselk (alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1019 %!assert (besselh (alpha,1,x), jx + I*yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1020 %!assert (besselh (alpha,2,x), jx - I*yx, 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1021 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1022 %!assert (besselj (alpha,x,1), jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1023 %!assert (bessely (alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1024 %!assert (besseli (alpha,x,1), ix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1025 %!assert (besselk (alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1026 %!assert (besselh (alpha,1,x,1), (jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1027 %!assert (besselh (alpha,2,x,1), (jx - I*yx)*exp (I*x), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1028 %!
17337
a31b54b5f84a Use only 1 space between '%!' and start of test/demo code.
Rik <rik@octave.org>
parents: 17336
diff changeset
1029 %! nix = 0.09822388691172060573913739253 - 0.7110230642207380127317227407*I;
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1030 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1031 %!assert (besselj (-alpha,x), yx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1032 %!assert (bessely (-alpha,x), -jx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1033 %!assert (besseli (-alpha,x), nix, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1034 %!assert (besselk (-alpha,x), kx, 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1035 %!assert (besselh (-alpha,1,x), -I*(jx + I*yx), 100*eps)
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1036 %!assert (besselh (-alpha,2,x), I*(jx - I*yx), 100*eps)
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1037 %!
28915
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1038 %!assert (besselj (-alpha,x,1), yx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1039 %!assert (bessely (-alpha,x,1), -jx*exp (-abs (imag (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1040 %!assert (besseli (-alpha,x,1), nix*exp (-abs (real (x))), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1041 %!assert (besselk (-alpha,x,1), kx*exp (x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1042 %!assert (besselh (-alpha,1,x,1), -I*(jx + I*yx)*exp (-I*x), 100*eps)
c40a367a84c0 maint: Use Octave convention of space after function name in libinterp/.
Rik <rik@octave.org>
parents: 27923
diff changeset
1043 %!assert (besselh (-alpha,2,x,1), I*(jx - I*yx)*exp (I*x), 100*eps)
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1044
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1045 Tests contributed by Robert T. Short.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1046 Tests are based on the properties and tables in A&S:
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1047 Abramowitz and Stegun, "Handbook of Mathematical Functions",
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1048 1972.
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1049
21751
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1050 For regular Bessel functions, there are 3 tests. These compare octave
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1051 results against Tables 9.1, 9.2, and 9.4 in A&S. Tables 9.1 and 9.2
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1052 are good to only a few decimal places, so any failures should be
21751
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1053 considered a broken implementation. Table 9.4 is an extended table
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1054 for larger orders and arguments. There are some differences between
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1055 Octave and Table 9.4, mostly in the last decimal place but in a very
21751
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1056 few instances the errors are in the last two places. The comparison
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1057 tolerance has been changed to reflect this.
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1058
21751
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1059 Similarly for modified Bessel functions, there are 3 tests. These
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1060 compare octave results against Tables 9.8, 9.9, and 9.11 in A&S.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1061 Tables 9.8 and 9.9 are good to only a few decimal places, so any
21751
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1062 failures should be considered a broken implementation. Table 9.11 is
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1063 an extended table for larger orders and arguments. There are some
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1064 differences between octave and Table 9.11, mostly in the last decimal
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1065 place but in a very few instances the errors are in the last two
21751
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1066 places. The comparison tolerance has been changed to reflect this.
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1067
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1068 For spherical Bessel functions, there are also three tests, comparing
21751
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1069 octave results to Tables 10.1, 10.2, and 10.4 in A&S. Very similar
b571fc85953f maint: Use two spaces after period to indicate sentence break.
Rik <rik@octave.org>
parents: 21724
diff changeset
1070 comments may be made here as in the previous lines. At this time,
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1071 modified spherical Bessel function tests are not included.
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1072
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1073 % Table 9.1 - J and Y for integer orders 0, 1, 2.
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1074 % Compare against excerpts of Table 9.1, Abramowitz and Stegun.
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1075 %!test
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1076 %! n = 0:2;
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1077 %! z = (0:2.5:17.5)';
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1078 %!
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1079 %! Jt = [[ 1.000000000000000, 0.0000000000, 0.0000000000];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1080 %! [-0.048383776468198, 0.4970941025, 0.4460590584];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1081 %! [-0.177596771314338, -0.3275791376, 0.0465651163];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1082 %! [ 0.266339657880378, 0.1352484276, -0.2302734105];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1083 %! [-0.245935764451348, 0.0434727462, 0.2546303137];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1084 %! [ 0.146884054700421, -0.1654838046, -0.1733614634];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1085 %! [-0.014224472826781, 0.2051040386, 0.0415716780];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1086 %! [-0.103110398228686, -0.1634199694, 0.0844338303]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1087 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1088 %! Yt = [[-Inf, -Inf, -Inf ];
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1089 %! [ 0.4980703596, 0.1459181380, -0.38133585 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1090 %! [-0.3085176252, 0.1478631434, 0.36766288 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1091 %! [ 0.1173132861, -0.2591285105, -0.18641422 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1092 %! [ 0.0556711673, 0.2490154242, -0.00586808 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1093 %! [-0.1712143068, -0.1538382565, 0.14660019 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1094 %! [ 0.2054642960, 0.0210736280, -0.20265448 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1095 %! [-0.1604111925, 0.0985727987, 0.17167666 ]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1096 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1097 %! J = besselj (n,z);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1098 %! Y = bessely (n,z);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1099 %! assert (Jt(:,1), J(:,1), 0.5e-10);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1100 %! assert (Yt(:,1), Y(:,1), 0.5e-10);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1101 %! assert (Jt(:,2:3), J(:,2:3), 0.5e-10);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1102
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1103 Table 9.2 - J and Y for integer orders 3-9.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1104
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1105 %!test
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1106 %! n = (3:9);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1107 %! z = (0:2:20).';
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1108 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1109 %! Jt = [[ 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00, 0.0000e+00];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1110 %! [ 1.2894e-01, 3.3996e-02, 7.0396e-03, 1.2024e-03, 1.7494e-04, 2.2180e-05, 2.4923e-06];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1111 %! [ 4.3017e-01, 2.8113e-01, 1.3209e-01, 4.9088e-02, 1.5176e-02, 4.0287e-03, 9.3860e-04];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1112 %! [ 1.1477e-01, 3.5764e-01, 3.6209e-01, 2.4584e-01, 1.2959e-01, 5.6532e-02, 2.1165e-02];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1113 %! [-2.9113e-01,-1.0536e-01, 1.8577e-01, 3.3758e-01, 3.2059e-01, 2.2345e-01, 1.2632e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1114 %! [ 5.8379e-02,-2.1960e-01,-2.3406e-01,-1.4459e-02, 2.1671e-01, 3.1785e-01, 2.9186e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1115 %! [ 1.9514e-01, 1.8250e-01,-7.3471e-02,-2.4372e-01,-1.7025e-01, 4.5095e-02, 2.3038e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1116 %! [-1.7681e-01, 7.6244e-02, 2.2038e-01, 8.1168e-02,-1.5080e-01,-2.3197e-01,-1.1431e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1117 %! [-4.3847e-02,-2.0264e-01,-5.7473e-02, 1.6672e-01, 1.8251e-01,-7.0211e-03,-1.8953e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1118 %! [ 1.8632e-01, 6.9640e-02,-1.5537e-01,-1.5596e-01, 5.1399e-02, 1.9593e-01, 1.2276e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1119 %! [-9.8901e-02, 1.3067e-01, 1.5117e-01,-5.5086e-02,-1.8422e-01,-7.3869e-02, 1.2513e-01]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1120 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1121 %! Yt = [[ -Inf, -Inf, -Inf, -Inf, -Inf, -Inf, -Inf];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1122 %! [-1.1278e+00,-2.7659e+00,-9.9360e+00,-4.6914e+01,-2.7155e+02,-1.8539e+03,-1.4560e+04];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1123 %! [-1.8202e-01,-4.8894e-01,-7.9585e-01,-1.5007e+00,-3.7062e+00,-1.1471e+01,-4.2178e+01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1124 %! [ 3.2825e-01, 9.8391e-02,-1.9706e-01,-4.2683e-01,-6.5659e-01,-1.1052e+00,-2.2907e+00];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1125 %! [ 2.6542e-02, 2.8294e-01, 2.5640e-01, 3.7558e-02,-2.0006e-01,-3.8767e-01,-5.7528e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1126 %! [-2.5136e-01,-1.4495e-01, 1.3540e-01, 2.8035e-01, 2.0102e-01, 1.0755e-03,-1.9930e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1127 %! [ 1.2901e-01,-1.5122e-01,-2.2982e-01,-4.0297e-02, 1.8952e-01, 2.6140e-01, 1.5902e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1128 %! [ 1.2350e-01, 2.0393e-01,-6.9717e-03,-2.0891e-01,-1.7209e-01, 3.6816e-02, 2.1417e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1129 %! [-1.9637e-01,-7.3222e-05, 1.9633e-01, 1.2278e-01,-1.0425e-01,-2.1399e-01,-1.0975e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1130 %! [ 3.3724e-02,-1.7722e-01,-1.1249e-01, 1.1472e-01, 1.8897e-01, 3.2253e-02,-1.6030e-01];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1131 %! [ 1.4967e-01, 1.2409e-01,-1.0004e-01,-1.7411e-01,-4.4312e-03, 1.7101e-01, 1.4124e-01]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1132 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1133 %! n = (3:9);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1134 %! z = (0:2:20).';
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1135 %! J = besselj (n,z);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1136 %! Y = bessely (n,z);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1137 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1138 %! assert (J(1,:), zeros (1, columns (J)));
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1139 %! assert (J(2:end,:), Jt(2:end,:), -5e-5);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1140 %! assert (Yt(1,:), Y(1,:));
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1141 %! assert (Y(2:end,:), Yt(2:end,:), -5e-5);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1142
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1143 Table 9.4 - J and Y for various integer orders and arguments.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1144
14450
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1145 %!test
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1146 %! Jt = [[ 7.651976866e-01, 2.238907791e-01, -1.775967713e-01, -2.459357645e-01, 5.581232767e-02, 1.998585030e-02];
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1147 %! [ 2.497577302e-04, 7.039629756e-03, 2.611405461e-01, -2.340615282e-01, -8.140024770e-02, -7.419573696e-02];
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1148 %! [ 2.630615124e-10, 2.515386283e-07, 1.467802647e-03, 2.074861066e-01, -1.138478491e-01, -5.473217694e-02];
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1149 %! [ 2.297531532e-17, 7.183016356e-13, 4.796743278e-07, 4.507973144e-03, -1.082255990e-01, 1.519812122e-02];
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1150 %! [ 3.873503009e-25, 3.918972805e-19, 2.770330052e-11, 1.151336925e-05, -1.167043528e-01, 6.221745850e-02];
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1151 %! [ 3.482869794e-42, 3.650256266e-33, 2.671177278e-21, 1.551096078e-12, 4.843425725e-02, 8.146012958e-02];
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1152 %! [ 1.107915851e-60, 1.196077458e-48, 8.702241617e-33, 6.030895312e-21, -1.381762812e-01, 7.270175482e-02];
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1153 %! [ 2.906004948e-80, 3.224095839e-65, 2.294247616e-45, 1.784513608e-30, 1.214090219e-01, -3.869833973e-02];
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1154 %! [ 8.431828790e-189, 1.060953112e-158, 6.267789396e-119, 6.597316064e-89, 1.115927368e-21, 9.636667330e-02]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1155 %!
14450
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1156 %! Yt = [[ 8.825696420e-02, 5.103756726e-01, -3.085176252e-01, 5.567116730e-02, -9.806499547e-02, -7.724431337e-02]
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1157 %! [-2.604058666e+02, -9.935989128e+00, -4.536948225e-01, 1.354030477e-01, -7.854841391e-02, -2.948019628e-02]
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1158 %! [-1.216180143e+08, -1.291845422e+05, -2.512911010e+01, -3.598141522e-01, 5.723897182e-03, 5.833157424e-02]
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1159 %! [-9.256973276e+14, -2.981023646e+10, -4.694049564e+04, -6.364745877e+00, 4.041280205e-02, 7.879068695e-02]
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1160 %! [-4.113970315e+22, -4.081651389e+16, -5.933965297e+08, -1.597483848e+03, 1.644263395e-02, 5.124797308e-02]
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1161 %! [-3.048128783e+39, -2.913223848e+30, -4.028568418e+18, -7.256142316e+09, -1.164572349e-01, 6.138839212e-03]
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1162 %! [-7.184874797e+57, -6.661541235e+45, -9.216816571e+29, -1.362803297e+18, -4.530801120e-02, 4.074685217e-02]
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1163 %! [-2.191142813e+77, -1.976150576e+62, -2.788837017e+42, -3.641066502e+27, -2.103165546e-01, 7.650526394e-02]
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1164 %! [-3.775287810e+185, -3.000826049e+155, -5.084863915e+115, -4.849148271e+85, -3.293800188e+18, -1.669214114e-01]];
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1165 %!
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1166 %! n = [(0:5:20).';30;40;50;100];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1167 %! z = [1,2,5,10,50,100];
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1168 %! J = besselj (n.', z.').';
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1169 %! Y = bessely (n.', z.').';
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1170 %! assert (J, Jt, -1e-9);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1171 %! assert (Y, Yt, -1e-9);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1172
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1173 Table 9.8 - I and K for integer orders 0, 1, 2.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1174
14450
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1175 %!test
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1176 %! n = 0:2;
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1177 %! z1 = [0.1;2.5;5.0];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1178 %! z2 = [7.5;10.0;15.0;20.0];
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1179 %! rtbl = [[ 0.9071009258 0.0452984468 0.1251041992 2.6823261023 10.890182683 1.995039646 ];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1180 %! [ 0.2700464416 0.2065846495 0.2042345837 0.7595486903 0.9001744239 0.759126289 ];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1181 %! [ 0.1835408126 0.1639722669 0.7002245988 0.5478075643 0.6002738588 0.132723593 ];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1182 %! [ 0.1483158301 0.1380412115 0.111504840 0.4505236991 0.4796689336 0.57843541 ];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1183 %! [ 0.1278333372 0.1212626814 0.103580801 0.3916319344 0.4107665704 0.47378525 ];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1184 %! [ 0.1038995314 0.1003741751 0.090516308 0.3210023535 0.3315348950 0.36520701 ];
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1185 %! [ 0.0897803119 0.0875062222 0.081029690 0.2785448768 0.2854254970 0.30708743 ]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1186 %!
14450
ab4676288414 Fix typos and errors in new bessel function tests
Mike Miller <mtmiller@ieee.org>
parents: 14426
diff changeset
1187 %! tbl = [besseli(n,z1,1), besselk(n,z1,1)];
14854
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1188 %! tbl(:,3) = tbl(:,3) .* (exp (z1) .* z1.^(-2));
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1189 %! tbl(:,6) = tbl(:,6) .* (exp (-z1) .* z1.^(2));
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1190 %! tbl = [tbl;[besseli(n,z2,1),besselk(n,z2,1)]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1191 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1192 %! assert (tbl, rtbl, -2e-8);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1193
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1194 Table 9.9 - I and K for orders 3-9.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1195
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1196 %!test
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1197 %! It = [[ 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00 0.0000e+00];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1198 %! [ 2.8791e-02 6.8654e-03 1.3298e-03 2.1656e-04 3.0402e-05 3.7487e-06 4.1199e-07];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1199 %! [ 6.1124e-02 2.5940e-02 9.2443e-03 2.8291e-03 7.5698e-04 1.7968e-04 3.8284e-05];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1200 %! [ 7.4736e-02 4.1238e-02 1.9752e-02 8.3181e-03 3.1156e-03 1.0484e-03 3.1978e-04];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1201 %! [ 7.9194e-02 5.0500e-02 2.8694e-02 1.4633e-02 6.7449e-03 2.8292e-03 1.0866e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1202 %! [ 7.9830e-02 5.5683e-02 3.5284e-02 2.0398e-02 1.0806e-02 5.2694e-03 2.3753e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1203 %! [ 7.8848e-02 5.8425e-02 3.9898e-02 2.5176e-02 1.4722e-02 8.0010e-03 4.0537e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1204 %! [ 7.7183e-02 5.9723e-02 4.3056e-02 2.8969e-02 1.8225e-02 1.0744e-02 5.9469e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1205 %! [ 7.5256e-02 6.0155e-02 4.5179e-02 3.1918e-02 2.1240e-02 1.3333e-02 7.9071e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1206 %! [ 7.3263e-02 6.0059e-02 4.6571e-02 3.4186e-02 2.3780e-02 1.5691e-02 9.8324e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1207 %! [ 7.1300e-02 5.9640e-02 4.7444e-02 3.5917e-02 2.5894e-02 1.7792e-02 1.1661e-02]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1208 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1209 %! Kt = [[ Inf Inf Inf Inf Inf Inf Inf];
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1210 %! [ 4.7836e+00 1.6226e+01 6.9687e+01 3.6466e+02 2.2576e+03 1.6168e+04 1.3160e+05];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1211 %! [ 1.6317e+00 3.3976e+00 8.4268e+00 2.4465e+01 8.1821e+01 3.1084e+02 1.3252e+03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1212 %! [ 9.9723e-01 1.6798e+00 3.2370e+00 7.0748e+00 1.7387e+01 4.7644e+01 1.4444e+02];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1213 %! [ 7.3935e-01 1.1069e+00 1.8463e+00 3.4148e+00 6.9684e+00 1.5610e+01 3.8188e+01];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1214 %! [ 6.0028e-01 8.3395e-01 1.2674e+00 2.1014e+00 3.7891e+00 7.4062e+00 1.5639e+01];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1215 %! [ 5.1294e-01 6.7680e-01 9.6415e-01 1.4803e+00 2.4444e+00 4.3321e+00 8.2205e+00];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1216 %! [ 4.5266e-01 5.7519e-01 7.8133e-01 1.1333e+00 1.7527e+00 2.8860e+00 5.0510e+00];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1217 %! [ 4.0829e-01 5.0414e-01 6.6036e-01 9.1686e-01 1.3480e+00 2.0964e+00 3.4444e+00];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1218 %! [ 3.7411e-01 4.5162e-01 5.7483e-01 7.7097e-01 1.0888e+00 1.6178e+00 2.5269e+00];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1219 %! [ 3.4684e-01 4.1114e-01 5.1130e-01 6.6679e-01 9.1137e-01 1.3048e+00 1.9552e+00]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1220 %!
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1221 %! n = (3:9);
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1222 %! z = (0:2:20).';
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1223 %! I = besseli (n,z,1);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1224 %! K = besselk (n,z,1);
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1225 %!
14854
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1226 %! assert (abs (I(1,:)), zeros (1, columns (I)));
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1227 %! assert (I(2:end,:), It(2:end,:), -5e-5);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1228 %! assert (Kt(1,:), K(1,:));
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1229 %! assert (K(2:end,:), Kt(2:end,:), -5e-5);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1230
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1231 Table 9.11 - I and K for various integer orders and arguments.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1232
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1233 %!test
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1234 %! It = [[ 1.266065878e+00 2.279585302e+00 2.723987182e+01 2.815716628e+03 2.93255378e+20 1.07375171e+42 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1235 %! [ 2.714631560e-04 9.825679323e-03 2.157974547e+00 7.771882864e+02 2.27854831e+20 9.47009387e+41 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1236 %! [ 2.752948040e-10 3.016963879e-07 4.580044419e-03 2.189170616e+01 1.07159716e+20 6.49897552e+41 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1237 %! [ 2.370463051e-17 8.139432531e-13 1.047977675e-06 1.043714907e-01 3.07376455e+19 3.47368638e+41 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1238 %! [ 3.966835986e-25 4.310560576e-19 5.024239358e-11 1.250799736e-04 5.44200840e+18 1.44834613e+41 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1239 %! [ 3.539500588e-42 3.893519664e-33 3.997844971e-21 7.787569783e-12 4.27499365e+16 1.20615487e+40 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1240 %! [ 1.121509741e-60 1.255869192e-48 1.180426980e-32 2.042123274e-20 6.00717897e+13 3.84170550e+38 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1241 %! [ 2.934635309e-80 3.353042830e-65 2.931469647e-45 4.756894561e-30 1.76508024e+10 4.82195809e+36 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1242 %! [ 8.473674008e-189 1.082171475e-158 7.093551489e-119 1.082344202e-88 2.72788795e-16 4.64153494e+21 ]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1243 %!
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1244 %! Kt = [[ 4.210244382e-01 1.138938727e-01 3.691098334e-03 1.778006232e-05 3.41016774e-23 4.65662823e-45 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1245 %! [ 3.609605896e+02 9.431049101e+00 3.270627371e-02 5.754184999e-05 4.36718224e-23 5.27325611e-45 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1246 %! [ 1.807132899e+08 1.624824040e+05 9.758562829e+00 1.614255300e-03 9.15098819e-23 7.65542797e-45 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1247 %! [ 1.403066801e+15 4.059213332e+10 3.016976630e+04 2.656563849e-01 3.11621117e-22 1.42348325e-44 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1248 %! [ 6.294369360e+22 5.770856853e+16 4.827000521e+08 1.787442782e+02 1.70614838e-21 3.38520541e-44 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1249 %! [ 4.706145527e+39 4.271125755e+30 4.112132063e+18 2.030247813e+09 2.00581681e-19 3.97060205e-43 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1250 %! [ 1.114220651e+58 9.940839886e+45 1.050756722e+30 5.938224681e+17 1.29986971e-16 1.20842080e-41 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1251 %! [ 3.406896854e+77 2.979981740e+62 3.394322243e+42 2.061373775e+27 4.00601349e-13 9.27452265e-40 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1252 %! [ 5.900333184e+185 4.619415978e+155 7.039860193e+115 4.596674084e+85 1.63940352e+13 7.61712963e-25 ]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1253 %!
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1254 %! n = [(0:5:20).';30;40;50;100];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1255 %! z = [1,2,5,10,50,100];
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1256 %! I = besseli (n.', z.').';
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1257 %! K = besselk (n.', z.').';
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1258 %! assert (I, It, -5e-9);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1259 %! assert (K, Kt, -5e-9);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1260
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1261 The next section checks that negative integer orders and positive
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1262 integer orders are appropriately related.
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1263
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1264 %!test
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1265 %! n = (0:2:20);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1266 %! assert (besselj (n,1), besselj (-n,1), 1e-8);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1267 %! assert (-besselj (n+1,1), besselj (-n-1,1), 1e-8);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1268
14854
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1269 besseli (n,z) = besseli (-n,z);
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1270
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1271 %!test
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1272 %! n = (0:2:20);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1273 %! assert (besseli (n,1), besseli (-n,1), 1e-8);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1274
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1275 Table 10.1 - j and y for integer orders 0, 1, 2.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1276 Compare against excerpts of Table 10.1, Abramowitz and Stegun.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1277
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1278 %!test
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1279 %! n = (0:2);
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1280 %! z = [0.1;(2.5:2.5:10.0).'];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1281 %!
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1282 %! jt = [[ 9.9833417e-01 3.33000119e-02 6.6619061e-04 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1283 %! [ 2.3938886e-01 4.16212989e-01 2.6006673e-01 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1284 %! [-1.9178485e-01 -9.50894081e-02 1.3473121e-01 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1285 %! [ 1.2507e-01 -2.9542e-02 -1.3688e-01 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1286 %! [ -5.4402e-02 7.8467e-02 7.7942e-02 ]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1287 %!
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1288 %! yt = [[-9.9500417e+00 -1.0049875e+02 -3.0050125e+03 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1289 %! [ 3.2045745e-01 -1.1120588e-01 -4.5390450e-01 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1290 %! [-5.6732437e-02 1.8043837e-01 1.6499546e-01 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1291 %! [ -4.6218e-02 -1.3123e-01 -6.2736e-03 ];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1292 %! [ 8.3907e-02 6.2793e-02 -6.5069e-02 ]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1293 %!
14854
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1294 %! j = sqrt ((pi/2)./z) .* besselj (n+1/2,z);
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1295 %! y = sqrt ((pi/2)./z) .* bessely (n+1/2,z);
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1296 %! assert (jt, j, -5e-5);
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1297 %! assert (yt, y, -5e-5);
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1298
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1299 Table 10.2 - j and y for orders 3-8.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1300 Compare against excerpts of Table 10.2, Abramowitzh and Stegun.
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1301
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1302 Important note: In A&S, y_4(0.1) = -1.0507e+7, but Octave returns
14854
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1303 y_4(0.1) = -1.0508e+07 (-10507503.75). If I compute the same term using
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1304 a series, the difference is in the eighth significant digit so I left
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1305 the Octave results in place.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1306
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1307 %!test
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1308 %! n = (3:8);
14854
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
1309 %! z = (0:2.5:10).'; z(1) = 0.1;
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1310 %!
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1311 %! jt = [[ 9.5185e-06 1.0577e-07 9.6163e-10 7.3975e-12 4.9319e-14 2.9012e-16];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1312 %! [ 1.0392e-01 3.0911e-02 7.3576e-03 1.4630e-03 2.5009e-04 3.7516e-05];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1313 %! [ 2.2982e-01 1.8702e-01 1.0681e-01 4.7967e-02 1.7903e-02 5.7414e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1314 %! [-6.1713e-02 7.9285e-02 1.5685e-01 1.5077e-01 1.0448e-01 5.8188e-02];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1315 %! [-3.9496e-02 -1.0559e-01 -5.5535e-02 4.4501e-02 1.1339e-01 1.2558e-01]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1316 %!
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1317 %! yt = [[-1.5015e+05 -1.0508e+07 -9.4553e+08 -1.0400e+11 -1.3519e+13 -2.0277e+15];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1318 %! [-7.9660e-01 -1.7766e+00 -5.5991e+00 -2.2859e+01 -1.1327e+02 -6.5676e+02];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1319 %! [-1.5443e-02 -1.8662e-01 -3.2047e-01 -5.1841e-01 -1.0274e+00 -2.5638e+00];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1320 %! [ 1.2705e-01 1.2485e-01 2.2774e-02 -9.1449e-02 -1.8129e-01 -2.7112e-01];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1321 %! [-9.5327e-02 -1.6599e-03 9.3834e-02 1.0488e-01 4.2506e-02 -4.1117e-02]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1322 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1323 %! j = sqrt ((pi/2)./z) .* besselj (n+1/2,z);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1324 %! y = sqrt ((pi/2)./z) .* bessely (n+1/2,z);
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1325 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1326 %! assert (jt, j, -5e-5);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1327 %! assert (yt, y, -5e-5);
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1328
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1329 Table 10.4 - j and y for various integer orders and arguments.
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1330
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1331 %!test
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1332 %! jt = [[ 8.414709848e-01 4.546487134e-01 -1.917848549e-01 -5.440211109e-02 -5.247497074e-03 -5.063656411e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1333 %! [ 9.256115861e-05 2.635169770e-03 1.068111615e-01 -5.553451162e-02 -2.004830056e-02 -9.290148935e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1334 %! [ 7.116552640e-11 6.825300865e-08 4.073442442e-04 6.460515449e-02 -1.503922146e-02 -1.956578597e-04];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1335 %! [ 5.132686115e-18 1.606982166e-13 1.084280182e-07 1.063542715e-03 -1.129084539e-02 7.877261748e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1336 %! [ 7.537795722e-26 7.632641101e-20 5.427726761e-12 2.308371961e-06 -1.578502990e-02 1.010767128e-02];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1337 %! [ 5.566831267e-43 5.836617888e-34 4.282730217e-22 2.512057385e-13 -1.494673454e-03 8.700628514e-03];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1338 %! [ 1.538210374e-61 1.660978779e-49 1.210347583e-33 8.435671634e-22 -2.606336952e-02 1.043410851e-02];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1339 %! [ 3.615274717e-81 4.011575290e-66 2.857479350e-46 2.230696023e-31 1.882910737e-02 5.797140882e-04];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1340 %! [7.444727742e-190 9.367832591e-160 5.535650303e-120 5.832040182e-90 1.019012263e-22 1.088047701e-02]];
14426
e28a1723d5a2 besselj: Style fixes on new tests due to Robert T. Short
Jordi Gutiérrez Hermoso <jordigh@octave.org>
parents: 14425
diff changeset
1341 %!
14425
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1342 %! yt = [[ -5.403023059e-01 2.080734183e-01 -5.673243709e-02 8.390715291e-02 -1.929932057e-02 -8.623188723e-03]
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1343 %! [ -9.994403434e+02 -1.859144531e+01 -3.204650467e-01 9.383354168e-02 -6.971131965e-04 3.720678486e-03]
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1344 %! [ -6.722150083e+08 -3.554147201e+05 -2.665611441e+01 -1.724536721e-01 1.352468751e-02 1.002577737e-02]
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1345 %! [ -6.298007233e+15 -1.012182944e+11 -6.288146513e+04 -3.992071745e+00 1.712319725e-02 6.258641510e-03]
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1346 %! [ -3.239592219e+23 -1.605436493e+17 -9.267951403e+08 -1.211210605e+03 1.375953130e-02 5.631729379e-05]
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1347 %! [ -2.946428547e+40 -1.407393871e+31 -7.760717570e+18 -6.908318646e+09 -2.241226812e-02 -5.412929349e-03]
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1348 %! [ -8.028450851e+58 -3.720929322e+46 -2.055758716e+30 -1.510304919e+18 4.978797221e-05 -7.048420407e-04]
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1349 %! [ -2.739192285e+78 -1.235021944e+63 -6.964109188e+42 -4.528227272e+27 -4.190000150e-02 1.074782297e-02]
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1350 %! [-6.683079463e+186 -2.655955830e+156 -1.799713983e+116 -8.573226309e+85 -1.125692891e+18 -2.298385049e-02]];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1351 %!
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1352 %! n = [(0:5:20).';30;40;50;100];
6646031e2450 besselj.cc: Added tests to cover a broad range of arguments and orders.
Robert T. Short <octave@phaselockedsystems.com>
parents: 14360
diff changeset
1353 %! z = [1,2,5,10,50,100];
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1354 %! j = sqrt ((pi/2)./z) .* besselj ((n+1/2).', z.').';
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1355 %! y = sqrt ((pi/2)./z) .* bessely ((n+1/2).', z.').';
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1356 %! assert (j, jt, -1e-9);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14450
diff changeset
1357 %! assert (y, yt, -1e-9);
8278
ab0674a8b345 fix scaling factor for negative alpha in zbesi,cbesi
Brian Gough <bjg@gnu.org>
parents: 7919
diff changeset
1358 */
29958
32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 29359
diff changeset
1359
32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 29359
diff changeset
1360 OCTAVE_NAMESPACE_END