octave: liboctave/array/CNDArray.cc annotate

annotate liboctave/array/CNDArray.cc @ 23443:3f1bf237908b

maint: Eliminate <cfloat.h> header from liboctave files. * build-aux/mk-opts.pl: Remove "#include <cfloat>" from script which creates *-opts.h * Array-util.cc, Array-util.h, Array.cc, Array.h, CColVector.h, CDiagMatrix.h, CMatrix.cc, CMatrix.h, CNDArray.cc, CNDArray.h, CRowVector.h, CSparse.cc, CSparse.h, MSparse.h, Matrix.h, Range.cc, Sparse.h, boolMatrix.h, boolNDArray.h, boolSparse.h, chMatrix.h, chNDArray.h, dColVector.h, dDiagMatrix.h, dMatrix.cc, dMatrix.h, dNDArray.cc, dNDArray.h, dRowVector.h, dSparse.cc, dSparse.h, dim-vector.cc, dim-vector.h, fCColVector.h, fCDiagMatrix.h, fCMatrix.cc, fCMatrix.h, fCNDArray.cc, fCNDArray.h, fCRowVector.h, fColVector.h, fDiagMatrix.h, fMatrix.cc, fMatrix.h, fNDArray.cc, fNDArray.h, fRowVector.h, int16NDArray.h, int32NDArray.h, int64NDArray.h, int8NDArray.h, intNDArray.h, uint16NDArray.h, uint32NDArray.h, uint64NDArray.h, uint8NDArray.h, CollocWt.cc, DASPK.cc, DASPK.h, DASRT.cc, DASRT.h, DASSL.cc, DASSL.h, LSODE.cc, LSODE.h, Quad.h, eigs-base.cc, eigs-base.h, lo-mappers.cc, lo-mappers.h, oct-norm.cc, lo-ieee.cc, lo-utils.cc, lo-utils.h, oct-refcount.h: Remove "#include <cfloat>". Add "#include <limits>" if necessary. Remove redundant #includes. Alphabetize #includes.

author	Rik <rik@octave.org>
date	Wed, 26 Apr 2017 09:57:28 -0700
parents	092078913d54
children	855122b993da

rev	line source
21562 6c2fd62db1f7 maint: Eliminate accidental double spaces in code. Rik <rik@octave.org> parents: 21301 diff changeset	1 // N-D Array manipulations.
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	2 /*
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	3
23219 3ac9f9ecfae5 maint: Update copyright dates. John W. Eaton <jwe@octave.org> parents: 23083 diff changeset	4 Copyright (C) 1996-2017 John W. Eaton
9601 a9b37bae1802 add a couple of missing copyright statements Jaroslav Hajek <highegg@gmail.com> parents: 9578 diff changeset	5 Copyright (C) 2009 VZLU Prague, a.s.
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	6
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	7 This file is part of Octave.
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	8
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	9 Octave is free software; you can redistribute it and/or modify it
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	10 under the terms of the GNU General Public License as published by
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	11 the Free Software Foundation; either version 3 of the License, or
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	12 (at your option) any later version.
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	13
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	14 Octave is distributed in the hope that it will be useful, but
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	15 WITHOUT ANY WARRANTY; without even the implied warranty of
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	16 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	17 GNU General Public License for more details.
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	18
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	19 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: 5775 diff changeset	20 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5775 diff changeset	21 <http://www.gnu.org/licenses/>.
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	22
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	23 */
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	24
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21723 diff changeset	25 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21202 diff changeset	26 # include "config.h"
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	27 #endif
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	28
23443 3f1bf237908b maint: Eliminate <cfloat.h> header from liboctave files. Rik <rik@octave.org> parents: 23220 diff changeset	29 #include <complex>
3f1bf237908b maint: Eliminate <cfloat.h> header from liboctave files. Rik <rik@octave.org> parents: 23220 diff changeset	30 #include <iostream>
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	31
4588 19bfd295f400 [project @ 2003-11-10 17:18:48 by jwe] jwe parents: 4584 diff changeset	32 #include "Array-util.h"
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	33 #include "CNDArray.h"
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	34 #include "f77-fcn.h"
7503 8c32f95c2639 convert mapper functions to new format David Bateman <dbateman@free.fr> parents: 7017 diff changeset	35 #include "functor.h"
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	36 #include "lo-ieee.h"
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	37 #include "lo-mappers.h"
9523 0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	38 #include "mx-base.h"
23443 3f1bf237908b maint: Eliminate <cfloat.h> header from liboctave files. Rik <rik@octave.org> parents: 23220 diff changeset	39 #include "mx-cnda-s.h"
8774 b756ce0002db split implementation and interface in mx-op-defs and MArray-defs Jaroslav Hajek <highegg@gmail.com> parents: 8751 diff changeset	40 #include "mx-op-defs.h"
4780 5eb5b8aaed8a [project @ 2004-02-17 15:32:13 by jwe] jwe parents: 4773 diff changeset	41 #include "oct-fftw.h"
9523 0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	42 #include "oct-locbuf.h"
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	43
9743 26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	44 #include "bsxfun-defs.cc"
26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	45
8956 d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	46 ComplexNDArray::ComplexNDArray (const charNDArray& a)
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10329 diff changeset	47 : MArray<Complex> (a.dims ())
8956 d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	48 {
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	49 octave_idx_type n = a.numel ();
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	50 for (octave_idx_type i = 0; i < n; i++)
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	51 xelem (i) = static_cast<unsigned char> (a(i));
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	52 }
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	53
9523 0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	54 #if defined (HAVE_FFTW)
0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	55
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	56 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	57 ComplexNDArray::fourier (int dim) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	58 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	59 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	60
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	61 if (dim > dv.ndims () \|\| dim < 0)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	62 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	63
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	64 octave_idx_type stride = 1;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	65 octave_idx_type n = dv(dim);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	66
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	67 for (int i = 0; i < dim; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	68 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	69
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	70 octave_idx_type howmany = numel () / dv(dim);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	71 howmany = (stride == 1 ? howmany : (howmany > stride ? stride : howmany));
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	72 octave_idx_type nloop = (stride == 1 ? 1 : numel () / dv(dim) / stride);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	73 octave_idx_type dist = (stride == 1 ? n : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	74
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	75 const Complex *in (fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	76 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	77 Complex *out (retval.fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	78
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	79 // Need to be careful here about the distance between fft's
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	80 for (octave_idx_type k = 0; k < nloop; k++)
23058 5839d16828cb move fftw classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	81 octave::fftw::fft (in + k * stride * n, out + k * stride * n,
5839d16828cb move fftw classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	82 n, howmany, stride, dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	83
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	84 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	85 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	86
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	87 ComplexNDArray
4816 4908c82cd1a1 [project @ 2004-03-04 18:13:46 by jwe] jwe parents: 4806 diff changeset	88 ComplexNDArray::ifourier (int dim) const
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	89 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	90 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	91
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	92 if (dim > dv.ndims () \|\| dim < 0)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	93 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	94
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	95 octave_idx_type stride = 1;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	96 octave_idx_type n = dv(dim);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	97
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	98 for (int i = 0; i < dim; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	99 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	100
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	101 octave_idx_type howmany = numel () / dv(dim);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	102 howmany = (stride == 1 ? howmany : (howmany > stride ? stride : howmany));
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	103 octave_idx_type nloop = (stride == 1 ? 1 : numel () / dv(dim) / stride);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	104 octave_idx_type dist = (stride == 1 ? n : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	105
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	106 const Complex *in (fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	107 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	108 Complex *out (retval.fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	109
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	110 // Need to be careful here about the distance between fft's
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	111 for (octave_idx_type k = 0; k < nloop; k++)
23058 5839d16828cb move fftw classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	112 octave::fftw::ifft (in + k * stride * n, out + k * stride * n,
5839d16828cb move fftw classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	113 n, howmany, stride, dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	114
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	115 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	116 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	117
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	118 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	119 ComplexNDArray::fourier2d (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	120 {
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14844 diff changeset	121 dim_vector dv = dims ();
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	122 if (dv.ndims () < 2)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	123 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	124
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	125 dim_vector dv2 (dv(0), dv(1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	126 const Complex *in = fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	127 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	128 Complex *out = retval.fortran_vec ();
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14844 diff changeset	129 octave_idx_type howmany = numel () / dv(0) / dv(1);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	130 octave_idx_type dist = dv(0) * dv(1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	131
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	132 for (octave_idx_type i=0; i < howmany; i++)
23058 5839d16828cb move fftw classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	133 octave::fftw::fftNd (in + idist, out + idist, 2, dv2);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	134
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	135 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	136 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	137
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	138 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	139 ComplexNDArray::ifourier2d (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	140 {
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14844 diff changeset	141 dim_vector dv = dims ();
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	142 if (dv.ndims () < 2)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	143 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	144
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	145 dim_vector dv2 (dv(0), dv(1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	146 const Complex *in = fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	147 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	148 Complex *out = retval.fortran_vec ();
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14844 diff changeset	149 octave_idx_type howmany = numel () / dv(0) / dv(1);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	150 octave_idx_type dist = dv(0) * dv(1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	151
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	152 for (octave_idx_type i=0; i < howmany; i++)
23058 5839d16828cb move fftw classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	153 octave::fftw::ifftNd (in + idist, out + idist, 2, dv2);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	154
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	155 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	156 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	157
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	158 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	159 ComplexNDArray::fourierNd (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	160 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	161 dim_vector dv = dims ();
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	162 int rank = dv.ndims ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	163
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	164 const Complex *in (fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	165 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	166 Complex *out (retval.fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	167
23058 5839d16828cb move fftw classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	168 octave::fftw::fftNd (in, out, rank, dv);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	169
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	170 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	171 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	172
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	173 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	174 ComplexNDArray::ifourierNd (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	175 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	176 dim_vector dv = dims ();
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	177 int rank = dv.ndims ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	178
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	179 const Complex *in (fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	180 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	181 Complex *out (retval.fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	182
23058 5839d16828cb move fftw classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22755 diff changeset	183 octave::fftw::ifftNd (in, out, rank, dv);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	184
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	185 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	186 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	187
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	188 #else
9523 0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	189
22322 93b3cdd36854 move most f77 function decls to separate header files John W. Eaton <jwe@octave.org> parents: 22274 diff changeset	190 #include "lo-fftpack-proto.h"
9523 0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	191
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	192 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	193 ComplexNDArray::fourier (int dim) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	194 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	195 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	196
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	197 if (dim > dv.ndims () \|\| dim < 0)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	198 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	199
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	200 ComplexNDArray retval (dv);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	201 octave_idx_type npts = dv(dim);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	202 octave_idx_type nn = 4*npts+15;
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	203 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	204 Complex *pwsave = wsave.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	205
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	206 OCTAVE_LOCAL_BUFFER (Complex, tmp, npts);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	207
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	208 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	209
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	210 for (int i = 0; i < dim; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	211 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	212
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	213 octave_idx_type howmany = numel () / npts;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	214 howmany = (stride == 1 ? howmany : (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	215 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	216 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	217
22261 77c4d43e06d1 Use new macros for passing complex arguments to Fortran FFTPACK subroutines Mike Miller <mtmiller@octave.org> parents: 22133 diff changeset	218 F77_FUNC (zffti, ZFFTI) (npts, F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	219
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	220 for (octave_idx_type k = 0; k < nloop; k++)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	221 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	222 for (octave_idx_type j = 0; j < howmany; j++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	223 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	224 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	225
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	226 for (octave_idx_type i = 0; i < npts; i++)
15018 3d8ace26c5b4 maint: Use Octave coding conventions for cuddled parentheses in liboctave/. Rik <rik@octave.org> parents: 14846 diff changeset	227 tmp[i] = elem ((i + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	228
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	229 F77_FUNC (zfftf, ZFFTF) (npts, F77_DBLE_CMPLX_ARG (tmp),
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	230 F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	231
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	232 for (octave_idx_type i = 0; i < npts; i++)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14557 diff changeset	233 retval((i + knpts)stride + j*dist) = tmp[i];
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	234 }
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	235 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	236
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	237 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	238 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	239
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	240 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	241 ComplexNDArray::ifourier (int dim) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	242 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	243 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	244
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	245 if (dim > dv.ndims () \|\| dim < 0)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	246 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	247
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	248 ComplexNDArray retval (dv);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	249 octave_idx_type npts = dv(dim);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	250 octave_idx_type nn = 4*npts+15;
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	251 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	252 Complex *pwsave = wsave.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	253
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	254 OCTAVE_LOCAL_BUFFER (Complex, tmp, npts);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	255
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	256 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	257
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	258 for (int i = 0; i < dim; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	259 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	260
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	261 octave_idx_type howmany = numel () / npts;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	262 howmany = (stride == 1 ? howmany : (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	263 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	264 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	265
22261 77c4d43e06d1 Use new macros for passing complex arguments to Fortran FFTPACK subroutines Mike Miller <mtmiller@octave.org> parents: 22133 diff changeset	266 F77_FUNC (zffti, ZFFTI) (npts, F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	267
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	268 for (octave_idx_type k = 0; k < nloop; k++)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	269 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	270 for (octave_idx_type j = 0; j < howmany; j++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	271 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	272 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	273
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	274 for (octave_idx_type i = 0; i < npts; i++)
15018 3d8ace26c5b4 maint: Use Octave coding conventions for cuddled parentheses in liboctave/. Rik <rik@octave.org> parents: 14846 diff changeset	275 tmp[i] = elem ((i + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	276
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	277 F77_FUNC (zfftb, ZFFTB) (npts, F77_DBLE_CMPLX_ARG (tmp),
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	278 F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	279
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	280 for (octave_idx_type i = 0; i < npts; i++)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14557 diff changeset	281 retval((i + knpts)stride + j*dist) = tmp[i] /
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	282 static_cast<double> (npts);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	283 }
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	284 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	285
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	286 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	287 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	288
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	289 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	290 ComplexNDArray::fourier2d (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	291 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	292 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	293 dim_vector dv2 (dv(0), dv(1));
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	294 int rank = 2;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	295 ComplexNDArray retval (*this);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	296 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	297
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	298 for (int i = 0; i < rank; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	299 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	300 octave_idx_type npts = dv2(i);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	301 octave_idx_type nn = 4*npts+15;
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	302 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	303 Complex *pwsave = wsave.fortran_vec ();
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	304 Array<Complex> row (dim_vector (npts, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	305 Complex *prow = row.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	306
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	307 octave_idx_type howmany = numel () / npts;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	308 howmany = (stride == 1 ? howmany :
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	309 (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	310 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	311 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	312
22261 77c4d43e06d1 Use new macros for passing complex arguments to Fortran FFTPACK subroutines Mike Miller <mtmiller@octave.org> parents: 22133 diff changeset	313 F77_FUNC (zffti, ZFFTI) (npts, F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	314
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	315 for (octave_idx_type k = 0; k < nloop; k++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	316 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	317 for (octave_idx_type j = 0; j < howmany; j++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	318 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	319 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	320
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	321 for (octave_idx_type l = 0; l < npts; l++)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14557 diff changeset	322 prow[l] = retval((l + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	323
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	324 F77_FUNC (zfftf, ZFFTF) (npts, F77_DBLE_CMPLX_ARG (prow),
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	325 F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	326
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	327 for (octave_idx_type l = 0; l < npts; l++)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14557 diff changeset	328 retval((l + knpts)stride + j*dist) = prow[l];
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	329 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	330 }
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	331
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	332 stride *= dv2(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	333 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	334
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	335 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	336 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	337
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	338 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	339 ComplexNDArray::ifourier2d (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	340 {
14846 460a3c6d8bf1 maint: Use Octave coding convention for cuddled parenthis in function calls with empty argument lists. Rik <octave@nomad.inbox5.com> parents: 14844 diff changeset	341 dim_vector dv = dims ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	342 dim_vector dv2 (dv(0), dv(1));
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	343 int rank = 2;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	344 ComplexNDArray retval (*this);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	345 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	346
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	347 for (int i = 0; i < rank; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	348 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	349 octave_idx_type npts = dv2(i);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	350 octave_idx_type nn = 4*npts+15;
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	351 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	352 Complex *pwsave = wsave.fortran_vec ();
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	353 Array<Complex> row (dim_vector (npts, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	354 Complex *prow = row.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	355
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	356 octave_idx_type howmany = numel () / npts;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	357 howmany = (stride == 1 ? howmany :
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	358 (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	359 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	360 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	361
22261 77c4d43e06d1 Use new macros for passing complex arguments to Fortran FFTPACK subroutines Mike Miller <mtmiller@octave.org> parents: 22133 diff changeset	362 F77_FUNC (zffti, ZFFTI) (npts, F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	363
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	364 for (octave_idx_type k = 0; k < nloop; k++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	365 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	366 for (octave_idx_type j = 0; j < howmany; j++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	367 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	368 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	369
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	370 for (octave_idx_type l = 0; l < npts; l++)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14557 diff changeset	371 prow[l] = retval((l + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	372
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	373 F77_FUNC (zfftb, ZFFTB) (npts, F77_DBLE_CMPLX_ARG (prow),
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	374 F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	375
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	376 for (octave_idx_type l = 0; l < npts; l++)
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	377 retval((l + knpts)stride + j*dist) =
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	378 prow[l] / static_cast<double> (npts);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	379 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	380 }
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	381
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	382 stride *= dv2(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	383 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	384
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	385 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	386 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	387
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	388 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	389 ComplexNDArray::fourierNd (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	390 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	391 dim_vector dv = dims ();
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	392 int rank = dv.ndims ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	393 ComplexNDArray retval (*this);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	394 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	395
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	396 for (int i = 0; i < rank; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	397 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	398 octave_idx_type npts = dv(i);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	399 octave_idx_type nn = 4*npts+15;
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	400 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	401 Complex *pwsave = wsave.fortran_vec ();
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	402 Array<Complex> row (dim_vector (npts, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	403 Complex *prow = row.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	404
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	405 octave_idx_type howmany = numel () / npts;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	406 howmany = (stride == 1 ? howmany :
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	407 (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	408 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	409 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	410
22261 77c4d43e06d1 Use new macros for passing complex arguments to Fortran FFTPACK subroutines Mike Miller <mtmiller@octave.org> parents: 22133 diff changeset	411 F77_FUNC (zffti, ZFFTI) (npts, F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	412
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	413 for (octave_idx_type k = 0; k < nloop; k++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	414 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	415 for (octave_idx_type j = 0; j < howmany; j++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	416 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	417 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	418
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	419 for (octave_idx_type l = 0; l < npts; l++)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14557 diff changeset	420 prow[l] = retval((l + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	421
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	422 F77_FUNC (zfftf, ZFFTF) (npts, F77_DBLE_CMPLX_ARG (prow),
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	423 F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	424
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	425 for (octave_idx_type l = 0; l < npts; l++)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14557 diff changeset	426 retval((l + knpts)stride + j*dist) = prow[l];
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	427 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	428 }
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	429
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	430 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	431 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	432
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	433 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	434 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	435
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	436 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	437 ComplexNDArray::ifourierNd (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	438 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	439 dim_vector dv = dims ();
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	440 int rank = dv.ndims ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	441 ComplexNDArray retval (*this);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	442 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	443
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	444 for (int i = 0; i < rank; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	445 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	446 octave_idx_type npts = dv(i);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	447 octave_idx_type nn = 4*npts+15;
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	448 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	449 Complex *pwsave = wsave.fortran_vec ();
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	450 Array<Complex> row (dim_vector (npts, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	451 Complex *prow = row.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	452
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	453 octave_idx_type howmany = numel () / npts;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	454 howmany = (stride == 1 ? howmany :
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	455 (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	456 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	457 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	458
22261 77c4d43e06d1 Use new macros for passing complex arguments to Fortran FFTPACK subroutines Mike Miller <mtmiller@octave.org> parents: 22133 diff changeset	459 F77_FUNC (zffti, ZFFTI) (npts, F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	460
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	461 for (octave_idx_type k = 0; k < nloop; k++)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	462 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	463 for (octave_idx_type j = 0; j < howmany; j++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	464 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	465 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	466
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	467 for (octave_idx_type l = 0; l < npts; l++)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14557 diff changeset	468 prow[l] = retval((l + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	469
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	470 F77_FUNC (zfftb, ZFFTB) (npts, F77_DBLE_CMPLX_ARG (prow),
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	471 F77_DBLE_CMPLX_ARG (pwsave));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	472
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	473 for (octave_idx_type l = 0; l < npts; l++)
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	474 retval((l + knpts)stride + j*dist) =
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	475 prow[l] / static_cast<double> (npts);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	476 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	477 }
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	478
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	479 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	480 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	481
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	482 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	483 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	484
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	485 #endif
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	486
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	487 // unary operations
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	488
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	489 boolNDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	490 ComplexNDArray::operator ! (void) const
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	491 {
11130 7c573eb981eb consistently give error for operator not applied to NaN values John W. Eaton <jwe@octave.org> parents: 11008 diff changeset	492 if (any_element_is_nan ())
22327 d0562b3159c7 move more classes inside octave namespace John W. Eaton <jwe@octave.org> parents: 22323 diff changeset	493 octave::err_nan_to_logical_conversion ();
11130 7c573eb981eb consistently give error for operator not applied to NaN values John W. Eaton <jwe@octave.org> parents: 11008 diff changeset	494
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	495 return do_mx_unary_op<bool, Complex> (*this, mx_inline_not);
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	496 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	497
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	498 // FIXME: this is not quite the right thing.
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	499
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	500 bool
7922 935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	501 ComplexNDArray::any_element_is_nan (void) const
935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	502 {
11008 3622db30ff05 simplify some array tests in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 10362 diff changeset	503 return do_mx_check<Complex> (*this, mx_inline_any_nan);
7922 935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	504 }
935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	505
935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	506 bool
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	507 ComplexNDArray::any_element_is_inf_or_nan (void) const
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	508 {
11008 3622db30ff05 simplify some array tests in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 10362 diff changeset	509 return ! do_mx_check<Complex> (*this, mx_inline_all_finite);
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	510 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	511
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	512 // Return true if no elements have imaginary components.
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	513
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	514 bool
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	515 ComplexNDArray::all_elements_are_real (void) const
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	516 {
11008 3622db30ff05 simplify some array tests in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 10362 diff changeset	517 return do_mx_check<Complex> (*this, mx_inline_all_real);
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	518 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	519
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	520 // Return nonzero if any element of CM has a non-integer real or
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	521 // imaginary part. Also extract the largest and smallest (real or
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	522 // imaginary) values and return them in MAX_VAL and MIN_VAL.
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	523
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	524 bool
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	525 ComplexNDArray::all_integers (double& max_val, double& min_val) const
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	526 {
20228 00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 20218 diff changeset	527 octave_idx_type nel = numel ();
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	528
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	529 if (nel > 0)
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	530 {
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	531 Complex val = elem (0);
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	532
22274 edcad35c364f don't rely on std:: overloads of arg, conj, imag, and real for real args John W. Eaton <jwe@octave.org> parents: 22261 diff changeset	533 double r_val = val.real ();
edcad35c364f don't rely on std:: overloads of arg, conj, imag, and real for real args John W. Eaton <jwe@octave.org> parents: 22261 diff changeset	534 double i_val = val.imag ();
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	535
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	536 max_val = r_val;
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	537 min_val = r_val;
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	538
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	539 if (i_val > max_val)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	540 max_val = i_val;
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	541
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	542 if (i_val < max_val)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	543 min_val = i_val;
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	544 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	545 else
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	546 return false;
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	547
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	548 for (octave_idx_type i = 0; i < nel; i++)
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	549 {
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	550 Complex val = elem (i);
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	551
22274 edcad35c364f don't rely on std:: overloads of arg, conj, imag, and real for real args John W. Eaton <jwe@octave.org> parents: 22261 diff changeset	552 double r_val = val.real ();
edcad35c364f don't rely on std:: overloads of arg, conj, imag, and real for real args John W. Eaton <jwe@octave.org> parents: 22261 diff changeset	553 double i_val = val.imag ();
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	554
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	555 if (r_val > max_val)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	556 max_val = r_val;
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	557
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	558 if (i_val > max_val)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	559 max_val = i_val;
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	560
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	561 if (r_val < min_val)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	562 min_val = r_val;
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	563
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	564 if (i_val < min_val)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	565 min_val = i_val;
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	566
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	567 if (octave::math::x_nint (r_val) != r_val
4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22327 diff changeset	568 \|\| octave::math::x_nint (i_val) != i_val)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	569 return false;
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	570 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	571
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	572 return true;
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	573 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	574
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	575 bool
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	576 ComplexNDArray::too_large_for_float (void) const
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	577 {
15212 4bbd3bbb8912 reduce code duplication in too_large_for_float array functions John W. Eaton <jwe@octave.org> parents: 15018 diff changeset	578 return test_any (xtoo_large_for_float);
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	579 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	580
4556 773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	581 boolNDArray
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	582 ComplexNDArray::all (int dim) const
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	583 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	584 return do_mx_red_op<bool, Complex> (*this, dim, mx_inline_all);
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	585 }
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	586
4556 773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	587 boolNDArray
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	588 ComplexNDArray::any (int dim) const
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	589 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	590 return do_mx_red_op<bool, Complex> (*this, dim, mx_inline_any);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	591 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	592
4584 f7697d703481 [project @ 2003-10-31 19:47:35 by jwe] jwe parents: 4569 diff changeset	593 ComplexNDArray
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	594 ComplexNDArray::cumprod (int dim) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	595 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	596 return do_mx_cum_op<Complex, Complex> (*this, dim, mx_inline_cumprod);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	597 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	598
4584 f7697d703481 [project @ 2003-10-31 19:47:35 by jwe] jwe parents: 4569 diff changeset	599 ComplexNDArray
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	600 ComplexNDArray::cumsum (int dim) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	601 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	602 return do_mx_cum_op<Complex, Complex> (*this, dim, mx_inline_cumsum);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	603 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	604
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	605 ComplexNDArray
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	606 ComplexNDArray::prod (int dim) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	607 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	608 return do_mx_red_op<Complex, Complex> (*this, dim, mx_inline_prod);
8736 53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	609 }
53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	610
53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	611 ComplexNDArray
53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	612 ComplexNDArray::sum (int dim) const
53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	613 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	614 return do_mx_red_op<Complex, Complex> (*this, dim, mx_inline_sum);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	615 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	616
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	617 ComplexNDArray
9721 192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9612 diff changeset	618 ComplexNDArray::xsum (int dim) const
192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9612 diff changeset	619 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	620 return do_mx_red_op<Complex, Complex> (*this, dim, mx_inline_xsum);
9721 192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9612 diff changeset	621 }
192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9612 diff changeset	622
192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9612 diff changeset	623 ComplexNDArray
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	624 ComplexNDArray::sumsq (int dim) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	625 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	626 return do_mx_red_op<double, Complex> (*this, dim, mx_inline_sumsq);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	627 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	628
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	629 ComplexNDArray
9513 9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	630 ComplexNDArray::diff (octave_idx_type order, int dim) const
9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	631 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	632 return do_mx_diff_op<Complex> (*this, dim, order, mx_inline_diff);
9513 9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	633 }
9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	634
9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	635 ComplexNDArray
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	636 ComplexNDArray::concat (const ComplexNDArray& rb,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	637 const Array<octave_idx_type>& ra_idx)
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	638 {
4940 954cc2ba6a49 [project @ 2004-08-24 19:30:31 by jwe] jwe parents: 4915 diff changeset	639 if (rb.numel () > 0)
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	640 insert (rb, ra_idx);
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	641 return *this;
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	642 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	643
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	644 ComplexNDArray
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	645 ComplexNDArray::concat (const NDArray& rb, const Array<octave_idx_type>& ra_idx)
4758 a308566c8b42 [project @ 2004-02-13 23:01:34 by jwe] jwe parents: 4687 diff changeset	646 {
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	647 ComplexNDArray tmp (rb);
4940 954cc2ba6a49 [project @ 2004-08-24 19:30:31 by jwe] jwe parents: 4915 diff changeset	648 if (rb.numel () > 0)
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	649 insert (tmp, ra_idx);
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	650 return *this;
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	651 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	652
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	653 ComplexNDArray
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	654 concat (NDArray& ra, ComplexNDArray& rb, const Array<octave_idx_type>& ra_idx)
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	655 {
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	656 ComplexNDArray retval (ra);
4940 954cc2ba6a49 [project @ 2004-08-24 19:30:31 by jwe] jwe parents: 4915 diff changeset	657 if (rb.numel () > 0)
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	658 retval.insert (rb, ra_idx);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	659 return retval;
4758 a308566c8b42 [project @ 2004-02-13 23:01:34 by jwe] jwe parents: 4687 diff changeset	660 }
a308566c8b42 [project @ 2004-02-13 23:01:34 by jwe] jwe parents: 4687 diff changeset	661
21723 bae585228161 use namespace for numeric_limits John W. Eaton <jwe@octave.org> parents: 21721 diff changeset	662 static const Complex Complex_NaN_result (octave::numeric_limits<double>::NaN (),
bae585228161 use namespace for numeric_limits John W. Eaton <jwe@octave.org> parents: 21721 diff changeset	663 octave::numeric_limits<double>::NaN ());
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	664
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	665 ComplexNDArray
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	666 ComplexNDArray::max (int dim) const
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	667 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	668 return do_mx_minmax_op<Complex> (*this, dim, mx_inline_max);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	669 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	670
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	671 ComplexNDArray
9732 b4fdfee405b5 remove ArrayN<T> + fix nonhom. diag-scalar ops Jaroslav Hajek <highegg@gmail.com> parents: 9721 diff changeset	672 ComplexNDArray::max (Array<octave_idx_type>& idx_arg, int dim) const
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	673 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	674 return do_mx_minmax_op<Complex> (*this, idx_arg, dim, mx_inline_max);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	675 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	676
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	677 ComplexNDArray
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	678 ComplexNDArray::min (int dim) const
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	679 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	680 return do_mx_minmax_op<Complex> (*this, dim, mx_inline_min);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	681 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	682
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	683 ComplexNDArray
9732 b4fdfee405b5 remove ArrayN<T> + fix nonhom. diag-scalar ops Jaroslav Hajek <highegg@gmail.com> parents: 9721 diff changeset	684 ComplexNDArray::min (Array<octave_idx_type>& idx_arg, int dim) const
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	685 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	686 return do_mx_minmax_op<Complex> (*this, idx_arg, dim, mx_inline_min);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	687 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	688
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	689 ComplexNDArray
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	690 ComplexNDArray::cummax (int dim) const
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	691 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	692 return do_mx_cumminmax_op<Complex> (*this, dim, mx_inline_cummax);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	693 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	694
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	695 ComplexNDArray
9732 b4fdfee405b5 remove ArrayN<T> + fix nonhom. diag-scalar ops Jaroslav Hajek <highegg@gmail.com> parents: 9721 diff changeset	696 ComplexNDArray::cummax (Array<octave_idx_type>& idx_arg, int dim) const
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	697 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	698 return do_mx_cumminmax_op<Complex> (*this, idx_arg, dim, mx_inline_cummax);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	699 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	700
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	701 ComplexNDArray
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	702 ComplexNDArray::cummin (int dim) const
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	703 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	704 return do_mx_cumminmax_op<Complex> (*this, dim, mx_inline_cummin);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	705 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	706
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	707 ComplexNDArray
9732 b4fdfee405b5 remove ArrayN<T> + fix nonhom. diag-scalar ops Jaroslav Hajek <highegg@gmail.com> parents: 9721 diff changeset	708 ComplexNDArray::cummin (Array<octave_idx_type>& idx_arg, int dim) const
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	709 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	710 return do_mx_cumminmax_op<Complex> (*this, idx_arg, dim, mx_inline_cummin);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	711 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	712
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	713 NDArray
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	714 ComplexNDArray::abs (void) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	715 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	716 return do_mx_unary_map<double, Complex, std::abs> (*this);
8650 a1ae2aae903e abs,real,imag,conj: use code from mx-inlines rather than the generic map Jaroslav Hajek <highegg@gmail.com> parents: 8377 diff changeset	717 }
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	718
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	719 boolNDArray
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	720 ComplexNDArray::isnan (void) const
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	721 {
21782 2aef506f3fec use namespace for lo-mappers.h functions John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	722 return do_mx_unary_map<bool, Complex, octave::math::isnan> (*this);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	723 }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	724
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	725 boolNDArray
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	726 ComplexNDArray::isinf (void) const
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	727 {
21782 2aef506f3fec use namespace for lo-mappers.h functions John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	728 return do_mx_unary_map<bool, Complex, octave::math::isinf> (*this);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	729 }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	730
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	731 boolNDArray
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	732 ComplexNDArray::isfinite (void) const
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	733 {
21782 2aef506f3fec use namespace for lo-mappers.h functions John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	734 return do_mx_unary_map<bool, Complex, octave::math::finite> (*this);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	735 }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	736
8650 a1ae2aae903e abs,real,imag,conj: use code from mx-inlines rather than the generic map Jaroslav Hajek <highegg@gmail.com> parents: 8377 diff changeset	737 ComplexNDArray
a1ae2aae903e abs,real,imag,conj: use code from mx-inlines rather than the generic map Jaroslav Hajek <highegg@gmail.com> parents: 8377 diff changeset	738 conj (const ComplexNDArray& a)
a1ae2aae903e abs,real,imag,conj: use code from mx-inlines rather than the generic map Jaroslav Hajek <highegg@gmail.com> parents: 8377 diff changeset	739 {
13107 353c71c76f22 maint: fix compilation problem with g++ -std=c++0x option Júlio Hoffimann <julio.hoffimann@gmail.com> parents: 11586 diff changeset	740 return do_mx_unary_map<Complex, Complex, std::conj<double> > (a);
4514 8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	741 }
8373b9b41eee [project @ 2003-09-19 21:47:10 by jwe] jwe parents: diff changeset	742
4765 e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	743 ComplexNDArray&
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	744 ComplexNDArray::insert (const NDArray& a, octave_idx_type r, octave_idx_type c)
4765 e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	745 {
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	746 dim_vector a_dv = a.dims ();
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	747
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	748 int n = a_dv.ndims ();
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	749
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	750 if (n != dimensions.ndims ())
4765 e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	751 (*current_liboctave_error_handler)
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	752 ("Array<T>::insert: invalid indexing operation");
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	753
21660 53728df3e4c9 maint: for clarity, call ndims() rather than length() on dimension vectors. Rik <rik@octave.org> parents: 21562 diff changeset	754 Array<octave_idx_type> a_ra_idx (dim_vector (a_dv.ndims (), 1), 0);
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	755
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	756 a_ra_idx.elem (0) = r;
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	757 a_ra_idx.elem (1) = c;
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	758
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	759 for (int i = 0; i < n; i++)
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	760 {
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	761 if (a_ra_idx(i) < 0 \|\| (a_ra_idx(i) + a_dv(i)) > dimensions(i))
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	762 (*current_liboctave_error_handler)
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	763 ("Array<T>::insert: range error for insert");
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	764 }
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	765
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	766 a_ra_idx.elem (0) = 0;
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	767 a_ra_idx.elem (1) = 0;
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	768
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	769 octave_idx_type n_elt = a.numel ();
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	770
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	771 // IS make_unique () NECESSARY HERE?
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	772
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	773 for (octave_idx_type i = 0; i < n_elt; i++)
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	774 {
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	775 Array<octave_idx_type> ra_idx = a_ra_idx;
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	776
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	777 ra_idx.elem (0) = a_ra_idx(0) + r;
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	778 ra_idx.elem (1) = a_ra_idx(1) + c;
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	779
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	780 elem (ra_idx) = a.elem (a_ra_idx);
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	781
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	782 increment_index (a_ra_idx, a_dv);
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	783 }
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	784
4765 e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	785 return *this;
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	786 }
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	787
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	788 ComplexNDArray&
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	789 ComplexNDArray::insert (const ComplexNDArray& a,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	790 octave_idx_type r, octave_idx_type c)
4765 e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	791 {
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	792 Array<Complex>::insert (a, r, c);
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	793 return *this;
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	794 }
e941e1470d7b [project @ 2004-02-16 05:56:50 by jwe] jwe parents: 4762 diff changeset	795
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	796 ComplexNDArray&
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	797 ComplexNDArray::insert (const ComplexNDArray& a,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	798 const Array<octave_idx_type>& ra_idx)
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	799 {
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	800 Array<Complex>::insert (a, ra_idx);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	801 return *this;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	802 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	803
4532 01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4514 diff changeset	804 void
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	805 ComplexNDArray::increment_index (Array<octave_idx_type>& ra_idx,
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	806 const dim_vector& dimensions,
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	807 int start_dimension)
4532 01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4514 diff changeset	808 {
01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4514 diff changeset	809 ::increment_index (ra_idx, dimensions, start_dimension);
01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4514 diff changeset	810 }
01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4514 diff changeset	811
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	812 octave_idx_type
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	813 ComplexNDArray::compute_index (Array<octave_idx_type>& ra_idx,
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	814 const dim_vector& dimensions)
4556 773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	815 {
773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	816 return ::compute_index (ra_idx, dimensions);
773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	817 }
773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	818
7620 36594d5bbe13 Move diag function into the octave_value class David Bateman <dbateman@free.fr> parents: 7600 diff changeset	819 ComplexNDArray
36594d5bbe13 Move diag function into the octave_value class David Bateman <dbateman@free.fr> parents: 7600 diff changeset	820 ComplexNDArray::diag (octave_idx_type k) const
36594d5bbe13 Move diag function into the octave_value class David Bateman <dbateman@free.fr> parents: 7600 diff changeset	821 {
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10329 diff changeset	822 return MArray<Complex>::diag (k);
7620 36594d5bbe13 Move diag function into the octave_value class David Bateman <dbateman@free.fr> parents: 7600 diff changeset	823 }
36594d5bbe13 Move diag function into the octave_value class David Bateman <dbateman@free.fr> parents: 7600 diff changeset	824
14557 e8e86ae3abbc make diag (x, m, n) return a proper diagonal matrix object (bug #36099) John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	825 ComplexNDArray
e8e86ae3abbc make diag (x, m, n) return a proper diagonal matrix object (bug #36099) John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	826 ComplexNDArray::diag (octave_idx_type m, octave_idx_type n) const
e8e86ae3abbc make diag (x, m, n) return a proper diagonal matrix object (bug #36099) John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	827 {
e8e86ae3abbc make diag (x, m, n) return a proper diagonal matrix object (bug #36099) John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	828 return MArray<Complex>::diag (m, n);
e8e86ae3abbc make diag (x, m, n) return a proper diagonal matrix object (bug #36099) John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	829 }
e8e86ae3abbc make diag (x, m, n) return a proper diagonal matrix object (bug #36099) John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	830
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	831 // This contains no information on the array structure !!!
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	832 std::ostream&
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	833 operator << (std::ostream& os, const ComplexNDArray& a)
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	834 {
20228 00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 20218 diff changeset	835 octave_idx_type nel = a.numel ();
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	836
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	837 for (octave_idx_type i = 0; i < nel; i++)
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	838 {
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	839 os << " ";
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	840 octave_write_complex (os, a.elem (i));
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	841 os << "\n";
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	842 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	843 return os;
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	844 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	845
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	846 std::istream&
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	847 operator >> (std::istream& is, ComplexNDArray& a)
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	848 {
20228 00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 20218 diff changeset	849 octave_idx_type nel = a.numel ();
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	850
8999 dc07bc4157b8 allow empty matrices in stream input operators Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	851 if (nel > 0)
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	852 {
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	853 Complex tmp;
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5260 diff changeset	854 for (octave_idx_type i = 0; i < nel; i++)
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	855 {
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	856 tmp = octave_read_value<Complex> (is);
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	857 if (is)
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	858 a.elem (i) = tmp;
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	859 else
21168 26f85aa072de maint: Replace instances of goto in liboctave where convenient. Rik <rik@octave.org> parents: 21136 diff changeset	860 return is;
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	861 }
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	862 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	863
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	864 return is;
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	865 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	866
10329 83fa590b8a09 simplify min/max definitions in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 10314 diff changeset	867 MINMAX_FCNS (ComplexNDArray, Complex)
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	868
9578 7dafdb8b062f refactor comparison ops implementations Jaroslav Hajek <highegg@gmail.com> parents: 9553 diff changeset	869 NDS_CMP_OPS (ComplexNDArray, Complex)
9550 3d6a9aea2aea refactor binary & bool ops in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 9546 diff changeset	870 NDS_BOOL_OPS (ComplexNDArray, Complex)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	871
9578 7dafdb8b062f refactor comparison ops implementations Jaroslav Hajek <highegg@gmail.com> parents: 9553 diff changeset	872 SND_CMP_OPS (Complex, ComplexNDArray)
9550 3d6a9aea2aea refactor binary & bool ops in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 9546 diff changeset	873 SND_BOOL_OPS (Complex, ComplexNDArray)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	874
9578 7dafdb8b062f refactor comparison ops implementations Jaroslav Hajek <highegg@gmail.com> parents: 9553 diff changeset	875 NDND_CMP_OPS (ComplexNDArray, ComplexNDArray)
9550 3d6a9aea2aea refactor binary & bool ops in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 9546 diff changeset	876 NDND_BOOL_OPS (ComplexNDArray, ComplexNDArray)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	877
9546 1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	878 ComplexNDArray& operator *= (ComplexNDArray& a, double s)
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	879 {
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	880 if (a.is_shared ())
9612 66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	881 a = a * s;
66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	882 else
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	883 do_ms_inplace_op<Complex, double> (a, s, mx_inline_mul2);
9546 1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	884 return a;
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	885 }
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	886
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	887 ComplexNDArray& operator /= (ComplexNDArray& a, double s)
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	888 {
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	889 if (a.is_shared ())
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	890 return a = a / s;
9612 66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	891 else
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	892 do_ms_inplace_op<Complex, double> (a, s, mx_inline_div2);
9546 1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	893 return a;
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	894 }
1beb23d2b892 optimize op= in common cases Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	895
9743 26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	896 BSXFUN_STDOP_DEFS_MXLOOP (ComplexNDArray)
26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	897 BSXFUN_STDREL_DEFS_MXLOOP (ComplexNDArray)
26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	898
9827 c15a5ed0da58 optimize bsxfun (@power, ...) Jaroslav Hajek <highegg@gmail.com> parents: 9825 diff changeset	899 BSXFUN_OP_DEF_MXLOOP (pow, ComplexNDArray, mx_inline_pow)

Mercurial > octave

annotate liboctave/array/CNDArray.cc @ 23443:3f1bf237908b