octave: liboctave/array/dNDArray.cc annotate

annotate liboctave/array/dNDArray.cc @ 21202:f7121e111991

maint: indent #ifdef blocks in liboctave and src directories. * Array-C.cc, Array-b.cc, Array-ch.cc, Array-d.cc, Array-f.cc, Array-fC.cc, Array-i.cc, Array-idx-vec.cc, Array-s.cc, Array-str.cc, Array-util.cc, Array-voidp.cc, Array.cc, CColVector.cc, CDiagMatrix.cc, CMatrix.cc, CNDArray.cc, CRowVector.cc, CSparse.cc, CSparse.h, DiagArray2.cc, MArray-C.cc, MArray-d.cc, MArray-f.cc, MArray-fC.cc, MArray-i.cc, MArray-s.cc, MArray.cc, MDiagArray2.cc, MSparse-C.cc, MSparse-d.cc, MSparse.h, MatrixType.cc, PermMatrix.cc, Range.cc, Sparse-C.cc, Sparse-b.cc, Sparse-d.cc, Sparse.cc, boolMatrix.cc, boolNDArray.cc, boolSparse.cc, chMatrix.cc, chNDArray.cc, dColVector.cc, dDiagMatrix.cc, dMatrix.cc, dNDArray.cc, dRowVector.cc, dSparse.cc, dSparse.h, dim-vector.cc, fCColVector.cc, fCDiagMatrix.cc, fCMatrix.cc, fCNDArray.cc, fCRowVector.cc, fColVector.cc, fDiagMatrix.cc, fMatrix.cc, fNDArray.cc, fRowVector.cc, idx-vector.cc, int16NDArray.cc, int32NDArray.cc, int64NDArray.cc, int8NDArray.cc, intNDArray.cc, uint16NDArray.cc, uint32NDArray.cc, uint64NDArray.cc, uint8NDArray.cc, blaswrap.c, cquit.c, f77-extern.cc, f77-fcn.c, f77-fcn.h, lo-error.c, quit.cc, quit.h, CmplxAEPBAL.cc, CmplxCHOL.cc, CmplxGEPBAL.cc, CmplxHESS.cc, CmplxLU.cc, CmplxQR.cc, CmplxQRP.cc, CmplxSCHUR.cc, CmplxSVD.cc, CollocWt.cc, DASPK.cc, DASRT.cc, DASSL.cc, EIG.cc, LSODE.cc, ODES.cc, Quad.cc, base-lu.cc, base-qr.cc, dbleAEPBAL.cc, dbleCHOL.cc, dbleGEPBAL.cc, dbleHESS.cc, dbleLU.cc, dbleQR.cc, dbleQRP.cc, dbleSCHUR.cc, dbleSVD.cc, eigs-base.cc, fCmplxAEPBAL.cc, fCmplxCHOL.cc, fCmplxGEPBAL.cc, fCmplxHESS.cc, fCmplxLU.cc, fCmplxQR.cc, fCmplxQRP.cc, fCmplxSCHUR.cc, fCmplxSVD.cc, fEIG.cc, floatAEPBAL.cc, floatCHOL.cc, floatGEPBAL.cc, floatHESS.cc, floatLU.cc, floatQR.cc, floatQRP.cc, floatSCHUR.cc, floatSVD.cc, lo-mappers.cc, lo-specfun.cc, oct-convn.cc, oct-fftw.cc, oct-fftw.h, oct-norm.cc, oct-rand.cc, oct-spparms.cc, randgamma.c, randmtzig.c, randpoisson.c, sparse-chol.cc, sparse-dmsolve.cc, sparse-lu.cc, sparse-qr.cc, mx-defs.h, dir-ops.cc, file-ops.cc, file-stat.cc, lo-sysdep.cc, mach-info.cc, oct-env.cc, oct-group.cc, oct-openmp.h, oct-passwd.cc, oct-syscalls.cc, oct-time.cc, oct-uname.cc, pathlen.h, sysdir.h, syswait.h, cmd-edit.cc, cmd-hist.cc, data-conv.cc, f2c-main.c, glob-match.cc, lo-array-errwarn.cc, lo-array-gripes.cc, lo-cutils.c, lo-cutils.h, lo-ieee.cc, lo-math.h, lo-regexp.cc, lo-utils.cc, oct-base64.cc, oct-glob.cc, oct-inttypes.cc, oct-inttypes.h, oct-locbuf.cc, oct-mutex.cc, oct-refcount.h, oct-rl-edit.c, oct-rl-hist.c, oct-shlib.cc, oct-sort.cc, pathsearch.cc, singleton-cleanup.cc, sparse-sort.cc, sparse-util.cc, statdefs.h, str-vec.cc, unwind-prot.cc, url-transfer.cc, display-available.h, main-cli.cc, main-gui.cc, main.in.cc, mkoctfile.in.cc, octave-config.in.cc, shared-fcns.h: indent #ifdef blocks in liboctave and src directories.

author	Rik <rik@octave.org>
date	Sat, 06 Feb 2016 06:40:13 -0800
parents	40051830f89b
children	40de9f8f23a6

rev	line source
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	1 // N-D Array manipulations.
4511 24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	2 /*
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	3
19697 4197fc428c7d maint: Update copyright notices for 2015. John W. Eaton <jwe@octave.org> parents: 19342 diff changeset	4 Copyright (C) 1996-2015 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.
4511 24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	6
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	7 This file is part of Octave.
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	8
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	9 Octave is free software; you can redistribute it and/or modify it
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	10 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6989 diff changeset	11 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6989 diff changeset	12 option) any later version.
4511 24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	13
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	14 Octave is distributed in the hope that it will be useful, but WITHOUT
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	15 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	16 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	17 for more details.
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	18
24af46b4ce84 [project @ 2003-09-12 16:46:04 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: 6989 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: 6989 diff changeset	21 <http://www.gnu.org/licenses/>.
4511 24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	22
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	23 */
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	24
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	25 #ifdef HAVE_CONFIG_H
21202 f7121e111991 maint: indent #ifdef blocks in liboctave and src directories. Rik <rik@octave.org> parents: 21169 diff changeset	26 # include <config.h>
4511 24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	27 #endif
24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	28
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	29 #include <cfloat>
5164 57077d0ddc8e [project @ 2005-02-25 19:55:24 by jwe] jwe parents: 5110 diff changeset	30
4780 5eb5b8aaed8a [project @ 2004-02-17 15:32:13 by jwe] jwe parents: 4773 diff changeset	31 #include <vector>
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	32
4588 19bfd295f400 [project @ 2003-11-10 17:18:48 by jwe] jwe parents: 4584 diff changeset	33 #include "Array-util.h"
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	34 #include "dNDArray.h"
9523 0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	35 #include "f77-fcn.h"
7503 8c32f95c2639 convert mapper functions to new format David Bateman <dbateman@free.fr> parents: 7017 diff changeset	36 #include "functor.h"
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	37 #include "lo-error.h"
4511 24af46b4ce84 [project @ 2003-09-12 16:46:04 by jwe] jwe parents: diff changeset	38 #include "lo-ieee.h"
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	39 #include "lo-mappers.h"
9523 0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	40 #include "mx-base.h"
8774 b756ce0002db split implementation and interface in mx-op-defs and MArray-defs Jaroslav Hajek <highegg@gmail.com> parents: 8751 diff changeset	41 #include "mx-op-defs.h"
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	42 #include "oct-fftw.h"
9523 0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	43 #include "oct-locbuf.h"
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	44
9743 26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	45 #include "bsxfun-defs.cc"
26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	46
7919 9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	47 NDArray::NDArray (const Array<octave_idx_type>& a, bool zero_based,
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	48 bool negative_to_nan)
7919 9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	49 {
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	50 const octave_idx_type *pa = a.fortran_vec ();
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	51 resize (a.dims ());
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	52 double *ptmp = fortran_vec ();
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	53 if (negative_to_nan)
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	54 {
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	55 double nan_val = lo_ieee_nan_value ();
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	56
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	57 if (zero_based)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	58 for (octave_idx_type i = 0; i < a.numel (); i++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	59 {
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	60 double val = static_cast<double>
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	61 (pa[i] + static_cast<octave_idx_type> (1));
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	62 if (val <= 0)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	63 ptmp[i] = nan_val;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	64 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	65 ptmp[i] = val;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	66 }
7919 9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	67 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	68 for (octave_idx_type i = 0; i < a.numel (); i++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	69 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	70 double val = static_cast<double> (pa[i]);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	71 if (val <= 0)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	72 ptmp[i] = nan_val;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	73 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	74 ptmp[i] = val;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	75 }
7919 9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	76 }
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	77 else
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	78 {
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	79 if (zero_based)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	80 for (octave_idx_type i = 0; i < a.numel (); i++)
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	81 ptmp[i] = static_cast<double>
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	82 (pa[i] + static_cast<octave_idx_type> (1));
7919 9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	83 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	84 for (octave_idx_type i = 0; i < a.numel (); i++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	85 ptmp[i] = static_cast<double> (pa[i]);
7919 9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	86 }
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	87 }
9d080df0c843 new NDArray constructor for ArrayN<octave_idx_type> David Bateman <dbateman@free.fr> parents: 7789 diff changeset	88
8956 d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	89 NDArray::NDArray (const charNDArray& a)
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10329 diff changeset	90 : MArray<double> (a.dims ())
8956 d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	91 {
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	92 octave_idx_type n = a.numel ();
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	93 for (octave_idx_type i = 0; i < n; i++)
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	94 xelem (i) = static_cast<unsigned char> (a(i));
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	95 }
d91fa4b20bbb ensure nonnegative char -> real conversion Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	96
9523 0ce82753dd72 more configure changes for libraries John W. Eaton <jwe@octave.org> parents: 9513 diff changeset	97 #if defined (HAVE_FFTW)
7941 f8cab9eeb128 Fix NDArray compilation/export John W. Eaton <jwe@octave.org> parents: 7922 diff changeset	98
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	99 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	100 NDArray::fourier (int dim) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	101 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	102 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	103
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	104 if (dim > dv.length () \|\| dim < 0)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	105 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	106
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	107 octave_idx_type stride = 1;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	108 octave_idx_type n = dv(dim);
4773 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 for (int i = 0; i < dim; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	111 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	112
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	113 octave_idx_type howmany = numel () / dv(dim);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	114 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	115 octave_idx_type nloop = (stride == 1 ? 1 : numel () / dv(dim) / stride);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	116 octave_idx_type dist = (stride == 1 ? n : 1);
4773 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 const double *in (fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	119 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	120 Complex *out (retval.fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	121
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	122 // Need to be careful here about the distance between fft's
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	123 for (octave_idx_type k = 0; k < nloop; k++)
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	124 octave_fftw::fft (in + k * stride * n, out + k * stride * n,
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	125 n, howmany, stride, dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	126
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	127 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	128 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	129
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	130 ComplexNDArray
4816 4908c82cd1a1 [project @ 2004-03-04 18:13:46 by jwe] jwe parents: 4806 diff changeset	131 NDArray::ifourier (int dim) const
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	132 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	133 dim_vector dv = dims ();
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 if (dim > dv.length () \|\| dim < 0)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	136 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	137
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	138 octave_idx_type stride = 1;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	139 octave_idx_type n = dv(dim);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	140
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	141 for (int i = 0; i < dim; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	142 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	143
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	144 octave_idx_type howmany = numel () / dv(dim);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	145 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	146 octave_idx_type nloop = (stride == 1 ? 1 : numel () / dv(dim) / stride);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	147 octave_idx_type dist = (stride == 1 ? n : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	148
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	149 ComplexNDArray retval (*this);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	150 Complex *out (retval.fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	151
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	152 // Need to be careful here about the distance between fft's
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	153 for (octave_idx_type k = 0; k < nloop; k++)
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	154 octave_fftw::ifft (out + k * stride * n, out + k * stride * n,
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	155 n, howmany, stride, dist);
4773 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 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	158 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	159
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	160 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	161 NDArray::fourier2d (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	162 {
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	163 dim_vector dv = dims ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	164 if (dv.length () < 2)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	165 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	166
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	167 dim_vector dv2 (dv(0), dv(1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	168 const double *in = fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	169 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	170 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	171 octave_idx_type howmany = numel () / dv(0) / dv(1);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	172 octave_idx_type dist = dv(0) * dv(1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	173
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	174 for (octave_idx_type i=0; i < howmany; i++)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	175 octave_fftw::fftNd (in + idist, out + idist, 2, dv2);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	176
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	177 return retval;
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
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	180 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	181 NDArray::ifourier2d (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	182 {
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	183 dim_vector dv = dims ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	184 if (dv.length () < 2)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	185 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	186
20218 b2100e1659ac maint: Use cuddled parentheses when indexing dimension_vectors. Rik <rik@octave.org> parents: 19697 diff changeset	187 dim_vector dv2 (dv(0), dv(1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	188 ComplexNDArray retval (*this);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	189 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	190 octave_idx_type howmany = numel () / dv(0) / dv(1);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	191 octave_idx_type dist = dv(0) * dv(1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	192
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	193 for (octave_idx_type i=0; i < howmany; i++)
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	194 octave_fftw::ifftNd (out + idist, out + idist, 2, dv2);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	195
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	196 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	197 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	198
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	199 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	200 NDArray::fourierNd (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	201 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	202 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	203 int rank = dv.length ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	204
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	205 const double *in (fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	206 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	207 Complex *out (retval.fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	208
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	209 octave_fftw::fftNd (in, out, rank, dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	210
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	211 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	212 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	213
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	214 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	215 NDArray::ifourierNd (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	216 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	217 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	218 int rank = dv.length ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	219
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	220 ComplexNDArray tmp (*this);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	221 Complex *in (tmp.fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	222 ComplexNDArray retval (dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	223 Complex *out (retval.fortran_vec ());
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	224
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	225 octave_fftw::ifftNd (in, out, rank, dv);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	226
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	227 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	228 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	229
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	230 #else
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	231
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	232 extern "C"
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	233 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	234 // Note that the original complex fft routines were not written for
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	235 // double complex arguments. They have been modified by adding an
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	236 // implicit double precision (a-h,o-z) statement at the beginning of
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	237 // each subroutine.
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 F77_RET_T
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7620 diff changeset	240 F77_FUNC (zffti, ZFFTI) (const octave_idx_type&, Complex*);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	241
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	242 F77_RET_T
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7620 diff changeset	243 F77_FUNC (zfftf, ZFFTF) (const octave_idx_type&, Complex, Complex);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	244
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	245 F77_RET_T
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7620 diff changeset	246 F77_FUNC (zfftb, ZFFTB) (const octave_idx_type&, Complex, Complex);
4773 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
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	249 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	250 NDArray::fourier (int dim) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	251 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	252 dim_vector dv = dims ();
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 if (dim > dv.length () \|\| dim < 0)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	255 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	256
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	257 ComplexNDArray retval (dv);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	258 octave_idx_type npts = dv(dim);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	259 octave_idx_type nn = 4*npts+15;
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	260 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	261 Complex *pwsave = wsave.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	262
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	263 OCTAVE_LOCAL_BUFFER (Complex, tmp, npts);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	264
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	265 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	266
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	267 for (int i = 0; i < dim; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	268 stride *= dv(i);
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: 5164 diff changeset	270 octave_idx_type howmany = numel () / npts;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	271 howmany = (stride == 1 ? howmany : (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	272 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	273 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	274
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7620 diff changeset	275 F77_FUNC (zffti, ZFFTI) (npts, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	276
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	277 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	278 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	279 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	280 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	281 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	282
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	283 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	284 tmp[i] = elem ((i + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	285
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	286 F77_FUNC (zfftf, ZFFTF) (npts, tmp, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	287
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	288 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	289 retval((i + knpts)stride + j*dist) = tmp[i];
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	290 }
4773 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
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	293 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	294 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	295
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	296 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	297 NDArray::ifourier (int dim) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	298 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	299 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	300
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	301 if (dim > dv.length () \|\| dim < 0)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	302 return ComplexNDArray ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	303
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	304 ComplexNDArray retval (dv);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	305 octave_idx_type npts = dv(dim);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	306 octave_idx_type nn = 4*npts+15;
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	307 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	308 Complex *pwsave = wsave.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	309
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	310 OCTAVE_LOCAL_BUFFER (Complex, tmp, npts);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	311
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	312 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	313
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	314 for (int i = 0; i < dim; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	315 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	316
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	317 octave_idx_type howmany = numel () / npts;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	318 howmany = (stride == 1 ? howmany : (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	319 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	320 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	321
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7620 diff changeset	322 F77_FUNC (zffti, ZFFTI) (npts, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	323
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	324 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	325 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	326 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	327 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	328 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	329
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	330 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	331 tmp[i] = elem ((i + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	332
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	333 F77_FUNC (zfftb, ZFFTB) (npts, tmp, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	334
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	335 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	336 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	337 static_cast<double> (npts);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	338 }
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	339 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	340
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	341 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	342 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	343
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	344 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	345 NDArray::fourier2d (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	346 {
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	347 dim_vector dv = dims ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	348 dim_vector dv2 (dv(0), dv(1));
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	349 int rank = 2;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	350 ComplexNDArray retval (*this);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	351 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	352
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	353 for (int i = 0; i < rank; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	354 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	355 octave_idx_type npts = dv2(i);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	356 octave_idx_type nn = 4*npts+15;
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	357 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	358 Complex *pwsave = wsave.fortran_vec ();
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	359 Array<Complex> row (dim_vector (npts, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	360 Complex *prow = row.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	361
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	362 octave_idx_type howmany = numel () / npts;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	363 howmany = (stride == 1 ? howmany :
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	364 (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	365 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	366 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	367
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7620 diff changeset	368 F77_FUNC (zffti, ZFFTI) (npts, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	369
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	370 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	371 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	372 for (octave_idx_type j = 0; j < howmany; j++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	373 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	374 octave_quit ();
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++)
14844 5bc9b9cb4362 maint: Use Octave coding conventions for cuddled parenthesis in retval assignments. Rik <octave@nomad.inbox5.com> parents: 14557 diff changeset	377 prow[l] = retval((l + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	378
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	379 F77_FUNC (zfftf, ZFFTF) (npts, prow, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	380
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	381 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	382 retval((l + knpts)stride + j*dist) = prow[l];
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	383 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	384 }
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	385
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	386 stride *= dv2(i);
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
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	389 return retval;
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
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	392 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	393 NDArray::ifourier2d (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	394 {
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	395 dim_vector dv = dims ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	396 dim_vector dv2 (dv(0), dv(1));
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	397 int rank = 2;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	398 ComplexNDArray retval (*this);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	399 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	400
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	401 for (int i = 0; i < rank; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	402 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	403 octave_idx_type npts = dv2(i);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	404 octave_idx_type nn = 4*npts+15;
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	405 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	406 Complex *pwsave = wsave.fortran_vec ();
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	407 Array<Complex> row (dim_vector (npts, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	408 Complex *prow = row.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	409
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	410 octave_idx_type howmany = numel () / npts;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	411 howmany = (stride == 1 ? howmany :
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	412 (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	413 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	414 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	415
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7620 diff changeset	416 F77_FUNC (zffti, ZFFTI) (npts, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	417
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	418 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	419 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	420 for (octave_idx_type j = 0; j < howmany; j++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	421 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	422 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	423
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	424 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	425 prow[l] = retval((l + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	426
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	427 F77_FUNC (zfftb, ZFFTB) (npts, prow, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	428
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	429 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	430 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	431 prow[l] / static_cast<double> (npts);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	432 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	433 }
4773 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 stride *= dv2(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	436 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	437
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	438 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	439 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	440
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	441 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	442 NDArray::fourierNd (void) const
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 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	445 int rank = dv.length ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	446 ComplexNDArray retval (*this);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	447 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	448
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	449 for (int i = 0; i < rank; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	450 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	451 octave_idx_type npts = dv(i);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	452 octave_idx_type nn = 4*npts+15;
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	453 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	454 Complex *pwsave = wsave.fortran_vec ();
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	455 Array<Complex> row (dim_vector (npts, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	456 Complex *prow = row.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	457
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	458 octave_idx_type howmany = numel () / npts;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	459 howmany = (stride == 1 ? howmany :
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	460 (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	461 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	462 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	463
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7620 diff changeset	464 F77_FUNC (zffti, ZFFTI) (npts, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	465
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	466 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	467 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	468 for (octave_idx_type j = 0; j < howmany; j++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	469 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	470 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	471
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	472 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	473 prow[l] = retval((l + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	474
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	475 F77_FUNC (zfftf, ZFFTF) (npts, prow, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	476
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	477 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	478 retval((l + knpts)stride + j*dist) = prow[l];
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	479 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	480 }
4773 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 stride *= dv(i);
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 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	486 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	487
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	488 ComplexNDArray
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	489 NDArray::ifourierNd (void) const
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	490 {
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	491 dim_vector dv = dims ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	492 int rank = dv.length ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	493 ComplexNDArray retval (*this);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	494 octave_idx_type stride = 1;
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	495
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	496 for (int i = 0; i < rank; i++)
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	497 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	498 octave_idx_type npts = dv(i);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	499 octave_idx_type nn = 4*npts+15;
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	500 Array<Complex> wsave (dim_vector (nn, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	501 Complex *pwsave = wsave.fortran_vec ();
19342 8f6a8422435d Replace instances of deprecated Array constructor in FFTPACK-related code. Rik <rik@octave.org> parents: 19320 diff changeset	502 Array<Complex> row (dim_vector (npts, 1));
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	503 Complex *prow = row.fortran_vec ();
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	504
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	505 octave_idx_type howmany = numel () / npts;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	506 howmany = (stride == 1 ? howmany :
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	507 (howmany > stride ? stride : howmany));
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	508 octave_idx_type nloop = (stride == 1 ? 1 : numel () / npts / stride);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	509 octave_idx_type dist = (stride == 1 ? npts : 1);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	510
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7620 diff changeset	511 F77_FUNC (zffti, ZFFTI) (npts, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	512
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	513 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	514 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	515 for (octave_idx_type j = 0; j < howmany; j++)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	516 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	517 octave_quit ();
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	518
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	519 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	520 prow[l] = retval((l + knpts)stride + j*dist);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	521
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	522 F77_FUNC (zfftb, ZFFTB) (npts, prow, pwsave);
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	523
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	524 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	525 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	526 prow[l] / static_cast<double> (npts);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	527 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	528 }
4773 ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	529
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	530 stride *= dv(i);
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	531 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	532
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	533 return retval;
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	534 }
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	535
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	536 #endif
ccfbd6047a54 [project @ 2004-02-16 19:02:32 by jwe] jwe parents: 4770 diff changeset	537
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	538 // unary operations
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	539
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	540 boolNDArray
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	541 NDArray::operator ! (void) const
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	542 {
11130 7c573eb981eb consistently give error for operator not applied to NaN values John W. Eaton <jwe@octave.org> parents: 11010 diff changeset	543 if (any_element_is_nan ())
21100 e39e05d90788 Switch gripe_XXX to either err_XXX or warn_XXX naming scheme. Rik <rik@octave.org> parents: 20974 diff changeset	544 err_nan_to_logical_conversion ();
11130 7c573eb981eb consistently give error for operator not applied to NaN values John W. Eaton <jwe@octave.org> parents: 11010 diff changeset	545
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	546 return do_mx_unary_op<bool, double> (*this, mx_inline_not);
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	547 }
79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	548
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	549 bool
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	550 NDArray::any_element_is_negative (bool neg_zero) const
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	551 {
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 11008 diff changeset	552 return (neg_zero ? test_all (xnegative_sign)
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 11008 diff changeset	553 : do_mx_check<double> (*this, mx_inline_any_negative));
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	554 }
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	555
7922 935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7919 diff changeset	556 bool
13756 6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13005 diff changeset	557 NDArray::any_element_is_positive (bool neg_zero) const
6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13005 diff changeset	558 {
6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13005 diff changeset	559 return (neg_zero ? test_all (xpositive_sign)
6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13005 diff changeset	560 : do_mx_check<double> (*this, mx_inline_any_positive));
6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13005 diff changeset	561 }
6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13005 diff changeset	562
6dfebfa334cb allow negative data log plots with OpenGL+FLTK graphics (bug #34232) John W. Eaton <jwe@octave.org> parents: 13005 diff changeset	563 bool
7922 935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7919 diff changeset	564 NDArray::any_element_is_nan (void) const
935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7919 diff changeset	565 {
11008 3622db30ff05 simplify some array tests in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 10362 diff changeset	566 return do_mx_check<double> (*this, mx_inline_any_nan);
7922 935be827eaf8 error for NaN values in & and \| expressions John W. Eaton <jwe@octave.org> parents: 7919 diff changeset	567 }
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	568
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	569 bool
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	570 NDArray::any_element_is_inf_or_nan (void) const
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	571 {
11008 3622db30ff05 simplify some array tests in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 10362 diff changeset	572 return ! do_mx_check<double> (*this, mx_inline_all_finite);
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	573 }
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	574
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	575 bool
5943 c038c2947ee1 [project @ 2006-08-18 21:27:03 by jwe] jwe parents: 5775 diff changeset	576 NDArray::any_element_not_one_or_zero (void) const
c038c2947ee1 [project @ 2006-08-18 21:27:03 by jwe] jwe parents: 5775 diff changeset	577 {
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 11008 diff changeset	578 return ! test_all (xis_one_or_zero);
5943 c038c2947ee1 [project @ 2006-08-18 21:27:03 by jwe] jwe parents: 5775 diff changeset	579 }
c038c2947ee1 [project @ 2006-08-18 21:27:03 by jwe] jwe parents: 5775 diff changeset	580
c038c2947ee1 [project @ 2006-08-18 21:27:03 by jwe] jwe parents: 5775 diff changeset	581 bool
6989 2d326000e09b [project @ 2007-10-09 20:32:42 by jwe] jwe parents: 5943 diff changeset	582 NDArray::all_elements_are_zero (void) const
2d326000e09b [project @ 2007-10-09 20:32:42 by jwe] jwe parents: 5943 diff changeset	583 {
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 11008 diff changeset	584 return test_all (xis_zero);
6989 2d326000e09b [project @ 2007-10-09 20:32:42 by jwe] jwe parents: 5943 diff changeset	585 }
2d326000e09b [project @ 2007-10-09 20:32:42 by jwe] jwe parents: 5943 diff changeset	586
2d326000e09b [project @ 2007-10-09 20:32:42 by jwe] jwe parents: 5943 diff changeset	587 bool
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	588 NDArray::all_elements_are_int_or_inf_or_nan (void) const
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	589 {
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 11008 diff changeset	590 return test_all (xis_int_or_inf_or_nan);
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	591 }
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	592
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	593 // Return nonzero if any element of M is not an integer. Also extract
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	594 // the largest and smallest values and return them in MAX_VAL and MIN_VAL.
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	595
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	596 bool
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	597 NDArray::all_integers (double& max_val, double& min_val) const
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	598 {
20228 00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 20218 diff changeset	599 octave_idx_type nel = numel ();
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	600
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	601 if (nel > 0)
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	602 {
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	603 max_val = elem (0);
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	604 min_val = elem (0);
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	605 }
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	606 else
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	607 return false;
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	608
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	609 for (octave_idx_type i = 0; i < nel; i++)
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	610 {
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	611 double val = elem (i);
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	612
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	613 if (val > max_val)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	614 max_val = val;
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	615
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	616 if (val < min_val)
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	617 min_val = val;
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	618
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 11008 diff changeset	619 if (! xisinteger (val))
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	620 return false;
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	621 }
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	622
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	623 return true;
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	624 }
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	625
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	626 bool
9827 c15a5ed0da58 optimize bsxfun (@power, ...) Jaroslav Hajek <highegg@gmail.com> parents: 9812 diff changeset	627 NDArray::all_integers (void) const
c15a5ed0da58 optimize bsxfun (@power, ...) Jaroslav Hajek <highegg@gmail.com> parents: 9812 diff changeset	628 {
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 11008 diff changeset	629 return test_all (xisinteger);
9827 c15a5ed0da58 optimize bsxfun (@power, ...) Jaroslav Hajek <highegg@gmail.com> parents: 9812 diff changeset	630 }
c15a5ed0da58 optimize bsxfun (@power, ...) Jaroslav Hajek <highegg@gmail.com> parents: 9812 diff changeset	631
c15a5ed0da58 optimize bsxfun (@power, ...) Jaroslav Hajek <highegg@gmail.com> parents: 9812 diff changeset	632 bool
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	633 NDArray::too_large_for_float (void) const
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	634 {
15212 4bbd3bbb8912 reduce code duplication in too_large_for_float array functions John W. Eaton <jwe@octave.org> parents: 15018 diff changeset	635 return test_any (xtoo_large_for_float);
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	636 }
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	637
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	638 // FIXME: this is not quite the right thing.
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	639
4556 773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	640 boolNDArray
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	641 NDArray::all (int dim) const
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	642 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	643 return do_mx_red_op<bool, double> (*this, dim, mx_inline_all);
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	644 }
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	645
4556 773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	646 boolNDArray
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	647 NDArray::any (int dim) const
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	648 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	649 return do_mx_red_op<bool, double> (*this, dim, mx_inline_any);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	650 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	651
4584 f7697d703481 [project @ 2003-10-31 19:47:35 by jwe] jwe parents: 4569 diff changeset	652 NDArray
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	653 NDArray::cumprod (int dim) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	654 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	655 return do_mx_cum_op<double, double> (*this, dim, mx_inline_cumprod);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	656 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	657
4584 f7697d703481 [project @ 2003-10-31 19:47:35 by jwe] jwe parents: 4569 diff changeset	658 NDArray
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	659 NDArray::cumsum (int dim) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	660 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	661 return do_mx_cum_op<double, double> (*this, dim, mx_inline_cumsum);
4513 508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	662 }
508238e65af7 [project @ 2003-09-19 21:40:57 by jwe] jwe parents: 4511 diff changeset	663
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	664 NDArray
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	665 NDArray::prod (int dim) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	666 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	667 return do_mx_red_op<double, double> (*this, dim, mx_inline_prod);
8736 53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	668 }
53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	669
53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	670 NDArray
53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	671 NDArray::sum (int dim) const
53b4fdeacc2e improve reduction functions Jaroslav Hajek <highegg@gmail.com> parents: 8650 diff changeset	672 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	673 return do_mx_red_op<double, double> (*this, dim, mx_inline_sum);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	674 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	675
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	676 NDArray
9721 192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	677 NDArray::xsum (int dim) const
192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	678 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	679 return do_mx_red_op<double, double> (*this, dim, mx_inline_xsum);
9721 192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	680 }
192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	681
192d94cff6c1 improve sum & implement the 'extra' option, refactor some code Jaroslav Hajek <highegg@gmail.com> parents: 9601 diff changeset	682 NDArray
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	683 NDArray::sumsq (int dim) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	684 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	685 return do_mx_red_op<double, double> (*this, dim, mx_inline_sumsq);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	686 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	687
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	688 NDArray
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	689 NDArray::max (int dim) const
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	690 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	691 return do_mx_minmax_op<double> (*this, dim, mx_inline_max);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	692 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	693
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	694 NDArray
9732 b4fdfee405b5 remove ArrayN<T> + fix nonhom. diag-scalar ops Jaroslav Hajek <highegg@gmail.com> parents: 9721 diff changeset	695 NDArray::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	696 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	697 return do_mx_minmax_op<double> (*this, idx_arg, dim, mx_inline_max);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	698 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	699
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	700 NDArray
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	701 NDArray::min (int dim) const
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	702 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	703 return do_mx_minmax_op<double> (*this, dim, mx_inline_min);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	704 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	705
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	706 NDArray
9732 b4fdfee405b5 remove ArrayN<T> + fix nonhom. diag-scalar ops Jaroslav Hajek <highegg@gmail.com> parents: 9721 diff changeset	707 NDArray::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	708 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	709 return do_mx_minmax_op<double> (*this, idx_arg, dim, mx_inline_min);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	710 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	711
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	712 NDArray
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	713 NDArray::cummax (int dim) const
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	714 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	715 return do_mx_cumminmax_op<double> (*this, dim, mx_inline_cummax);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	716 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	717
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	718 NDArray
9732 b4fdfee405b5 remove ArrayN<T> + fix nonhom. diag-scalar ops Jaroslav Hajek <highegg@gmail.com> parents: 9721 diff changeset	719 NDArray::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	720 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	721 return do_mx_cumminmax_op<double> (*this, idx_arg, dim, mx_inline_cummax);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	722 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	723
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	724 NDArray
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	725 NDArray::cummin (int dim) const
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	726 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	727 return do_mx_cumminmax_op<double> (*this, dim, mx_inline_cummin);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	728 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	729
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	730 NDArray
9732 b4fdfee405b5 remove ArrayN<T> + fix nonhom. diag-scalar ops Jaroslav Hajek <highegg@gmail.com> parents: 9721 diff changeset	731 NDArray::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	732 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	733 return do_mx_cumminmax_op<double> (*this, idx_arg, dim, mx_inline_cummin);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	734 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	735
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 8774 diff changeset	736 NDArray
9513 9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	737 NDArray::diff (octave_idx_type order, int dim) const
9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	738 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	739 return do_mx_diff_op<double> (*this, dim, order, mx_inline_diff);
9513 9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	740 }
9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	741
9f870f73ab7d implement built-in diff Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	742 NDArray
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	743 NDArray::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	744 {
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	745 if (rb.numel () > 0)
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	746 insert (rb, ra_idx);
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	747 return *this;
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	748 }
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	749
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	750 ComplexNDArray
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	751 NDArray::concat (const ComplexNDArray& rb, const Array<octave_idx_type>& ra_idx)
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	752 {
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	753 ComplexNDArray retval (*this);
4940 954cc2ba6a49 [project @ 2004-08-24 19:30:31 by jwe] jwe parents: 4915 diff changeset	754 if (rb.numel () > 0)
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	755 retval.insert (rb, ra_idx);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	756 return retval;
4758 a308566c8b42 [project @ 2004-02-13 23:01:34 by jwe] jwe parents: 4687 diff changeset	757 }
a308566c8b42 [project @ 2004-02-13 23:01:34 by jwe] jwe parents: 4687 diff changeset	758
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	759 charNDArray
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	760 NDArray::concat (const charNDArray& rb, const Array<octave_idx_type>& ra_idx)
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	761 {
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	762 charNDArray retval (dims ());
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	763 octave_idx_type nel = numel ();
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	764
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	765 for (octave_idx_type i = 0; i < nel; i++)
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	766 {
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	767 double d = elem (i);
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	768
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	769 if (xisnan (d))
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	770 (*current_liboctave_error_handler)
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	771 ("invalid conversion from NaN to character");
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 octave_idx_type ival = NINTbig (d);
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	774
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	775 if (ival < 0 \|\| ival > std::numeric_limits<unsigned char>::max ())
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	776 // FIXME: is there something better to do? Should we warn the user?
7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	777 ival = 0;
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	778
21136 7cac4e7458f2 maint: clean up code around calls to current_liboctave_error_handler. Rik <rik@octave.org> parents: 21100 diff changeset	779 retval.elem (i) = static_cast<char>(ival);
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	780 }
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	781
20974 1edf15793cac maint: Use is_empty rather than "numel () == 0" for clarity. Rik <rik@octave.org> parents: 20228 diff changeset	782 if (rb.is_empty ())
5073 e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	783 return retval;
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	784
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	785 retval.insert (rb, ra_idx);
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	786 return retval;
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	787 }
e71be9c548f2 [project @ 2004-11-09 18:31:25 by jwe] jwe parents: 4940 diff changeset	788
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	789 NDArray
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	790 real (const ComplexNDArray& a)
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	791 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	792 return do_mx_unary_op<double, Complex> (a, mx_inline_real);
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	793 }
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	794
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	795 NDArray
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	796 imag (const ComplexNDArray& a)
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	797 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	798 return do_mx_unary_op<double, Complex> (a, mx_inline_imag);
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	799 }
79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	800
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	801 NDArray&
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	802 NDArray::insert (const NDArray& a, octave_idx_type r, octave_idx_type c)
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	803 {
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	804 Array<double>::insert (a, r, c);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	805 return *this;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	806 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	807
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	808 NDArray&
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	809 NDArray::insert (const NDArray& a, const Array<octave_idx_type>& ra_idx)
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	810 {
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	811 Array<double>::insert (a, ra_idx);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	812 return *this;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	813 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4902 diff changeset	814
4634 79fe96966ca0 [project @ 2003-11-19 21:22:39 by jwe] jwe parents: 4588 diff changeset	815 NDArray
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	816 NDArray::abs (void) const
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	817 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	818 return do_mx_unary_map<double, double, std::abs> (*this);
4569 01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	819 }
01e4957409a4 [project @ 2003-10-30 05:48:17 by jwe] jwe parents: 4563 diff changeset	820
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	821 boolNDArray
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	822 NDArray::isnan (void) const
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	823 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	824 return do_mx_unary_map<bool, double, xisnan> (*this);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	825 }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	826
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	827 boolNDArray
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	828 NDArray::isinf (void) const
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	829 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	830 return do_mx_unary_map<bool, double, xisinf> (*this);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	831 }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	832
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	833 boolNDArray
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	834 NDArray::isfinite (void) const
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	835 {
10362 b47ab50a6aa8 simplify appliers in mx-inlines.cc Jaroslav Hajek <highegg@gmail.com> parents: 10352 diff changeset	836 return do_mx_unary_map<bool, double, xfinite> (*this);
8998 a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	837 }
a48fba01e4ac optimize isnan/isinf/isfinite mappers Jaroslav Hajek <highegg@gmail.com> parents: 8981 diff changeset	838
4532 01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4513 diff changeset	839 void
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	840 NDArray::increment_index (Array<octave_idx_type>& ra_idx,
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	841 const dim_vector& dimensions,
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	842 int start_dimension)
4532 01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4513 diff changeset	843 {
01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4513 diff changeset	844 ::increment_index (ra_idx, dimensions, start_dimension);
01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4513 diff changeset	845 }
01ee68d18069 [project @ 2003-10-10 00:49:54 by jwe] jwe parents: 4513 diff changeset	846
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	847 octave_idx_type
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	848 NDArray::compute_index (Array<octave_idx_type>& ra_idx,
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	849 const dim_vector& dimensions)
4556 773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	850 {
773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	851 return ::compute_index (ra_idx, dimensions);
773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	852 }
773a21e4fce8 [project @ 2003-10-28 05:02:10 by jwe] jwe parents: 4543 diff changeset	853
7503 8c32f95c2639 convert mapper functions to new format David Bateman <dbateman@free.fr> parents: 7017 diff changeset	854 NDArray
7620 36594d5bbe13 Move diag function into the octave_value class David Bateman <dbateman@free.fr> parents: 7600 diff changeset	855 NDArray::diag (octave_idx_type k) const
36594d5bbe13 Move diag function into the octave_value class David Bateman <dbateman@free.fr> parents: 7600 diff changeset	856 {
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10329 diff changeset	857 return MArray<double>::diag (k);
7620 36594d5bbe13 Move diag function into the octave_value class David Bateman <dbateman@free.fr> parents: 7600 diff changeset	858 }
36594d5bbe13 Move diag function into the octave_value class David Bateman <dbateman@free.fr> parents: 7600 diff changeset	859
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	860 NDArray
e8e86ae3abbc make diag (x, m, n) return a proper diagonal matrix object (bug #36099) John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	861 NDArray::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	862 {
e8e86ae3abbc make diag (x, m, n) return a proper diagonal matrix object (bug #36099) John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	863 return MArray<double>::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	864 }
e8e86ae3abbc make diag (x, m, n) return a proper diagonal matrix object (bug #36099) John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	865
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	866 // This contains no information on the array structure !!!
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	867 std::ostream&
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	868 operator << (std::ostream& os, const NDArray& a)
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	869 {
20228 00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 20218 diff changeset	870 octave_idx_type nel = a.numel ();
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	871
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	872 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	873 {
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	874 os << " ";
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	875 octave_write_double (os, a.elem (i));
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	876 os << "\n";
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	877 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	878 return os;
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	879 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	880
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	881 std::istream&
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	882 operator >> (std::istream& is, NDArray& a)
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	883 {
20228 00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel(). Carnë Draug <carandraug@octave.org> parents: 20218 diff changeset	884 octave_idx_type nel = a.numel ();
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	885
8999 dc07bc4157b8 allow empty matrices in stream input operators Jaroslav Hajek <highegg@gmail.com> parents: 8998 diff changeset	886 if (nel > 0)
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	887 {
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	888 double tmp;
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5164 diff changeset	889 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	890 {
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	891 tmp = octave_read_value<double> (is);
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	892 if (is)
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	893 a.elem (i) = tmp;
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	894 else
21169 40051830f89b maint: Replace more instances of goto in liboctave with other constructs. Rik <rik@octave.org> parents: 21136 diff changeset	895 return is;
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	896 }
4687 e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	897 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	898
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	899 return is;
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	900 }
e95c86d48732 [project @ 2004-01-06 21:53:34 by jwe] jwe parents: 4634 diff changeset	901
10329 83fa590b8a09 simplify min/max definitions in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 10314 diff changeset	902 MINMAX_FCNS (NDArray, double)
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4816 diff changeset	903
9578 7dafdb8b062f refactor comparison ops implementations Jaroslav Hajek <highegg@gmail.com> parents: 9553 diff changeset	904 NDS_CMP_OPS (NDArray, double)
9550 3d6a9aea2aea refactor binary & bool ops in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	905 NDS_BOOL_OPS (NDArray, double)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	906
9578 7dafdb8b062f refactor comparison ops implementations Jaroslav Hajek <highegg@gmail.com> parents: 9553 diff changeset	907 SND_CMP_OPS (double, NDArray)
9550 3d6a9aea2aea refactor binary & bool ops in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	908 SND_BOOL_OPS (double, NDArray)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	909
9578 7dafdb8b062f refactor comparison ops implementations Jaroslav Hajek <highegg@gmail.com> parents: 9553 diff changeset	910 NDND_CMP_OPS (NDArray, NDArray)
9550 3d6a9aea2aea refactor binary & bool ops in liboctave Jaroslav Hajek <highegg@gmail.com> parents: 9523 diff changeset	911 NDND_BOOL_OPS (NDArray, NDArray)
4543 79df15d4470c [project @ 2003-10-18 03:53:52 by jwe] jwe parents: 4532 diff changeset	912
9743 26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	913 BSXFUN_STDOP_DEFS_MXLOOP (NDArray)
26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	914 BSXFUN_STDREL_DEFS_MXLOOP (NDArray)
26abff55f6fe optimize bsxfun for common built-in operations Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	915
9827 c15a5ed0da58 optimize bsxfun (@power, ...) Jaroslav Hajek <highegg@gmail.com> parents: 9812 diff changeset	916 BSXFUN_OP_DEF_MXLOOP (pow, NDArray, mx_inline_pow)
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	917 BSXFUN_OP2_DEF_MXLOOP (pow, ComplexNDArray, ComplexNDArray,
9827 c15a5ed0da58 optimize bsxfun (@power, ...) Jaroslav Hajek <highegg@gmail.com> parents: 9812 diff changeset	918 NDArray, mx_inline_pow)
13005 4061106b1c4b Enable automatic bsxfun for power operators Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 11586 diff changeset	919 BSXFUN_OP2_DEF_MXLOOP (pow, ComplexNDArray, NDArray,
4061106b1c4b Enable automatic bsxfun for power operators Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 11586 diff changeset	920 ComplexNDArray, mx_inline_pow)

Mercurial > octave

annotate liboctave/array/dNDArray.cc @ 21202:f7121e111991