Mercurial > forge

annotate extra/NaN/src/histo_mex.cpp @ 6551:95ad0310becb octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
(none)
author schloegl
date Sun, 10 Jan 2010 22:22:46 +0000
parents ea89f0ae2495
children 8b1ab99c0e1c
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
1 //-------------------------------------------------------------------
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
2 #pragma hdrstop
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
3 //-------------------------------------------------------------------
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
4 // C-MEX implementation of SUMSKIPNAN - this function is part of the NaN-toolbox.
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
5 //
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
6 //
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
7 // This program is free software; you can redistribute it and/or modify
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
8 // it under the terms of the GNU General Public License as published by
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
9 // the Free Software Foundation; either version 3 of the License, or
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
10 // (at your option) any later version.
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
11 //
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
12 // This program is distributed in the hope that it will be useful,
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
13 // but WITHOUT ANY WARRANTY; without even the implied warranty of
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
15 // GNU General Public License for more details.
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
16 //
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
17 // You should have received a copy of the GNU General Public License
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
18 // along with this program; if not, see <http://www.gnu.org/licenses/>.
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
19 //
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
20 //
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
21 // histo_mex: computes histogram
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
22 //
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
23 // Input:
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
24 // - data matrix
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
25 // - flag for row-wise histogram
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
26 //
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
27 // Output:
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
28 // - histogram
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
29 // HIS.X
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
30 // HIS.H
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
31 //
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
32 // $Id$
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
33 // Copyright (C) 2009 Alois Schloegl <a.schloegl@ieee.org>
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
34 // This function is part of the NaN-toolbox
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
35 // http://hci.tugraz.at/~schloegl/matlab/NaN/
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
36 //
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
37 //-------------------------------------------------------------------
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
38
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
39 /* TODO:
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
40 speed: its slower than the m-functions histo2/3/4
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
41 |-> use a more efficient sorting function
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
42 resembling of histo3 for multicolumn data.
5894
d14e3342cb28 fix h.X(256) for uint8
schloegl
parents: 5761
diff changeset
43 support of complex data and char-strings
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
44 */
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
45
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
46 #include <inttypes.h>
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
47 #include <math.h>
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
48 #include <string.h>
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
49 #include "mex.h"
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
50
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
51
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
52 #ifdef tmwtypes_h
6121
ea89f0ae2495 fix version check for defining mwSize; set default ACC_LEVEL 0
schloegl
parents: 5900
diff changeset
53 #if (MX_API_VER<=0x07020000)
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
54 typedef int mwSize;
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
55 #endif
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
56 #endif
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
57
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
58 struct sort_t {
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
59 uint8_t *Table; // data table
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
60 size_t Size; // sizeof elements e.g. 4 for single
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
61 size_t Stride; // for multicolumn data
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
62 size_t N; // number of rows
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
63 mxClassID Type; // data type
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
64 } Sort;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
65
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
66 //inline int compare(const sqize_t *a, const size_t *b) {
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
67 int compare(const void *a, const void *b) {
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
68 int z = 0;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
69 size_t i = 0;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
70 size_t ix1 = *(size_t*)a;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
71 size_t ix2 = *(size_t*)b;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
72
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
73 while ((i<Sort.N) && !z) {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
74 switch (Sort.Type) {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
75 case mxCHAR_CLASS:
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
76 z = memcmp(Sort.Table+ix1*Sort.Size,Sort.Table+ix2*Sort.Size,Sort.Size);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
77 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
78 case mxINT32_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
79 int32_t f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
80 f1 = ((int32_t*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
81 f2 = ((int32_t*)Sort.Table)[ix2];
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
82 if (f1<f2) z = -1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
83 else if (f1>f2) z = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
84 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
85 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
86 case mxUINT32_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
87 uint32_t f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
88 f1 = ((uint32_t*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
89 f2 = ((uint32_t*)Sort.Table)[ix2];
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
90 if (f1<f2) z = -1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
91 else if (f1>f2) z = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
92 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
93 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
94 case mxINT64_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
95 int64_t f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
96 f1 = ((int64_t*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
97 f2 = ((int64_t*)Sort.Table)[ix2];
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
98 if (f1<f2) z = -1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
99 else if (f1>f2) z = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
100 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
101 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
102 case mxUINT64_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
103 uint64_t f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
104 f1 = ((uint64_t*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
105 f2 = ((uint64_t*)Sort.Table)[ix2];
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
106 if (f1<f2) z = -1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
107 else if (f1>f2) z = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
108 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
109 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
110 case mxSINGLE_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
111 float f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
112 f1 = ((float*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
113 f2 = ((float*)Sort.Table)[ix2];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
114 switch (isnan(f1) + 2*isnan(f2)) {
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
115 case 0:
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
116 if (f1<f2) z = -1;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
117 else if (f1>f2) z = 1;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
118 case 3:
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
119 break;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
120 case 1: z = 1; break;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
121 case 2: z = -1; break;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
122 }
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
123 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
124 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
125 case mxDOUBLE_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
126 double f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
127 f1 = ((double*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
128 f2 = ((double*)Sort.Table)[ix2];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
129 switch (isnan(f1) + 2*isnan(f2)) {
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
130 case 0:
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
131 if (f1<f2) z = -1;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
132 else if (f1>f2) z = 1;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
133 case 3:
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
134 break;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
135 case 1: z = 1; break;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
136 case 2: z = -1; break;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
137 }
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
138 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
139 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
140 case mxINT16_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
141 int16_t f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
142 f1 = ((int16_t*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
143 f2 = ((int16_t*)Sort.Table)[ix2];
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
144 if (f1<f2) z = -1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
145 else if (f1>f2) z = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
146 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
147 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
148 case mxUINT16_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
149 uint16_t f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
150 f1 = ((uint16_t*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
151 f2 = ((uint16_t*)Sort.Table)[ix2];
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
152 if (f1<f2) z = -1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
153 else if (f1>f2) z = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
154 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
155 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
156 case mxINT8_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
157 int8_t f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
158 f1 = ((int8_t*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
159 f2 = ((int8_t*)Sort.Table)[ix2];
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
160 if (f1<f2) z = -1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
161 else if (f1>f2) z = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
162 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
163 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
164 case mxUINT8_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
165 uint8_t f1,f2;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
166 f1 = ((uint8_t*)Sort.Table)[ix1];
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
167 f2 = ((uint8_t*)Sort.Table)[ix2];
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
168 if (f1<f2) z = -1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
169 else if (f1>f2) z = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
170 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
171 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
172 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
173 i++;
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
174 ix1 += Sort.Stride;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
175 ix2 += Sort.Stride;
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
176 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
177 return(z);
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
178 }
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
179
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
180
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
181 void mexFunction(int POutputCount, mxArray* POutput[], int PInputCount, const mxArray *PInputs[])
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
182 {
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
183
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
184 const mwSize *SZ;
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
185 char flag_rows = 0;
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
186 char done = 0;
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
187 mwSize j, k; // running indices
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
188 ssize_t l;
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
189 const mxArray *W = NULL;
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
190 double *w = NULL;
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
191
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
192 // check for proper number of input and output arguments
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
193 if ((PInputCount <= 0) || (PInputCount > 3)) {
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
194 mexPrintf("HISTO_MEX computes histogram from vector or column matrices\n\n");
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
195 mexPrintf("usage:\tHIS = histo_mex(Y)\n\t\tComputes histogram from each column\n");
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
196 mexPrintf("\t[HIS,tix] = histo_mex(Y,'rows')\n\t\tComputes row-wise histogram, tix is useful for data compression.\n\t\t Y = HIS.X(tix,:); \n\n");
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
197
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
198 mexPrintf("see also: HISTO2, HISTO3, HISTO4\n\n");
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
199 mexErrMsgTxt("HISTO_MEX requires 1 or 2 input arguments\n");
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
200 }
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
201 if (POutputCount > 2)
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
202 mexErrMsgTxt("histo.MEX has 1 output arguments.");
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
203
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
204 // get 1st argument
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
205 if (mxIsComplex(PInputs[0]))
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
206 mexErrMsgTxt("complex argument not supported (yet). ");
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
207 // TODO: support complex argument!
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
208
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
209 if (PInputCount==1)
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
210 ; // histo_mex(X)
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
211 else if (mxIsChar(PInputs[1])) {
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
212 // histo_mex(X,'rows')
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
213 char *t = mxArrayToString(PInputs[1]);
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
214 flag_rows = !strcmp(t,"rows");
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
215 mxFree(t);
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
216 // histo_mex(X,'rows',W)
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
217 if ((PInputCount>2) && mxIsDouble(PInputs[2])) W = PInputs[2];
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
218 }
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
219 // histo_mex(X,W)
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
220 else if (mxIsDouble(PInputs[1])) {
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
221 W = PInputs[1];
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
222 }
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
223 else
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
224 mexErrMsgTxt("Weight vector must be REAL/DOUBLE.");
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
225
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
226 if (W != NULL) {
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
227 if (mxGetM(PInputs[0])==mxGetM(W) )
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
228 w = (double*)mxGetData(W);
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
229 else
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
230 mexErrMsgTxt("number of rows in X and W do not match.");
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
231
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
232 for (k=0; (k<mxGetM(W)) && (w[k]>=0.0); k++);
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
233 if (k<mxGetM(W))
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
234 mexWarnMsgTxt("weight vector contains also non-negative numbers or NaN.");
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
235 }
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
236
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
237 // get size
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
238 mwSize ND = mxGetNumberOfDimensions(PInputs[0]);
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
239 // NN = mxGetNumberOfElements(PInputs[0]);
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
240 SZ = mxGetDimensions(PInputs[0]);
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
241
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
242 if (ND>2)
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
243 mexErrMsgTxt("Error HISTO.MEX: input must be vector or matrix (no more than two dimensions)");
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
244
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
245 size_t n = SZ[0];
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
246 size_t sz = 1;
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
247 char flag = 0;
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
248
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
249 const char *fnames[] = {"datatype","X","H"};
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
250 mxArray *HIS = mxCreateStructMatrix(1, 1, 3, fnames);
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
251 mxSetField(HIS,0,"datatype",mxCreateString("HISTOGRAM"));
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
252
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
253 if (flag_rows || (SZ[1]==1)) {
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
254
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
255 ///***** SORT each column: initialize sorting algorithm
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
256 int (*compar)(const void*, const void*);
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
257 size_t *idx = NULL;
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
258 idx = (size_t*) mxMalloc(SZ[0]*sizeof(size_t));
6551
schloegl
parents: 6121
diff changeset
259 for (n=0; n<SZ[0]; n++) {
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
260 idx[n]=n;
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
261 }
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
262 Sort.Type = mxGetClassID(PInputs[0]);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
263 Sort.Table = (uint8_t*) mxGetData(PInputs[0]);
5900
264f3fd4b27f support of logical data added
schloegl
parents: 5894
diff changeset
264
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
265 switch (mxGetClassID(PInputs[0])) {
5900
264f3fd4b27f support of logical data added
schloegl
parents: 5894
diff changeset
266 case mxLOGICAL_CLASS:
264f3fd4b27f support of logical data added
schloegl
parents: 5894
diff changeset
267 Sort.Size = sizeof(mxLogical);
264f3fd4b27f support of logical data added
schloegl
parents: 5894
diff changeset
268 break;
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
269 case mxINT8_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
270 case mxUINT8_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
271 Sort.Size = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
272 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
273 case mxINT16_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
274 case mxUINT16_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
275 Sort.Size = 2;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
276 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
277 case mxINT32_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
278 case mxUINT32_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
279 case mxSINGLE_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
280 Sort.Size = 4;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
281 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
282 case mxINT64_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
283 case mxUINT64_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
284 case mxDOUBLE_CLASS:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
285 Sort.Size = 8;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
286 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
287 default:
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
288 mexErrMsgTxt("unsupported input type");
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
289 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
290 Sort.N = flag_rows ? SZ[1] : 1;
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
291 qsort(idx,SZ[0],sizeof(*idx),compare);
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
292
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
293 // count number of different elements
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
294 n = SZ[0] ? 1 : 0;
6551
schloegl
parents: 6121
diff changeset
295 size_t k;
schloegl
parents: 6121
diff changeset
296 for (k=1; k<SZ[0]; k++) {
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
297 if (compare(idx+k-1,idx+k)) n++;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
298 }
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
299
5730
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
300 uint64_t *tix = NULL;
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
301 if (POutputCount>1) {
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
302 POutput[1] = mxCreateNumericMatrix(SZ[0], 1, mxUINT64_CLASS,mxREAL);
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
303 tix = (uint64_t*)mxGetData(POutput[1]);
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
304 }
c45fb9ded279 fix NaN support in HISTO_MEX; support of tix output
schloegl
parents: 5729
diff changeset
305
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
306 // fill HIS.H and HIS.X
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
307 mxArray *H = mxCreateNumericMatrix(n, 1, mxDOUBLE_CLASS,mxREAL);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
308 mxArray *X = mxCreateNumericMatrix(n, SZ[1], mxGetClassID(PInputs[0]),mxREAL);
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
309 mxSetField(HIS,0,"H",H);
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
310 mxSetField(HIS,0,"X",X);
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
311 double *h = (double*)mxGetData(H);
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
312 uint8_t *x = (uint8_t*)mxGetData(X);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
313
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
314 l = -1;
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
315 if (tix) tix[idx[0]] = 1;
6551
schloegl
parents: 6121
diff changeset
316 for (k=0; k<SZ[0]; k++) {
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
317 if ((k==0) || compare(&idx[k-1], &idx[k])) {
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
318 l++;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
319 for (j=0; j<SZ[1]; j++) {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
320 memcpy(x + (l+j*n)*Sort.Size, Sort.Table+(idx[k] + j*Sort.Stride)*Sort.Size, Sort.Size);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
321 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
322 }
5760
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
323 if (tix) tix[idx[k]] = l+1;
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
324 if (w != NULL)
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
325 h[l]+=w[idx[k]];
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
326 else
0b1b67486b3a support for data weigthening added
schloegl
parents: 5730
diff changeset
327 h[l]+=1.0;
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
328 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
329 mxFree(idx);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
330 done = 1;
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
331 }
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
332
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
333 else {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
334 switch (mxGetClassID(PInputs[0])) {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
335 /*
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
336 case mxCHAR_CLASS: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
337 mxArray *H = mxCreateNumericMatrix(256, SZ[1], mxUINT64_CLASS,mxREAL);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
338 mxArray *X = mxCreateNumericMatrix(256, 1, mxINT8_CLASS,mxREAL);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
339 mxSetField(HIS,0,"H",H);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
340 mxSetField(HIS,0,"X",X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
341
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
342 char *x;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
343 x = (char*)mxGetData(X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
344 qsort(x,n,sz,strcpy);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
345 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
346 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
347 */
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
348 case mxINT8_CLASS: {
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
349 mxArray *H = mxCreateNumericMatrix(256, SZ[1], mxDOUBLE_CLASS,mxREAL);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
350 mxArray *X = mxCreateNumericMatrix(256, 1, mxINT8_CLASS,mxREAL);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
351 mxSetField(HIS,0,"H",H);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
352 mxSetField(HIS,0,"X",X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
353
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
354 int8_t *x;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
355 x = (int8_t*)mxGetData(X);
5894
d14e3342cb28 fix h.X(256) for uint8
schloegl
parents: 5761
diff changeset
356 for (k=0; k<0x0100; k++)
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
357 x[k]=k-128;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
358
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
359 x = (int8_t*)mxGetData(PInputs[0]);
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
360 double *h = (double*)mxGetData(H);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
361 for (k=0; k<SZ[0]*SZ[1]; k++)
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
362 h[x[k]+128+(k/SZ[0]<<8)] += (w!=NULL ? w[k%SZ[0]] : 1.0);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
363
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
364 done = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
365 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
366 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
367
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
368 case mxUINT8_CLASS: {
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
369 mxArray *H = mxCreateNumericMatrix(256, SZ[1], mxDOUBLE_CLASS,mxREAL);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
370 mxArray *X = mxCreateNumericMatrix(256, 1, mxUINT8_CLASS,mxREAL);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
371 mxSetField(HIS,0,"H",H);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
372 mxSetField(HIS,0,"X",X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
373
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
374 uint8_t *x = (uint8_t*)mxGetData(X);
5894
d14e3342cb28 fix h.X(256) for uint8
schloegl
parents: 5761
diff changeset
375 for (k=0; k<0x0100; k++) x[k]=k;
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
376
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
377 x = (uint8_t*)mxGetData(PInputs[0]);
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
378 double *h = (double*)mxGetData(H);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
379 for (k=0; k<SZ[0]*SZ[1]; k++)
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
380 h[x[k]+(k/SZ[0]<<8)] += (w!=NULL ? w[k%SZ[0]] : 1.0);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
381
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
382 done = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
383 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
384 }
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
385
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
386 case mxINT16_CLASS: {
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
387 mxArray *H = mxCreateNumericMatrix(0x10000, SZ[1], mxDOUBLE_CLASS,mxREAL);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
388 mxArray *X = mxCreateNumericMatrix(0x10000, 1, mxINT16_CLASS,mxREAL);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
389 mxSetField(HIS,0,"H",H);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
390 mxSetField(HIS,0,"X",X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
391
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
392 double *h = (double*)mxGetData(H);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
393 int16_t *x = (int16_t*)mxGetData(X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
394 for (k=0; k<0x10000; k++)
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
395 x[k]=k-0x8000;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
396
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
397 x = (int16_t*)mxGetData(PInputs[0]);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
398 for (k=0; k<SZ[0]*SZ[1]; k++)
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
399 h[x[k]+0x8000+(k/SZ[0]<<16)] += (w!=NULL ? w[k%SZ[0]] : 1.0);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
400
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
401 done = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
402 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
403 }
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
404
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
405 case mxUINT16_CLASS: {
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
406 mxArray *H = mxCreateNumericMatrix(0x10000, SZ[1], mxDOUBLE_CLASS,mxREAL);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
407 mxArray *X = mxCreateNumericMatrix(0x10000, 1, mxINT16_CLASS,mxREAL);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
408 mxSetField(HIS,0,"H",H);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
409 mxSetField(HIS,0,"X",X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
410
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
411 double *h = (double*)mxGetData(H);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
412 int16_t *x = (int16_t*)mxGetData(X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
413 for (k=0; k<0x10000; k++)
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
414 x[k]=k-0x8000;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
415
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
416 uint16_t *x16 = (uint16_t*)mxGetData(PInputs[0]);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
417 for (k=0; k<SZ[0]*SZ[1]; k++)
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
418 h[x16[k]+(k/SZ[0]<<16)] += (w!=NULL ? w[k%SZ[0]] : 1.0);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
419 done = 1;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
420 break;
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
421 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
422
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
423 default: {
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
424 /* FIXME */
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
425 mexErrMsgTxt("multicolumns with int32 or larger not supported!");
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
426
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
427 mxArray *H = mxCreateNumericMatrix(0x10000, SZ[1], mxDOUBLE_CLASS,mxREAL);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
428 mxArray *X = mxCreateNumericMatrix(0x10000, SZ[1], mxGetClassID(PInputs[0]),mxREAL);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
429 mxSetField(HIS,0,"H",H);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
430 mxSetField(HIS,0,"X",X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
431
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
432 double *h = (double*)mxGetData(H);
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
433 int16_t *x = (int16_t*)mxGetData(X);
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
434
6551
schloegl
parents: 6121
diff changeset
435 for (n=0; n<SZ[1]; n++) {
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
436 }
5761
1dbdb1ab56ba support for data weigthening added - fix for (u)int8/16 data
schloegl
parents: 5760
diff changeset
437
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
438 }
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
439 } // end switch
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
440 }
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
441
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
442
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
443 /******* output *******/
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
444 if (done) POutput[0] = HIS;
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
445
5729
49cec4fd294b flag_rows done, single vector, and (u)int8/16 done;
schloegl
parents: 5726
diff changeset
446 return;
5726
5beb367bca41 add histo_mex, currently only (u)int8/16 are supported
schloegl
parents:
diff changeset
447 }