Mercurial > octave

annotate libinterp/corefcn/fft.cc @ 33617:ec2635a02328 bytecode-interpreter tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: Merge default to bytecode-interpreter.
author Markus Mützel <markus.muetzel@gmx.de>
date Tue, 21 May 2024 18:29:03 +0200
parents 273cf7c2302a
children
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 //
32632
2e484f9f1f18 maint: update Octave Project Developers copyright for the new year
John W. Eaton <jwe@octave.org>
parents: 32368
diff changeset
3 // Copyright (C) 1996-2024 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 ////////////////////////////////////////////////////////////////////////
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
25
21724
aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef"
John W. Eaton <jwe@octave.org>
parents: 21660
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"
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
28 #endif
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
29
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
30 #include "lo-mappers.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
31
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14854
diff changeset
32 #include "defun.h"
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
33 #include "error.h"
21100
e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme.
Rik <rik@octave.org>
parents: 20940
diff changeset
34 #include "errwarn.h"
20940
48b2ad5ee801 maint: Rename oct-obj.[cc|h] to ovl.[cc|h] for clarity.
Rik <rik@octave.org>
parents: 20918
diff changeset
35 #include "ovl.h"
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
36 #include "utils.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
37
31605
e88a07dec498 maint: Use macros to begin/end C++ namespaces.
Rik <rik@octave.org>
parents: 30888
diff changeset
38 OCTAVE_BEGIN_NAMESPACE(octave)
29958
32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 29688
diff changeset
39
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
40 static octave_value
23433
c9fab0bc983e maint: Use convention 'int& x' for naming references.
Rik <rik@octave.org>
parents: 23220
diff changeset
41 do_fft (const octave_value_list& args, const char *fcn, int type)
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
42 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
43 int nargin = args.length ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
44
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
45 if (nargin < 1 || nargin > 3)
20802
8bb38ba1bad6 eliminate return statements after calls to print_usage
John W. Eaton <jwe@octave.org>
parents: 20555
diff changeset
46 print_usage ();
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
47
20918
6f0bd96f93c0 maint: Use new C++ archetype in more files.
Rik <rik@octave.org>
parents: 20853
diff changeset
48 octave_value retval;
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
49 octave_value arg = args(0);
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
50 octave_idx_type n_points = -1;
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
51 dim_vector dims = arg.dims ();
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
52 int ndims = dims.ndims ();
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
53 int dim = -1;
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11553
diff changeset
54
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
55 if (nargin > 1)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
56 {
23577
80c42f4cca13 maint: Deprecate is_empty and replace with isempty.
Rik <rik@octave.org>
parents: 23433
diff changeset
57 if (! args(1).isempty ())
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
58 {
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
59 double dval = args(1).double_value ();
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
60 if (math::isnan (dval))
11553
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents: 11523
diff changeset
61 error ("%s: number of points (N) cannot be NaN", fcn);
20918
6f0bd96f93c0 maint: Use new C++ archetype in more files.
Rik <rik@octave.org>
parents: 20853
diff changeset
62
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
63 n_points = math::nint_big (dval);
20918
6f0bd96f93c0 maint: Use new C++ archetype in more files.
Rik <rik@octave.org>
parents: 20853
diff changeset
64 if (n_points < 0)
6f0bd96f93c0 maint: Use new C++ archetype in more files.
Rik <rik@octave.org>
parents: 20853
diff changeset
65 error ("%s: number of points (N) must be greater than zero", fcn);
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
66 }
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
67 }
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
68
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
69 if (nargin > 2)
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
70 {
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
71 double dval = args(2).double_value ();
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
72 if (math::isnan (dval))
11553
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents: 11523
diff changeset
73 error ("%s: DIM cannot be NaN", fcn);
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
74 else if (dval < 1 || dval > ndims)
17787
175b392e91fe Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents: 17744
diff changeset
75 error ("%s: DIM must be a valid dimension along which to perform FFT",
175b392e91fe Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents: 17744
diff changeset
76 fcn);
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
77 else
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
78 // to be safe, cast it back to int since dim is an int
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
79 dim = math::nint (dval) - 1;
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
80 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
81
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
82 if (dim < 0)
4869
b92d59213e63 [project @ 2004-04-21 17:03:02 by jwe]
jwe
parents: 4782
diff changeset
83 {
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
84 dim = dims.first_non_singleton ();
4869
b92d59213e63 [project @ 2004-04-21 17:03:02 by jwe]
jwe
parents: 4782
diff changeset
85
b92d59213e63 [project @ 2004-04-21 17:03:02 by jwe]
jwe
parents: 4782
diff changeset
86 // And if the first argument is scalar?
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
87 if (dim == ndims)
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
88 dim = 1;
4869
b92d59213e63 [project @ 2004-04-21 17:03:02 by jwe]
jwe
parents: 4782
diff changeset
89 }
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
90
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
91 if (n_points < 0)
20218
b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors.
Rik <rik@octave.org>
parents: 19697
diff changeset
92 n_points = dims(dim);
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
93 else
20218
b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors.
Rik <rik@octave.org>
parents: 19697
diff changeset
94 dims(dim) = n_points;
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
95
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
96 if (n_points == 0 || dims.any_zero ())
7924
4976f66d469b miscellaneous cleanup
John W. Eaton <jwe@octave.org>
parents: 7814
diff changeset
97 {
4976f66d469b miscellaneous cleanup
John W. Eaton <jwe@octave.org>
parents: 7814
diff changeset
98 if (arg.is_single_type ())
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
99 return octave_value (FloatNDArray (dims));
7924
4976f66d469b miscellaneous cleanup
John W. Eaton <jwe@octave.org>
parents: 7814
diff changeset
100 else
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
101 return octave_value (NDArray (dims));
7924
4976f66d469b miscellaneous cleanup
John W. Eaton <jwe@octave.org>
parents: 7814
diff changeset
102 }
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
103
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
104 if (n_points == 1)
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
105 {
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
106 octave_value_list idx (ndims);
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
107 for (octave_idx_type i = 0; i < ndims; i++)
29961
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
108 idx(i) = idx_vector::colon;
7d6709900da7 eliminate octave:: namespace tags in DEFUN and DEFMETHOD and more
John W. Eaton <jwe@octave.org>
parents: 29958
diff changeset
109 idx(dim) = idx_vector (0);
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
110
28631
70cdf8de553d move non-member octave_value operator functions to octave namespace
John W. Eaton <jwe@octave.org>
parents: 27923
diff changeset
111 return arg.index_op (idx);
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
112 }
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
113
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7562
diff changeset
114 if (arg.is_single_type ())
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
115 {
23582
0cc2011d800e maint: Deprecate is_real_type and replace with isreal.
Rik <rik@octave.org>
parents: 23581
diff changeset
116 if (arg.isreal ())
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
117 {
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
118 FloatNDArray nda = arg.float_array_value ();
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
119
20555
f90c8372b7ba eliminate many more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20218
diff changeset
120 nda.resize (dims, 0.0);
f90c8372b7ba eliminate many more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20218
diff changeset
121 retval = (type != 0 ? nda.ifourier (dim) : nda.fourier (dim));
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
122 }
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7562
diff changeset
123 else
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
124 {
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
125 FloatComplexNDArray cnda = arg.float_complex_array_value ();
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
126
20555
f90c8372b7ba eliminate many more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20218
diff changeset
127 cnda.resize (dims, 0.0);
f90c8372b7ba eliminate many more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20218
diff changeset
128 retval = (type != 0 ? cnda.ifourier (dim) : cnda.fourier (dim));
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
129 }
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
130 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
131 else
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
132 {
23582
0cc2011d800e maint: Deprecate is_real_type and replace with isreal.
Rik <rik@octave.org>
parents: 23581
diff changeset
133 if (arg.isreal ())
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
134 {
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
135 NDArray nda = arg.array_value ();
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7562
diff changeset
136
20555
f90c8372b7ba eliminate many more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20218
diff changeset
137 nda.resize (dims, 0.0);
f90c8372b7ba eliminate many more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20218
diff changeset
138 retval = (type != 0 ? nda.ifourier (dim) : nda.fourier (dim));
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
139 }
23581
c3075ae020e1 maint: Deprecate is_complex_type and replace with iscomplex.
Rik <rik@octave.org>
parents: 23577
diff changeset
140 else if (arg.iscomplex ())
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
141 {
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
142 ComplexNDArray cnda = arg.complex_array_value ();
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7562
diff changeset
143
20555
f90c8372b7ba eliminate many more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20218
diff changeset
144 cnda.resize (dims, 0.0);
f90c8372b7ba eliminate many more simple uses of error_state
John W. Eaton <jwe@octave.org>
parents: 20218
diff changeset
145 retval = (type != 0 ? cnda.ifourier (dim) : cnda.fourier (dim));
10154
40dfc0c99116 DLD-FUNCTIONS/*.cc: untabify
John W. Eaton <jwe@octave.org>
parents: 9523
diff changeset
146 }
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7562
diff changeset
147 else
21100
e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme.
Rik <rik@octave.org>
parents: 20940
diff changeset
148 err_wrong_type_arg (fcn, arg);
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
149 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
150
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
151 return retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
152 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
153
6620
bf4bdc21dc8d [project @ 2007-05-14 17:35:46 by jwe]
jwe
parents: 5823
diff changeset
154 /*
25586
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
155 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
156 %! assert (fft ([]), [])
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
157 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
158 %! assert (fft (zeros (10,0)), zeros (10,0))
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
159 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
160 %! assert (fft (zeros (0,10)), zeros (0,10))
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
161 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
162 %! assert (fft (0), 0)
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
163 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
164 %! assert (fft (1), 1)
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
165 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
166 %! assert (fft (ones (2,2)), [2,2; 0,0])
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
167 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
168 %! assert (fft (eye (2,2)), [1,1; 1,-1])
6620
bf4bdc21dc8d [project @ 2007-05-14 17:35:46 by jwe]
jwe
parents: 5823
diff changeset
169
25586
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
170 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
171 %! assert (fft (single ([])), single ([]))
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
172 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
173 %! assert (fft (zeros (10,0,"single")), zeros (10,0,"single"))
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
174 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
175 %! assert (fft (zeros (0,10,"single")), zeros (0,10,"single"))
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
176 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
177 %! assert (fft (single (0)), single (0))
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
178 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
179 %! assert (fft (single (1)), single (1))
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
180 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
181 %! assert (fft (ones (2,2,"single")), single ([2,2; 0,0]))
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
182 %!testif HAVE_FFTW
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
183 %! assert (fft (eye (2,2,"single")), single ([1,1; 1,-1]))
6620
bf4bdc21dc8d [project @ 2007-05-14 17:35:46 by jwe]
jwe
parents: 5823
diff changeset
184
28930
2883b3d08b7e maint: Remove unnecessary parentheses bracketing %!error BIST tests.
Rik <rik@octave.org>
parents: 28631
diff changeset
185 %!error fft ()
6620
bf4bdc21dc8d [project @ 2007-05-14 17:35:46 by jwe]
jwe
parents: 5823
diff changeset
186 */
bf4bdc21dc8d [project @ 2007-05-14 17:35:46 by jwe]
jwe
parents: 5823
diff changeset
187
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
188
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14854
diff changeset
189 DEFUN (fft, args, ,
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
190 doc: /* -*- texinfo -*-
30888
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
191 @deftypefn {} {@var{y} =} fft (@var{x})
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
192 @deftypefnx {} {@var{y} =} fft (@var{x}, @var{n})
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
193 @deftypefnx {} {@var{y} =} fft (@var{x}, @var{n}, @var{dim})
24575
93714122ed59 doc: Fixes for signal, image, audio, and OOP chapters (bug #52844).
Rik <rik@octave.org>
parents: 23219
diff changeset
194 Compute the discrete Fourier transform of @var{x} using
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
195 a Fast Fourier Transform (FFT) algorithm.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
196
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
197 The FFT is calculated along the first non-singleton dimension of the
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
198 array. Thus if @var{x} is a matrix, @code{fft (@var{x})} computes the
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
199 FFT for each column of @var{x}.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
200
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
201 If called with two arguments, @var{n} is expected to be an integer
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
202 specifying the number of elements of @var{x} to use, or an empty
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
203 matrix to specify that its value should be ignored. If @var{n} is
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
204 larger than the dimension along which the FFT is calculated, then
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
205 @var{x} is resized and padded with zeros. Otherwise, if @var{n} is
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
206 smaller than the dimension along which the FFT is calculated, then
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
207 @var{x} is truncated.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
208
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
209 If called with three arguments, @var{dim} is an integer specifying the
24575
93714122ed59 doc: Fixes for signal, image, audio, and OOP chapters (bug #52844).
Rik <rik@octave.org>
parents: 23219
diff changeset
210 dimension of the matrix along which the FFT is performed.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
211 @seealso{ifft, fft2, fftn, fftw}
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
212 @end deftypefn */)
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
213 {
4782
3dfdf6f36854 [project @ 2004-02-17 16:23:22 by jwe]
jwe
parents: 4773
diff changeset
214 return do_fft (args, "fft", 0);
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
215 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
216
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
217
15039
e753177cde93 maint: Move non-dynamically linked functions from DLD-FUNCTIONS/ to corefcn/ directory
Rik <rik@octave.org>
parents: 14854
diff changeset
218 DEFUN (ifft, args, ,
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
219 doc: /* -*- texinfo -*-
30888
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
220 @deftypefn {} {@var{x} =} ifft (@var{y})
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
221 @deftypefnx {} {@var{x} =} ifft (@var{y}, @var{n})
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
222 @deftypefnx {} {@var{x} =} ifft (@var{y}, @var{n}, @var{dim})
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
223 Compute the inverse discrete Fourier transform of @var{y}
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
224 using a Fast Fourier Transform (FFT) algorithm.
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
225
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
226 The inverse FFT is calculated along the first non-singleton dimension
30888
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
227 of the array. Thus if @var{y} is a matrix, @code{ifft (@var{y})} computes
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
228 the inverse FFT for each column of @var{y}.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
229
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
230 If called with two arguments, @var{n} is expected to be an integer
30888
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
231 specifying the number of elements of @var{y} to use, or an empty
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
232 matrix to specify that its value should be ignored. If @var{n} is
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
233 larger than the dimension along which the inverse FFT is calculated, then
30888
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
234 @var{y} is resized and padded with zeros. Otherwise, if @var{n} is
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
235 smaller than the dimension along which the inverse FFT is calculated,
30888
32d2b6604a9f doc: Ensure documentation lists output argument when it exists for functions in libinterp/
Rik <rik@octave.org>
parents: 30564
diff changeset
236 then @var{y} is truncated.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
237
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
238 If called with three arguments, @var{dim} is an integer specifying the
24575
93714122ed59 doc: Fixes for signal, image, audio, and OOP chapters (bug #52844).
Rik <rik@octave.org>
parents: 23219
diff changeset
239 dimension of the matrix along which the inverse FFT is performed.
21966
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
240 @seealso{fft, ifft2, ifftn, fftw}
112b20240c87 move docstrings in C++ files out of C strings and into comments
John W. Eaton <jwe@octave.org>
parents: 21782
diff changeset
241 @end deftypefn */)
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
242 {
4782
3dfdf6f36854 [project @ 2004-02-17 16:23:22 by jwe]
jwe
parents: 4773
diff changeset
243 return do_fft (args, "ifft", 1);
4773
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
244 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe]
jwe
parents: 4233
diff changeset
245
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
246 /*
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
247 ## Author: David Billinghurst (David.Billinghurst@riotinto.com.au)
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
248 ## Comalco Research and Technology
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
249 ## 02 May 2000
25586
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
250 %!testif HAVE_FFTW
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
251 %! N = 64;
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
252 %! n = 4;
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
253 %! t = 2*pi*(0:1:N-1)/N;
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
254 %! s = cos (n*t);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
255 %! S = fft (s);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
256 %!
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
257 %! answer = zeros (size (t));
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
258 %! answer(n+1) = N/2;
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
259 %! answer(N-n+1) = N/2;
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
260 %!
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
261 %! assert (S, answer, 4*N*eps);
7562
c827f5673321 move tests to individual source files
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
262
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
263 ## Author: David Billinghurst (David.Billinghurst@riotinto.com.au)
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
264 ## Comalco Research and Technology
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
265 ## 02 May 2000
25586
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
266 %!testif HAVE_FFTW
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
267 %! N = 64;
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
268 %! n = 7;
7562
c827f5673321 move tests to individual source files
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
269 %! t = 2*pi*(0:1:N-1)/N;
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
270 %! s = cos (n*t);
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11553
diff changeset
271 %!
14854
5ae9f0f77635 maint: Use Octave coding conventions for coddling parenthis is DLD-FUNCTIONS directory
Rik <octave@nomad.inbox5.com>
parents: 14501
diff changeset
272 %! S = zeros (size (t));
7562
c827f5673321 move tests to individual source files
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
273 %! S(n+1) = N/2;
c827f5673321 move tests to individual source files
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
274 %! S(N-n+1) = N/2;
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11553
diff changeset
275 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
276 %! assert (ifft (S), s, 4*N*eps);
7562
c827f5673321 move tests to individual source files
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
277
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
278 ## Author: David Billinghurst (David.Billinghurst@riotinto.com.au)
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
279 ## Comalco Research and Technology
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
280 ## 02 May 2000
25586
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
281 %!testif HAVE_FFTW
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
282 %! N = 64;
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
283 %! n = 4;
7814
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
284 %! t = single (2*pi*(0:1:N-1)/N);
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
285 %! s = cos (n*t);
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
286 %! S = fft (s);
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11553
diff changeset
287 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
288 %! answer = zeros (size (t), "single");
7814
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
289 %! answer(n+1) = N/2;
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
290 %! answer(N-n+1) = N/2;
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11553
diff changeset
291 %!
31607
aac27ad79be6 maint: Re-indent code after switch to using namespace macros.
Rik <rik@octave.org>
parents: 31605
diff changeset
292 %! assert (S, answer, 4*N* eps ("single"));
7562
c827f5673321 move tests to individual source files
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
293
22662
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
294 ## Author: David Billinghurst (David.Billinghurst@riotinto.com.au)
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
295 ## Comalco Research and Technology
655157b34a9f Fix segfault when calling fft on last dimension (bug #49026).
Rik <rik@octave.org>
parents: 22407
diff changeset
296 ## 02 May 2000
25586
59f4fda5c548 remove fftpack sources
John W. Eaton <jwe@octave.org>
parents: 25054
diff changeset
297 %!testif HAVE_FFTW
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
298 %! N = 64;
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
299 %! n = 7;
7814
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
300 %! t = 2*pi*(0:1:N-1)/N;
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
301 %! s = cos (n*t);
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11553
diff changeset
302 %!
14501
60e5cf354d80 Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
303 %! S = zeros (size (t), "single");
7814
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
304 %! S(n+1) = N/2;
87865ed7405f Second set of single precision test code and fix of resulting bugs
David Bateman <dbateman@free.fr>
parents: 7789
diff changeset
305 %! S(N-n+1) = N/2;
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11553
diff changeset
306 %!
31607
aac27ad79be6 maint: Re-indent code after switch to using namespace macros.
Rik <rik@octave.org>
parents: 31605
diff changeset
307 %! assert (ifft (S), s, 4*N* eps ("single"));
32365
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
308
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
309 %!testif HAVE_FFTW <*64729>
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
310 %! x = rand (5, 5, 5, 5);
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
311 %! old_planner = fftw ('planner', 'measure');
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
312 %! unwind_protect
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
313 %! y = fft (x);
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
314 %! y = fft (x); # second invocation might crash Octave
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
315 %! assert (size (y), [5, 5, 5, 5]);
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
316 %! unwind_protect_cleanup
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
317 %! fftw ('planner', old_planner);
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
318 %! end_unwind_protect
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
319
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
320 %!testif HAVE_FFTW <*64729>
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
321 %! x = rand (5, 5, 5, 5, 'single');
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
322 %! old_planner = fftw ('planner', 'measure');
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
323 %! unwind_protect
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
324 %! y = fft (x);
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
325 %! y = fft (x); # second invocation might crash Octave
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
326 %! assert (size (y), [5, 5, 5, 5]);
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
327 %! unwind_protect_cleanup
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
328 %! fftw ('planner', old_planner);
65e27afd86a2 fft: Avoid segmentation fault with ND-arrays (bug #64729).
Markus Mützel <markus.muetzel@gmx.de>
parents: 31706
diff changeset
329 %! end_unwind_protect
32368
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
330
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
331 %!testif HAVE_FFTW <*64733>
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
332 %! old_planner = fftw ('planner', 'measure');
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
333 %! unwind_protect
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
334 %! assert (ifft ([2, 4, 6, 8]), [5, -1-1i, -1, -1+1i]);
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
335 %! unwind_protect_cleanup
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
336 %! fftw ('planner', old_planner);
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
337 %! end_unwind_protect
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
338
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
339 %!testif HAVE_FFTW <*64733>
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
340 %! old_planner = fftw ('planner', 'measure');
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
341 %! unwind_protect
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
342 %! assert (ifft (single ([2, 4, 6, 8])), single ([5, -1-1i, -1, -1+1i]));
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
343 %! unwind_protect_cleanup
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
344 %! fftw ('planner', old_planner);
903c9100178b fft: Avoid overwriting input for inplace operations (bug #64733).
Markus Mützel <markus.muetzel@gmx.de>
parents: 32365
diff changeset
345 %! end_unwind_protect
7562
c827f5673321 move tests to individual source files
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
346 */
29958
32c3a5805893 move DEFUN and DEFMETHOD functions inside octave namespace
John W. Eaton <jwe@octave.org>
parents: 29688
diff changeset
347
31605
e88a07dec498 maint: Use macros to begin/end C++ namespaces.
Rik <rik@octave.org>
parents: 30888
diff changeset
348 OCTAVE_END_NAMESPACE(octave)