Mercurial > octave-libgccjit

annotate src/DLD-FUNCTIONS/besselj.cc @ 7919:9d080df0c843

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
new NDArray constructor for ArrayN<octave_idx_type>
author David Bateman <dbateman@free.fr>
date Mon, 30 Jun 2008 15:51:31 +0200
parents 82be108cc558
children ab0674a8b345
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
1 /*
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
2
7017
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe]
jwe
parents: 7016
diff changeset
3 Copyright (C) 1997, 1998, 1999, 2000, 2001, 2002, 2004, 2005, 2006,
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe]
jwe
parents: 7016
diff changeset
4 2007 John W. Eaton
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
5
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
6 This file is part of Octave.
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
7
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
8 Octave is free software; you can redistribute it and/or modify it
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
9 under the terms of the GNU General Public License as published by the
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5823
diff changeset
10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5823
diff changeset
11 option) any later version.
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
12
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
13 Octave is distributed in the hope that it will be useful, but WITHOUT
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
16 for more details.
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
17
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
18 You should have received a copy of the GNU General Public License
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5823
diff changeset
19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5823
diff changeset
20 <http://www.gnu.org/licenses/>.
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
21
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
22 */
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
23
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
24 #ifdef HAVE_CONFIG_H
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
25 #include <config.h>
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
26 #endif
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
27
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
28 #include "lo-specfun.h"
4153
6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe]
jwe
parents: 3810
diff changeset
29 #include "quit.h"
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
30
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
31 #include "defun-dld.h"
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
32 #include "error.h"
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
33 #include "gripes.h"
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
34 #include "oct-obj.h"
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
35 #include "utils.h"
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
36
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
37 enum bessel_type
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
38 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
39 BESSEL_J,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
40 BESSEL_Y,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
41 BESSEL_I,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
42 BESSEL_K,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
43 BESSEL_H1,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
44 BESSEL_H2
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
45 };
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
46
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
47 #define DO_BESSEL(type, alpha, x, scaled, ierr, result) \
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
48 do \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
49 { \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
50 switch (type) \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
51 { \
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
52 case BESSEL_J: \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
53 result = besselj (alpha, x, scaled, ierr); \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
54 break; \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
55 \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
56 case BESSEL_Y: \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
57 result = bessely (alpha, x, scaled, ierr); \
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
58 break; \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
59 \
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
60 case BESSEL_I: \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
61 result = besseli (alpha, x, scaled, ierr); \
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
62 break; \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
63 \
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
64 case BESSEL_K: \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
65 result = besselk (alpha, x, scaled, ierr); \
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
66 break; \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
67 \
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
68 case BESSEL_H1: \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
69 result = besselh1 (alpha, x, scaled, ierr); \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
70 break; \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
71 \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
72 case BESSEL_H2: \
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
73 result = besselh2 (alpha, x, scaled, ierr); \
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
74 break; \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
75 \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
76 default: \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
77 break; \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
78 } \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
79 } \
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
80 while (0)
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
81
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
82 static void
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
83 gripe_bessel_arg (const char *fn, const char *arg)
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
84 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
85 error ("%s: expecting scalar or matrix as %s argument", fn, arg);
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
86 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
87
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
88 octave_value_list
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
89 do_bessel (enum bessel_type type, const char *fn,
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
90 const octave_value_list& args, int nargout)
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
91 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
92 octave_value_list retval;
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
93
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
94 int nargin = args.length ();
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
95
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
96 if (nargin == 2 || nargin == 3)
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
97 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
98 bool scaled = (nargin == 3);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
99
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
100 octave_value alpha_arg = args(0);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
101 octave_value x_arg = args(1);
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
102
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
103 if (alpha_arg.is_single_type () || x_arg.is_single_type ())
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
104 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
105 if (alpha_arg.is_scalar_type ())
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
106 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
107 float alpha = args(0).float_value ();
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
108
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
109 if (! error_state)
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
110 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
111 if (x_arg.is_scalar_type ())
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
112 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
113 FloatComplex x = x_arg.float_complex_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
114
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
115 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
116 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
117 octave_idx_type ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
118 octave_value result;
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
119
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
120 DO_BESSEL (type, alpha, x, scaled, ierr, result);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
121
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
122 if (nargout > 1)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
123 retval(1) = static_cast<float> (ierr);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
124
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
125 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
126 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
127 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
128 gripe_bessel_arg (fn, "second");
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
129 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
130 else
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
131 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
132 FloatComplexNDArray x = x_arg.float_complex_array_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
133
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
134 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
135 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
136 ArrayN<octave_idx_type> ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
137 octave_value result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
138
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
139 DO_BESSEL (type, alpha, x, scaled, ierr, result);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
140
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
141 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
142 retval(1) = NDArray (ierr);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
143
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
144 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
145 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
146 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
147 gripe_bessel_arg (fn, "second");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
148 }
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
149 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
150 else
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
151 gripe_bessel_arg (fn, "first");
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
152 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
153 else
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
154 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
155 dim_vector dv0 = args(0).dims ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
156 dim_vector dv1 = args(1).dims ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
157
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
158 bool args0_is_row_vector = (dv0 (1) == dv0.numel ());
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
159 bool args1_is_col_vector = (dv1 (0) == dv1.numel ());
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
160
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
161 if (args0_is_row_vector && args1_is_col_vector)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
162 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
163 FloatRowVector ralpha = args(0).float_row_vector_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
164
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
165 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
166 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
167 FloatComplexColumnVector cx =
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
168 x_arg.float_complex_column_vector_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
169
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
170 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
171 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
172 Array2<octave_idx_type> ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
173 octave_value result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
174
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
175 DO_BESSEL (type, ralpha, cx, scaled, ierr, result);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
176
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
177 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
178 retval(1) = NDArray (ierr);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
179
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
180 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
181 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
182 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
183 gripe_bessel_arg (fn, "second");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
184 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
185 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
186 gripe_bessel_arg (fn, "first");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
187 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
188 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
189 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
190 FloatNDArray alpha = args(0).float_array_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
191
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
192 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
193 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
194 if (x_arg.is_scalar_type ())
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
195 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
196 FloatComplex x = x_arg.float_complex_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
197
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
198 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
199 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
200 ArrayN<octave_idx_type> ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
201 octave_value result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
202
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
203 DO_BESSEL (type, alpha, x, scaled, ierr, result);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
204
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
205 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
206 retval(1) = NDArray (ierr);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
207
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
208 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
209 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
210 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
211 gripe_bessel_arg (fn, "second");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
212 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
213 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
214 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
215 FloatComplexNDArray x = x_arg.float_complex_array_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
216
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
217 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
218 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
219 ArrayN<octave_idx_type> ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
220 octave_value result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
221
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
222 DO_BESSEL (type, alpha, x, scaled, ierr, result);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
223
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
224 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
225 retval(1) = NDArray (ierr);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
226
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
227 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
228 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
229 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
230 gripe_bessel_arg (fn, "second");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
231 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
232 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
233 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
234 gripe_bessel_arg (fn, "first");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
235 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
236 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
237 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
238 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
239 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
240 if (alpha_arg.is_scalar_type ())
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
241 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
242 double alpha = args(0).double_value ();
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
243
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
244 if (! error_state)
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
245 {
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
246 if (x_arg.is_scalar_type ())
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
247 {
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
248 Complex x = x_arg.complex_value ();
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
249
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
250 if (! error_state)
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
251 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
252 octave_idx_type ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
253 octave_value result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
254
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
255 DO_BESSEL (type, alpha, x, scaled, ierr, result);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
256
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
257 if (nargout > 1)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
258 retval(1) = static_cast<double> (ierr);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
259
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
260 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
261 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
262 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
263 gripe_bessel_arg (fn, "second");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
264 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
265 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
266 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
267 ComplexNDArray x = x_arg.complex_array_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
268
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
269 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
270 {
5275
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe]
jwe
parents: 4844
diff changeset
271 ArrayN<octave_idx_type> ierr;
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
272 octave_value result;
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
273
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
274 DO_BESSEL (type, alpha, x, scaled, ierr, result);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
275
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
276 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
277 retval(1) = NDArray (ierr);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
278
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
279 retval(0) = result;
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
280 }
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
281 else
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
282 gripe_bessel_arg (fn, "second");
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
283 }
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
284 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
285 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
286 gripe_bessel_arg (fn, "first");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
287 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
288 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
289 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
290 dim_vector dv0 = args(0).dims ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
291 dim_vector dv1 = args(1).dims ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
292
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
293 bool args0_is_row_vector = (dv0 (1) == dv0.numel ());
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
294 bool args1_is_col_vector = (dv1 (0) == dv1.numel ());
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
295
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
296 if (args0_is_row_vector && args1_is_col_vector)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
297 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
298 RowVector ralpha = args(0).row_vector_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
299
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
300 if (! error_state)
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
301 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
302 ComplexColumnVector cx =
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
303 x_arg.complex_column_vector_value ();
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
304
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
305 if (! error_state)
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
306 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
307 Array2<octave_idx_type> ierr;
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
308 octave_value result;
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
309
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
310 DO_BESSEL (type, ralpha, cx, scaled, ierr, result);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
311
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
312 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
313 retval(1) = NDArray (ierr);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
314
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
315 retval(0) = result;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
316 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
317 else
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
318 gripe_bessel_arg (fn, "second");
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
319 }
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
320 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
321 gripe_bessel_arg (fn, "first");
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
322 }
4844
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe]
jwe
parents: 4153
diff changeset
323 else
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
324 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
325 NDArray alpha = args(0).array_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
326
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
327 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
328 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
329 if (x_arg.is_scalar_type ())
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
330 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
331 Complex x = x_arg.complex_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
332
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
333 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
334 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
335 ArrayN<octave_idx_type> ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
336 octave_value result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
337
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
338 DO_BESSEL (type, alpha, x, scaled, ierr, result);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
339
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
340 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
341 retval(1) = NDArray (ierr);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
342
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
343 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
344 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
345 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
346 gripe_bessel_arg (fn, "second");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
347 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
348 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
349 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
350 ComplexNDArray x = x_arg.complex_array_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
351
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
352 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
353 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
354 ArrayN<octave_idx_type> ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
355 octave_value result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
356
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
357 DO_BESSEL (type, alpha, x, scaled, ierr, result);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
358
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
359 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
360 retval(1) = NDArray (ierr);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
361
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
362 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
363 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
364 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
365 gripe_bessel_arg (fn, "second");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
366 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
367 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
368 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
369 gripe_bessel_arg (fn, "first");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
370 }
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
371 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
372 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
373 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
374 else
5823
080c08b192d8 [project @ 2006-05-19 05:32:17 by jwe]
jwe
parents: 5760
diff changeset
375 print_usage ();
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
376
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
377 return retval;
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
378 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
379
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
380 DEFUN_DLD (besselj, args, nargout,
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
381 "-*- texinfo -*-\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
382 @deftypefn {Loadable Function} {[@var{j}, @var{ierr}] =} besselj (@var{alpha}, @var{x}, @var{opt})\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
383 @deftypefnx {Loadable Function} {[@var{y}, @var{ierr}] =} bessely (@var{alpha}, @var{x}, @var{opt})\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
384 @deftypefnx {Loadable Function} {[@var{i}, @var{ierr}] =} besseli (@var{alpha}, @var{x}, @var{opt})\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
385 @deftypefnx {Loadable Function} {[@var{k}, @var{ierr}] =} besselk (@var{alpha}, @var{x}, @var{opt})\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
386 @deftypefnx {Loadable Function} {[@var{h}, @var{ierr}] =} besselh (@var{alpha}, @var{k}, @var{x}, @var{opt})\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
387 Compute Bessel or Hankel functions of various kinds:\n\
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
388 \n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
389 @table @code\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
390 @item besselj\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
391 Bessel functions of the first kind.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
392 @item bessely\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
393 Bessel functions of the second kind.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
394 @item besseli\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
395 Modified Bessel functions of the first kind.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
396 @item besselk\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
397 Modified Bessel functions of the second kind.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
398 @item besselh\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
399 Compute Hankel functions of the first (@var{k} = 1) or second (@var{k}\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
400 = 2) kind.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
401 @end table\n\
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
402 \n\
3568
58549b931ad5 [project @ 2000-02-03 10:21:54 by jwe]
jwe
parents: 3549
diff changeset
403 If the argument @var{opt} is supplied, the result is scaled by the\n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
404 @code{exp (-I*@var{x})} for @var{k} = 1 or @code{exp (I*@var{x})} for\n\
3499
3e3e14ad5149 [project @ 2000-01-31 05:18:07 by jwe]
jwe
parents: 3459
diff changeset
405 @var{k} = 2.\n\
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
406 \n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
407 If @var{alpha} is a scalar, the result is the same size as @var{x}.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
408 If @var{x} is a scalar, the result is the same size as @var{alpha}.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
409 If @var{alpha} is a row vector and @var{x} is a column vector, the\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
410 result is a matrix with @code{length (@var{x})} rows and\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
411 @code{length (@var{alpha})} columns. Otherwise, @var{alpha} and\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
412 @var{x} must conform and the result will be the same size.\n\
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
413 \n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
414 The value of @var{alpha} must be real. The value of @var{x} may be\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
415 complex.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
416 \n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
417 If requested, @var{ierr} contains the following status information\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
418 and is the same size as the result.\n\
3548
ab7fa5a8f23f [project @ 2000-02-03 01:17:15 by jwe]
jwe
parents: 3499
diff changeset
419 \n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
420 @enumerate 0\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
421 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
422 Normal return.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
423 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
424 Input error, return @code{NaN}.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
425 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
426 Overflow, return @code{Inf}.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
427 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
428 Loss of significance by argument reduction results in less than\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
429 half of machine accuracy.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
430 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
431 Complete loss of significance by argument reduction, return @code{NaN}.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
432 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
433 Error---no computation, algorithm termination condition not met,\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
434 return @code{NaN}.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
435 @end enumerate\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
436 @end deftypefn")
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
437 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
438 return do_bessel (BESSEL_J, "besselj", args, nargout);
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
439 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
440
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
441 DEFUN_DLD (bessely, args, nargout,
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
442 "-*- texinfo -*-\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
443 @deftypefn {Loadable Function} {[@var{y}, @var{ierr}] =} bessely (@var{alpha}, @var{x}, @var{opt})\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
444 See besselj.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
445 @end deftypefn")
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
446 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
447 return do_bessel (BESSEL_Y, "bessely", args, nargout);
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
448 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
449
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
450 DEFUN_DLD (besseli, args, nargout,
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
451 "-*- texinfo -*-\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
452 @deftypefn {Loadable Function} {[@var{i}, @var{ierr}] =} besseli (@var{alpha}, @var{x}, @var{opt})\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
453 See besselj.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
454 @end deftypefn")
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
455 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
456 return do_bessel (BESSEL_I, "besseli", args, nargout);
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
457 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
458
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
459 DEFUN_DLD (besselk, args, nargout,
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
460 "-*- texinfo -*-\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
461 @deftypefn {Loadable Function} {[@var{k}, @var{ierr}] =} besselk (@var{alpha}, @var{x}, @var{opt})\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
462 See besselj.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
463 @end deftypefn")
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
464 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
465 return do_bessel (BESSEL_K, "besselk", args, nargout);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
466 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
467
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
468 DEFUN_DLD (besselh, args, nargout,
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
469 "-*- texinfo -*-\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
470 @deftypefn {Loadable Function} {[@var{h}, @var{ierr}] =} besselh (@var{alpha}, @var{k}, @var{x}, @var{opt})\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
471 See besselj.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
472 @end deftypefn")
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
473 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
474 octave_value_list retval;
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
475
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
476 int nargin = args.length ();
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
477
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
478 if (nargin == 2)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
479 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
480 retval = do_bessel (BESSEL_H1, "besselh", args, nargout);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
481 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
482 else if (nargin == 3 || nargin == 4)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
483 {
5275
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe]
jwe
parents: 4844
diff changeset
484 octave_idx_type kind = args(1).int_value ();
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
485
3810
f19f00339363 [project @ 2001-03-29 21:03:01 by jwe]
jwe
parents: 3568
diff changeset
486 if (! error_state)
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
487 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
488 octave_value_list tmp_args;
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
489
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
490 if (nargin == 4)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
491 tmp_args(2) = args(3);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
492
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
493 tmp_args(1) = args(2);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
494 tmp_args(0) = args(0);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
495
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
496 if (kind == 1)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
497 retval = do_bessel (BESSEL_H1, "besselh", tmp_args, nargout);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
498 else if (kind == 2)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
499 retval = do_bessel (BESSEL_H2, "besselh", tmp_args, nargout);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
500 else
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
501 error ("besselh: expecting K = 1 or 2");
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
502 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
503 else
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
504 error ("besselh: invalid value of K");
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
505 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
506 else
5823
080c08b192d8 [project @ 2006-05-19 05:32:17 by jwe]
jwe
parents: 5760
diff changeset
507 print_usage ();
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
508
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
509 return retval;
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
510 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
511
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
512 DEFUN_DLD (airy, args, nargout,
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
513 "-*- texinfo -*-\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
514 @deftypefn {Loadable Function} {[@var{a}, @var{ierr}] =} airy (@var{k}, @var{z}, @var{opt})\n\
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
515 Compute Airy functions of the first and second kind, and their\n\
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
516 derivatives.\n\
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
517 \n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
518 @example\n\
7031
120f3135952f [project @ 2007-10-15 15:30:03 by jwe]
jwe
parents: 7017
diff changeset
519 K Function Scale factor (if 'opt' is supplied)\n\
120f3135952f [project @ 2007-10-15 15:30:03 by jwe]
jwe
parents: 7017
diff changeset
520 --- -------- ---------------------------------------\n\
120f3135952f [project @ 2007-10-15 15:30:03 by jwe]
jwe
parents: 7017
diff changeset
521 0 Ai (Z) exp ((2/3) * Z * sqrt (Z))\n\
120f3135952f [project @ 2007-10-15 15:30:03 by jwe]
jwe
parents: 7017
diff changeset
522 1 dAi(Z)/dZ exp ((2/3) * Z * sqrt (Z))\n\
120f3135952f [project @ 2007-10-15 15:30:03 by jwe]
jwe
parents: 7017
diff changeset
523 2 Bi (Z) exp (-abs (real ((2/3) * Z *sqrt (Z))))\n\
120f3135952f [project @ 2007-10-15 15:30:03 by jwe]
jwe
parents: 7017
diff changeset
524 3 dBi(Z)/dZ exp (-abs (real ((2/3) * Z *sqrt (Z))))\n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
525 @end example\n\
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
526 \n\
3549
03025e79d8b9 [project @ 2000-02-03 01:28:24 by jwe]
jwe
parents: 3548
diff changeset
527 The function call @code{airy (@var{z})} is equivalent to\n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
528 @code{airy (0, @var{z})}.\n\
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
529 \n\
3549
03025e79d8b9 [project @ 2000-02-03 01:28:24 by jwe]
jwe
parents: 3548
diff changeset
530 The result is the same size as @var{z}.\n\
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
531 \n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
532 If requested, @var{ierr} contains the following status information and\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
533 is the same size as the result.\n\
3548
ab7fa5a8f23f [project @ 2000-02-03 01:17:15 by jwe]
jwe
parents: 3499
diff changeset
534 \n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
535 @enumerate 0\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
536 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
537 Normal return.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
538 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
539 Input error, return @code{NaN}.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
540 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
541 Overflow, return @code{Inf}.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
542 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
543 Loss of significance by argument reduction results in less than half\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
544 of machine accuracy.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
545 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
546 Complete loss of significance by argument reduction, return @code{NaN}.\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
547 @item\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
548 Error---no computation, algorithm termination condition not met,\n\
5448
ebe5d7d15522 [project @ 2005-09-14 18:55:04 by jwe]
jwe
parents: 5307
diff changeset
549 return @code{NaN}.\n\
3459
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
550 @end enumerate\n\
8e36c45e3a61 [project @ 2000-01-19 10:26:18 by jwe]
jwe
parents: 3325
diff changeset
551 @end deftypefn")
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
552 {
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
553 octave_value_list retval;
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
554
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
555 int nargin = args.length ();
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
556
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
557 if (nargin > 0 && nargin < 4)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
558 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
559 bool scale = (nargin == 3);
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
560
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
561 int kind = 0;
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
562
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
563 if (nargin > 1)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
564 {
5760
8d7162924bd3 [project @ 2006-04-14 04:01:37 by jwe]
jwe
parents: 5448
diff changeset
565 kind = args(0).int_value ();
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
566
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
567 if (! error_state)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
568 {
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
569 if (kind < 0 || kind > 3)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
570 error ("airy: expecting K = 0, 1, 2, or 3");
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
571 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
572 else
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
573 error ("airy: expecting integer value for K");
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
574 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
575
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
576 if (! error_state)
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
577 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
578 int idx = nargin == 1 ? 0 : 1;
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
579
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
580 if (args (idx).is_single_type ())
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
581 {
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
582 FloatComplexNDArray z = args(idx).float_complex_array_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
583
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
584 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
585 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
586 ArrayN<octave_idx_type> ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
587 octave_value result;
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
588
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
589 if (kind > 1)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
590 result = biry (z, kind == 3, scale, ierr);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
591 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
592 result = airy (z, kind == 1, scale, ierr);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
593
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
594 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
595 retval(1) = NDArray (ierr);
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
596
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
597 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
598 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
599 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
600 error ("airy: expecting complex matrix for Z");
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
601 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
602 else
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
603 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
604 ComplexNDArray z = args(idx).complex_array_value ();
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
605
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
606 if (! error_state)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
607 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
608 ArrayN<octave_idx_type> ierr;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
609 octave_value result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
610
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
611 if (kind > 1)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
612 result = biry (z, kind == 3, scale, ierr);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
613 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
614 result = airy (z, kind == 1, scale, ierr);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
615
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
616 if (nargout > 1)
7919
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type>
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
617 retval(1) = NDArray (ierr);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
618
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
619 retval(0) = result;
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
620 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
621 else
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
622 error ("airy: expecting complex matrix for Z");
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7031
diff changeset
623 }
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
624 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
625 }
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
626 else
5823
080c08b192d8 [project @ 2006-05-19 05:32:17 by jwe]
jwe
parents: 5760
diff changeset
627 print_usage ();
3220
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
628
3deb1105fbc1 [project @ 1998-11-19 00:06:30 by jwe]
jwe
parents: 3187
diff changeset
629 return retval;
3155
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
630 }
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
631
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
632 /*
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
633 ;;; Local Variables: ***
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
634 ;;; mode: C++ ***
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
635 ;;; End: ***
1016520a9d38 [project @ 1998-02-18 21:51:50 by jwe]
jwe
parents:
diff changeset
636 */