File Coverage

lib/PDL/Image2D-pp-pnpoly_pp.c
Criterion Covered Total %
statement 40 184 21.7
branch 35 562 6.2
condition n/a
subroutine n/a
pod n/a
total 75 746 10.0


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/Image2D.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_Image2D
21             extern Core* PDL; /* Structure hold core C functions */
22              
23             #line 1857 "lib/PDL/PP.pm"
24             pdl_error pdl_pnpoly_pp_readdata(pdl_trans *__privtrans) {
25             pdl_error PDL_err = {0, NULL, 0};
26             #line 27 "lib/PDL/Image2D-pp-pnpoly_pp.c"
27 2           register PDL_Indx __l_size = __privtrans->ind_sizes[1];
28 2           register PDL_Indx __m_size = __privtrans->ind_sizes[2];
29 2           register PDL_Indx __n_size = __privtrans->ind_sizes[3];
30 2 50         if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in pnpoly_pp:" "broadcast.incs NULL");
31             /* broadcastloop declarations */
32             int __brcloopval;
33             register PDL_Indx __tind0,__tind1; /* counters along dim */
34 2           register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
35             /* dims here are how many steps along those dims */
36 2           register PDL_Indx __tinc0_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
37 2           register PDL_Indx __tinc0_ps = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
38 2           register PDL_Indx __tinc0_msk = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,0);
39 2           register PDL_Indx __tinc1_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
40 2           register PDL_Indx __tinc1_ps = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
41 2           register PDL_Indx __tinc1_msk = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,1);
42             #define PDL_BROADCASTLOOP_START_pnpoly_pp_readdata PDL_BROADCASTLOOP_START( \
43             readdata, \
44             __privtrans->broadcast, \
45             __privtrans->vtable, \
46             a_datap += __offsp[0]; \
47             ps_datap += __offsp[1]; \
48             msk_datap += __offsp[2]; \
49             , \
50             ( ,a_datap += __tinc1_a - __tinc0_a * __tdims0 \
51             ,ps_datap += __tinc1_ps - __tinc0_ps * __tdims0 \
52             ,msk_datap += __tinc1_msk - __tinc0_msk * __tdims0 \
53             ), \
54             ( ,a_datap += __tinc0_a \
55             ,ps_datap += __tinc0_ps \
56             ,msk_datap += __tinc0_msk \
57             ) \
58             )
59             #define PDL_BROADCASTLOOP_END_pnpoly_pp_readdata PDL_BROADCASTLOOP_END( \
60             __privtrans->broadcast, \
61             a_datap -= __tinc1_a * __tdims1 + __offsp[0]; \
62             ps_datap -= __tinc1_ps * __tdims1 + __offsp[1]; \
63             msk_datap -= __tinc1_msk * __tdims1 + __offsp[2]; \
64             )
65 2           register PDL_Indx __inc_a_m = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,0,0)]; (void)__inc_a_m;register PDL_Indx __inc_a_n = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,0,1)]; (void)__inc_a_n;
66 2           register PDL_Indx __inc_msk_m = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,2,0)]; (void)__inc_msk_m;register PDL_Indx __inc_msk_n = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,2,1)]; (void)__inc_msk_n;
67 2           register PDL_Indx __inc_ps_k = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,1,0)]; (void)__inc_ps_k;register PDL_Indx __inc_ps_l = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,1,1)]; (void)__inc_ps_l;
68             #ifndef PDL_DECLARE_PARAMS_pnpoly_pp_1
69             #define PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_TYPE_OP,PDL_PPSYM_OP,PDL_TYPE_PARAM_msk,PDL_PPSYM_PARAM_msk) \
70             PDL_DECLARE_PARAMETER(PDL_TYPE_OP, a, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
71             PDL_DECLARE_PARAMETER(PDL_TYPE_OP, ps, (__privtrans->pdls[1]), 1, PDL_PPSYM_OP) \
72             PDL_DECLARE_PARAMETER(PDL_TYPE_PARAM_msk, msk, (__privtrans->pdls[2]), 1, PDL_PPSYM_PARAM_msk)
73             #endif
74             #define PDL_IF_BAD(t,f) f
75 2           switch (__privtrans->__datatype) { /* Start generic switch */
76 0           case PDL_SB: {
77 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_SByte,A,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
78             {
79             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
80             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
81              
82 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
83 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
84 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
85 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
86 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
87 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
88 0           c = !c;
89 0           j = l;
90             }} /* Close l */
91 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
92             }} /* Close n */}} /* Close m */
93 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
94              
95             #undef VERTX
96             #undef VERTY
97             }
98 0           } break;
99 0           case PDL_B: {
100 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_Byte,B,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
101             {
102             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
103             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
104              
105 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
106 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
107 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
108 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
109 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
110 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
111 0           c = !c;
112 0           j = l;
113             }} /* Close l */
114 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
115             }} /* Close n */}} /* Close m */
116 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
117              
118             #undef VERTX
119             #undef VERTY
120             }
121 0           } break;
122 0           case PDL_S: {
123 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_Short,S,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
124             {
125             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
126             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
127              
128 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
129 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
130 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
131 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
132 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
133 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
134 0           c = !c;
135 0           j = l;
136             }} /* Close l */
137 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
138             }} /* Close n */}} /* Close m */
139 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
140              
141             #undef VERTX
142             #undef VERTY
143             }
144 0           } break;
145 0           case PDL_US: {
146 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_Ushort,U,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
147             {
148             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
149             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
150              
151 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
152 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
153 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
154 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
155 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
156 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
157 0           c = !c;
158 0           j = l;
159             }} /* Close l */
160 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
161             }} /* Close n */}} /* Close m */
162 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
163              
164             #undef VERTX
165             #undef VERTY
166             }
167 0           } break;
168 0           case PDL_L: {
169 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_Long,L,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
170             {
171             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
172             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
173              
174 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
175 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
176 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
177 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
178 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
179 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
180 0           c = !c;
181 0           j = l;
182             }} /* Close l */
183 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
184             }} /* Close n */}} /* Close m */
185 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
186              
187             #undef VERTX
188             #undef VERTY
189             }
190 0           } break;
191 0           case PDL_UL: {
192 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_ULong,K,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
193             {
194             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
195             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
196              
197 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
198 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
199 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
200 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
201 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
202 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
203 0           c = !c;
204 0           j = l;
205             }} /* Close l */
206 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
207             }} /* Close n */}} /* Close m */
208 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
209              
210             #undef VERTX
211             #undef VERTY
212             }
213 0           } break;
214 0           case PDL_IND: {
215 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_Indx,N,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
216             {
217             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
218             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
219              
220 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
221 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
222 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
223 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
224 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
225 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
226 0           c = !c;
227 0           j = l;
228             }} /* Close l */
229 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
230             }} /* Close n */}} /* Close m */
231 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
232              
233             #undef VERTX
234             #undef VERTY
235             }
236 0           } break;
237 0           case PDL_ULL: {
238 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_ULongLong,P,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
239             {
240             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
241             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
242              
243 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
244 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
245 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
246 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
247 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
248 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
249 0           c = !c;
250 0           j = l;
251             }} /* Close l */
252 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
253             }} /* Close n */}} /* Close m */
254 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
255              
256             #undef VERTX
257             #undef VERTY
258             }
259 0           } break;
260 0           case PDL_LL: {
261 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_LongLong,Q,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
262             {
263             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
264             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
265              
266 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
267 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
268 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
269 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
270 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
271 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
272 0           c = !c;
273 0           j = l;
274             }} /* Close l */
275 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
276             }} /* Close n */}} /* Close m */
277 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
278              
279             #undef VERTX
280             #undef VERTY
281             }
282 0           } break;
283 0           case PDL_F: {
284 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_Float,F,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
285             {
286             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
287             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
288              
289 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
290 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
291 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
292 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
293 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
294 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
295 0           c = !c;
296 0           j = l;
297             }} /* Close l */
298 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
299             }} /* Close n */}} /* Close m */
300 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
301              
302             #undef VERTX
303             #undef VERTY
304             }
305 0           } break;
306 2           case PDL_D: {
307 2 50         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_Double,D,PDL_Long,L)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
308             {
309             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
310             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
311              
312 8 50         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    50          
    50          
    50          
    50          
    100          
    100          
313 62 100         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    100          
314 50           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
315 200 100         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
316 150 100         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
317 80 100         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
318 22           c = !c;
319 150           j = l;
320             }} /* Close l */
321 50           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
322             }} /* Close n */}} /* Close m */
323 2 50         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    50          
324              
325             #undef VERTX
326             #undef VERTY
327             }
328 2           } break;
329 0           case PDL_LD: {
330 0 0         PDL_DECLARE_PARAMS_pnpoly_pp_1(PDL_LDouble,E,PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
331             {
332             #define VERTX(q) (ps_datap)[0+(__inc_ps_k*(0))+(__inc_ps_l*(q))]
333             #define VERTY(q) (ps_datap)[0+(__inc_ps_k*(1))+(__inc_ps_l*(q))]
334              
335 0 0         PDL_BROADCASTLOOP_START_pnpoly_pp_readdata
    0          
    0          
    0          
    0          
    0          
    0          
336 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
    0          
337 0           PDL_Indx j = __privtrans->ind_sizes[1]-1, c = 0;
338 0 0         {/* Open l */ PDL_EXPAND2(register PDL_Indx l=0, __l_stop=(__l_size)); for(; l<__l_stop; l+=1) {
339 0 0         if ( ((VERTY(l)>n) != (VERTY(j)>n)) &&
340 0 0         (m < (VERTX(j)-VERTX(l)) * (n-VERTY(l)) / (VERTY(j)-VERTY(l)) + VERTX(l)) )
341 0           c = !c;
342 0           j = l;
343             }} /* Close l */
344 0           (msk_datap)[0+(__inc_msk_m*(m))+(__inc_msk_n*(n))] = c;
345             }} /* Close n */}} /* Close m */
346 0 0         PDL_BROADCASTLOOP_END_pnpoly_pp_readdata
    0          
347              
348             #undef VERTX
349             #undef VERTY
350             }
351 0           } break;
352 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in pnpoly_pp: unhandled datatype(%d), only handles (ABSULKNPQFDE)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
353             }
354             #undef PDL_IF_BAD
355 2           return PDL_err;
356             }
357              
358             static pdl_datatypes pdl_pnpoly_pp_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, -1 };
359             static PDL_Indx pdl_pnpoly_pp_vtable_realdims[] = { 2, 2, 2 };
360             static char *pdl_pnpoly_pp_vtable_parnames[] = { "a","ps","msk" };
361             static short pdl_pnpoly_pp_vtable_parflags[] = {
362             0,
363             0,
364             PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISTYPED|PDL_PARAM_ISWRITE
365             };
366             static pdl_datatypes pdl_pnpoly_pp_vtable_partypes[] = { -1, -1, PDL_L };
367             static PDL_Indx pdl_pnpoly_pp_vtable_realdims_starts[] = { 0, 2, 4 };
368             static PDL_Indx pdl_pnpoly_pp_vtable_realdims_ind_ids[] = { 2, 3, 0, 1, 2, 3 };
369             static char *pdl_pnpoly_pp_vtable_indnames[] = { "k","l","m","n" };
370             pdl_transvtable pdl_pnpoly_pp_vtable = {
371             PDL_TRANS_DO_BROADCAST|PDL_TRANS_BADIGNORE, 0, pdl_pnpoly_pp_vtable_gentypes, 2, 3, NULL /*CORE21*/,
372             pdl_pnpoly_pp_vtable_realdims, pdl_pnpoly_pp_vtable_parnames,
373             pdl_pnpoly_pp_vtable_parflags, pdl_pnpoly_pp_vtable_partypes,
374             pdl_pnpoly_pp_vtable_realdims_starts, pdl_pnpoly_pp_vtable_realdims_ind_ids, 6,
375             4, pdl_pnpoly_pp_vtable_indnames,
376             NULL, pdl_pnpoly_pp_readdata, NULL,
377             NULL,
378             0,"PDL::Image2D::pnpoly_pp"
379             };
380              
381              
382 2           pdl_error pdl_run_pnpoly_pp(pdl *a,pdl *ps,pdl *msk) {
383 2           pdl_error PDL_err = {0, NULL, 0};
384 2 50         if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
385 2           pdl_trans *__privtrans = PDL->create_trans(&pdl_pnpoly_pp_vtable);
386 2 50         if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
387 2           __privtrans->pdls[0] = a;
388 2           __privtrans->pdls[1] = ps;
389 2           __privtrans->pdls[2] = msk;
390 2 50         PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
391 2 50         PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
392 2           return PDL_err;
393             }