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

annotate src/DLD-FUNCTIONS/typecast.cc @ 9671:e4b383a0effa

merge reinterpret.cc into typecast.cc

author	Jaroslav Hajek <highegg@gmail.com>
date	Tue, 29 Sep 2009 09:13:24 +0200
parents	58604c45ca74
children	34d6f005db4b

rev	line source
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	1 /*
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	2
8920 eb63fbe60fab update copyright notices John W. Eaton <jwe@octave.org> parents: 8715 diff changeset	3 Copyright (C) 2007, 2008, 2009 David Bateman
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	4 Copyright (C) 2009 VZLU Prague
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	5
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	6 This file is part of Octave.
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	7
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	8 Octave is free software; you can redistribute it and/or modify it
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	9 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: 6948 diff changeset	10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6948 diff changeset	11 option) any later version.
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	12
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	13 Octave is distributed in the hope that it will be useful, but WITHOUT
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	16 for more details.
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	17
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	18 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: 6948 diff changeset	19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6948 diff changeset	20 <http://www.gnu.org/licenses/>.
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	21
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	22 */
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	23
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	24 #ifdef HAVE_CONFIG_H
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	25 #include <config.h>
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	26 #endif
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	27
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	28 #include "mx-base.h"
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	29
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	30 #include "defun-dld.h"
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	31 #include "error.h"
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	32 #include "gripes.h"
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	33 #include "oct-obj.h"
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	34 #include "unwind-prot.h"
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	35
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	36 static dim_vector
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	37 get_vec_dims (const dim_vector& old_dims, octave_idx_type n)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	38 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	39 if (old_dims.length () == 2 && old_dims(0) == 1)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	40 return dim_vector (1, n);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	41 else if (old_dims.length () == 2 && old_dims (0) == 0 && old_dims (1) == 0)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	42 return dim_vector ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	43 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	44 return dim_vector (n, 1);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	45 }
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	46
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	47 template <class ArrayType>
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	48 static void
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	49 get_data_and_bytesize (const ArrayType& array,
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	50 const void *& data,
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	51 octave_idx_type& byte_size,
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	52 dim_vector& old_dims)
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	53 {
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	54 // The array given may be a temporary, constructed from a scalar or sparse
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	55 // array. This will ensure the data will be deallocated after we exit.
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	56 unwind_protect::add_delete (new ArrayType (array));
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	57
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	58 data = reinterpret_cast<const void *> (array.data ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	59 byte_size = array.byte_size ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	60
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	61 old_dims = array.dims ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	62 }
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	63
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	64 template <class ArrayType>
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	65 static ArrayType
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	66 reinterpret_copy (const void *data, octave_idx_type byte_size,
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	67 const dim_vector& old_dims)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	68 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	69 typedef typename ArrayType::element_type T;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	70 octave_idx_type n = byte_size / sizeof (T);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	71
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	72 if (n * static_cast<int> (sizeof (T)) == byte_size)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	73 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	74 ArrayType retval (get_vec_dims (old_dims, n));
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	75 T *dest = retval.fortran_vec ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	76 std::memcpy (dest, data, n * sizeof (T));
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	77
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	78 return retval;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	79 }
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	80 else
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	81 {
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	82 error ("typecast: incorrect number of input values to make output value");
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	83 return ArrayType ();
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	84 }
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	85 }
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	86
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	87
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	88 DEFUN_DLD (typecast, args, ,
8715 954b6f69f51d doc fixes John W. Eaton <jwe@octave.org> parents: 8100 diff changeset	89 "-- texinfo --\n\
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	90 @deftypefn {Loadable Function} {} typecast (@var{x}, @var{class})\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	91 Returns a new array @var{y} resulting from interpreting the data of\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	92 @var{x} in memory as data of the numeric class @var{class}. Both the class\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	93 of @var{x} and @var{class} must be one of the built-in numeric classes:\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	94 \n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	95 @example\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	96 \"logical\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	97 \"char\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	98 \"int8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	99 \"int16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	100 \"int32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	101 \"int64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	102 \"uint8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	103 \"uint16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	104 \"uint32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	105 \"uint64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	106 \"double\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	107 \"single\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	108 \"double complex\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	109 \"single complex\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	110 @end example\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	111 \n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	112 the last two are reserved for @var{class}; they indicate that a complex-valued result\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	113 is requested. Complex arrays are stored in memory as consecutive pairs of real numbers.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	114 The sizes of integer types are given by their bit counts. Both logical and char are typically\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	115 one byte wide; however, this is not guaranteed by C++. If your system is IEEE conformant,\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	116 single and double should be 4 bytes and 8 bytes wide, respectively.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	117 \"logical\" is not allowed for @var{class}.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	118 If the input is a row vector, the return value is a row vector, otherwise it is a column vector.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	119 If the bit length of @var{x} is not divisible by that of @var{class}, an error occurs.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	120 @seealso{bitunpack,bitpack}\n\
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	121 \n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	122 An example of the use of typecast on a little-endian machine is\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	123 \n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	124 @example\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	125 @group\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	126 @var{x} = uint16 ([1, 65535]);\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	127 typecast (@var{x}, 'uint8')\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	128 @result{} [ 0, 1, 255, 255]\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	129 @end group\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	130 @end example\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	131 @seealso{cast, swapbytes}\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	132 @end deftypefn")
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	133 {
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	134 octave_value retval;
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	135
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	136 if (args.length () == 2)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	137 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	138 const void *data = 0;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	139 octave_idx_type byte_size = 0;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	140 dim_vector old_dims;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	141
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	142 octave_value array = args(0);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	143
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	144 if (array.is_bool_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	145 get_data_and_bytesize (array.bool_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	146 else if (array.is_string ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	147 get_data_and_bytesize (array.char_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	148 else if (array.is_integer_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	149 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	150 if (array.is_int8_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	151 get_data_and_bytesize (array.int8_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	152 else if (array.is_int16_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	153 get_data_and_bytesize (array.int16_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	154 else if (array.is_int32_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	155 get_data_and_bytesize (array.int32_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	156 else if (array.is_int64_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	157 get_data_and_bytesize (array.int64_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	158 else if (array.is_uint8_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	159 get_data_and_bytesize (array.uint8_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	160 else if (array.is_uint16_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	161 get_data_and_bytesize (array.uint16_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	162 else if (array.is_uint32_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	163 get_data_and_bytesize (array.uint32_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	164 else if (array.is_uint64_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	165 get_data_and_bytesize (array.uint64_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	166 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	167 assert (0);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	168 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	169 else if (array.is_complex_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	170 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	171 if (array.is_single_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	172 get_data_and_bytesize (array.float_complex_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	173 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	174 get_data_and_bytesize (array.complex_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	175 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	176 else if (array.is_real_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	177 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	178 if (array.is_single_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	179 get_data_and_bytesize (array.float_array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	180 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	181 get_data_and_bytesize (array.array_value (), data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	182 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	183 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	184 error ("typecast: invalid input class: %s", array.class_name ().c_str ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	185
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	186 std::string numclass = args(1).string_value ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	187
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	188 if (error_state \|\| numclass.size () == 0)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	189 ;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	190 else if (numclass == "char")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	191 retval = octave_value (reinterpret_copy<charNDArray> (data, byte_size, old_dims), true);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	192 else if (numclass[0] == 'i')
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	193 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	194 if (numclass == "int8")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	195 retval = reinterpret_copy<int8NDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	196 else if (numclass == "int16")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	197 retval = reinterpret_copy<int16NDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	198 else if (numclass == "int32")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	199 retval = reinterpret_copy<int32NDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	200 else if (numclass == "int64")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	201 retval = reinterpret_copy<int64NDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	202 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	203 else if (numclass[0] == 'u')
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	204 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	205 if (numclass == "uint8")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	206 retval = reinterpret_copy<uint8NDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	207 else if (numclass == "uint16")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	208 retval = reinterpret_copy<uint16NDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	209 else if (numclass == "uint32")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	210 retval = reinterpret_copy<uint32NDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	211 else if (numclass == "uint64")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	212 retval = reinterpret_copy<uint64NDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	213 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	214 else if (numclass == "single")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	215 retval = reinterpret_copy<FloatNDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	216 else if (numclass == "double")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	217 retval = reinterpret_copy<NDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	218 else if (numclass == "single complex")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	219 retval = reinterpret_copy<FloatComplexNDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	220 else if (numclass == "double complex")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	221 retval = reinterpret_copy<ComplexNDArray> (data, byte_size, old_dims);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	222
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	223 if (! error_state && retval.is_undefined ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	224 error ("typecast: cannot convert to %s class", numclass.c_str ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	225 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	226 else
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	227 print_usage ();
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	228
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	229 return retval;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	230 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	231
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	232 template <class ArrayType>
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	233 ArrayType
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	234 do_bitpack (const boolNDArray& bitp)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	235 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	236 typedef typename ArrayType::element_type T;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	237 octave_idx_type n = bitp.numel () / (sizeof (T) * CHAR_BIT);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	238
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	239 if (n * static_cast<int> (sizeof (T)) * CHAR_BIT == bitp.numel ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	240 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	241
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	242 ArrayType retval (get_vec_dims (bitp.dims (), n));
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	243
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	244 const bool *bits = bitp.fortran_vec ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	245 char packed = reinterpret_cast<char > (retval.fortran_vec ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	246
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	247 octave_idx_type m = n * sizeof (T);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	248
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	249 for (octave_idx_type i = 0; i < m; i++)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	250 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	251 char c = bits[0];
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	252 for (int j = 1; j < CHAR_BIT; j++)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	253 c \|= bits[j] << j;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	254
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	255 packed[i] = c;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	256 bits += CHAR_BIT;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	257 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	258
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	259 return retval;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	260 }
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	261 else
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	262 {
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	263 error ("bitpack: incorrect number of bits to make up output value");
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	264 return ArrayType ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	265 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	266 }
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	267
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	268 DEFUN_DLD (bitpack, args, ,
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	269 "-- texinfo --\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	270 @deftypefn {Loadable Function} {@var{y} =} bitpack (@var{x}, @var{class})\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	271 Returns a new array @var{y} resulting from interpreting a logical array @var{x}\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	272 as raw bit pattern for data of the numeric class @var{class}. @var{class} must be\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	273 one of the built-in numeric classes:\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	274 \n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	275 @example\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	276 \"char\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	277 \"int8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	278 \"int16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	279 \"int32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	280 \"int64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	281 \"uint8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	282 \"uint16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	283 \"uint32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	284 \"uint64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	285 \"double\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	286 \"single\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	287 @end example\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	288 \n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	289 The number of elements of @var{x} should be divisible by the bit length of @var{class}.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	290 If it is not, excess bits are discarded. Bits come in increasing order of significance, i.e.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	291 @code{x(1)} is bit 0, @code{x(2)} is bit 1, etc.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	292 The result is a row vector if @var{x} is a row vector, otherwise it is a column vector.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	293 @seealso{bitunpack,typecast,bitget,bitset}\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	294 @end deftypefn")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	295 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	296 octave_value retval;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	297
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	298 if (args.length () == 2 && args(0).is_bool_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	299 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	300 boolNDArray bitp = args(0).bool_array_value ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	301
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	302 std::string numclass = args(1).string_value ();
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	303
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	304 if (error_state \|\| numclass.size () == 0)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	305 ;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	306 else if (numclass == "char")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	307 retval = octave_value (do_bitpack<charNDArray> (bitp), true);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	308 else if (numclass[0] == 'i')
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	309 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	310 if (numclass == "int8")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	311 retval = do_bitpack<int8NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	312 else if (numclass == "int16")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	313 retval = do_bitpack<int16NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	314 else if (numclass == "int32")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	315 retval = do_bitpack<int32NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	316 else if (numclass == "int64")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	317 retval = do_bitpack<int64NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	318 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	319 else if (numclass[0] == 'u')
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	320 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	321 if (numclass == "uint8")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	322 retval = do_bitpack<uint8NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	323 else if (numclass == "uint16")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	324 retval = do_bitpack<uint16NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	325 else if (numclass == "uint32")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	326 retval = do_bitpack<uint32NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	327 else if (numclass == "uint64")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	328 retval = do_bitpack<uint64NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	329 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	330 else if (numclass == "single")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	331 retval = do_bitpack<FloatNDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	332 else if (numclass == "double")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	333 retval = do_bitpack<NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	334 else if (numclass == "single complex")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	335 retval = do_bitpack<FloatComplexNDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	336 else if (numclass == "double complex")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	337 retval = do_bitpack<ComplexNDArray> (bitp);
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	338
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	339 if (! error_state && retval.is_undefined ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	340 error ("bitpack: cannot pack to %s class", numclass.c_str ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	341 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	342 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	343 print_usage ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	344
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	345 return retval;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	346 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	347
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	348 template <class ArrayType>
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	349 boolNDArray
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	350 do_bitunpack (const ArrayType& array)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	351 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	352 typedef typename ArrayType::element_type T;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	353 octave_idx_type n = array.numel () * sizeof (T) * CHAR_BIT;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	354
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	355 boolNDArray retval (get_vec_dims (array.dims (), n));
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	356
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	357 const char packed = reinterpret_cast<const char > (array.fortran_vec ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	358 bool *bits = retval.fortran_vec ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	359
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	360 octave_idx_type m = n / CHAR_BIT;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	361
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	362 for (octave_idx_type i = 0; i < m; i++)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	363 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	364 char c = packed[i];
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	365 bits[0] = c & 1;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	366 for (int j = 1; j < CHAR_BIT; j++)
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	367 bits[j] = (c >>= 1) & 1;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	368 bits += CHAR_BIT;
6869 f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	369 }
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	370
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	371 return retval;
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman] dbateman parents: diff changeset	372 }
9671 e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	373
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	374 DEFUN_DLD (bitunpack, args, ,
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	375 "-- texinfo --\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	376 @deftypefn {Loadable Function} {@var{y} =} bitpack (@var{x})\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	377 Returns a logical array @var{y} corresponding to the raw bit pattern of @var{x}.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	378 @var{x} must belong to one of the built-in numeric classes:\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	379 \n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	380 @example\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	381 \"char\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	382 \"int8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	383 \"int16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	384 \"int32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	385 \"int64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	386 \"uint8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	387 \"uint16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	388 \"uint32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	389 \"uint64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	390 \"double\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	391 \"single\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	392 @end example\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	393 \n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	394 The result is a row vector if @var{x} is a row vector, otherwise it is a column vector.\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	395 @seealso{bitpack,typecast,bitget,bitset}\n\
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	396 @end deftypefn")
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	397 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	398 octave_value retval;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	399
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	400 if (args.length () == 1 && (args(0).is_numeric_type () \|\| args(0).is_string ()))
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	401 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	402 octave_value array = args(0);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	403
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	404 if (array.is_string ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	405 retval = do_bitunpack (array.char_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	406 else if (array.is_integer_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	407 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	408 if (array.is_int8_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	409 retval = do_bitunpack (array.int8_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	410 else if (array.is_int16_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	411 retval = do_bitunpack (array.int16_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	412 else if (array.is_int32_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	413 retval = do_bitunpack (array.int32_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	414 else if (array.is_int64_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	415 retval = do_bitunpack (array.int64_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	416 else if (array.is_uint8_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	417 retval = do_bitunpack (array.uint8_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	418 else if (array.is_uint16_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	419 retval = do_bitunpack (array.uint16_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	420 else if (array.is_uint32_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	421 retval = do_bitunpack (array.uint32_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	422 else if (array.is_uint64_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	423 retval = do_bitunpack (array.uint64_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	424 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	425 assert (0);
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	426 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	427 else if (array.is_complex_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	428 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	429 if (array.is_single_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	430 retval = do_bitunpack (array.float_complex_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	431 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	432 retval = do_bitunpack (array.complex_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	433 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	434 else if (array.is_real_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	435 {
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	436 if (array.is_single_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	437 retval = do_bitunpack (array.float_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	438 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	439 retval = do_bitunpack (array.array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	440 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	441 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	442 error ("bitunpack: invalid input class: %s", array.class_name ().c_str ());
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	443 }
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	444 else
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	445 print_usage ();
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	446
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	447 return retval;
e4b383a0effa merge reinterpret.cc into typecast.cc Jaroslav Hajek <highegg@gmail.com> parents: 9036 diff changeset	448 }

Mercurial > octave

annotate src/DLD-FUNCTIONS/typecast.cc @ 9671:e4b383a0effa