File Coverage

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