octave: src/DLD-FUNCTIONS/max.cc annotate

annotate src/DLD-FUNCTIONS/max.cc @ 9791:6e425f6be618

simplify cummin, cummax docs

author	Jaroslav Hajek <highegg@gmail.com>
date	Tue, 10 Nov 2009 08:58:26 +0100
parents	a5035bc7fbfb
children	b22a2f4b34aa

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
7017 a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe] jwe parents: 7016 diff changeset	3 Copyright (C) 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004,
8920 eb63fbe60fab update copyright notices John W. Eaton <jwe@octave.org> parents: 8780 diff changeset	4 2005, 2006, 2007, 2008 John W. Eaton
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	5 Copyright (C) 2009 VZLU Prague
2928 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 This file is part of Octave.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	8
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	9 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	10 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6711 diff changeset	11 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6711 diff changeset	12 option) any later version.
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	13
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	14 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	15 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	16 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	17 for more details.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	18
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	19 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6711 diff changeset	20 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6711 diff changeset	21 <http://www.gnu.org/licenses/>.
2928 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 */
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	24
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	25 #ifdef HAVE_CONFIG_H
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	26 #include <config.h>
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	27 #endif
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	28
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	29 #include "lo-ieee.h"
3248 68259f410026 [project @ 1999-07-13 03:34:54 by jwe] jwe parents: 3243 diff changeset	30 #include "lo-mappers.h"
7231 2eb392d058bb [project @ 2007-11-30 18:53:29 by jwe] jwe parents: 7189 diff changeset	31 #include "lo-math.h"
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4522 diff changeset	32 #include "dNDArray.h"
9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4522 diff changeset	33 #include "CNDArray.h"
4153 6b96ce9f5743 [project @ 2002-11-06 20:38:49 by jwe] jwe parents: 4102 diff changeset	34 #include "quit.h"
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	35
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	36 #include "defun-dld.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	37 #include "error.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	38 #include "gripes.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
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4522 diff changeset	41 #include "ov-cx-mat.h"
7505 f5005d9510f4 Remove dispatched sparse functions and treat in the generic versions of the functions David Bateman <dbateman@free.fr> parents: 7231 diff changeset	42 #include "ov-re-sparse.h"
f5005d9510f4 Remove dispatched sparse functions and treat in the generic versions of the functions David Bateman <dbateman@free.fr> parents: 7231 diff changeset	43 #include "ov-cx-sparse.h"
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4522 diff changeset	44
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	45 template <class ArrayType>
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	46 static octave_value_list
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	47 do_minmax_red_op (const octave_value& arg,
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	48 int nargout, int dim, bool ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	49 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	50 octave_value_list retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	51 ArrayType array = octave_value_extract<ArrayType> (arg);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	52
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	53 if (error_state)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	54 return retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	55
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	56 if (nargout == 2)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	57 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	58 retval.resize (2);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	59 Array<octave_idx_type> idx;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	60 if (ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	61 retval(0) = array.min (idx, dim);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	62 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	63 retval(0) = array.max (idx, dim);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	64
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	65 retval(1) = octave_value (idx, true, true);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	66 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	67 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	68 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	69 if (ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	70 retval(0) = array.min (dim);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	71 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	72 retval(0) = array.max (dim);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	73 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	74
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	75 return retval;
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	76 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	77
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	78 template <class ArrayType>
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	79 static octave_value
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	80 do_minmax_bin_op (const octave_value& argx, const octave_value& argy,
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	81 bool ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	82 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	83 typedef typename ArrayType::element_type ScalarType;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	84
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	85 octave_value retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	86
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	87 if (argx.is_scalar_type () == 1)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	88 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	89 ScalarType x = octave_value_extract<ScalarType> (argx);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	90 ArrayType y = octave_value_extract<ArrayType> (argy);
7189 e8d953d03f6a [project @ 2007-11-26 20:42:09 by dbateman] dbateman parents: 7017 diff changeset	91
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	92 if (error_state)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	93 ;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	94 else if (ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	95 retval = min (x, y);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	96 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	97 retval = max (x, y);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	98 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	99 else if (argy.is_scalar_type () == 1)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	100 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	101 ArrayType x = octave_value_extract<ArrayType> (argx);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	102 ScalarType y = octave_value_extract<ScalarType> (argy);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	103
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	104 if (error_state)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	105 ;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	106 else if (ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	107 retval = min (x, y);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	108 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	109 retval = max (x, y);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	110 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	111 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	112 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	113 ArrayType x = octave_value_extract<ArrayType> (argx);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	114 ArrayType y = octave_value_extract<ArrayType> (argy);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	115
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	116 if (error_state)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	117 ;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	118 else if (ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	119 retval = min (x, y);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	120 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	121 retval = max (x, y);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	122 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	123
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	124 return retval;
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	125 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	126
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	127 static octave_value_list
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	128 do_minmax_body (const octave_value_list& args,
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	129 int nargout, bool ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	130 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	131 octave_value_list retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	132
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	133 const char *func = ismin ? "min" : "max";
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	134
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	135 int nargin = args.length ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7662 diff changeset	136
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	137 if (nargin == 3 \|\| nargin == 1)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	138 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	139 octave_value arg = args(0);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	140 int dim = -1;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	141 if (nargin == 3)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	142 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	143 dim = args(2).int_value (true) - 1;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	144 if (error_state \|\| dim < 0)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	145 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	146 error ("%s: invalid dimension", func);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	147 return retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	148 }
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	149
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	150 if (! args(1).is_empty ())
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	151 warning ("%s: second argument is ignored");
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	152 }
7189 e8d953d03f6a [project @ 2007-11-26 20:42:09 by dbateman] dbateman parents: 7017 diff changeset	153
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	154 switch (arg.builtin_type ())
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	155 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	156 case btyp_double:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	157 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	158 if (arg.is_range () && (dim == -1 \|\| dim == 1))
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	159 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	160 Range range = arg.range_value ();
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	161 if (range.nelem () == 0)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	162 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	163 retval(0) = arg;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	164 if (nargout > 1)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	165 retval(1) = arg;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	166 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	167 else if (ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	168 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	169 retval(0) = range.min ();
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	170 if (nargout > 1)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	171 retval(1) = static_cast<double> (range.inc () < 0 ? range.nelem () : 1);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	172 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	173 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	174 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	175 retval(0) = range.max ();
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	176 if (nargout > 1)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	177 retval(1) = static_cast<double> (range.inc () >= 0 ? range.nelem () : 1);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	178 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	179 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	180 else if (arg.is_sparse_type ())
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	181 retval = do_minmax_red_op<SparseMatrix> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	182 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	183 retval = do_minmax_red_op<NDArray> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	184 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	185 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	186 case btyp_complex:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	187 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	188 if (arg.is_sparse_type ())
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	189 retval = do_minmax_red_op<SparseComplexMatrix> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	190 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	191 retval = do_minmax_red_op<ComplexNDArray> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	192 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	193 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	194 case btyp_float:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	195 retval = do_minmax_red_op<FloatNDArray> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	196 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	197 case btyp_float_complex:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	198 retval = do_minmax_red_op<FloatComplexNDArray> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	199 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	200 #define MAKE_INT_BRANCH(X) \
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	201 case btyp_ ## X: \
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	202 retval = do_minmax_red_op<X ## NDArray> (arg, nargout, dim, ismin); \
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	203 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	204 MAKE_INT_BRANCH (int8);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	205 MAKE_INT_BRANCH (int16);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	206 MAKE_INT_BRANCH (int32);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	207 MAKE_INT_BRANCH (int64);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	208 MAKE_INT_BRANCH (uint8);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	209 MAKE_INT_BRANCH (uint16);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	210 MAKE_INT_BRANCH (uint32);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	211 MAKE_INT_BRANCH (uint64);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	212 #undef MAKE_INT_BRANCH
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	213 default:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	214 gripe_wrong_type_arg (func, arg);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	215 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	216 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	217 else if (nargin == 2)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	218 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	219 octave_value argx = args(0), argy = args(1);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	220 builtin_type_t xtyp = argx.builtin_type (), ytyp = argy.builtin_type ();
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	221 builtin_type_t rtyp = btyp_mixed_numeric (xtyp, ytyp);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	222
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	223 switch (rtyp)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	224 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	225 case btyp_double:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	226 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	227 if ((argx.is_sparse_type ()
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	228 && (argy.is_sparse_type () \|\| argy.is_scalar_type ()))
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	229 \|\| (argy.is_sparse_type () && argx.is_scalar_type ()))
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	230 retval = do_minmax_bin_op<SparseMatrix> (argx, argy, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	231 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	232 retval = do_minmax_bin_op<NDArray> (argx, argy, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	233 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	234 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	235 case btyp_complex:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	236 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	237 if ((argx.is_sparse_type ()
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	238 && (argy.is_sparse_type () \|\| argy.is_scalar_type ()))
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	239 \|\| (argy.is_sparse_type () && argx.is_scalar_type ()))
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	240 retval = do_minmax_bin_op<SparseComplexMatrix> (argx, argy, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	241 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	242 retval = do_minmax_bin_op<ComplexNDArray> (argx, argy, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	243 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	244 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	245 case btyp_float:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	246 retval = do_minmax_bin_op<FloatNDArray> (argx, argy, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	247 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	248 case btyp_float_complex:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	249 retval = do_minmax_bin_op<FloatComplexNDArray> (argx, argy, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	250 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	251 #define MAKE_INT_BRANCH(X) \
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	252 case btyp_ ## X: \
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	253 retval = do_minmax_bin_op<X ## NDArray> (argx, argy, ismin); \
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	254 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	255 MAKE_INT_BRANCH (int8);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	256 MAKE_INT_BRANCH (int16);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	257 MAKE_INT_BRANCH (int32);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	258 MAKE_INT_BRANCH (int64);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	259 MAKE_INT_BRANCH (uint8);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	260 MAKE_INT_BRANCH (uint16);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	261 MAKE_INT_BRANCH (uint32);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	262 MAKE_INT_BRANCH (uint64);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	263 #undef MAKE_INT_BRANCH
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	264 default:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	265 error ("%s: cannot compute %s (%s, %s)", func, func,
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	266 argx.type_name ().c_str (), argy.type_name ().c_str ());
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	267 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	268 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	269 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	270 print_usage ();
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	271
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	272 return retval;
7505 f5005d9510f4 Remove dispatched sparse functions and treat in the generic versions of the functions David Bateman <dbateman@free.fr> parents: 7231 diff changeset	273 }
f5005d9510f4 Remove dispatched sparse functions and treat in the generic versions of the functions David Bateman <dbateman@free.fr> parents: 7231 diff changeset	274
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	275 DEFUN_DLD (min, args, nargout,
3443 bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	276 "-- texinfo --\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	277 @deftypefn {Loadable Function} {} min (@var{x})\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	278 @deftypefnx {Loadable Function} {} min (@var{x}, @var{y})\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	279 @deftypefnx {Loadable Function} {} min (@var{x}, @var{y}, @var{dim})\n\
8780 ea76466605ba support native cumsum, gripe on overflow in sum/cumsum Jaroslav Hajek <highegg@gmail.com> parents: 8778 diff changeset	280 @deftypefnx {Loadable Function} {[@var{w}, @var{iw}] =} min (@var{x})\n\
3443 bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	281 For a vector argument, return the minimum value. For a matrix\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	282 argument, return the minimum value from each column, as a row\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8949 diff changeset	283 vector, or over the dimension @var{dim} if defined. For two matrices\n\
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4522 diff changeset	284 (or a matrix and scalar), return the pair-wise minimum.\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	285 Thus,\n\
3443 bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	286 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	287 @example\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	288 min (min (@var{x}))\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	289 @end example\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	290 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	291 @noindent\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	292 returns the smallest element of @var{x}, and\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	293 \n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	294 @example\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	295 @group\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	296 min (2:5, pi)\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	297 @result{} 2.0000 3.0000 3.1416 3.1416\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	298 @end group\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	299 @end example\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	300 @noindent\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	301 compares each element of the range @code{2:5} with @code{pi}, and\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	302 returns a row vector of the minimum values.\n\
3443 bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	303 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	304 For complex arguments, the magnitude of the elements are used for\n\
3657 a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe] jwe parents: 3443 diff changeset	305 comparison.\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe] jwe parents: 3443 diff changeset	306 \n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	307 If called with one input and two output arguments,\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	308 @code{min} also returns the first index of the\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8949 diff changeset	309 minimum value(s). Thus,\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	310 \n\
3775 13905c3a24af [project @ 2001-02-06 01:56:59 by jwe] jwe parents: 3747 diff changeset	311 @example\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	312 @group\n\
9165 8c71a86c4bf4 Update section 17.5 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9141 diff changeset	313 [x, ix] = min ([1, 3, 0, 2, 0])\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	314 @result{} x = 0\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	315 ix = 3\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	316 @end group\n\
3657 a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe] jwe parents: 3443 diff changeset	317 @end example\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	318 @seealso{max, cummin, cummax}\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	319 @end deftypefn")
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	320 {
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	321 return do_minmax_body (args, nargout, true);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	322 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	323
7600 24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	324 /*
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	325
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	326 %% test/octave.test/arith/min-1.m
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	327 %!assert (min ([1, 4, 2, 3]) == 1);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	328 %!assert (min ([1; -10; 5; -2]) == -10);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	329
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	330 %% test/octave.test/arith/min-2.m
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	331 %!assert(all (min ([4, i; -2, 2]) == [-2, i]));
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	332
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	333 %% test/octave.test/arith/min-3.m
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	334 %!error <Invalid call to min.*> min ();
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	335
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	336 %% test/octave.test/arith/min-4.m
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	337 %!error <Invalid call to min.*> min (1, 2, 3, 4);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	338
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	339 %!test
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	340 %! x = reshape (1:8,[2,2,2]);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	341 %! assert (max (x,[],1), reshape ([2, 4, 6, 8], [1,2,2]));
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	342 %! assert (max (x,[],2), reshape ([3, 4, 7, 8], [2,1,2]));
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	343 %! [y, i ] = max (x, [], 3);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	344 %! assert (y, [5, 7; 6, 8]);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	345 %! assert (ndims(y), 2);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	346 %! assert (i, [2, 2; 2, 2]);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	347 %! assert (ndims(i), 2);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	348
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	349 */
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	350
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	351 DEFUN_DLD (max, args, nargout,
3443 bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	352 "-- texinfo --\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	353 @deftypefn {Loadable Function} {} max (@var{x})\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	354 @deftypefnx {Loadable Function} {} max (@var{x}, @var{y})\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	355 @deftypefnx {Loadable Function} {} max (@var{x}, @var{y}, @var{dim})\n\
8780 ea76466605ba support native cumsum, gripe on overflow in sum/cumsum Jaroslav Hajek <highegg@gmail.com> parents: 8778 diff changeset	356 @deftypefnx {Loadable Function} {[@var{w}, @var{iw}] =} max (@var{x})\n\
3443 bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	357 For a vector argument, return the maximum value. For a matrix\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	358 argument, return the maximum value from each column, as a row\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8949 diff changeset	359 vector, or over the dimension @var{dim} if defined. For two matrices\n\
4844 9f7ef92b50b0 [project @ 2004-04-02 17:26:53 by jwe] jwe parents: 4522 diff changeset	360 (or a matrix and scalar), return the pair-wise maximum.\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	361 Thus,\n\
3443 bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	362 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	363 @example\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	364 max (max (@var{x}))\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	365 @end example\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	366 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	367 @noindent\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	368 returns the largest element of the matrix @var{x}, and\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	369 \n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	370 @example\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	371 @group\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	372 max (2:5, pi)\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	373 @result{} 3.1416 3.1416 4.0000 5.0000\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	374 @end group\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	375 @end example\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	376 @noindent\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	377 compares each element of the range @code{2:5} with @code{pi}, and\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	378 returns a row vector of the maximum values.\n\
3443 bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	379 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe] jwe parents: 3418 diff changeset	380 For complex arguments, the magnitude of the elements are used for\n\
3775 13905c3a24af [project @ 2001-02-06 01:56:59 by jwe] jwe parents: 3747 diff changeset	381 comparison.\n\
3657 a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe] jwe parents: 3443 diff changeset	382 \n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	383 If called with one input and two output arguments,\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	384 @code{max} also returns the first index of the\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8949 diff changeset	385 maximum value(s). Thus,\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	386 \n\
3775 13905c3a24af [project @ 2001-02-06 01:56:59 by jwe] jwe parents: 3747 diff changeset	387 @example\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	388 @group\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	389 [x, ix] = max ([1, 3, 5, 2, 5])\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	390 @result{} x = 5\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	391 ix = 3\n\
abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	392 @end group\n\
3657 a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe] jwe parents: 3443 diff changeset	393 @end example\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	394 @seealso{min, cummax, cummin}\n\
4522 abdcb14e598d [project @ 2003-09-28 22:04:01 by jwe] jwe parents: 4309 diff changeset	395 @end deftypefn")
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	396 {
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	397 return do_minmax_body (args, nargout, false);
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	398 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	399
7600 24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	400 /*
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	401
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	402 %% test/octave.test/arith/max-1.m
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	403 %!assert (max ([1, 4, 2, 3]) == 4);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	404 %!assert (max ([1; -10; 5; -2]) == 5);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	405
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	406 %% test/octave.test/arith/max-2.m
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	407 %!assert(all (max ([4, i 4.999; -2, 2, 3+4i]) == [4, 2, 3+4i]));
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	408
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	409 %% test/octave.test/arith/max-3.m
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	410 %!error <Invalid call to max.*> max ();
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	411
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	412 %% test/octave.test/arith/max-4.m
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	413 %!error <Invalid call to max.*> max (1, 2, 3, 4);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	414
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	415 %!test
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	416 %! x = reshape (1:8,[2,2,2]);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	417 %! assert (min (x,[],1), reshape ([1, 3, 5, 7], [1,2,2]));
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	418 %! assert (min (x,[],2), reshape ([1, 2, 5, 6], [2,1,2]));
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	419 %! [y, i ] = min (x, [], 3);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	420 %! assert (y, [1, 3; 2, 4]);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	421 %! assert (ndims(y), 2);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	422 %! assert (i, [1, 1; 1, 1]);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	423 %! assert (ndims(i), 2);
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	424
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	425
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	426 */
24abf5a702d9 Chop trailing singletons in min/max functions David Bateman <dbateman@free.fr> parents: 7505 diff changeset	427
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	428 template <class ArrayType>
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	429 static octave_value_list
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	430 do_cumminmax_red_op (const octave_value& arg,
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	431 int nargout, int dim, bool ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	432 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	433 octave_value_list retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	434 ArrayType array = octave_value_extract<ArrayType> (arg);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	435
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	436 if (error_state)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	437 return retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	438
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	439 if (nargout == 2)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	440 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	441 retval.resize (2);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	442 Array<octave_idx_type> idx;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	443 if (ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	444 retval(0) = array.cummin (idx, dim);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	445 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	446 retval(0) = array.cummax (idx, dim);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	447
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	448 retval(1) = octave_value (idx, true, true);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	449 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	450 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	451 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	452 if (ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	453 retval(0) = array.cummin (dim);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	454 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	455 retval(0) = array.cummax (dim);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	456 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	457
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	458 return retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	459 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	460
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	461 static octave_value_list
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	462 do_cumminmax_body (const octave_value_list& args,
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	463 int nargout, bool ismin)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	464 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	465 octave_value_list retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	466
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	467 const char *func = ismin ? "cummin" : "cummax";
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	468
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	469 int nargin = args.length ();
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	470
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	471 if (nargin == 1 \|\| nargin == 2)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	472 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	473 octave_value arg = args(0);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	474 int dim = -1;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	475 if (nargin == 2)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	476 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	477 dim = args(1).int_value (true) - 1;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	478 if (error_state \|\| dim < 0)
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	479 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	480 error ("%s: invalid dimension", func);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	481 return retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	482 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	483 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	484
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	485 switch (arg.builtin_type ())
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	486 {
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	487 case btyp_double:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	488 retval = do_cumminmax_red_op<NDArray> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	489 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	490 case btyp_complex:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	491 retval = do_cumminmax_red_op<ComplexNDArray> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	492 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	493 case btyp_float:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	494 retval = do_cumminmax_red_op<FloatNDArray> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	495 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	496 case btyp_float_complex:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	497 retval = do_cumminmax_red_op<FloatComplexNDArray> (arg, nargout, dim, ismin);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	498 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	499 #define MAKE_INT_BRANCH(X) \
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	500 case btyp_ ## X: \
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	501 retval = do_cumminmax_red_op<X ## NDArray> (arg, nargout, dim, ismin); \
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	502 break;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	503 MAKE_INT_BRANCH (int8);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	504 MAKE_INT_BRANCH (int16);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	505 MAKE_INT_BRANCH (int32);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	506 MAKE_INT_BRANCH (int64);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	507 MAKE_INT_BRANCH (uint8);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	508 MAKE_INT_BRANCH (uint16);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	509 MAKE_INT_BRANCH (uint32);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	510 MAKE_INT_BRANCH (uint64);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	511 #undef MAKE_INT_BRANCH
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	512 default:
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	513 gripe_wrong_type_arg (func, arg);
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	514 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	515 }
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	516 else
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	517 print_usage ();
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	518
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	519 return retval;
a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	520 }
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	521
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	522 DEFUN_DLD (cummin, args, nargout,
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	523 "-- texinfo --\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	524 @deftypefn {Loadable Function} {} cummin (@var{x})\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	525 @deftypefnx {Loadable Function} {} cummin (@var{x}, @var{dim})\n\
8780 ea76466605ba support native cumsum, gripe on overflow in sum/cumsum Jaroslav Hajek <highegg@gmail.com> parents: 8778 diff changeset	526 @deftypefnx {Loadable Function} {[@var{w}, @var{iw}] =} cummin (@var{x})\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	527 Return the cumulative minimum values along dimension @var{dim}. If @var{dim}\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	528 is unspecified it defaults to column-wise operation. For example,\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	529 \n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	530 @example\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	531 @group\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	532 cummin ([5 4 6 2 3 1])\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	533 @result{} 5 4 4 2 2 1\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	534 @end group\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	535 @end example\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	536 \n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	537 \n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	538 The call\n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	539 @example\n\
9791 6e425f6be618 simplify cummin, cummax docs Jaroslav Hajek <highegg@gmail.com> parents: 9790 diff changeset	540 [w, iw] = cummin (x)\n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	541 @end example\n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	542 \n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	543 @noindent\n\
9791 6e425f6be618 simplify cummin, cummax docs Jaroslav Hajek <highegg@gmail.com> parents: 9790 diff changeset	544 with @code{x} a vector, is equivalent to the following code:\n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	545 @example\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8949 diff changeset	546 @group\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	547 w = iw = zeros (size (x));\n\
9791 6e425f6be618 simplify cummin, cummax docs Jaroslav Hajek <highegg@gmail.com> parents: 9790 diff changeset	548 for i = 1:length (x)\n\
6e425f6be618 simplify cummin, cummax docs Jaroslav Hajek <highegg@gmail.com> parents: 9790 diff changeset	549 [w(i), iw(i)] = max (x(1:i));\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	550 endfor\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8949 diff changeset	551 @end group\n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	552 @end example\n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	553 \n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	554 @noindent\n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	555 but computed in a much faster manner.\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	556 @seealso{cummax, min, max}\n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	557 @end deftypefn")
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	558 {
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	559 return do_cumminmax_body (args, nargout, true);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	560 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	561
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	562 DEFUN_DLD (cummax, args, nargout,
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	563 "-- texinfo --\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	564 @deftypefn {Loadable Function} {} cummax (@var{x})\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	565 @deftypefnx {Loadable Function} {} cummax (@var{x}, @var{dim})\n\
8780 ea76466605ba support native cumsum, gripe on overflow in sum/cumsum Jaroslav Hajek <highegg@gmail.com> parents: 8778 diff changeset	566 @deftypefnx {Loadable Function} {[@var{w}, @var{iw}] =} cummax (@var{x})\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	567 Return the cumulative maximum values along dimension @var{dim}. If @var{dim}\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	568 is unspecified it defaults to column-wise operation. For example,\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	569 \n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	570 @example\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	571 @group\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	572 cummax ([1 3 2 6 4 5])\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	573 @result{} 1 3 3 6 6 6\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	574 @end group\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	575 @end example\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	576 \n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	577 The call\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	578 @example\n\
c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	579 [w, iw] = cummax (x, dim)\n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	580 @end example\n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	581 \n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	582 @noindent\n\
9791 6e425f6be618 simplify cummin, cummax docs Jaroslav Hajek <highegg@gmail.com> parents: 9790 diff changeset	583 with @code{x} a vector, is equivalent to the following code:\n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	584 @example\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8949 diff changeset	585 @group\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	586 w = iw = zeros (size (x));\n\
9791 6e425f6be618 simplify cummin, cummax docs Jaroslav Hajek <highegg@gmail.com> parents: 9790 diff changeset	587 for i = 1:length (x)\n\
6e425f6be618 simplify cummin, cummax docs Jaroslav Hajek <highegg@gmail.com> parents: 9790 diff changeset	588 [w(i), iw(i)] = max (x(1:i));\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	589 endfor\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8949 diff changeset	590 @end group\n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	591 @end example\n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	592 \n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	593 @noindent\n\
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	594 but computed in a much faster manner.\n\
9141 c1fff751b5a8 Update section 17.1 (Utility Functions) of arith.txi Rik <rdrider0-list@yahoo.com> parents: 9086 diff changeset	595 @seealso{cummin, max, min}\n\
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	596 @end deftypefn")
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	597 {
9790 a5035bc7fbfb rewrite dispatch part & slightly improve min,max,cummin,cummax Jaroslav Hajek <highegg@gmail.com> parents: 9732 diff changeset	598 return do_cumminmax_body (args, nargout, false);
8777 724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	599 }
724c0f46d9d4 implement cummin/cummax functions Jaroslav Hajek <highegg@gmail.com> parents: 7919 diff changeset	600
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	601 /*
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	602 ;;; Local Variables: ***
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	603 ;;; mode: C++ ***
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	604 ;;; End: ***
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	605 */

Mercurial > octave

annotate src/DLD-FUNCTIONS/max.cc @ 9791:6e425f6be618