File Coverage

lib/PDL/Ops-pp-power.c
Criterion Covered Total %
statement 32 61 52.4
branch 49 372 13.1
condition n/a
subroutine n/a
pod n/a
total 81 433 18.7


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/Ops.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_Ops
21             extern Core* PDL; /* Structure hold core C functions */
22             #line 23 "lib/PDL/Ops-pp-power.c"
23              
24             #include
25              
26             #define MOD(X,N) (((N) == 0) ? 0 : ( (X) - (PDL_ABS(N)) * ((long long)((X)/(PDL_ABS(N))) + ( ( ((N) * ((long long)((X)/(N)))) != (X) ) ? ( ( ((N)<0) ? 1 : 0 ) + ( (((X)<0) ? -1 : 0))) : 0 ))))
27             #define BU_MOD(X,N)(((N) == 0) ? 0 : ( (X)-(N)*((uint64_t)((X)/(N))) ))
28             #define SPACE(A,B) ( ((A)<(B)) ? -1 : ((A)!=(B)) )
29              
30             #line 83 "lib/PDL/Ops.pd"
31             #define PDL_BADVAL_WARN_X(datatype, ctype, ppsym, ...) \
32             bad_anyval.type = datatype; bad_anyval.value.ppsym = PDL->bvals.ppsym;
33             #define PDL_BADVAL_WARN(var) \
34             { \
35             PDL_Anyval bad_anyval = { PDL_INVALID, {0} }; \
36             if (!(var->has_badvalue && var->badvalue.type != var->datatype)) { \
37             if (var->has_badvalue) \
38             bad_anyval = var->badvalue; \
39             else { \
40             PDL_GENERICSWITCH(PDL_TYPELIST_ALL, var->datatype, PDL_BADVAL_WARN_X, ) \
41             } \
42             } \
43             if (bad_anyval.type < 0) \
44             barf("Error getting badvalue, type=%d", bad_anyval.type); \
45             complex double bad_c; \
46             ANYVAL_TO_CTYPE(bad_c, complex double, bad_anyval); \
47             if( bad_c == 0 || bad_c == 1 ) \
48             warn(#var " badvalue is set to 0 or 1. This will cause data loss when using badvalues for comparison operators."); \
49             }
50              
51             #line 1857 "lib/PDL/PP.pm"
52             pdl_error pdl_power_readdata(pdl_trans *__privtrans) {
53             pdl_error PDL_err = {0, NULL, 0};
54             #line 55 "lib/PDL/Ops-pp-power.c"
55 440 50         if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in power:" "broadcast.incs NULL");
56             /* broadcastloop declarations */
57             int __brcloopval;
58             register PDL_Indx __tind0,__tind1; /* counters along dim */
59 440           register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
60             /* dims here are how many steps along those dims */
61 440           register PDL_Indx __tinc0_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
62 440           register PDL_Indx __tinc0_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
63 440           register PDL_Indx __tinc0_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,0);
64 440           register PDL_Indx __tinc1_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
65 440           register PDL_Indx __tinc1_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
66 440           register PDL_Indx __tinc1_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,1);
67             #define PDL_BROADCASTLOOP_START_power_readdata PDL_BROADCASTLOOP_START( \
68             readdata, \
69             __privtrans->broadcast, \
70             __privtrans->vtable, \
71             a_datap += __offsp[0]; \
72             b_datap += __offsp[1]; \
73             c_datap += __offsp[2]; \
74             , \
75             ( ,a_datap += __tinc1_a - __tinc0_a * __tdims0 \
76             ,b_datap += __tinc1_b - __tinc0_b * __tdims0 \
77             ,c_datap += __tinc1_c - __tinc0_c * __tdims0 \
78             ), \
79             ( ,a_datap += __tinc0_a \
80             ,b_datap += __tinc0_b \
81             ,c_datap += __tinc0_c \
82             ) \
83             )
84             #define PDL_BROADCASTLOOP_END_power_readdata PDL_BROADCASTLOOP_END( \
85             __privtrans->broadcast, \
86             a_datap -= __tinc1_a * __tdims1 + __offsp[0]; \
87             b_datap -= __tinc1_b * __tdims1 + __offsp[1]; \
88             c_datap -= __tinc1_c * __tdims1 + __offsp[2]; \
89             )
90             #ifndef PDL_DECLARE_PARAMS_power_1
91             #define PDL_DECLARE_PARAMS_power_1(PDL_TYPE_OP,PDL_PPSYM_OP) \
92             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, a, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
93             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, b, (__privtrans->pdls[1]), 1, PDL_PPSYM_OP) \
94             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, c, (__privtrans->pdls[2]), 1, PDL_PPSYM_OP)
95             #endif
96 440 100         if ( __privtrans->bvalflag ) { /* ** do 'bad' Code ** */
97             #define PDL_BAD_CODE
98             #define PDL_IF_BAD(t,f) t
99 2           switch (__privtrans->__datatype) { /* Start generic switch */
100 0           case PDL_F: {
101 0 0         PDL_DECLARE_PARAMS_power_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
102             {
103             #line 199 "lib/PDL/Ops.pd"
104             PDL_IF_BAD(char anybad = 0;,)
105             PDL_BROADCASTLOOP_START_power_readdata
106             #line 201 "lib/PDL/Ops.pd"
107             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,F,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,F,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
108             (c_datap)[0] = (PDL_Float)pow((a_datap)[0],(b_datap)[0]);
109             }
110             PDL_BROADCASTLOOP_END_power_readdata
111             #line 205 "lib/PDL/Ops.pd"
112             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
113             #line 114 "lib/PDL/Ops-pp-power.c"
114             }
115 0           } break;
116 2           case PDL_D: {
117 2 50         PDL_DECLARE_PARAMS_power_1(PDL_Double,D)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
118             {
119             #line 199 "lib/PDL/Ops.pd"
120             PDL_IF_BAD(char anybad = 0;,)
121             PDL_BROADCASTLOOP_START_power_readdata
122             #line 201 "lib/PDL/Ops.pd"
123             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,D,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,D,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
124             (c_datap)[0] = (PDL_Double)pow((a_datap)[0],(b_datap)[0]);
125             }
126             PDL_BROADCASTLOOP_END_power_readdata
127             #line 205 "lib/PDL/Ops.pd"
128             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
129             #line 130 "lib/PDL/Ops-pp-power.c"
130             }
131 2           } break;
132 0           case PDL_LD: {
133 0 0         PDL_DECLARE_PARAMS_power_1(PDL_LDouble,E)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
134             {
135             #line 199 "lib/PDL/Ops.pd"
136             PDL_IF_BAD(char anybad = 0;,)
137             PDL_BROADCASTLOOP_START_power_readdata
138             #line 201 "lib/PDL/Ops.pd"
139             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,E,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,E,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
140             (c_datap)[0] = (PDL_LDouble)pow((a_datap)[0],(b_datap)[0]);
141             }
142             PDL_BROADCASTLOOP_END_power_readdata
143             #line 205 "lib/PDL/Ops.pd"
144             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
145             #line 146 "lib/PDL/Ops-pp-power.c"
146             }
147 0           } break;
148 0           case PDL_CF: {
149 0 0         PDL_DECLARE_PARAMS_power_1(PDL_CFloat,G)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
150             {
151             #line 199 "lib/PDL/Ops.pd"
152             PDL_IF_BAD(char anybad = 0;,)
153             PDL_BROADCASTLOOP_START_power_readdata
154             #line 201 "lib/PDL/Ops.pd"
155             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,G,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,G,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
156             (c_datap)[0] = (PDL_CFloat)pow((a_datap)[0],(b_datap)[0]);
157             }
158             PDL_BROADCASTLOOP_END_power_readdata
159             #line 205 "lib/PDL/Ops.pd"
160             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
161             #line 162 "lib/PDL/Ops-pp-power.c"
162             }
163 0           } break;
164 0           case PDL_CD: {
165 0 0         PDL_DECLARE_PARAMS_power_1(PDL_CDouble,C)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
166             {
167             #line 199 "lib/PDL/Ops.pd"
168             PDL_IF_BAD(char anybad = 0;,)
169             PDL_BROADCASTLOOP_START_power_readdata
170             #line 201 "lib/PDL/Ops.pd"
171             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,C,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,C,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
172             (c_datap)[0] = (PDL_CDouble)pow((a_datap)[0],(b_datap)[0]);
173             }
174             PDL_BROADCASTLOOP_END_power_readdata
175             #line 205 "lib/PDL/Ops.pd"
176             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
177             #line 178 "lib/PDL/Ops-pp-power.c"
178             }
179 0           } break;
180 0           case PDL_CLD: {
181 0 0         PDL_DECLARE_PARAMS_power_1(PDL_CLDouble,H)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
182             {
183             #line 199 "lib/PDL/Ops.pd"
184             PDL_IF_BAD(char anybad = 0;,)
185             PDL_BROADCASTLOOP_START_power_readdata
186             #line 201 "lib/PDL/Ops.pd"
187             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,H,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,H,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
188             (c_datap)[0] = (PDL_CLDouble)pow((a_datap)[0],(b_datap)[0]);
189             }
190             PDL_BROADCASTLOOP_END_power_readdata
191             #line 205 "lib/PDL/Ops.pd"
192             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
193             #line 194 "lib/PDL/Ops-pp-power.c"
194             }
195 0           } break;
196 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in power: unhandled datatype(%d), only handles (FDEGCH)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
197             }
198             #undef PDL_BAD_CODE
199             #undef PDL_IF_BAD
200             } else { /* ** else do 'good' Code ** */
201             #define PDL_IF_BAD(t,f) f
202 438           switch (__privtrans->__datatype) { /* Start generic switch */
203 0           case PDL_F: {
204 0 0         PDL_DECLARE_PARAMS_power_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
205             {
206             #line 199 "lib/PDL/Ops.pd"
207             PDL_IF_BAD(char anybad = 0;,)
208             PDL_BROADCASTLOOP_START_power_readdata
209             #line 201 "lib/PDL/Ops.pd"
210             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,F,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,F,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
211             (c_datap)[0] = (PDL_Float)pow((a_datap)[0],(b_datap)[0]);
212             }
213             PDL_BROADCASTLOOP_END_power_readdata
214             #line 205 "lib/PDL/Ops.pd"
215             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
216             #line 217 "lib/PDL/Ops-pp-power.c"
217             }
218 0           } break;
219 432           case PDL_D: {
220 432 50         PDL_DECLARE_PARAMS_power_1(PDL_Double,D)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
221             {
222             #line 199 "lib/PDL/Ops.pd"
223             PDL_IF_BAD(char anybad = 0;,)
224             PDL_BROADCASTLOOP_START_power_readdata
225             #line 201 "lib/PDL/Ops.pd"
226             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,D,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,D,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
227             (c_datap)[0] = (PDL_Double)pow((a_datap)[0],(b_datap)[0]);
228             }
229             PDL_BROADCASTLOOP_END_power_readdata
230             #line 205 "lib/PDL/Ops.pd"
231             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
232             #line 233 "lib/PDL/Ops-pp-power.c"
233             }
234 392           } break;
235 0           case PDL_LD: {
236 0 0         PDL_DECLARE_PARAMS_power_1(PDL_LDouble,E)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
237             {
238             #line 199 "lib/PDL/Ops.pd"
239             PDL_IF_BAD(char anybad = 0;,)
240             PDL_BROADCASTLOOP_START_power_readdata
241             #line 201 "lib/PDL/Ops.pd"
242             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,E,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,E,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
243             (c_datap)[0] = (PDL_LDouble)pow((a_datap)[0],(b_datap)[0]);
244             }
245             PDL_BROADCASTLOOP_END_power_readdata
246             #line 205 "lib/PDL/Ops.pd"
247             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
248             #line 249 "lib/PDL/Ops-pp-power.c"
249             }
250 0           } break;
251 0           case PDL_CF: {
252 0 0         PDL_DECLARE_PARAMS_power_1(PDL_CFloat,G)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
253             {
254             #line 199 "lib/PDL/Ops.pd"
255             PDL_IF_BAD(char anybad = 0;,)
256             PDL_BROADCASTLOOP_START_power_readdata
257             #line 201 "lib/PDL/Ops.pd"
258             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,G,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,G,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
259             (c_datap)[0] = (PDL_CFloat)pow((a_datap)[0],(b_datap)[0]);
260             }
261             PDL_BROADCASTLOOP_END_power_readdata
262             #line 205 "lib/PDL/Ops.pd"
263             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
264             #line 265 "lib/PDL/Ops-pp-power.c"
265             }
266 0           } break;
267 6           case PDL_CD: {
268 6 50         PDL_DECLARE_PARAMS_power_1(PDL_CDouble,C)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
269             {
270             #line 199 "lib/PDL/Ops.pd"
271             PDL_IF_BAD(char anybad = 0;,)
272             PDL_BROADCASTLOOP_START_power_readdata
273             #line 201 "lib/PDL/Ops.pd"
274             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,C,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,C,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
275             (c_datap)[0] = (PDL_CDouble)pow((a_datap)[0],(b_datap)[0]);
276             }
277             PDL_BROADCASTLOOP_END_power_readdata
278             #line 205 "lib/PDL/Ops.pd"
279             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
280             #line 281 "lib/PDL/Ops-pp-power.c"
281             }
282 6           } break;
283 0           case PDL_CLD: {
284 0 0         PDL_DECLARE_PARAMS_power_1(PDL_CLDouble,H)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
285             {
286             #line 199 "lib/PDL/Ops.pd"
287             PDL_IF_BAD(char anybad = 0;,)
288             PDL_BROADCASTLOOP_START_power_readdata
289             #line 201 "lib/PDL/Ops.pd"
290             PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,H,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,H,b_badval_isnan) ) { anybad = 1; (c_datap)[0]=c_badval; } else ,) {
291             (c_datap)[0] = (PDL_CLDouble)pow((a_datap)[0],(b_datap)[0]);
292             }
293             PDL_BROADCASTLOOP_END_power_readdata
294             #line 205 "lib/PDL/Ops.pd"
295             PDL_IF_BAD(if (anybad) __privtrans->pdls[2]->state |= PDL_BADVAL;,)
296             #line 297 "lib/PDL/Ops-pp-power.c"
297             }
298 0           } break;
299 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in power: unhandled datatype(%d), only handles (FDEGCH)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
300             }
301             #undef PDL_IF_BAD
302             }
303 400           return PDL_err;
304             }
305              
306             static pdl_datatypes pdl_power_vtable_gentypes[] = { PDL_CF, PDL_CD, PDL_CLD, PDL_F, PDL_LD, PDL_D, -1 };
307             static PDL_Indx pdl_power_vtable_realdims[] = { 0, 0, 0 };
308             static char *pdl_power_vtable_parnames[] = { "a","b","c" };
309             static short pdl_power_vtable_parflags[] = {
310             0,
311             0,
312             PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
313             };
314             static pdl_datatypes pdl_power_vtable_partypes[] = { -1, -1, -1 };
315             static PDL_Indx pdl_power_vtable_realdims_starts[] = { 0, 0, 0 };
316             static PDL_Indx pdl_power_vtable_realdims_ind_ids[] = { 0 };
317             static char *pdl_power_vtable_indnames[] = { "" };
318             pdl_transvtable pdl_power_vtable = {
319             PDL_TRANS_DO_BROADCAST|PDL_TRANS_BADPROCESS, 0, pdl_power_vtable_gentypes, 2, 3, NULL /*CORE21*/,
320             pdl_power_vtable_realdims, pdl_power_vtable_parnames,
321             pdl_power_vtable_parflags, pdl_power_vtable_partypes,
322             pdl_power_vtable_realdims_starts, pdl_power_vtable_realdims_ind_ids, 0,
323             0, pdl_power_vtable_indnames,
324             NULL, pdl_power_readdata, NULL,
325             NULL,
326             0,"PDL::Ops::power"
327             };
328              
329              
330 80           pdl_error pdl_run_power(pdl *a,pdl *b,pdl *c) {
331 80           pdl_error PDL_err = {0, NULL, 0};
332 80 50         if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
333 80           pdl_trans *__privtrans = PDL->create_trans(&pdl_power_vtable);
334 80 50         if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
335 80           __privtrans->pdls[0] = a;
336 80           __privtrans->pdls[1] = b;
337 80           __privtrans->pdls[2] = c;
338 80 50         PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
339 80 50         PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
340 80           return PDL_err;
341             }