octave: src/bitfcns.cc annotate

annotate src/bitfcns.cc @ 11586:12df7854fa7c

strip trailing whitespace from source files

author	John W. Eaton <jwe@octave.org>
date	Thu, 20 Jan 2011 17:24:59 -0500
parents	fd0a3ac60b0e
children	7a5aacf65f81

rev	line source
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	1 /*
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	2
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 10897 diff changeset	3 Copyright (C) 2004-2011 John W. Eaton
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	4
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	5 This file is part of Octave.
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	6
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	7 Octave is free software; you can redistribute it and/or modify it
66af4b9f451e [project @ 2004-06-29 19:20:00 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: 6678 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: 6678 diff changeset	10 option) any later version.
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	11
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	12 Octave is distributed in the hope that it will be useful, but WITHOUT
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	15 for more details.
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	16
66af4b9f451e [project @ 2004-06-29 19:20:00 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: 6678 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: 6678 diff changeset	19 <http://www.gnu.org/licenses/>.
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	20
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	21 */
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	22
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	23 #ifdef HAVE_CONFIG_H
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	24 #include <config.h>
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	25 #endif
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	26
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	27 #include "str-vec.h"
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	28 #include "quit.h"
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	29
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	30 #include "defun.h"
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	31 #include "error.h"
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	32 #include "ov.h"
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	33 #include "ov-uint64.h"
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	34 #include "ov-uint32.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	35 #include "ov-uint16.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	36 #include "ov-uint8.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	37 #include "ov-int64.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	38 #include "ov-int32.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	39 #include "ov-int16.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	40 #include "ov-int8.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	41 #include "ov-scalar.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	42 #include "ov-re-mat.h"
7763 0c6b4c7d7117 Treat bool as a scalar in the bit functions David Bateman <dbateman@free.fr> parents: 7097 diff changeset	43 #include "ov-bool.h"
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	44
5775 ace8d8d26933 [project @ 2006-04-24 19:13:06 by jwe] jwe parents: 5642 diff changeset	45 // FIXME -- could probably eliminate some code duplication by
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	46 // clever use of templates.
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	47
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	48 #define BITOPX(OP, FNAME, RET) \
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	49 { \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	50 int nelx = x.numel (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	51 int nely = y.numel (); \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	52 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	53 bool is_scalar_op = (nelx == 1 \|\| nely == 1); \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	54 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	55 dim_vector dvx = x.dims (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	56 dim_vector dvy = y.dims (); \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	57 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	58 bool is_array_op = (dvx == dvy); \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	59 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	60 if (is_array_op \|\| is_scalar_op) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	61 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	62 RET result; \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	63 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	64 if (nelx != 1) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	65 result.resize (dvx); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	66 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	67 result.resize (dvy); \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	68 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	69 for (int i = 0; i < nelx; i++) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	70 if (is_scalar_op) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	71 for (int k = 0; k < nely; k++) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	72 result(i+k) = x(i) OP y(k); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	73 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	74 result(i) = x(i) OP y(i); \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	75 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	76 retval = result; \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	77 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	78 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	79 error ("%s: size of x and y must match, or one operand must be a scalar", FNAME); \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	80 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	81
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	82 #define BITOP(OP, FNAME) \
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	83 \
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	84 octave_value retval; \
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	85 \
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	86 int nargin = args.length (); \
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	87 \
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	88 if (nargin == 2) \
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	89 { \
4952 bfd57b466752 [project @ 2004-09-01 00:49:05 by jwe] jwe parents: 4920 diff changeset	90 if ((args(0).class_name () == octave_scalar::static_class_name ()) \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	91 \|\| (args(0).class_name () == octave_bool::static_class_name ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	92 \|\| (args(1).class_name () == octave_scalar::static_class_name ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	93 \|\| (args(1).class_name () == octave_bool::static_class_name ())) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	94 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	95 bool arg0_is_int = (args(0).class_name () != \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	96 octave_scalar::static_class_name () && \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	97 args(0).class_name () != \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	98 octave_bool::static_class_name ()); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	99 bool arg1_is_int = (args(1).class_name () != \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	100 octave_scalar::static_class_name () && \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	101 args(1).class_name () != \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	102 octave_bool::static_class_name ()); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	103 \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	104 if (! (arg0_is_int \|\| arg1_is_int)) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	105 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	106 uint64NDArray x (args(0).array_value ()); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	107 uint64NDArray y (args(1).array_value ()); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	108 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	109 BITOPX (OP, FNAME, uint64NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	110 retval = retval.array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	111 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	112 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	113 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	114 int p = (arg0_is_int ? 1 : 0); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	115 int q = (arg0_is_int ? 0 : 1); \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	116 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	117 NDArray dx = args(p).array_value (); \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	118 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	119 if (args(q).type_id () == octave_uint64_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	120 \|\| args(q).type_id () == octave_uint64_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	121 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	122 uint64NDArray x (dx); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	123 uint64NDArray y = args(q).uint64_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	124 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	125 BITOPX (OP, FNAME, uint64NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	126 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	127 else if (args(q).type_id () == octave_uint32_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	128 \|\| args(q).type_id () == octave_uint32_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	129 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	130 uint32NDArray x (dx); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	131 uint32NDArray y = args(q).uint32_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	132 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	133 BITOPX (OP, FNAME, uint32NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	134 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	135 else if (args(q).type_id () == octave_uint16_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	136 \|\| args(q).type_id () == octave_uint16_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	137 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	138 uint16NDArray x (dx); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	139 uint16NDArray y = args(q).uint16_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	140 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	141 BITOPX (OP, FNAME, uint16NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	142 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	143 else if (args(q).type_id () == octave_uint8_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	144 \|\| args(q).type_id () == octave_uint8_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	145 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	146 uint8NDArray x (dx); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	147 uint8NDArray y = args(q).uint8_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	148 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	149 BITOPX (OP, FNAME, uint8NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	150 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	151 else if (args(q).type_id () == octave_int64_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	152 \|\| args(q).type_id () == octave_int64_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	153 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	154 int64NDArray x (dx); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	155 int64NDArray y = args(q).int64_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	156 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	157 BITOPX (OP, FNAME, int64NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	158 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	159 else if (args(q).type_id () == octave_int32_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	160 \|\| args(q).type_id () == octave_int32_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	161 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	162 int32NDArray x (dx); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	163 int32NDArray y = args(q).int32_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	164 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	165 BITOPX (OP, FNAME, int32NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	166 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	167 else if (args(q).type_id () == octave_int16_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	168 \|\| args(q).type_id () == octave_int16_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	169 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	170 int16NDArray x (dx); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	171 int16NDArray y = args(q).int16_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	172 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	173 BITOPX (OP, FNAME, int16NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	174 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	175 else if (args(q).type_id () == octave_int8_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	176 \|\| args(q).type_id () == octave_int8_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	177 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	178 int8NDArray x (dx); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	179 int8NDArray y = args(q).int8_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	180 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	181 BITOPX (OP, FNAME, int8NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	182 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	183 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	184 error ("%s: invalid operand type", FNAME); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	185 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	186 } \
4952 bfd57b466752 [project @ 2004-09-01 00:49:05 by jwe] jwe parents: 4920 diff changeset	187 else if (args(0).class_name () == args(1).class_name ()) \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	188 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	189 if (args(0).type_id () == octave_uint64_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	190 \|\| args(0).type_id () == octave_uint64_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	191 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	192 uint64NDArray x = args(0).uint64_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	193 uint64NDArray y = args(1).uint64_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	194 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	195 BITOPX (OP, FNAME, uint64NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	196 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	197 else if (args(0).type_id () == octave_uint32_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	198 \|\| args(0).type_id () == octave_uint32_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	199 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	200 uint32NDArray x = args(0).uint32_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	201 uint32NDArray y = args(1).uint32_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	202 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	203 BITOPX (OP, FNAME, uint32NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	204 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	205 else if (args(0).type_id () == octave_uint16_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	206 \|\| args(0).type_id () == octave_uint16_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	207 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	208 uint16NDArray x = args(0).uint16_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	209 uint16NDArray y = args(1).uint16_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	210 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	211 BITOPX (OP, FNAME, uint16NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	212 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	213 else if (args(0).type_id () == octave_uint8_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	214 \|\| args(0).type_id () == octave_uint8_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	215 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	216 uint8NDArray x = args(0).uint8_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	217 uint8NDArray y = args(1).uint8_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	218 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	219 BITOPX (OP, FNAME, uint8NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	220 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	221 else if (args(0).type_id () == octave_int64_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	222 \|\| args(0).type_id () == octave_int64_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	223 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	224 int64NDArray x = args(0).int64_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	225 int64NDArray y = args(1).int64_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	226 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	227 BITOPX (OP, FNAME, int64NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	228 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	229 else if (args(0).type_id () == octave_int32_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	230 \|\| args(0).type_id () == octave_int32_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	231 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	232 int32NDArray x = args(0).int32_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	233 int32NDArray y = args(1).int32_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	234 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	235 BITOPX (OP, FNAME, int32NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	236 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	237 else if (args(0).type_id () == octave_int16_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	238 \|\| args(0).type_id () == octave_int16_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	239 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	240 int16NDArray x = args(0).int16_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	241 int16NDArray y = args(1).int16_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	242 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	243 BITOPX (OP, FNAME, int16NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	244 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	245 else if (args(0).type_id () == octave_int8_matrix::static_type_id () \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	246 \|\| args(0).type_id () == octave_int8_scalar::static_type_id ()) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	247 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	248 int8NDArray x = args(0).int8_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	249 int8NDArray y = args(1).int8_array_value (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	250 if (! error_state) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	251 BITOPX (OP, FNAME, int8NDArray); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	252 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	253 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	254 error ("%s: invalid operand type", FNAME); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	255 } \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	256 else \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	257 error ("%s: must have matching operand types", FNAME); \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	258 } \
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	259 else \
5823 080c08b192d8 [project @ 2006-05-19 05:32:17 by jwe] jwe parents: 5775 diff changeset	260 print_usage (); \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	261 \
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	262 return retval
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	263
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	264 DEFUN (bitand, args, ,
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	265 "-- texinfo --\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	266 @deftypefn {Built-in Function} {} bitand (@var{x}, @var{y})\n\
9209 923c7cb7f13f Simplify TeXinfo files by eliminating redundant @iftex followed by @tex construction. Rik <rdrider0-list@yahoo.com> parents: 9064 diff changeset	267 Return the bitwise AND of non-negative integers.\n\
8828 8463d1a2e544 Doc fixes. Brian Gough <bjg@network-theory.co.uk> parents: 8492 diff changeset	268 @var{x}, @var{y} must be in the range [0,bitmax]\n\
5642 2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	269 @seealso{bitor, bitxor, bitset, bitget, bitcmp, bitshift, bitmax}\n\
2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	270 @end deftypefn")
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	271 {
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	272 BITOP (&, "bitand");
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	273 }
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	274
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	275 DEFUN (bitor, args, ,
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	276 "-- texinfo --\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	277 @deftypefn {Built-in Function} {} bitor (@var{x}, @var{y})\n\
9209 923c7cb7f13f Simplify TeXinfo files by eliminating redundant @iftex followed by @tex construction. Rik <rdrider0-list@yahoo.com> parents: 9064 diff changeset	278 Return the bitwise OR of non-negative integers.\n\
8828 8463d1a2e544 Doc fixes. Brian Gough <bjg@network-theory.co.uk> parents: 8492 diff changeset	279 @var{x}, @var{y} must be in the range [0,bitmax]\n\
5642 2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	280 @seealso{bitor, bitxor, bitset, bitget, bitcmp, bitshift, bitmax}\n\
2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	281 @end deftypefn")
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	282 {
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	283 BITOP (\|, "bitor");
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	284 }
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	285
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	286 DEFUN (bitxor, args, ,
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	287 "-- texinfo --\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	288 @deftypefn {Built-in Function} {} bitxor (@var{x}, @var{y})\n\
9209 923c7cb7f13f Simplify TeXinfo files by eliminating redundant @iftex followed by @tex construction. Rik <rdrider0-list@yahoo.com> parents: 9064 diff changeset	289 Return the bitwise XOR of non-negative integers.\n\
8828 8463d1a2e544 Doc fixes. Brian Gough <bjg@network-theory.co.uk> parents: 8492 diff changeset	290 @var{x}, @var{y} must be in the range [0,bitmax]\n\
5642 2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	291 @seealso{bitand, bitor, bitset, bitget, bitcmp, bitshift, bitmax}\n\
2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	292 @end deftypefn")
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	293 {
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	294 BITOP (^, "bitxor");
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	295 }
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	296
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	297 static int64_t
22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	298 bitshift (double a, int n, int64_t mask)
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	299 {
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5828 diff changeset	300 // In the name of bug-for-bug compatibility.
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5828 diff changeset	301 if (a < 0)
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5828 diff changeset	302 return -bitshift (-a, n, mask);
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5828 diff changeset	303
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	304 if (n > 0)
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	305 return (static_cast<int64_t> (a) << n) & mask;
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	306 else if (n < 0)
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	307 return (static_cast<int64_t> (a) >> -n) & mask;
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	308 else
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	309 return static_cast<int64_t> (a) & mask;
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	310 }
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	311
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	312 static int64_t
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	313 bitshift (float a, int n, int64_t mask)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	314 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	315 // In the name of bug-for-bug compatibility.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	316 if (a < 0)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	317 return -bitshift (-a, n, mask);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	318
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	319 if (n > 0)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	320 return (static_cast<int64_t> (a) << n) & mask;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	321 else if (n < 0)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	322 return (static_cast<int64_t> (a) >> -n) & mask;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	323 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	324 return static_cast<int64_t> (a) & mask;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	325 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7763 diff changeset	326
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	327 // Note that the bitshift operators are undefined if shifted by more
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	328 // bits than in the type, so we need to test for the size of the
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	329 // shift.
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	330
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	331 #define DO_BITSHIFT(T) \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	332 if (! error_state) \
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	333 { \
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	334 double d1, d2; \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	335 \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	336 if (n.all_integers (d1, d2)) \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	337 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	338 int m_nel = m.numel (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	339 int n_nel = n.numel (); \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	340 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	341 bool is_scalar_op = (m_nel == 1 \|\| n_nel == 1); \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	342 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	343 dim_vector m_dv = m.dims (); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	344 dim_vector n_dv = n.dims (); \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	345 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	346 bool is_array_op = (m_dv == n_dv); \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	347 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	348 if (is_array_op \|\| is_scalar_op) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	349 { \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	350 T ## NDArray result; \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	351 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	352 if (m_nel != 1) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	353 result.resize (m_dv); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	354 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	355 result.resize (n_dv); \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	356 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	357 for (int i = 0; i < m_nel; i++) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	358 if (is_scalar_op) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	359 for (int k = 0; k < n_nel; k++) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	360 if (static_cast<int> (n(k)) >= bits_in_type) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	361 result(i+k) = 0; \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	362 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	363 result(i+k) = bitshift (m(i), static_cast<int> (n(k)), mask); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	364 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	365 if (static_cast<int> (n(i)) >= bits_in_type) \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	366 result(i) = 0; \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	367 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	368 result(i) = bitshift (m(i), static_cast<int> (n(i)), mask); \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	369 \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	370 retval = result; \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	371 } \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	372 else \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	373 error ("bitshift: size of A and N must match, or one operand must be a scalar"); \
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	374 } \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	375 else \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	376 error ("bitshift: expecting second argument to be integer"); \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	377 }
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	378
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	379 #define DO_UBITSHIFT(T, N) \
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	380 do \
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	381 { \
4920 b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	382 int bits_in_type = octave_ ## T :: nbits (); \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	383 T ## NDArray m = m_arg.T ## _array_value (); \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	384 octave_ ## T mask = octave_ ## T::max (); \
4920 b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	385 if ((N) < bits_in_type) \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	386 mask = bitshift (mask, (N) - bits_in_type); \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	387 else if ((N) < 1) \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	388 mask = 0; \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	389 DO_BITSHIFT (T); \
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	390 } \
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	391 while (0)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	392
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	393 #define DO_SBITSHIFT(T, N) \
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	394 do \
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	395 { \
4920 b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	396 int bits_in_type = octave_ ## T :: nbits (); \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	397 T ## NDArray m = m_arg.T ## _array_value (); \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	398 octave_ ## T mask = octave_ ## T::max (); \
4920 b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	399 if ((N) < bits_in_type) \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	400 mask = bitshift (mask, (N) - bits_in_type); \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	401 else if ((N) < 1) \
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	402 mask = 0; \
9440 357cff83985d fix signed integer shift Jaroslav Hajek <highegg@gmail.com> parents: 9209 diff changeset	403 mask = mask \| octave_ ## T :: min (); /* FIXME: 2's complement only? */ \
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	404 DO_BITSHIFT (T); \
a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	405 } \
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	406 while (0)
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	407
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	408 DEFUN (bitshift, args, ,
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	409 "-- texinfo --\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	410 @deftypefn {Built-in Function} {} bitshift (@var{a}, @var{k})\n\
6678 49724abe1236 [project @ 2007-05-31 19:44:45 by jwe] jwe parents: 6439 diff changeset	411 @deftypefnx {Built-in Function} {} bitshift (@var{a}, @var{k}, @var{n})\n\
8492 9e18aff77e9e [docs] remove Brian Gough <bjg@gnu.org> parents: 8347 diff changeset	412 Return a @var{k} bit shift of @var{n}-digit unsigned\n\
4920 b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	413 integers in @var{a}. A positive @var{k} leads to a left shift.\n\
b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	414 A negative value to a right shift. If @var{n} is omitted it defaults\n\
b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	415 to log2(bitmax)+1.\n\
8347 fa78cb8d8a5c corrections for typos Brian Gough<bjg@network-theory.co.uk> parents: 8325 diff changeset	416 @var{n} must be in the range [1,log2(bitmax)+1] usually [1,33]\n\
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	417 \n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	418 @example\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8920 diff changeset	419 @group\n\
7097 e8c94e473c68 [project @ 2007-11-02 17:24:23 by jwe] jwe parents: 7017 diff changeset	420 bitshift (eye (3), 1)\n\
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	421 @result{}\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	422 @group\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	423 2 0 0\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	424 0 2 0\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	425 0 0 2\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	426 @end group\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	427 \n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	428 bitshift (10, [-2, -1, 0, 1, 2])\n\
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	429 @result{} 2 5 10 20 40\n\
6439 a37c427ae9d1 [project @ 2007-03-23 15:04:01 by jwe] jwe parents: 6108 diff changeset	430 @c FIXME -- restore this example when third arg is allowed to be an array.\n\
a37c427ae9d1 [project @ 2007-03-23 15:04:01 by jwe] jwe parents: 6108 diff changeset	431 @c \n\
a37c427ae9d1 [project @ 2007-03-23 15:04:01 by jwe] jwe parents: 6108 diff changeset	432 @c \n\
a37c427ae9d1 [project @ 2007-03-23 15:04:01 by jwe] jwe parents: 6108 diff changeset	433 @c bitshift ([1, 10], 2, [3,4])\n\
a37c427ae9d1 [project @ 2007-03-23 15:04:01 by jwe] jwe parents: 6108 diff changeset	434 @c @result{} 4 8\n\
9064 7c02ec148a3c Check grammar on all .cc files Rik <rdrider0-list@yahoo.com> parents: 8920 diff changeset	435 @end group\n\
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	436 @end example\n\
5642 2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	437 @seealso{bitand, bitor, bitxor, bitset, bitget, bitcmp, bitmax}\n\
2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	438 @end deftypefn")
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	439 {
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	440 octave_value retval;
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	441
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	442 int nargin = args.length ();
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	443
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	444 if (nargin == 2 \|\| nargin == 3)
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	445 {
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	446 int nbits = 64;
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	447
4920 b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	448 NDArray n = args(1).array_value ();
b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	449
b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	450 if (error_state)
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	451 error ("bitshift: expecting integer as second argument");
4920 b22a7a1db0d5 [project @ 2004-07-27 18:15:25 by jwe] jwe parents: 4919 diff changeset	452 else
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	453 {
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	454 if (nargin == 3)
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	455 {
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	456 // FIXME -- for compatibility, we should accept an array
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	457 // or a scalar as the third argument.
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	458 if (args(2).numel () > 1)
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	459 error ("bitshift: expecting scalar integer as third argument");
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	460 else
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	461 {
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	462 nbits = args(2).int_value ();
11586 12df7854fa7c strip trailing whitespace from source files John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	463
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	464 if (error_state)
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	465 error ("bitshift: expecting integer as third argument");
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	466 else if (nbits < 0)
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	467 error ("bitshift: number of bits to mask must be positive");
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	468 }
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	469 }
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	470 }
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	471
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	472 if (error_state)
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	473 return retval;
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	474
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	475 octave_value m_arg = args(0);
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	476 std::string cname = m_arg.class_name ();
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	477
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	478 if (cname == "uint8")
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	479 DO_UBITSHIFT (uint8, nbits < 8 ? nbits : 8);
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	480 else if (cname == "uint16")
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	481 DO_UBITSHIFT (uint16, nbits < 16 ? nbits : 16);
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	482 else if (cname == "uint32")
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	483 DO_UBITSHIFT (uint32, nbits < 32 ? nbits : 32);
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	484 else if (cname == "uint64")
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	485 DO_UBITSHIFT (uint64, nbits < 64 ? nbits : 64);
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	486 else if (cname == "int8")
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	487 DO_SBITSHIFT (int8, nbits < 8 ? nbits : 8);
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	488 else if (cname == "int16")
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	489 DO_SBITSHIFT (int16, nbits < 16 ? nbits : 16);
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	490 else if (cname == "int32")
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	491 DO_SBITSHIFT (int32, nbits < 32 ? nbits : 32);
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	492 else if (cname == "int64")
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	493 DO_SBITSHIFT (int64, nbits < 64 ? nbits : 64);
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	494 else if (cname == "double")
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	495 {
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	496 nbits = (nbits < 53 ? nbits : 53);
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	497 int64_t mask = 0x1FFFFFFFFFFFFFLL;
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	498 if (nbits < 53)
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	499 mask = mask >> (53 - nbits);
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	500 else if (nbits < 1)
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	501 mask = 0;
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	502 int bits_in_type = 64;
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	503 NDArray m = m_arg.array_value ();
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	504 DO_BITSHIFT ( );
57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	505 }
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	506 else
10315 57a59eae83cc untabify src C++ source files John W. Eaton <jwe@octave.org> parents: 10160 diff changeset	507 error ("bitshift: not defined for %s objects", cname.c_str ());
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	508 }
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	509 else
5823 080c08b192d8 [project @ 2006-05-19 05:32:17 by jwe] jwe parents: 5775 diff changeset	510 print_usage ();
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	511
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	512 return retval;
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	513 }
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	514
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	515 DEFUN (bitmax, args, ,
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	516 "-- texinfo --\n\
10897 dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	517 @deftypefn {Built-in Function} {} bitmax ()\n\
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	518 @deftypefnx {Built-in Function} {} bitmax (\"double\")\n\
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	519 @deftypefnx {Built-in Function} {} bitmax (\"single\")\n\
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	520 Return the largest integer that can be represented within a floating point\n\
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	521 value. The default class is \"double\", but \"single\" is a valid option. \n\
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	522 On IEEE-754 compatible systems, @code{bitmax} is @w{@math{2^{53} - 1}}.\n\
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	523 @end deftypefn")
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	524 {
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	525 octave_value retval;
10897 dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	526 std::string cname = "double";
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	527 int nargin = args.length ();
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	528
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	529 if (nargin == 1 && args(0).is_string ())
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	530 cname = args(0).string_value ();
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	531 else if (nargin != 0)
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	532 {
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	533 print_usage ();
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	534 return retval;
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	535 }
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	536
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	537 if (cname == "double")
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	538 retval = (static_cast<double> (0x1FFFFFFFFFFFFFLL));
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	539 else if (cname == "single")
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	540 retval = (static_cast<double> (0xFFFFFFL));
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	541 else
10897 dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	542 error ("bitmax: not defined for class '%s'", cname.c_str ());
dbec9d590756 Add CLASS argument to bitmax Rik <octave@nomad.inbox5.com> parents: 10895 diff changeset	543
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	544 return retval;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	545 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	546
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	547 DEFUN (intmax, args, ,
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	548 "-- texinfo --\n\
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	549 @deftypefn {Built-in Function} {} intmax (@var{type})\n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	550 Return the largest integer that can be represented in an integer type.\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	551 The variable @var{type} can be\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	552 \n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	553 @table @code\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	554 @item int8\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	555 signed 8-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	556 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	557 @item int16\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	558 signed 16-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	559 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	560 @item int32\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	561 signed 32-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	562 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	563 @item int64\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	564 signed 64-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	565 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	566 @item uint8\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	567 unsigned 8-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	568 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	569 @item uint16\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	570 unsigned 16-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	571 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	572 @item uint32\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	573 unsigned 32-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	574 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	575 @item uint64\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	576 unsigned 64-bit integer.\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	577 @end table\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	578 \n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	579 The default for @var{type} is @code{uint32}.\n\
5642 2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	580 @seealso{intmin, bitmax}\n\
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	581 @end deftypefn")
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	582 {
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	583 octave_value retval;
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	584 std::string cname = "int32";
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	585 int nargin = args.length ();
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	586
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	587 if (nargin == 1 && args(0).is_string ())
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	588 cname = args(0).string_value ();
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	589 else if (nargin != 0)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	590 {
5823 080c08b192d8 [project @ 2006-05-19 05:32:17 by jwe] jwe parents: 5775 diff changeset	591 print_usage ();
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	592 return retval;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	593 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	594
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	595 if (cname == "uint8")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	596 retval = octave_uint8 (std::numeric_limits<uint8_t>::max ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	597 else if (cname == "uint16")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	598 retval = octave_uint16 (std::numeric_limits<uint16_t>::max ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	599 else if (cname == "uint32")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	600 retval = octave_uint32 (std::numeric_limits<uint32_t>::max ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	601 else if (cname == "uint64")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	602 retval = octave_uint64 (std::numeric_limits<uint64_t>::max ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	603 else if (cname == "int8")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	604 retval = octave_int8 (std::numeric_limits<int8_t>::max ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	605 else if (cname == "int16")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	606 retval = octave_int16 (std::numeric_limits<int16_t>::max ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	607 else if (cname == "int32")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	608 retval = octave_int32 (std::numeric_limits<int32_t>::max ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	609 else if (cname == "int64")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	610 retval = octave_int64 (std::numeric_limits<int64_t>::max ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	611 else
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	612 error ("intmax: not defined for '%s' objects", cname.c_str ());
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	613
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	614 return retval;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	615 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	616
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	617 DEFUN (intmin, args, ,
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	618 "-- texinfo --\n\
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	619 @deftypefn {Built-in Function} {} intmin (@var{type})\n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	620 Return the smallest integer that can be represented in an integer type.\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	621 The variable @var{type} can be\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	622 \n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	623 @table @code\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	624 @item int8\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	625 signed 8-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	626 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	627 @item int16\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	628 signed 16-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	629 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	630 @item int32\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	631 signed 32-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	632 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	633 @item int64\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	634 signed 64-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	635 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	636 @item uint8\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	637 unsigned 8-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	638 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	639 @item uint16\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	640 unsigned 16-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	641 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	642 @item uint32\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	643 unsigned 32-bit integer.\n\
10840 89f4d7e294cc Grammarcheck .cc files Rik <octave@nomad.inbox5.com> parents: 10830 diff changeset	644 \n\
5040 b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	645 @item uint64\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	646 unsigned 64-bit integer.\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	647 @end table\n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	648 \n\
b6ff3a0d3d17 [project @ 2004-09-27 13:50:01 by jwe] jwe parents: 4952 diff changeset	649 The default for @var{type} is @code{uint32}.\n\
5642 2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe] jwe parents: 5307 diff changeset	650 @seealso{intmax, bitmax}\n\
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	651 @end deftypefn")
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	652 {
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	653 octave_value retval;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	654 std::string cname = "int32";
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	655 int nargin = args.length ();
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	656
4919 a1073eef650c [project @ 2004-07-27 15:22:50 by jwe] jwe parents: 4915 diff changeset	657 if (nargin == 1 && args(0).is_string ())
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	658 cname = args(0).string_value ();
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	659 else if (nargin != 0)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	660 {
5823 080c08b192d8 [project @ 2006-05-19 05:32:17 by jwe] jwe parents: 5775 diff changeset	661 print_usage ();
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	662 return retval;
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	663 }
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	664
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	665 if (cname == "uint8")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	666 retval = octave_uint8 (std::numeric_limits<uint8_t>::min ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	667 else if (cname == "uint16")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	668 retval = octave_uint16 (std::numeric_limits<uint16_t>::min());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	669 else if (cname == "uint32")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	670 retval = octave_uint32 (std::numeric_limits<uint32_t>::min ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	671 else if (cname == "uint64")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	672 retval = octave_uint64 (std::numeric_limits<uint64_t>::min ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	673 else if (cname == "int8")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	674 retval = octave_int8 (std::numeric_limits<int8_t>::min ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	675 else if (cname == "int16")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	676 retval = octave_int16 (std::numeric_limits<int16_t>::min ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	677 else if (cname == "int32")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	678 retval = octave_int32 (std::numeric_limits<int32_t>::min ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	679 else if (cname == "int64")
5828 22e23bee74c8 [project @ 2006-05-23 06:05:14 by jwe] jwe parents: 5823 diff changeset	680 retval = octave_int64 (std::numeric_limits<int64_t>::min ());
4915 c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	681 else
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	682 error ("intmin: not defined for '%s' objects", cname.c_str ());
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe] jwe parents: 4908 diff changeset	683
4908 66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	684 return retval;
66af4b9f451e [project @ 2004-06-29 19:20:00 by jwe] jwe parents: diff changeset	685 }
10811 e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	686
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	687 DEFUN (sizemax, args, ,
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	688 "-- texinfo --\n\
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	689 @deftypefn {Built-in Function} {} sizemax ()\n\
10895 4176c5c62138 Imrove documentation string for sizemax Rik <octave@nomad.inbox5.com> parents: 10840 diff changeset	690 Return the largest value allowed for the size of an array.\n\
10811 e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	691 If Octave is compiled with 64-bit indexing, the result is of class int64,\n\
10895 4176c5c62138 Imrove documentation string for sizemax Rik <octave@nomad.inbox5.com> parents: 10840 diff changeset	692 otherwise it is of class int32. The maximum array size is slightly\n\
4176c5c62138 Imrove documentation string for sizemax Rik <octave@nomad.inbox5.com> parents: 10840 diff changeset	693 smaller than the maximum value allowable for the relevant class as reported\n\
4176c5c62138 Imrove documentation string for sizemax Rik <octave@nomad.inbox5.com> parents: 10840 diff changeset	694 by @code{intmax}.\n\
10811 e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	695 @seealso{intmax}\n\
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	696 @end deftypefn")
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	697 {
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	698 octave_value retval;
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	699
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	700 if (args.length () == 0)
10830 b4ebfd675321 avoid static initialization disaster in dim_vector Jaroslav Hajek <highegg@gmail.com> parents: 10811 diff changeset	701 retval = octave_int<octave_idx_type> (dim_vector::dim_max ());
10811 e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	702 else
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	703 print_usage ();
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	704
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	705 return retval;
e38c071bbc41 allow user query the maximum array size Jaroslav Hajek <highegg@gmail.com> parents: 10315 diff changeset	706 }

Mercurial > octave

annotate src/bitfcns.cc @ 11586:12df7854fa7c