octave-nkf: src/DLD-FUNCTIONS/filter.cc annotate

annotate src/DLD-FUNCTIONS/filter.cc @ 12706:7c8b50cd4fb4 stable

Backout 12700:9843b3b055e0 from stable branch

author	Jordi Gutiérrez Hermoso <jordigh@gmail.com>
date	Mon, 06 Jun 2011 12:38:15 -0500
parents	9843b3b055e0
children	a5f4dad4ec27

rev	line source
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	1 /*
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	2
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 10846 diff changeset	3 Copyright (C) 1996-2011 John W. Eaton
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	4
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	5 This file is part of Octave.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	6
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	7 Octave is free software; you can redistribute it and/or modify it
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	8 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: 6476 diff changeset	9 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6476 diff changeset	10 option) any later version.
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	11
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	12 Octave is distributed in the hope that it will be useful, but WITHOUT
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	15 for more details.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	16
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	17 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: 6476 diff changeset	18 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6476 diff changeset	19 <http://www.gnu.org/licenses/>.
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	20
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	21 */
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	22
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	23 // Based on Tony Richardson's filter.m.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	24 //
5428 2a16423e4aa0 [project @ 2005-08-23 18:38:27 by jwe] jwe parents: 5307 diff changeset	25 // Originally translated to C++ by KH (Kurt.Hornik@wu-wien.ac.at)
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	26 // with help from Fritz Leisch and Andreas Weingessel on Oct 20, 1994.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	27 //
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	28 // Rewritten to use templates to handle both real and complex cases by
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	29 // jwe, Wed Nov 1 19:15:29 1995.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	30
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	31 #ifdef HAVE_CONFIG_H
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	32 #include <config.h>
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	33 #endif
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	34
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3419 diff changeset	35 #include "quit.h"
6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3419 diff changeset	36
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	37 #include "defun-dld.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	38 #include "error.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	39 #include "oct-obj.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	40
3109 fcb2f5063e7f [project @ 1997-11-20 08:52:42 by jwe] jwe parents: 3107 diff changeset	41 #if !defined (CXX_NEW_FRIEND_TEMPLATE_DECL)
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	42 extern MArray<double>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	43 filter (MArray<double>&, MArray<double>&, MArray<double>&, int dim);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	44
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	45 extern MArray<Complex>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	46 filter (MArray<Complex>&, MArray<Complex>&, MArray<Complex>&, int dim);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	47
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	48 extern MArray<float>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	49 filter (MArray<float>&, MArray<float>&, MArray<float>&, int dim);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	50
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	51 extern MArray<FloatComplex>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	52 filter (MArray<FloatComplex>&, MArray<FloatComplex>&, MArray<FloatComplex>&, int dim);
3107 a8dcfbf87ea3 [project @ 1997-11-19 22:34:59 by jwe] jwe parents: 3097 diff changeset	53 #endif
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	54
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	55 template <class T>
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	56 MArray<T>
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11574 diff changeset	57 filter (MArray<T>& b, MArray<T>& a, MArray<T>& x, MArray<T>& si,
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	58 int dim = 0)
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	59 {
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	60 MArray<T> y;
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	61
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	62 octave_idx_type a_len = a.length ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	63 octave_idx_type b_len = b.length ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	64
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	65 octave_idx_type ab_len = a_len > b_len ? a_len : b_len;
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	66
12637 d02798e9eeae filter.cc: Stop test errors in residue.m (bug #33164) Rik <octave@nomad.inbox5.com> parents: 12566 diff changeset	67 // FIXME: The two lines below should be unecessary because
d02798e9eeae filter.cc: Stop test errors in residue.m (bug #33164) Rik <octave@nomad.inbox5.com> parents: 12566 diff changeset	68 // this template is called with a and b as column vectors
d02798e9eeae filter.cc: Stop test errors in residue.m (bug #33164) Rik <octave@nomad.inbox5.com> parents: 12566 diff changeset	69 // already. However the a.resize line is currently (2011/04/26)
d02798e9eeae filter.cc: Stop test errors in residue.m (bug #33164) Rik <octave@nomad.inbox5.com> parents: 12566 diff changeset	70 // necessary to stop bug #33164.
d02798e9eeae filter.cc: Stop test errors in residue.m (bug #33164) Rik <octave@nomad.inbox5.com> parents: 12566 diff changeset	71 b.resize (dim_vector (ab_len, 1), 0.0);
d02798e9eeae filter.cc: Stop test errors in residue.m (bug #33164) Rik <octave@nomad.inbox5.com> parents: 12566 diff changeset	72 if (a_len > 1)
d02798e9eeae filter.cc: Stop test errors in residue.m (bug #33164) Rik <octave@nomad.inbox5.com> parents: 12566 diff changeset	73 a.resize (dim_vector (ab_len, 1), 0.0);
d02798e9eeae filter.cc: Stop test errors in residue.m (bug #33164) Rik <octave@nomad.inbox5.com> parents: 12566 diff changeset	74
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	75 T norm = a (0);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	76
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	77 if (norm == static_cast<T>(0.0))
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	78 {
11553 01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	79 error ("filter: the first element of A must be non-zero");
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	80 return y;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	81 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	82
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	83 dim_vector x_dims = x.dims ();
5056 9379d8716407 [project @ 2004-10-22 14:58:05 by jwe] jwe parents: 5055 diff changeset	84 if (dim < 0 \|\| dim > x_dims.length ())
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	85 {
12483 7a5aacf65f81 Rewrite error strings in src/ to use variables named in documentation. Rik <octave@nomad.inbox5.com> parents: 11586 diff changeset	86 error ("filter: DIM must be a valid dimension");
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	87 return y;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	88 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	89
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	90 octave_idx_type x_len = x_dims(dim);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	91
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	92 dim_vector si_dims = si.dims ();
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	93 octave_idx_type si_len = si_dims(0);
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	94
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	95 if (si_len != ab_len - 1)
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	96 {
11553 01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	97 error ("filter: first dimension of SI must be of length max (length (a), length (b)) - 1");
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	98 return y;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	99 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	100
5059 38dd53c8935e [project @ 2004-11-02 02:42:25 by jwe] jwe parents: 5056 diff changeset	101 if (si_dims.length () != x_dims.length ())
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	102 {
11553 01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	103 error ("filter: dimensionality of SI and X must agree");
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	104 return y;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	105 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	106
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	107 octave_idx_type si_dim = 0;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	108 for (octave_idx_type i = 0; i < x_dims.length (); i++)
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	109 {
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	110 if (i == dim)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	111 continue;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11574 diff changeset	112
5055 51a4406317e9 [project @ 2004-10-22 14:46:52 by jwe] jwe parents: 4901 diff changeset	113 if (x_dims(i) == 1)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	114 continue;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11574 diff changeset	115
5059 38dd53c8935e [project @ 2004-11-02 02:42:25 by jwe] jwe parents: 5056 diff changeset	116 if (si_dims(++si_dim) != x_dims(i))
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	117 {
11553 01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	118 error ("filter: dimensionality of SI and X must agree");
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	119 return y;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	120 }
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	121 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	122
10241 a277ba5da4dc filter.cc: return early if X is empty John W. Eaton <jwe@octave.org> parents: 10155 diff changeset	123 if (x_len == 0)
a277ba5da4dc filter.cc: return early if X is empty John W. Eaton <jwe@octave.org> parents: 10155 diff changeset	124 return x;
a277ba5da4dc filter.cc: return early if X is empty John W. Eaton <jwe@octave.org> parents: 10155 diff changeset	125
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	126 if (norm != static_cast<T>(1.0))
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	127 {
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	128 a = a / norm;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	129 b = b / norm;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	130 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	131
5056 9379d8716407 [project @ 2004-10-22 14:58:05 by jwe] jwe parents: 5055 diff changeset	132 if (a_len <= 1 && si_len <= 0)
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	133 return b(0) * x;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	134
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	135 y.resize (x_dims, 0.0);
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	136
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	137 int x_stride = 1;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	138 for (int i = 0; i < dim; i++)
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	139 x_stride *= x_dims(i);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	140
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	141 octave_idx_type x_num = x_dims.numel () / x_len;
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	142 for (octave_idx_type num = 0; num < x_num; num++)
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	143 {
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	144 octave_idx_type x_offset;
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	145 if (x_stride == 1)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	146 x_offset = num * x_len;
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	147 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	148 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	149 octave_idx_type x_offset2 = 0;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	150 x_offset = num;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	151 while (x_offset >= x_stride)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	152 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	153 x_offset -= x_stride;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	154 x_offset2++;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	155 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	156 x_offset += x_offset2 * x_stride * x_len;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	157 }
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	158 octave_idx_type si_offset = num * si_len;
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	159
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	160 if (a_len > 1)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	161 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	162 T *py = y.fortran_vec ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	163 T *psi = si.fortran_vec ();
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	164
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	165 const T *pa = a.data ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	166 const T *pb = b.data ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	167 const T *px = x.data ();
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	168
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	169 psi += si_offset;
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	170
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	171 for (octave_idx_type i = 0, idx = x_offset; i < x_len; i++, idx += x_stride)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	172 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	173 py[idx] = psi[0] + pb[0] * px[idx];
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	174
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	175 if (si_len > 0)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	176 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	177 for (octave_idx_type j = 0; j < si_len - 1; j++)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	178 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	179 OCTAVE_QUIT;
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	180
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	181 psi[j] = psi[j+1] - pa[j+1] * py[idx] + pb[j+1] * px[idx];
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	182 }
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	183
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	184 psi[si_len-1] = pb[si_len] * px[idx] - pa[si_len] * py[idx];
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	185 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	186 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	187 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	188 OCTAVE_QUIT;
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	189
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	190 psi[0] = pb[si_len] * px[idx] - pa[si_len] * py[idx];
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	191 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	192 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	193 }
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	194 else if (si_len > 0)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	195 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	196 T *py = y.fortran_vec ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	197 T *psi = si.fortran_vec ();
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	198
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	199 const T *pb = b.data ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	200 const T *px = x.data ();
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	201
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	202 psi += si_offset;
5114 cda8c0a823c5 [project @ 2005-01-21 17:02:21 by jwe] jwe parents: 5059 diff changeset	203
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	204 for (octave_idx_type i = 0, idx = x_offset; i < x_len; i++, idx += x_stride)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	205 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	206 py[idx] = psi[0] + pb[0] * px[idx];
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	207
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	208 if (si_len > 1)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	209 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	210 for (octave_idx_type j = 0; j < si_len - 1; j++)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	211 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	212 OCTAVE_QUIT;
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 3419 diff changeset	213
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	214 psi[j] = psi[j+1] + pb[j+1] * px[idx];
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	215 }
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	216
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	217 psi[si_len-1] = pb[si_len] * px[idx];
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	218 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	219 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	220 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	221 OCTAVE_QUIT;
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	222
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	223 psi[0] = pb[1] * px[idx];
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	224 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	225 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	226 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	227 }
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11574 diff changeset	228
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	229 return y;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	230 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	231
3109 fcb2f5063e7f [project @ 1997-11-20 08:52:42 by jwe] jwe parents: 3107 diff changeset	232 #if !defined (CXX_NEW_FRIEND_TEMPLATE_DECL)
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	233 extern MArray<double>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	234 filter (MArray<double>&, MArray<double>&, MArray<double>&,
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	235 MArray<double>&, int dim);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	236
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	237 extern MArray<Complex>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	238 filter (MArray<Complex>&, MArray<Complex>&, MArray<Complex>&,
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	239 MArray<Complex>&, int dim);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	240
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	241 extern MArray<float>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	242 filter (MArray<float>&, MArray<float>&, MArray<float>&,
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	243 MArray<float>&, int dim);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	244
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	245 extern MArray<FloatComplex>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	246 filter (MArray<FloatComplex>&, MArray<FloatComplex>&, MArray<FloatComplex>&,
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	247 MArray<FloatComplex>&, int dim);
3107 a8dcfbf87ea3 [project @ 1997-11-19 22:34:59 by jwe] jwe parents: 3097 diff changeset	248 #endif
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	249
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	250 template <class T>
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	251 MArray<T>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	252 filter (MArray<T>& b, MArray<T>& a, MArray<T>& x, int dim = -1)
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	253 {
5059 38dd53c8935e [project @ 2004-11-02 02:42:25 by jwe] jwe parents: 5056 diff changeset	254 dim_vector x_dims = x.dims();
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	255
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	256 if (dim < 0)
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	257 {
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	258 // Find first non-singleton dimension
5059 38dd53c8935e [project @ 2004-11-02 02:42:25 by jwe] jwe parents: 5056 diff changeset	259 while (dim < x_dims.length () && x_dims(dim) <= 1)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	260 dim++;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11574 diff changeset	261
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	262 // All dimensions singleton, pick first dimension
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	263 if (dim == x_dims.length ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	264 dim = 0;
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	265 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	266 else
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	267 if (dim < 0 \|\| dim > x_dims.length ())
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	268 {
12483 7a5aacf65f81 Rewrite error strings in src/ to use variables named in documentation. Rik <octave@nomad.inbox5.com> parents: 11586 diff changeset	269 error ("filter: DIM must be a valid dimension");
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	270 return MArray<T> ();
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	271 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	272
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	273 octave_idx_type a_len = a.length ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	274 octave_idx_type b_len = b.length ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	275
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5114 diff changeset	276 octave_idx_type si_len = (a_len > b_len ? a_len : b_len) - 1;
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	277 dim_vector si_dims = x.dims ();
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	278 for (int i = dim; i > 0; i--)
5059 38dd53c8935e [project @ 2004-11-02 02:42:25 by jwe] jwe parents: 5056 diff changeset	279 si_dims(i) = si_dims(i-1);
38dd53c8935e [project @ 2004-11-02 02:42:25 by jwe] jwe parents: 5056 diff changeset	280 si_dims(0) = si_len;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11574 diff changeset	281
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	282 MArray<T> si (si_dims, T (0.0));
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	283
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	284 return filter (b, a, x, si, dim);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	285 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	286
3097 d12c312d1cfb [project @ 1997-11-14 08:03:47 by jwe] jwe parents: 3088 diff changeset	287 DEFUN_DLD (filter, args, nargout,
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	288 "-- texinfo --\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10350 diff changeset	289 @deftypefn {Loadable Function} {y =} filter (@var{b}, @var{a}, @var{x})\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	290 @deftypefnx {Loadable Function} {[@var{y}, @var{sf}] =} filter (@var{b}, @var{a}, @var{x}, @var{si})\n\
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	291 @deftypefnx {Loadable Function} {[@var{y}, @var{sf}] =} filter (@var{b}, @var{a}, @var{x}, [], @var{dim})\n\
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	292 @deftypefnx {Loadable Function} {[@var{y}, @var{sf}] =} filter (@var{b}, @var{a}, @var{x}, @var{si}, @var{dim})\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	293 Return the solution to the following linear, time-invariant difference\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	294 equation:\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	295 @tex\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	296 $$\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	297 \\sum_{k=0}^N a_{k+1} y_{n-k} = \\sum_{k=0}^M b_{k+1} x_{n-k}, \\qquad\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	298 1 \\le n \\le P\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	299 $$\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	300 @end tex\n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	301 @ifnottex\n\
10846 a4f482e66b65 Grammarcheck more of the documentation. Rik <octave@nomad.inbox5.com> parents: 10840 diff changeset	302 @c Set example in small font to prevent overfull line\n\
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	303 \n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	304 @smallexample\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10350 diff changeset	305 @group\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	306 N M\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	307 SUM a(k+1) y(n-k) = SUM b(k+1) x(n-k) for 1<=n<=length(x)\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	308 k=0 k=0\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10350 diff changeset	309 @end group\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	310 @end smallexample\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10350 diff changeset	311 \n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	312 @end ifnottex\n\
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	313 \n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	314 @noindent\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	315 where\n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	316 @ifnottex\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	317 N=length(a)-1 and M=length(b)-1.\n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	318 @end ifnottex\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	319 @tex\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	320 $a \\in \\Re^{N-1}$, $b \\in \\Re^{M-1}$, and $x \\in \\Re^P$.\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	321 @end tex\n\
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	322 over the first non-singleton dimension of @var{x} or over @var{dim} if\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8920 diff changeset	323 supplied. An equivalent form of this equation is:\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	324 @tex\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	325 $$\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	326 y_n = -\\sum_{k=1}^N c_{k+1} y_{n-k} + \\sum_{k=0}^M d_{k+1} x_{n-k}, \\qquad\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	327 1 \\le n \\le P\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	328 $$\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	329 @end tex\n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	330 @ifnottex\n\
10846 a4f482e66b65 Grammarcheck more of the documentation. Rik <octave@nomad.inbox5.com> parents: 10840 diff changeset	331 @c Set example in small font to prevent overfull line\n\
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	332 \n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	333 @smallexample\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10350 diff changeset	334 @group\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	335 N M\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	336 y(n) = - SUM c(k+1) y(n-k) + SUM d(k+1) x(n-k) for 1<=n<=length(x)\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	337 k=1 k=0\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10350 diff changeset	338 @end group\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	339 @end smallexample\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10350 diff changeset	340 \n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	341 @end ifnottex\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	342 \n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	343 @noindent\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	344 where\n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	345 @ifnottex\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	346 c = a/a(1) and d = b/a(1).\n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	347 @end ifnottex\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	348 @tex\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	349 $c = a/a_1$ and $d = b/a_1$.\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	350 @end tex\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	351 \n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	352 If the fourth argument @var{si} is provided, it is taken as the\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	353 initial state of the system and the final state is returned as\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	354 @var{sf}. The state vector is a column vector whose length is\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	355 equal to the length of the longest coefficient vector minus one.\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	356 If @var{si} is not supplied, the initial state vector is set to all\n\
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	357 zeros.\n\
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	358 \n\
11553 01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	359 In terms of the Z Transform, y is the result of passing the discrete-\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	360 time signal x through a system characterized by the following rational\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	361 system function:\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	362 @tex\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	363 $$\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	364 H(z) = {\\displaystyle\\sum_{k=0}^M d_{k+1} z^{-k}\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	365 \\over 1 + \\displaystyle\\sum_{k+1}^N c_{k+1} z^{-k}}\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	366 $$\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	367 @end tex\n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	368 @ifnottex\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	369 \n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	370 @example\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8920 diff changeset	371 @group\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	372 M\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	373 SUM d(k+1) z^(-k)\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	374 k=0\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	375 H(z) = ----------------------\n\
0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	376 N\n\
6476 2533c8c41aad [project @ 2007-03-29 13:33:37 by jwe] jwe parents: 5823 diff changeset	377 1 + SUM c(k+1) z^(-k)\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	378 k=1\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8920 diff changeset	379 @end group\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	380 @end example\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10350 diff changeset	381 \n\
8517 81d6ab3ac93c Allow documentation tobe built for other formats than tex and info sh@sh-laptop parents: 7789 diff changeset	382 @end ifnottex\n\
11572 7d6d8c1e471f Grammarcheck Texinfo for files in src directory. Rik <octave@nomad.inbox5.com> parents: 11553 diff changeset	383 @seealso{filter2, fftfilt, freqz}\n\
3367 0748b03c3510 [project @ 1999-11-20 14:52:38 by jwe] jwe parents: 3109 diff changeset	384 @end deftypefn")
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	385 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	386 octave_value_list retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	387
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	388 int nargin = args.length ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	389
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	390 if (nargin < 3 \|\| nargin > 5)
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	391 {
5823 080c08b192d8 [project @ 2006-05-19 05:32:17 by jwe] jwe parents: 5428 diff changeset	392 print_usage ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	393 return retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	394 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	395
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	396 const char *errmsg = "filter: arguments a and b must be vectors";
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	397
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	398 int dim;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	399 dim_vector x_dims = args(2).dims ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	400
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	401 if (nargin == 5)
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	402 {
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	403 dim = args(4).nint_value() - 1;
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	404 if (dim < 0 \|\| dim >= x_dims.length ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	405 {
12483 7a5aacf65f81 Rewrite error strings in src/ to use variables named in documentation. Rik <octave@nomad.inbox5.com> parents: 11586 diff changeset	406 error ("filter: DIM must be a valid dimension");
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	407 return retval;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	408 }
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	409 }
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	410 else
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	411 {
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	412 // Find first non-singleton dimension
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	413 dim = 0;
5059 38dd53c8935e [project @ 2004-11-02 02:42:25 by jwe] jwe parents: 5056 diff changeset	414 while (dim < x_dims.length () && x_dims(dim) <= 1)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	415 dim++;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11574 diff changeset	416
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	417 // All dimensions singleton, pick first dimension
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	418 if (dim == x_dims.length ())
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	419 dim = 0;
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	420 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	421
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	422 bool isfloat = (args(0).is_single_type ()
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	423 \|\| args(1).is_single_type ()
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	424 \|\| args(2).is_single_type ()
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	425 \|\| (nargin >= 4 && args(3).is_single_type ()));
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11574 diff changeset	426
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	427 if (args(0).is_complex_type ()
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	428 \|\| args(1).is_complex_type ()
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	429 \|\| args(2).is_complex_type ()
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4153 diff changeset	430 \|\| (nargin >= 4 && args(3).is_complex_type ()))
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	431 {
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	432 if (isfloat)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	433 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	434 FloatComplexColumnVector b (args(0).float_complex_vector_value ());
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	435 FloatComplexColumnVector a (args(1).float_complex_vector_value ());
3088 f33738510eb2 [project @ 1997-10-02 06:38:36 by jwe] jwe parents: 3014 diff changeset	436
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	437 FloatComplexNDArray x (args(2).float_complex_array_value ());
3088 f33738510eb2 [project @ 1997-10-02 06:38:36 by jwe] jwe parents: 3014 diff changeset	438
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	439 if (! error_state)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	440 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	441 FloatComplexNDArray si;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	442
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	443 if (nargin == 3 \|\| args(3).is_empty ())
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	444 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	445 octave_idx_type a_len = a.length ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	446 octave_idx_type b_len = b.length ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	447
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	448 octave_idx_type si_len = (a_len > b_len ? a_len : b_len) - 1;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	449
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	450 dim_vector si_dims = x.dims ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	451 for (int i = dim; i > 0; i--)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	452 si_dims(i) = si_dims(i-1);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	453 si_dims(0) = si_len;
3088 f33738510eb2 [project @ 1997-10-02 06:38:36 by jwe] jwe parents: 3014 diff changeset	454
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	455 si.resize (si_dims, 0.0);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	456 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	457 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	458 {
12706 7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	459 dim_vector si_dims = args (3).dims ();
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	460 bool si_is_vector = true;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	461 for (int i = 0; i < si_dims.length (); i++)
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	462 if (si_dims(i) != 1 && si_dims(i) < si_dims.numel ())
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	463 {
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	464 si_is_vector = false;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	465 break;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	466 }
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	467
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	468 si = args(3).float_complex_array_value ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	469
12706 7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	470 if (si_is_vector)
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	471 si = si.reshape (dim_vector (1, si.numel ()));
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	472 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	473
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	474 if (! error_state)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	475 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	476 FloatComplexNDArray y (filter (b, a, x, si, dim));
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	477
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	478 if (nargout == 2)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	479 retval(1) = si;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	480
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	481 retval(0) = y;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	482 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	483 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	484 error (errmsg);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	485 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	486 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	487 error (errmsg);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	488 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	489 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	490 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	491 ComplexColumnVector b (args(0).complex_vector_value ());
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	492 ComplexColumnVector a (args(1).complex_vector_value ());
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	493
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	494 ComplexNDArray x (args(2).complex_array_value ());
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	495
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	496 if (! error_state)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	497 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	498 ComplexNDArray si;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	499
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	500 if (nargin == 3 \|\| args(3).is_empty ())
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	501 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	502 octave_idx_type a_len = a.length ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	503 octave_idx_type b_len = b.length ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	504
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	505 octave_idx_type si_len = (a_len > b_len ? a_len : b_len) - 1;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	506
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	507 dim_vector si_dims = x.dims ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	508 for (int i = dim; i > 0; i--)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	509 si_dims(i) = si_dims(i-1);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	510 si_dims(0) = si_len;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	511
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	512 si.resize (si_dims, 0.0);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	513 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	514 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	515 {
12706 7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	516 dim_vector si_dims = args (3).dims ();
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	517 bool si_is_vector = true;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	518 for (int i = 0; i < si_dims.length (); i++)
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	519 if (si_dims(i) != 1 && si_dims(i) < si_dims.numel ())
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	520 {
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	521 si_is_vector = false;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	522 break;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	523 }
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	524
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	525 si = args(3).complex_array_value ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	526
12706 7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	527 if (si_is_vector)
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	528 si = si.reshape (dim_vector (1, si.numel ()));
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	529 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	530
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	531 if (! error_state)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	532 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	533 ComplexNDArray y (filter (b, a, x, si, dim));
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	534
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	535 if (nargout == 2)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	536 retval(1) = si;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	537
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	538 retval(0) = y;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	539 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	540 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	541 error (errmsg);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	542 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	543 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	544 error (errmsg);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	545 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	546 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	547 else
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	548 {
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	549 if (isfloat)
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	550 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	551 FloatColumnVector b (args(0).float_vector_value ());
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	552 FloatColumnVector a (args(1).float_vector_value ());
3088 f33738510eb2 [project @ 1997-10-02 06:38:36 by jwe] jwe parents: 3014 diff changeset	553
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	554 FloatNDArray x (args(2).float_array_value ());
3088 f33738510eb2 [project @ 1997-10-02 06:38:36 by jwe] jwe parents: 3014 diff changeset	555
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	556 if (! error_state)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	557 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	558 FloatNDArray si;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	559
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	560 if (nargin == 3 \|\| args(3).is_empty ())
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	561 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	562 octave_idx_type a_len = a.length ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	563 octave_idx_type b_len = b.length ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	564
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	565 octave_idx_type si_len = (a_len > b_len ? a_len : b_len) - 1;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	566
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	567 dim_vector si_dims = x.dims ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	568 for (int i = dim; i > 0; i--)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	569 si_dims(i) = si_dims(i-1);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	570 si_dims(0) = si_len;
3088 f33738510eb2 [project @ 1997-10-02 06:38:36 by jwe] jwe parents: 3014 diff changeset	571
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	572 si.resize (si_dims, 0.0);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	573 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	574 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	575 {
12706 7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	576 dim_vector si_dims = args (3).dims ();
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	577 bool si_is_vector = true;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	578 for (int i = 0; i < si_dims.length (); i++)
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	579 if (si_dims(i) != 1 && si_dims(i) < si_dims.numel ())
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	580 {
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	581 si_is_vector = false;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	582 break;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	583 }
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	584
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	585 si = args(3).float_array_value ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	586
12706 7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	587 if (si_is_vector)
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	588 si = si.reshape (dim_vector (1, si.numel ()));
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	589 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	590
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	591 if (! error_state)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	592 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	593 FloatNDArray y (filter (b, a, x, si, dim));
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	594
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	595 if (nargout == 2)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	596 retval(1) = si;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	597
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	598 retval(0) = y;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	599 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	600 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	601 error (errmsg);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	602 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	603 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	604 error (errmsg);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	605 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	606 else
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	607 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	608 ColumnVector b (args(0).vector_value ());
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	609 ColumnVector a (args(1).vector_value ());
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	610
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	611 NDArray x (args(2).array_value ());
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	612
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	613 if (! error_state)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	614 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	615 NDArray si;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	616
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	617 if (nargin == 3 \|\| args(3).is_empty ())
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	618 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	619 octave_idx_type a_len = a.length ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	620 octave_idx_type b_len = b.length ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	621
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	622 octave_idx_type si_len = (a_len > b_len ? a_len : b_len) - 1;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	623
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	624 dim_vector si_dims = x.dims ();
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	625 for (int i = dim; i > 0; i--)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	626 si_dims(i) = si_dims(i-1);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	627 si_dims(0) = si_len;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	628
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	629 si.resize (si_dims, 0.0);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	630 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	631 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	632 {
12706 7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	633 dim_vector si_dims = args (3).dims ();
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	634 bool si_is_vector = true;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	635 for (int i = 0; i < si_dims.length (); i++)
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	636 if (si_dims(i) != 1 && si_dims(i) < si_dims.numel ())
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	637 {
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	638 si_is_vector = false;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	639 break;
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	640 }
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	641
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	642 si = args(3).array_value ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	643
12706 7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	644 if (si_is_vector)
7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	645 si = si.reshape (dim_vector (1, si.numel ()));
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	646 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	647
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	648 if (! error_state)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	649 {
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	650 NDArray y (filter (b, a, x, si, dim));
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	651
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	652 if (nargout == 2)
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	653 retval(1) = si;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	654
10154 40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	655 retval(0) = y;
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	656 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	657 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	658 error (errmsg);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	659 }
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	660 else
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	661 error (errmsg);
40dfc0c99116 DLD-FUNCTIONS/.cc: untabify John W. Eaton <jwe@octave.org>* parents: 9758 diff changeset	662 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	663 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	664
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	665 return retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	666 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	667
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	668 template MArray<double>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	669 filter (MArray<double>&, MArray<double>&, MArray<double>&,
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	670 MArray<double>&, int dim);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	671
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	672 template MArray<double>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	673 filter (MArray<double>&, MArray<double>&, MArray<double>&, int dim);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	674
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	675 template MArray<Complex>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	676 filter (MArray<Complex>&, MArray<Complex>&, MArray<Complex>&,
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	677 MArray<Complex>&, int dim);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	678
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	679 template MArray<Complex>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	680 filter (MArray<Complex>&, MArray<Complex>&, MArray<Complex>&, int dim);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	681
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	682 template MArray<float>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	683 filter (MArray<float>&, MArray<float>&, MArray<float>&,
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	684 MArray<float>&, int dim);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	685
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	686 template MArray<float>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	687 filter (MArray<float>&, MArray<float>&, MArray<float>&, int dim);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	688
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	689 template MArray<FloatComplex>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	690 filter (MArray<FloatComplex>&, MArray<FloatComplex>&, MArray<FloatComplex>&,
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	691 MArray<FloatComplex>&, int dim);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	692
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	693 template MArray<FloatComplex>
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10247 diff changeset	694 filter (MArray<FloatComplex>&, MArray<FloatComplex>&, MArray<FloatComplex>&, int dim);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7017 diff changeset	695
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	696 /*
9432 8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	697 %!shared a, b, x, r
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	698 %!test
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	699 %! a = [1 1];
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	700 %! b = [1 1];
12566 3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	701 %! x = zeros (1,10); x(1) = 1;
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	702 %! assert(filter(b, [1], x ), [1 1 0 0 0 0 0 0 0 0]);
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	703 %! assert(filter(b, [1], x.'), [1 1 0 0 0 0 0 0 0 0].');
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	704 %! assert(filter(b.', [1], x ), [1 1 0 0 0 0 0 0 0 0] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	705 %! assert(filter(b.', [1], x.'), [1 1 0 0 0 0 0 0 0 0].');
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	706 %! assert(filter([1], a, x ), [+1 -1 +1 -1 +1 -1 +1 -1 +1 -1] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	707 %! assert(filter([1], a, x.'), [+1 -1 +1 -1 +1 -1 +1 -1 +1 -1].');
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	708 %! assert(filter([1], a.', x ), [+1 -1 +1 -1 +1 -1 +1 -1 +1 -1] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	709 %! assert(filter([1], a.', x.'), [+1 -1 +1 -1 +1 -1 +1 -1 +1 -1].');
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	710 %! assert(filter(b, a, x ), [1 0 0 0 0 0 0 0 0 0] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	711 %! assert(filter(b.', a, x ), [1 0 0 0 0 0 0 0 0 0] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	712 %! assert(filter(b, a.', x ), [1 0 0 0 0 0 0 0 0 0] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	713 %! assert(filter(b.', a, x ), [1 0 0 0 0 0 0 0 0 0] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	714 %! assert(filter(b, a, x.'), [1 0 0 0 0 0 0 0 0 0].');
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	715 %! assert(filter(b.', a, x.'), [1 0 0 0 0 0 0 0 0 0].');
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	716 %! assert(filter(b, a.', x.'), [1 0 0 0 0 0 0 0 0 0].');
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	717 %! assert(filter(b.', a, x.'), [1 0 0 0 0 0 0 0 0 0].');
9432 8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	718 %!
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	719 %!test
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	720 %! r = sqrt(1/2)*(1+i);
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	721 %! a = a*r;
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	722 %! b = b*r;
12566 3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	723 %! assert(filter(b, [1], x ), r*[1 1 0 0 0 0 0 0 0 0] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	724 %! assert(filter(b, [1], rx ), rr*[1 1 0 0 0 0 0 0 0 0] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	725 %! assert(filter(b, [1], x.' ), r*[1 1 0 0 0 0 0 0 0 0].' );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	726 %! assert(filter(b, a, x ), [1 0 0 0 0 0 0 0 0 0] );
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	727 %! assert(filter(b, a, rx ), r[1 0 0 0 0 0 0 0 0 0] );
9432 8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	728 %!
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	729 %!shared a, b, x, y, so
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	730 %!test
12566 3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	731 %! a = [1,1]; b = [1,1];
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	732 %! x = zeros (1,10); x(1) = 1;
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	733 %! [y, so] = filter (b, [1], x, [-1]);
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	734 %! assert(y, [0 1 0 0 0 0 0 0 0 0]);
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	735 %! assert(so,0);
9432 8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	736 %!
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	737 %!test
12566 3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	738 %! x = zeros (10,3); x(1,1)=-1; x(1,2)=1;
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	739 %! y0 = zeros (10,3); y0(1:2,1)=-1; y0(1:2,2)=1;
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	740 %! y = filter (b, [1], x);
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	741 %! assert(y,y0);
9432 8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	742 %!
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	743 %!test
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	744 %! a = [1,1]; b=[1,1];
12566 3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	745 %! x = zeros (4,4,2); x(1,1:4,1) = +1; x(1,1:4,2) = -1;
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	746 %! y0 = zeros (4,4,2); y0(1:2,1:4,1) = +1; y0(1:2,1:4,2) = -1;
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	747 %! y = filter (b, [1], x);
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	748 %! assert(y, y0);
10247 edfceeab945a Add test case for empty argument in filter() Thomas Weber <thomas.weber.mail@gmail.com> parents: 10241 diff changeset	749 %!
12566 3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	750 %!assert(filter (1, ones(10,1)/10, []), []);
3b0eb443c519 Fix orientation of initial conditions vector for filter (bug #32741). Rik <octave@nomad.inbox5.com> parents: 12483 diff changeset	751 %!assert(filter (1, ones(10,1)/10, zeros(0,10)), zeros(0,10));
12706 7c8b50cd4fb4 Backout 12700:9843b3b055e0 from stable branch Jordi Gutiérrez Hermoso <jordigh@gmail.com> parents: 12700 diff changeset	752 %!assert(filter([1, 3], [1], [1 2; 3 4; 5 6], [4, 5]), [5 7; 6 10; 14 18]);
12700 9843b3b055e0 Fix bug with initial conditions in filter (Bug #32741). Rik <octave@nomad.inbox5.com> parents: 12637 diff changeset	753 %!assert(filter (1, ones(10,1)/10, single (1:5)), repmat (single (10), 1, 5));
9432 8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	754 %% Should put some tests of the "DIM" parameter in here.
8cc2d087f3c1 src/DLD-FUNCTIONS/filter.cc: new tests Robert T. Short <octave@phaselockedsystems.com> parents: 9153 diff changeset	755
10155 d0ce5e973937 DLD-FUNCTIONS/.cc: delete Emacs local-variables settings John W. Eaton <jwe@octave.org>* parents: 10154 diff changeset	756 */

Mercurial > octave-nkf

annotate src/DLD-FUNCTIONS/filter.cc @ 12706:7c8b50cd4fb4 stable