File Coverage

lib/PDL/Slices-pp-rle.c
Criterion Covered Total %
statement 34 105 32.3
branch 26 552 4.7
condition n/a
subroutine n/a
pod n/a
total 60 657 9.1


line stmt bran cond sub pod time code
1              
2             #line 453 "lib/PDL/PP.pm"
3             /*
4             * THIS FILE WAS GENERATED BY PDL::PP from lib/PDL/Slices.pd! Do not modify!
5             */
6              
7             #define PDL_FREE_CODE(trans, destroy, comp_free_code, ntpriv_free_code) \
8             if (destroy) { \
9             comp_free_code \
10             } \
11             if ((trans)->dims_redone) { \
12             ntpriv_free_code \
13             }
14              
15             #include "EXTERN.h"
16             #include "perl.h"
17             #include "XSUB.h"
18             #include "pdl.h"
19             #include "pdlcore.h"
20             #define PDL PDL_Slices
21             extern Core* PDL; /* Structure hold core C functions */
22              
23             #line 1857 "lib/PDL/PP.pm"
24             pdl_error pdl_rle_redodims(pdl_trans *__privtrans) {
25             pdl_error PDL_err = {0, NULL, 0};
26             #line 27 "lib/PDL/Slices-pp-rle.c"
27 5           __privtrans->ind_sizes[0] = __privtrans->ind_sizes[1];
28 5 50         PDL_RETERROR(PDL_err, PDL->redodims_default(__privtrans));
29 5           return PDL_err;
30             }
31              
32              
33             #line 1857 "lib/PDL/PP.pm"
34             pdl_error pdl_rle_readdata(pdl_trans *__privtrans) {
35             pdl_error PDL_err = {0, NULL, 0};
36             #line 37 "lib/PDL/Slices-pp-rle.c"
37 5           register PDL_Indx __m_size = __privtrans->ind_sizes[0];
38 5           register PDL_Indx __n_size = __privtrans->ind_sizes[1];
39 5 50         if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in rle:" "broadcast.incs NULL");
40             /* broadcastloop declarations */
41             int __brcloopval;
42             register PDL_Indx __tind0,__tind1; /* counters along dim */
43 5           register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
44             /* dims here are how many steps along those dims */
45 5           register PDL_Indx __tinc0_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
46 5           register PDL_Indx __tinc0_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
47 5           register PDL_Indx __tinc0_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,0);
48 5           register PDL_Indx __tinc1_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
49 5           register PDL_Indx __tinc1_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
50 5           register PDL_Indx __tinc1_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,1);
51             #define PDL_BROADCASTLOOP_START_rle_readdata PDL_BROADCASTLOOP_START( \
52             readdata, \
53             __privtrans->broadcast, \
54             __privtrans->vtable, \
55             c_datap += __offsp[0]; \
56             a_datap += __offsp[1]; \
57             b_datap += __offsp[2]; \
58             , \
59             ( ,c_datap += __tinc1_c - __tinc0_c * __tdims0 \
60             ,a_datap += __tinc1_a - __tinc0_a * __tdims0 \
61             ,b_datap += __tinc1_b - __tinc0_b * __tdims0 \
62             ), \
63             ( ,c_datap += __tinc0_c \
64             ,a_datap += __tinc0_a \
65             ,b_datap += __tinc0_b \
66             ) \
67             )
68             #define PDL_BROADCASTLOOP_END_rle_readdata PDL_BROADCASTLOOP_END( \
69             __privtrans->broadcast, \
70             c_datap -= __tinc1_c * __tdims1 + __offsp[0]; \
71             a_datap -= __tinc1_a * __tdims1 + __offsp[1]; \
72             b_datap -= __tinc1_b * __tdims1 + __offsp[2]; \
73             )
74 5           register PDL_Indx __inc_a_m = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,1,0)]; (void)__inc_a_m;
75 5           register PDL_Indx __inc_b_m = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,2,0)]; (void)__inc_b_m;
76 5           register PDL_Indx __inc_c_n = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,0,0)]; (void)__inc_c_n;
77             #ifndef PDL_DECLARE_PARAMS_rle_1
78             #define PDL_DECLARE_PARAMS_rle_1(PDL_TYPE_OP,PDL_PPSYM_OP,PDL_TYPE_PARAM_a,PDL_PPSYM_PARAM_a) \
79             PDL_DECLARE_PARAMETER(PDL_TYPE_OP, c, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
80             PDL_DECLARE_PARAMETER(PDL_TYPE_PARAM_a, a, (__privtrans->pdls[1]), 1, PDL_PPSYM_PARAM_a) \
81             PDL_DECLARE_PARAMETER(PDL_TYPE_OP, b, (__privtrans->pdls[2]), 1, PDL_PPSYM_OP)
82             #endif
83             #define PDL_IF_BAD(t,f) f
84 5           switch (__privtrans->__datatype) { /* Start generic switch */
85 0           case PDL_SB: {
86 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_SByte,A,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
87 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
88             #line 1072 "lib/PDL/Slices.pd"
89             PDL_Indx j=0;
90             PDL_SByte cv, clv = (c_datap)[0+(__inc_c_n*(0))];
91             (b_datap)[0+(__inc_b_m*(0))] = clv;
92             (a_datap)[0+(__inc_a_m*(0))] = 0;
93             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
94             #line 1077 "lib/PDL/Slices.pd"
95             cv = (c_datap)[0+(__inc_c_n*(n))];
96             if (cv == clv) {
97             (a_datap)[0+(__inc_a_m*(j))]++;
98             } else {
99             j++;
100             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
101             (a_datap)[0+(__inc_a_m*(j))] = 1;
102             }
103             }} /* Close n */
104             #line 1086 "lib/PDL/Slices.pd"
105             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
106             #line 1087 "lib/PDL/Slices.pd"
107             (a_datap)[0+(__inc_a_m*(m))] = 0;
108             (b_datap)[0+(__inc_b_m*(m))] = 0;
109             }} /* Close m=j+1 */
110             #line 111 "lib/PDL/Slices-pp-rle.c"
111 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
112 0           } break;
113 0           case PDL_B: {
114 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_Byte,B,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
115 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
116             #line 1072 "lib/PDL/Slices.pd"
117             PDL_Indx j=0;
118             PDL_Byte cv, clv = (c_datap)[0+(__inc_c_n*(0))];
119             (b_datap)[0+(__inc_b_m*(0))] = clv;
120             (a_datap)[0+(__inc_a_m*(0))] = 0;
121             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
122             #line 1077 "lib/PDL/Slices.pd"
123             cv = (c_datap)[0+(__inc_c_n*(n))];
124             if (cv == clv) {
125             (a_datap)[0+(__inc_a_m*(j))]++;
126             } else {
127             j++;
128             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
129             (a_datap)[0+(__inc_a_m*(j))] = 1;
130             }
131             }} /* Close n */
132             #line 1086 "lib/PDL/Slices.pd"
133             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
134             #line 1087 "lib/PDL/Slices.pd"
135             (a_datap)[0+(__inc_a_m*(m))] = 0;
136             (b_datap)[0+(__inc_b_m*(m))] = 0;
137             }} /* Close m=j+1 */
138             #line 139 "lib/PDL/Slices-pp-rle.c"
139 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
140 0           } break;
141 0           case PDL_S: {
142 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_Short,S,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
143 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
144             #line 1072 "lib/PDL/Slices.pd"
145             PDL_Indx j=0;
146             PDL_Short cv, clv = (c_datap)[0+(__inc_c_n*(0))];
147             (b_datap)[0+(__inc_b_m*(0))] = clv;
148             (a_datap)[0+(__inc_a_m*(0))] = 0;
149             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
150             #line 1077 "lib/PDL/Slices.pd"
151             cv = (c_datap)[0+(__inc_c_n*(n))];
152             if (cv == clv) {
153             (a_datap)[0+(__inc_a_m*(j))]++;
154             } else {
155             j++;
156             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
157             (a_datap)[0+(__inc_a_m*(j))] = 1;
158             }
159             }} /* Close n */
160             #line 1086 "lib/PDL/Slices.pd"
161             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
162             #line 1087 "lib/PDL/Slices.pd"
163             (a_datap)[0+(__inc_a_m*(m))] = 0;
164             (b_datap)[0+(__inc_b_m*(m))] = 0;
165             }} /* Close m=j+1 */
166             #line 167 "lib/PDL/Slices-pp-rle.c"
167 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
168 0           } break;
169 0           case PDL_US: {
170 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_Ushort,U,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
171 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
172             #line 1072 "lib/PDL/Slices.pd"
173             PDL_Indx j=0;
174             PDL_Ushort cv, clv = (c_datap)[0+(__inc_c_n*(0))];
175             (b_datap)[0+(__inc_b_m*(0))] = clv;
176             (a_datap)[0+(__inc_a_m*(0))] = 0;
177             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
178             #line 1077 "lib/PDL/Slices.pd"
179             cv = (c_datap)[0+(__inc_c_n*(n))];
180             if (cv == clv) {
181             (a_datap)[0+(__inc_a_m*(j))]++;
182             } else {
183             j++;
184             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
185             (a_datap)[0+(__inc_a_m*(j))] = 1;
186             }
187             }} /* Close n */
188             #line 1086 "lib/PDL/Slices.pd"
189             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
190             #line 1087 "lib/PDL/Slices.pd"
191             (a_datap)[0+(__inc_a_m*(m))] = 0;
192             (b_datap)[0+(__inc_b_m*(m))] = 0;
193             }} /* Close m=j+1 */
194             #line 195 "lib/PDL/Slices-pp-rle.c"
195 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
196 0           } break;
197 0           case PDL_L: {
198 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_Long,L,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
199 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
200             #line 1072 "lib/PDL/Slices.pd"
201             PDL_Indx j=0;
202             PDL_Long cv, clv = (c_datap)[0+(__inc_c_n*(0))];
203             (b_datap)[0+(__inc_b_m*(0))] = clv;
204             (a_datap)[0+(__inc_a_m*(0))] = 0;
205             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
206             #line 1077 "lib/PDL/Slices.pd"
207             cv = (c_datap)[0+(__inc_c_n*(n))];
208             if (cv == clv) {
209             (a_datap)[0+(__inc_a_m*(j))]++;
210             } else {
211             j++;
212             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
213             (a_datap)[0+(__inc_a_m*(j))] = 1;
214             }
215             }} /* Close n */
216             #line 1086 "lib/PDL/Slices.pd"
217             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
218             #line 1087 "lib/PDL/Slices.pd"
219             (a_datap)[0+(__inc_a_m*(m))] = 0;
220             (b_datap)[0+(__inc_b_m*(m))] = 0;
221             }} /* Close m=j+1 */
222             #line 223 "lib/PDL/Slices-pp-rle.c"
223 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
224 0           } break;
225 0           case PDL_UL: {
226 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_ULong,K,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
227 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
228             #line 1072 "lib/PDL/Slices.pd"
229             PDL_Indx j=0;
230             PDL_ULong cv, clv = (c_datap)[0+(__inc_c_n*(0))];
231             (b_datap)[0+(__inc_b_m*(0))] = clv;
232             (a_datap)[0+(__inc_a_m*(0))] = 0;
233             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
234             #line 1077 "lib/PDL/Slices.pd"
235             cv = (c_datap)[0+(__inc_c_n*(n))];
236             if (cv == clv) {
237             (a_datap)[0+(__inc_a_m*(j))]++;
238             } else {
239             j++;
240             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
241             (a_datap)[0+(__inc_a_m*(j))] = 1;
242             }
243             }} /* Close n */
244             #line 1086 "lib/PDL/Slices.pd"
245             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
246             #line 1087 "lib/PDL/Slices.pd"
247             (a_datap)[0+(__inc_a_m*(m))] = 0;
248             (b_datap)[0+(__inc_b_m*(m))] = 0;
249             }} /* Close m=j+1 */
250             #line 251 "lib/PDL/Slices-pp-rle.c"
251 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
252 0           } break;
253 0           case PDL_IND: {
254 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_Indx,N,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
255 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
256             #line 1072 "lib/PDL/Slices.pd"
257             PDL_Indx j=0;
258             PDL_Indx cv, clv = (c_datap)[0+(__inc_c_n*(0))];
259             (b_datap)[0+(__inc_b_m*(0))] = clv;
260             (a_datap)[0+(__inc_a_m*(0))] = 0;
261             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
262             #line 1077 "lib/PDL/Slices.pd"
263             cv = (c_datap)[0+(__inc_c_n*(n))];
264             if (cv == clv) {
265             (a_datap)[0+(__inc_a_m*(j))]++;
266             } else {
267             j++;
268             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
269             (a_datap)[0+(__inc_a_m*(j))] = 1;
270             }
271             }} /* Close n */
272             #line 1086 "lib/PDL/Slices.pd"
273             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
274             #line 1087 "lib/PDL/Slices.pd"
275             (a_datap)[0+(__inc_a_m*(m))] = 0;
276             (b_datap)[0+(__inc_b_m*(m))] = 0;
277             }} /* Close m=j+1 */
278             #line 279 "lib/PDL/Slices-pp-rle.c"
279 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
280 0           } break;
281 0           case PDL_ULL: {
282 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_ULongLong,P,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
283 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
284             #line 1072 "lib/PDL/Slices.pd"
285             PDL_Indx j=0;
286             PDL_ULongLong cv, clv = (c_datap)[0+(__inc_c_n*(0))];
287             (b_datap)[0+(__inc_b_m*(0))] = clv;
288             (a_datap)[0+(__inc_a_m*(0))] = 0;
289             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
290             #line 1077 "lib/PDL/Slices.pd"
291             cv = (c_datap)[0+(__inc_c_n*(n))];
292             if (cv == clv) {
293             (a_datap)[0+(__inc_a_m*(j))]++;
294             } else {
295             j++;
296             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
297             (a_datap)[0+(__inc_a_m*(j))] = 1;
298             }
299             }} /* Close n */
300             #line 1086 "lib/PDL/Slices.pd"
301             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
302             #line 1087 "lib/PDL/Slices.pd"
303             (a_datap)[0+(__inc_a_m*(m))] = 0;
304             (b_datap)[0+(__inc_b_m*(m))] = 0;
305             }} /* Close m=j+1 */
306             #line 307 "lib/PDL/Slices-pp-rle.c"
307 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
308 0           } break;
309 0           case PDL_LL: {
310 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_LongLong,Q,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
311 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
312             #line 1072 "lib/PDL/Slices.pd"
313             PDL_Indx j=0;
314             PDL_LongLong cv, clv = (c_datap)[0+(__inc_c_n*(0))];
315             (b_datap)[0+(__inc_b_m*(0))] = clv;
316             (a_datap)[0+(__inc_a_m*(0))] = 0;
317             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
318             #line 1077 "lib/PDL/Slices.pd"
319             cv = (c_datap)[0+(__inc_c_n*(n))];
320             if (cv == clv) {
321             (a_datap)[0+(__inc_a_m*(j))]++;
322             } else {
323             j++;
324             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
325             (a_datap)[0+(__inc_a_m*(j))] = 1;
326             }
327             }} /* Close n */
328             #line 1086 "lib/PDL/Slices.pd"
329             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
330             #line 1087 "lib/PDL/Slices.pd"
331             (a_datap)[0+(__inc_a_m*(m))] = 0;
332             (b_datap)[0+(__inc_b_m*(m))] = 0;
333             }} /* Close m=j+1 */
334             #line 335 "lib/PDL/Slices-pp-rle.c"
335 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
336 0           } break;
337 0           case PDL_F: {
338 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_Float,F,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
339 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
340             #line 1072 "lib/PDL/Slices.pd"
341             PDL_Indx j=0;
342             PDL_Float cv, clv = (c_datap)[0+(__inc_c_n*(0))];
343             (b_datap)[0+(__inc_b_m*(0))] = clv;
344             (a_datap)[0+(__inc_a_m*(0))] = 0;
345             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
346             #line 1077 "lib/PDL/Slices.pd"
347             cv = (c_datap)[0+(__inc_c_n*(n))];
348             if (cv == clv) {
349             (a_datap)[0+(__inc_a_m*(j))]++;
350             } else {
351             j++;
352             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
353             (a_datap)[0+(__inc_a_m*(j))] = 1;
354             }
355             }} /* Close n */
356             #line 1086 "lib/PDL/Slices.pd"
357             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
358             #line 1087 "lib/PDL/Slices.pd"
359             (a_datap)[0+(__inc_a_m*(m))] = 0;
360             (b_datap)[0+(__inc_b_m*(m))] = 0;
361             }} /* Close m=j+1 */
362             #line 363 "lib/PDL/Slices-pp-rle.c"
363 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
364 0           } break;
365 5           case PDL_D: {
366 5 50         PDL_DECLARE_PARAMS_rle_1(PDL_Double,D,PDL_Indx,N)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
367 28 50         PDL_BROADCASTLOOP_START_rle_readdata {
    50          
    50          
    50          
    50          
    100          
    100          
368             #line 1072 "lib/PDL/Slices.pd"
369             PDL_Indx j=0;
370             PDL_Double cv, clv = (c_datap)[0+(__inc_c_n*(0))];
371             (b_datap)[0+(__inc_b_m*(0))] = clv;
372             (a_datap)[0+(__inc_a_m*(0))] = 0;
373             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
374             #line 1077 "lib/PDL/Slices.pd"
375             cv = (c_datap)[0+(__inc_c_n*(n))];
376             if (cv == clv) {
377             (a_datap)[0+(__inc_a_m*(j))]++;
378             } else {
379             j++;
380             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
381             (a_datap)[0+(__inc_a_m*(j))] = 1;
382             }
383             }} /* Close n */
384             #line 1086 "lib/PDL/Slices.pd"
385             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
386             #line 1087 "lib/PDL/Slices.pd"
387             (a_datap)[0+(__inc_a_m*(m))] = 0;
388             (b_datap)[0+(__inc_b_m*(m))] = 0;
389             }} /* Close m=j+1 */
390             #line 391 "lib/PDL/Slices-pp-rle.c"
391 5 50         }PDL_BROADCASTLOOP_END_rle_readdata
    50          
392 5           } break;
393 0           case PDL_LD: {
394 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_LDouble,E,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
395 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
396             #line 1072 "lib/PDL/Slices.pd"
397             PDL_Indx j=0;
398             PDL_LDouble cv, clv = (c_datap)[0+(__inc_c_n*(0))];
399             (b_datap)[0+(__inc_b_m*(0))] = clv;
400             (a_datap)[0+(__inc_a_m*(0))] = 0;
401             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
402             #line 1077 "lib/PDL/Slices.pd"
403             cv = (c_datap)[0+(__inc_c_n*(n))];
404             if (cv == clv) {
405             (a_datap)[0+(__inc_a_m*(j))]++;
406             } else {
407             j++;
408             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
409             (a_datap)[0+(__inc_a_m*(j))] = 1;
410             }
411             }} /* Close n */
412             #line 1086 "lib/PDL/Slices.pd"
413             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
414             #line 1087 "lib/PDL/Slices.pd"
415             (a_datap)[0+(__inc_a_m*(m))] = 0;
416             (b_datap)[0+(__inc_b_m*(m))] = 0;
417             }} /* Close m=j+1 */
418             #line 419 "lib/PDL/Slices-pp-rle.c"
419 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
420 0           } break;
421 0           case PDL_CF: {
422 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_CFloat,G,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
423 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
424             #line 1072 "lib/PDL/Slices.pd"
425             PDL_Indx j=0;
426             PDL_CFloat cv, clv = (c_datap)[0+(__inc_c_n*(0))];
427             (b_datap)[0+(__inc_b_m*(0))] = clv;
428             (a_datap)[0+(__inc_a_m*(0))] = 0;
429             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
430             #line 1077 "lib/PDL/Slices.pd"
431             cv = (c_datap)[0+(__inc_c_n*(n))];
432             if (cv == clv) {
433             (a_datap)[0+(__inc_a_m*(j))]++;
434             } else {
435             j++;
436             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
437             (a_datap)[0+(__inc_a_m*(j))] = 1;
438             }
439             }} /* Close n */
440             #line 1086 "lib/PDL/Slices.pd"
441             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
442             #line 1087 "lib/PDL/Slices.pd"
443             (a_datap)[0+(__inc_a_m*(m))] = 0;
444             (b_datap)[0+(__inc_b_m*(m))] = 0;
445             }} /* Close m=j+1 */
446             #line 447 "lib/PDL/Slices-pp-rle.c"
447 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
448 0           } break;
449 0           case PDL_CD: {
450 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_CDouble,C,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
451 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
452             #line 1072 "lib/PDL/Slices.pd"
453             PDL_Indx j=0;
454             PDL_CDouble cv, clv = (c_datap)[0+(__inc_c_n*(0))];
455             (b_datap)[0+(__inc_b_m*(0))] = clv;
456             (a_datap)[0+(__inc_a_m*(0))] = 0;
457             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
458             #line 1077 "lib/PDL/Slices.pd"
459             cv = (c_datap)[0+(__inc_c_n*(n))];
460             if (cv == clv) {
461             (a_datap)[0+(__inc_a_m*(j))]++;
462             } else {
463             j++;
464             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
465             (a_datap)[0+(__inc_a_m*(j))] = 1;
466             }
467             }} /* Close n */
468             #line 1086 "lib/PDL/Slices.pd"
469             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
470             #line 1087 "lib/PDL/Slices.pd"
471             (a_datap)[0+(__inc_a_m*(m))] = 0;
472             (b_datap)[0+(__inc_b_m*(m))] = 0;
473             }} /* Close m=j+1 */
474             #line 475 "lib/PDL/Slices-pp-rle.c"
475 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
476 0           } break;
477 0           case PDL_CLD: {
478 0 0         PDL_DECLARE_PARAMS_rle_1(PDL_CLDouble,H,PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
479 0 0         PDL_BROADCASTLOOP_START_rle_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
480             #line 1072 "lib/PDL/Slices.pd"
481             PDL_Indx j=0;
482             PDL_CLDouble cv, clv = (c_datap)[0+(__inc_c_n*(0))];
483             (b_datap)[0+(__inc_b_m*(0))] = clv;
484             (a_datap)[0+(__inc_a_m*(0))] = 0;
485             {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
486             #line 1077 "lib/PDL/Slices.pd"
487             cv = (c_datap)[0+(__inc_c_n*(n))];
488             if (cv == clv) {
489             (a_datap)[0+(__inc_a_m*(j))]++;
490             } else {
491             j++;
492             (b_datap)[0+(__inc_b_m*(j))] = clv = cv;
493             (a_datap)[0+(__inc_a_m*(j))] = 1;
494             }
495             }} /* Close n */
496             #line 1086 "lib/PDL/Slices.pd"
497             {/* Open m=j+1 */ PDL_EXPAND2(register PDL_Indx m=PDLMAX((j+1),0), __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
498             #line 1087 "lib/PDL/Slices.pd"
499             (a_datap)[0+(__inc_a_m*(m))] = 0;
500             (b_datap)[0+(__inc_b_m*(m))] = 0;
501             }} /* Close m=j+1 */
502             #line 503 "lib/PDL/Slices-pp-rle.c"
503 0 0         }PDL_BROADCASTLOOP_END_rle_readdata
    0          
504 0           } break;
505 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in rle: unhandled datatype(%d), only handles (ABSULKNPQFDEGCH)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
506             }
507             #undef PDL_IF_BAD
508 5           return PDL_err;
509             }
510              
511             static pdl_datatypes pdl_rle_vtable_gentypes[] = { PDL_SB, PDL_B, PDL_S, PDL_US, PDL_L, PDL_UL, PDL_IND, PDL_ULL, PDL_LL, PDL_F, PDL_D, PDL_LD, PDL_CF, PDL_CD, PDL_CLD, -1 };
512             static PDL_Indx pdl_rle_vtable_realdims[] = { 1, 1, 1 };
513             static char *pdl_rle_vtable_parnames[] = { "c","a","b" };
514             static short pdl_rle_vtable_parflags[] = {
515             0,
516             PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISTYPED|PDL_PARAM_ISWRITE,
517             PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
518             };
519             static pdl_datatypes pdl_rle_vtable_partypes[] = { -1, PDL_IND, -1 };
520             static PDL_Indx pdl_rle_vtable_realdims_starts[] = { 0, 1, 2 };
521             static PDL_Indx pdl_rle_vtable_realdims_ind_ids[] = { 1, 0, 0 };
522             static char *pdl_rle_vtable_indnames[] = { "m","n" };
523             pdl_transvtable pdl_rle_vtable = {
524             PDL_TRANS_DO_BROADCAST, 0, pdl_rle_vtable_gentypes, 1, 3, NULL /*CORE21*/,
525             pdl_rle_vtable_realdims, pdl_rle_vtable_parnames,
526             pdl_rle_vtable_parflags, pdl_rle_vtable_partypes,
527             pdl_rle_vtable_realdims_starts, pdl_rle_vtable_realdims_ind_ids, 3,
528             2, pdl_rle_vtable_indnames,
529             pdl_rle_redodims, pdl_rle_readdata, NULL,
530             NULL,
531             0,"PDL::Slices::rle"
532             };
533              
534              
535 5           pdl_error pdl_run_rle(pdl *c,pdl *a,pdl *b) {
536 5           pdl_error PDL_err = {0, NULL, 0};
537 5 50         if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
538 5           pdl_trans *__privtrans = PDL->create_trans(&pdl_rle_vtable);
539 5 50         if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
540 5           __privtrans->pdls[0] = c;
541 5           __privtrans->pdls[1] = a;
542 5           __privtrans->pdls[2] = b;
543 5 50         PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
544 5 50         PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
545 5           return PDL_err;
546             }