octave-nkf: liboctave/CNDArray.cc annotate

annotate liboctave/CNDArray.cc @ 11008:3622db30ff05

simplify some array tests in liboctave

author	Jaroslav Hajek <highegg@gmail.com>
date	Tue, 21 Sep 2010 13:09:12 +0200
parents	b47ab50a6aa8
children	7c573eb981eb

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

Mercurial > hg > octave-nkf

annotate liboctave/CNDArray.cc @ 11008:3622db30ff05