annotate src/DLD-FUNCTIONS/typecast.cc @ 11572:7d6d8c1e471f

Grammarcheck Texinfo for files in src directory.
author Rik <octave@nomad.inbox5.com>
date Wed, 19 Jan 2011 17:51:22 -0800
parents 01f703952eff
children 12df7854fa7c
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
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
11523
fd0a3ac60b0e update copyright notices
John W. Eaton <jwe@octave.org>
parents: 10846
diff changeset
3 Copyright (C) 2007-2011 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,
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
52 dim_vector& old_dims,
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
53 unwind_protect& frame)
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
54 {
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
55 // 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
56 // array. This will ensure the data will be deallocated after we exit.
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
57 frame.add_delete (new ArrayType (array));
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
58
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
59 data = reinterpret_cast<const void *> (array.data ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
60 byte_size = array.byte_size ();
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
61
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
62 old_dims = array.dims ();
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
63 }
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
64
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
65 template <class ArrayType>
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
66 static ArrayType
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
67 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
68 const dim_vector& old_dims)
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
69 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
70 typedef typename ArrayType::element_type T;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
71 octave_idx_type n = byte_size / sizeof (T);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
72
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
73 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
74 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
75 ArrayType retval (get_vec_dims (old_dims, n));
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
76 T *dest = retval.fortran_vec ();
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
77 std::memcpy (dest, data, n * sizeof (T));
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
78
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
79 return retval;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
80 }
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
81 else
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
82 {
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
83 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
84 return ArrayType ();
6869
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
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
89 DEFUN_DLD (typecast, args, ,
8715
954b6f69f51d doc fixes
John W. Eaton <jwe@octave.org>
parents: 8100
diff changeset
90 "-*- texinfo -*-\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
91 @deftypefn {Loadable Function} {} typecast (@var{x}, @var{class})\n\
11553
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents: 11523
diff changeset
92 Return a new array @var{y} resulting from interpreting the data of\n\
9758
09da0bd91412 Periodic grammar check of Octave documentation files to ensure common format
Rik <rdrider0-list@yahoo.com>
parents: 9689
diff changeset
93 @var{x} in memory as data of the numeric class @var{class}. Both the class\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
94 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
95 \n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
96 @example\n\
9758
09da0bd91412 Periodic grammar check of Octave documentation files to ensure common format
Rik <rdrider0-list@yahoo.com>
parents: 9689
diff changeset
97 @group\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
98 \"logical\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
99 \"char\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
100 \"int8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
101 \"int16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
102 \"int32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
103 \"int64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
104 \"uint8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
105 \"uint16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
106 \"uint32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
107 \"uint64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
108 \"double\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
109 \"single\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
110 \"double complex\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
111 \"single complex\"\n\
9758
09da0bd91412 Periodic grammar check of Octave documentation files to ensure common format
Rik <rdrider0-list@yahoo.com>
parents: 9689
diff changeset
112 @end group\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
113 @end example\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
114 \n\
10846
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
115 @noindent\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
116 the last two are reserved for @var{class}; they indicate that a\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
117 complex-valued result is requested. Complex arrays are stored in memory as\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
118 consecutive pairs of real numbers. The sizes of integer types are given by\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
119 their bit counts. Both logical and char are typically one byte wide;\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
120 however, this is not guaranteed by C++. If your system is IEEE conformant,\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
121 single and double should be 4 bytes and 8 bytes wide, respectively. \n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
122 \"logical\" is not allowed for @var{class}. If the input is a row vector,\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
123 the return value is a row vector, otherwise it is a column vector. If the\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
124 bit length of @var{x} is not divisible by that of @var{class}, an error\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
125 occurs.\n\
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
126 \n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
127 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
128 \n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
129 @example\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
130 @group\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
131 @var{x} = uint16 ([1, 65535]);\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
132 typecast (@var{x}, 'uint8')\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
133 @result{} [ 0, 1, 255, 255]\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
134 @end group\n\
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
135 @end example\n\
10807
bdde9c59f77a typecast: Combine two @seealso macro uses to one.
Rik <octave@nomad.inbox5.com>
parents: 10066
diff changeset
136 @seealso{cast, bitunpack, bitpack, swapbytes}\n\
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
137 @end deftypefn")
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
138 {
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
139 octave_value retval;
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
140
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
141 if (args.length () == 2)
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
142 {
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
143 unwind_protect frame;
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
144 const void *data = 0;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
145 octave_idx_type byte_size = 0;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
146 dim_vector old_dims;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
147
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
148 octave_value array = args(0);
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_bool_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
151 get_data_and_bytesize (array.bool_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
152 else if (array.is_string ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
153 get_data_and_bytesize (array.char_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
154 else if (array.is_integer_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
155 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
156 if (array.is_int8_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
157 get_data_and_bytesize (array.int8_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
158 else if (array.is_int16_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
159 get_data_and_bytesize (array.int16_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
160 else if (array.is_int32_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
161 get_data_and_bytesize (array.int32_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
162 else if (array.is_int64_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
163 get_data_and_bytesize (array.int64_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
164 else if (array.is_uint8_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
165 get_data_and_bytesize (array.uint8_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
166 else if (array.is_uint16_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
167 get_data_and_bytesize (array.uint16_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
168 else if (array.is_uint32_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
169 get_data_and_bytesize (array.uint32_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
170 else if (array.is_uint64_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
171 get_data_and_bytesize (array.uint64_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
172 else
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
173 assert (0);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
174 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
175 else if (array.is_complex_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
176 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
177 if (array.is_single_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
178 get_data_and_bytesize (array.float_complex_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
179 else
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
180 get_data_and_bytesize (array.complex_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
181 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
182 else if (array.is_real_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
183 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
184 if (array.is_single_type ())
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
185 get_data_and_bytesize (array.float_array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
186 else
10066
2cd940306a06 make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents: 9758
diff changeset
187 get_data_and_bytesize (array.array_value (), data, byte_size, old_dims, frame);
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
188 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
189 else
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
190 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
191
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
192 std::string numclass = args(1).string_value ();
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 (error_state || numclass.size () == 0)
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
195 ;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
196 else if (numclass == "char")
9689
34d6f005db4b eliminate is_string argument from octave_value character array constructors
John W. Eaton <jwe@octave.org>
parents: 9671
diff changeset
197 retval = octave_value (reinterpret_copy<charNDArray> (data, byte_size, old_dims), array.is_dq_string () ? '"' : '\'');
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
198 else if (numclass[0] == 'i')
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
199 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
200 if (numclass == "int8")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
201 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
202 else if (numclass == "int16")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
203 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
204 else if (numclass == "int32")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
205 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
206 else if (numclass == "int64")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
207 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
208 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
209 else if (numclass[0] == 'u')
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
210 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
211 if (numclass == "uint8")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
212 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
213 else if (numclass == "uint16")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
214 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
215 else if (numclass == "uint32")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
216 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
217 else if (numclass == "uint64")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
218 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
219 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
220 else if (numclass == "single")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
221 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
222 else if (numclass == "double")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
223 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
224 else if (numclass == "single complex")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
225 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
226 else if (numclass == "double complex")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
227 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
228
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
229 if (! error_state && retval.is_undefined ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
230 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
231 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
232 else
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
233 print_usage ();
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
234
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
235 return retval;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
236 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
237
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
238 template <class ArrayType>
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
239 ArrayType
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
240 do_bitpack (const boolNDArray& bitp)
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 typedef typename ArrayType::element_type T;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
243 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
244
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
245 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
246 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
247
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
248 ArrayType retval (get_vec_dims (bitp.dims (), n));
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
249
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
250 const bool *bits = bitp.fortran_vec ();
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
251 char *packed = reinterpret_cast<char *> (retval.fortran_vec ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
252
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
253 octave_idx_type m = n * sizeof (T);
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 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
256 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
257 char c = bits[0];
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
258 for (int j = 1; j < CHAR_BIT; j++)
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
259 c |= bits[j] << j;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
260
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
261 packed[i] = c;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
262 bits += CHAR_BIT;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
263 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
264
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
265 return retval;
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 else
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
268 {
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
269 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
270 return ArrayType ();
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
271 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
272 }
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
273
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
274 DEFUN_DLD (bitpack, args, ,
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
275 "-*- texinfo -*-\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
276 @deftypefn {Loadable Function} {@var{y} =} bitpack (@var{x}, @var{class})\n\
11553
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents: 11523
diff changeset
277 Return a new array @var{y} resulting from interpreting an array\n\
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents: 11523
diff changeset
278 @var{x} as raw bit patterns for data of the numeric class @var{class}.\n\
10846
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
279 @var{class} must be one of the built-in numeric classes:\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
280 \n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
281 @example\n\
9758
09da0bd91412 Periodic grammar check of Octave documentation files to ensure common format
Rik <rdrider0-list@yahoo.com>
parents: 9689
diff changeset
282 @group\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
283 \"char\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
284 \"int8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
285 \"int16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
286 \"int32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
287 \"int64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
288 \"uint8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
289 \"uint16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
290 \"uint32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
291 \"uint64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
292 \"double\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
293 \"single\"\n\
9758
09da0bd91412 Periodic grammar check of Octave documentation files to ensure common format
Rik <rdrider0-list@yahoo.com>
parents: 9689
diff changeset
294 @end group\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
295 @end example\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
296 \n\
10840
89f4d7e294cc Grammarcheck .cc files
Rik <octave@nomad.inbox5.com>
parents: 10807
diff changeset
297 The number of elements of @var{x} should be divisible by the bit length of\n\
10846
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
298 @var{class}. If it is not, excess bits are discarded. Bits come in\n\
11553
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents: 11523
diff changeset
299 increasing order of significance, i.e., @code{x(1)} is bit 0, @code{x(2)} is\n\
10846
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
300 bit 1, etc. The result is a row vector if @var{x} is a row vector, otherwise\n\
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
301 it is a column vector.\n\
11572
7d6d8c1e471f Grammarcheck Texinfo for files in src directory.
Rik <octave@nomad.inbox5.com>
parents: 11553
diff changeset
302 @seealso{bitunpack, typecast}\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
303 @end deftypefn")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
304 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
305 octave_value retval;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
306
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
307 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
308 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
309 boolNDArray bitp = args(0).bool_array_value ();
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
310
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
311 std::string numclass = args(1).string_value ();
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
312
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
313 if (error_state || numclass.size () == 0)
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
314 ;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
315 else if (numclass == "char")
9689
34d6f005db4b eliminate is_string argument from octave_value character array constructors
John W. Eaton <jwe@octave.org>
parents: 9671
diff changeset
316 retval = octave_value (do_bitpack<charNDArray> (bitp), '\'');
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
317 else if (numclass[0] == 'i')
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 if (numclass == "int8")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
320 retval = do_bitpack<int8NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
321 else if (numclass == "int16")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
322 retval = do_bitpack<int16NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
323 else if (numclass == "int32")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
324 retval = do_bitpack<int32NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
325 else if (numclass == "int64")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
326 retval = do_bitpack<int64NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
327 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
328 else if (numclass[0] == 'u')
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 if (numclass == "uint8")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
331 retval = do_bitpack<uint8NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
332 else if (numclass == "uint16")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
333 retval = do_bitpack<uint16NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
334 else if (numclass == "uint32")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
335 retval = do_bitpack<uint32NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
336 else if (numclass == "uint64")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
337 retval = do_bitpack<uint64NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
338 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
339 else if (numclass == "single")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
340 retval = do_bitpack<FloatNDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
341 else if (numclass == "double")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
342 retval = do_bitpack<NDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
343 else if (numclass == "single complex")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
344 retval = do_bitpack<FloatComplexNDArray> (bitp);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
345 else if (numclass == "double complex")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
346 retval = do_bitpack<ComplexNDArray> (bitp);
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
347
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
348 if (! error_state && retval.is_undefined ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
349 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
350 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
351 else
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
352 print_usage ();
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
353
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
354 return retval;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
355 }
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 template <class ArrayType>
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
358 boolNDArray
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
359 do_bitunpack (const ArrayType& array)
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
360 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
361 typedef typename ArrayType::element_type T;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
362 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
363
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
364 boolNDArray retval (get_vec_dims (array.dims (), n));
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
365
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
366 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
367 bool *bits = retval.fortran_vec ();
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
368
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
369 octave_idx_type m = n / CHAR_BIT;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
370
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
371 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
372 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
373 char c = packed[i];
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
374 bits[0] = c & 1;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
375 for (int j = 1; j < CHAR_BIT; j++)
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
376 bits[j] = (c >>= 1) & 1;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
377 bits += CHAR_BIT;
6869
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
378 }
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
379
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
380 return retval;
f9c893831e68 [project @ 2007-09-06 16:38:44 by dbateman]
dbateman
parents:
diff changeset
381 }
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
382
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
383 DEFUN_DLD (bitunpack, args, ,
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
384 "-*- texinfo -*-\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
385 @deftypefn {Loadable Function} {@var{y} =} bitpack (@var{x})\n\
11553
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents: 11523
diff changeset
386 Return an array @var{y} corresponding to the raw bit patterns of\n\
10846
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
387 @var{x}. @var{x} must belong to one of the built-in numeric classes:\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
388 \n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
389 @example\n\
9758
09da0bd91412 Periodic grammar check of Octave documentation files to ensure common format
Rik <rdrider0-list@yahoo.com>
parents: 9689
diff changeset
390 @group\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
391 \"char\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
392 \"int8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
393 \"int16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
394 \"int32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
395 \"int64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
396 \"uint8\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
397 \"uint16\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
398 \"uint32\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
399 \"uint64\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
400 \"double\"\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
401 \"single\"\n\
9758
09da0bd91412 Periodic grammar check of Octave documentation files to ensure common format
Rik <rdrider0-list@yahoo.com>
parents: 9689
diff changeset
402 @end group\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
403 @end example\n\
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
404 \n\
11553
01f703952eff Improve docstrings for functions in DLD-FUNCTIONS directory.
Rik <octave@nomad.inbox5.com>
parents: 11523
diff changeset
405 The result is a row vector if @var{x} is a row vector; otherwise, it is a\n\
10846
a4f482e66b65 Grammarcheck more of the documentation.
Rik <octave@nomad.inbox5.com>
parents: 10840
diff changeset
406 column vector.\n\
11572
7d6d8c1e471f Grammarcheck Texinfo for files in src directory.
Rik <octave@nomad.inbox5.com>
parents: 11553
diff changeset
407 @seealso{bitpack, typecast}\n\
9671
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
408 @end deftypefn")
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
409 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
410 octave_value retval;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
411
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
412 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
413 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
414 octave_value array = args(0);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
415
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
416 if (array.is_string ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
417 retval = do_bitunpack (array.char_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
418 else if (array.is_integer_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
419 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
420 if (array.is_int8_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
421 retval = do_bitunpack (array.int8_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
422 else if (array.is_int16_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
423 retval = do_bitunpack (array.int16_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
424 else if (array.is_int32_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
425 retval = do_bitunpack (array.int32_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
426 else if (array.is_int64_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
427 retval = do_bitunpack (array.int64_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
428 else if (array.is_uint8_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
429 retval = do_bitunpack (array.uint8_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
430 else if (array.is_uint16_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
431 retval = do_bitunpack (array.uint16_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
432 else if (array.is_uint32_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
433 retval = do_bitunpack (array.uint32_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
434 else if (array.is_uint64_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
435 retval = do_bitunpack (array.uint64_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
436 else
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
437 assert (0);
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
438 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
439 else if (array.is_complex_type ())
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 if (array.is_single_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
442 retval = do_bitunpack (array.float_complex_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
443 else
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
444 retval = do_bitunpack (array.complex_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
445 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
446 else if (array.is_real_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
447 {
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
448 if (array.is_single_type ())
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
449 retval = do_bitunpack (array.float_array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
450 else
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
451 retval = do_bitunpack (array.array_value ());
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
452 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
453 else
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
454 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
455 }
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
456 else
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
457 print_usage ();
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
458
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
459 return retval;
e4b383a0effa merge reinterpret.cc into typecast.cc
Jaroslav Hajek <highegg@gmail.com>
parents: 9036
diff changeset
460 }