File Coverage

lib/PDL/Ufunc-pp-magnover.c
Criterion Covered Total %
statement 46 138 33.3
branch 75 540 13.8
condition n/a
subroutine n/a
pod n/a
total 121 678 17.8


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_magnover_readdata(pdl_trans *__privtrans) {
25             pdl_error PDL_err = {0, NULL, 0};
26             #line 27 "lib/PDL/Ufunc-pp-magnover.c"
27 12           register PDL_Indx __n_size = __privtrans->ind_sizes[0];
28 12 50         if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in magnover:" "broadcast.incs NULL");
29             /* broadcastloop declarations */
30             int __brcloopval;
31             register PDL_Indx __tind0,__tind1; /* counters along dim */
32 12           register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
33             /* dims here are how many steps along those dims */
34 12           register PDL_Indx __tinc0_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
35 12           register PDL_Indx __tinc0_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
36 12           register PDL_Indx __tinc1_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
37 12           register PDL_Indx __tinc1_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
38             #define PDL_BROADCASTLOOP_START_magnover_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_magnover_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 12           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_magnover_1
59             #define PDL_DECLARE_PARAMS_magnover_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 12 100         if ( __privtrans->bvalflag ) { /* ** do 'bad' Code ** */
64             #define PDL_BAD_CODE
65             #define PDL_IF_BAD(t,f) t
66 1           switch (__privtrans->__datatype) { /* Start generic switch */
67 1           case PDL_F: {
68             #define PDL_IF_GENTYPE_INTEGER(t,f) f
69             #define PDL_IF_GENTYPE_REAL(t,f) t
70             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
71 1 50         PDL_DECLARE_PARAMS_magnover_1(PDL_Float,F,PDL_Float,F)
    50          
    0          
    50          
    50          
    50          
    50          
    50          
72 4 50         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    50          
    50          
    50          
    50          
    100          
    100          
73 1           PDL_IF_BAD(int flag = 0;,)
74 1 50         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
75 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          
76 0           sum += PDL_IF_GENTYPE_REAL(
77             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
78             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
79             );
80             }} /* Close n */
81 1 50         PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
82 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
83 1 50         }PDL_BROADCASTLOOP_END_magnover_readdata
    50          
84             #undef PDL_IF_GENTYPE_INTEGER
85             #undef PDL_IF_GENTYPE_REAL
86             #undef PDL_IF_GENTYPE_UNSIGNED
87 1           } break;
88 0           case PDL_D: {
89             #define PDL_IF_GENTYPE_INTEGER(t,f) f
90             #define PDL_IF_GENTYPE_REAL(t,f) t
91             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
92 0 0         PDL_DECLARE_PARAMS_magnover_1(PDL_Double,D,PDL_Double,D)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
93 0 0         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    0          
    0          
    0          
    0          
    0          
    0          
94 0           PDL_IF_BAD(int flag = 0;,)
95 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
96 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          
97 0           sum += PDL_IF_GENTYPE_REAL(
98             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
99             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
100             );
101             }} /* Close n */
102 0 0         PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
103 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
104 0 0         }PDL_BROADCASTLOOP_END_magnover_readdata
    0          
105             #undef PDL_IF_GENTYPE_INTEGER
106             #undef PDL_IF_GENTYPE_REAL
107             #undef PDL_IF_GENTYPE_UNSIGNED
108 0           } break;
109 0           case PDL_LD: {
110             #define PDL_IF_GENTYPE_INTEGER(t,f) f
111             #define PDL_IF_GENTYPE_REAL(t,f) t
112             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
113 0 0         PDL_DECLARE_PARAMS_magnover_1(PDL_LDouble,E,PDL_LDouble,E)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
114 0 0         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    0          
    0          
    0          
    0          
    0          
    0          
115 0           PDL_IF_BAD(int flag = 0;,)
116 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
117 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          
118 0           sum += PDL_IF_GENTYPE_REAL(
119             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
120             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
121             );
122             }} /* Close n */
123 0 0         PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
124 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
125 0 0         }PDL_BROADCASTLOOP_END_magnover_readdata
    0          
126             #undef PDL_IF_GENTYPE_INTEGER
127             #undef PDL_IF_GENTYPE_REAL
128             #undef PDL_IF_GENTYPE_UNSIGNED
129 0           } break;
130 0           case PDL_CF: {
131             #define PDL_IF_GENTYPE_INTEGER(t,f) f
132             #define PDL_IF_GENTYPE_REAL(t,f) f
133             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
134 0 0         PDL_DECLARE_PARAMS_magnover_1(PDL_CFloat,G,PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
135 0 0         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    0          
    0          
    0          
    0          
    0          
    0          
136 0           PDL_IF_BAD(int flag = 0;,)
137 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
138 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n*(n))],a_badval,G,a_badval_isnan)) continue; flag = 1;,)
    0          
    0          
    0          
139 0           sum += PDL_IF_GENTYPE_REAL(
140             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
141             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
142             );
143             }} /* Close n */
144 0 0         PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
145 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
146 0 0         }PDL_BROADCASTLOOP_END_magnover_readdata
    0          
147             #undef PDL_IF_GENTYPE_INTEGER
148             #undef PDL_IF_GENTYPE_REAL
149             #undef PDL_IF_GENTYPE_UNSIGNED
150 0           } break;
151 0           case PDL_CD: {
152             #define PDL_IF_GENTYPE_INTEGER(t,f) f
153             #define PDL_IF_GENTYPE_REAL(t,f) f
154             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
155 0 0         PDL_DECLARE_PARAMS_magnover_1(PDL_CDouble,C,PDL_Double,D)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
156 0 0         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    0          
    0          
    0          
    0          
    0          
    0          
157 0           PDL_IF_BAD(int flag = 0;,)
158 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
159 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n*(n))],a_badval,C,a_badval_isnan)) continue; flag = 1;,)
    0          
    0          
    0          
160 0           sum += PDL_IF_GENTYPE_REAL(
161             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
162             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
163             );
164             }} /* Close n */
165 0 0         PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
166 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
167 0 0         }PDL_BROADCASTLOOP_END_magnover_readdata
    0          
168             #undef PDL_IF_GENTYPE_INTEGER
169             #undef PDL_IF_GENTYPE_REAL
170             #undef PDL_IF_GENTYPE_UNSIGNED
171 0           } break;
172 0           case PDL_CLD: {
173             #define PDL_IF_GENTYPE_INTEGER(t,f) f
174             #define PDL_IF_GENTYPE_REAL(t,f) f
175             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
176 0 0         PDL_DECLARE_PARAMS_magnover_1(PDL_CLDouble,H,PDL_LDouble,E)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
177 0 0         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    0          
    0          
    0          
    0          
    0          
    0          
178 0           PDL_IF_BAD(int flag = 0;,)
179 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
180 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n*(n))],a_badval,H,a_badval_isnan)) continue; flag = 1;,)
    0          
    0          
    0          
181 0           sum += PDL_IF_GENTYPE_REAL(
182             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
183             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
184             );
185             }} /* Close n */
186 0 0         PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
187 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
188 0 0         }PDL_BROADCASTLOOP_END_magnover_readdata
    0          
189             #undef PDL_IF_GENTYPE_INTEGER
190             #undef PDL_IF_GENTYPE_REAL
191             #undef PDL_IF_GENTYPE_UNSIGNED
192 0           } break;
193 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in magnover: unhandled datatype(%d), only handles (FDEGCH)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
194             }
195             #undef PDL_BAD_CODE
196             #undef PDL_IF_BAD
197             } else { /* ** else do 'good' Code ** */
198             #define PDL_IF_BAD(t,f) f
199 11           switch (__privtrans->__datatype) { /* Start generic switch */
200 0           case PDL_F: {
201             #define PDL_IF_GENTYPE_INTEGER(t,f) f
202             #define PDL_IF_GENTYPE_REAL(t,f) t
203             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
204 0 0         PDL_DECLARE_PARAMS_magnover_1(PDL_Float,F,PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
205 0 0         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    0          
    0          
    0          
    0          
    0          
    0          
206             PDL_IF_BAD(int flag = 0;,)
207 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
208             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n*(n))],a_badval,F,a_badval_isnan)) continue; flag = 1;,)
209 0           sum += PDL_IF_GENTYPE_REAL(
210             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
211             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
212             );
213             }} /* Close n */
214             PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
215 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
216 0 0         }PDL_BROADCASTLOOP_END_magnover_readdata
    0          
217             #undef PDL_IF_GENTYPE_INTEGER
218             #undef PDL_IF_GENTYPE_REAL
219             #undef PDL_IF_GENTYPE_UNSIGNED
220 0           } break;
221 10           case PDL_D: {
222             #define PDL_IF_GENTYPE_INTEGER(t,f) f
223             #define PDL_IF_GENTYPE_REAL(t,f) t
224             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
225 10 100         PDL_DECLARE_PARAMS_magnover_1(PDL_Double,D,PDL_Double,D)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
226 904 50         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    50          
    50          
    50          
    50          
    100          
    100          
227             PDL_IF_BAD(int flag = 0;,)
228 2410 100         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
229             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n*(n))],a_badval,D,a_badval_isnan)) continue; flag = 1;,)
230 1608           sum += PDL_IF_GENTYPE_REAL(
231             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
232             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
233             );
234             }} /* Close n */
235             PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
236 802 100         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
237 10 50         }PDL_BROADCASTLOOP_END_magnover_readdata
    50          
238             #undef PDL_IF_GENTYPE_INTEGER
239             #undef PDL_IF_GENTYPE_REAL
240             #undef PDL_IF_GENTYPE_UNSIGNED
241 10           } break;
242 0           case PDL_LD: {
243             #define PDL_IF_GENTYPE_INTEGER(t,f) f
244             #define PDL_IF_GENTYPE_REAL(t,f) t
245             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
246 0 0         PDL_DECLARE_PARAMS_magnover_1(PDL_LDouble,E,PDL_LDouble,E)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
247 0 0         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    0          
    0          
    0          
    0          
    0          
    0          
248             PDL_IF_BAD(int flag = 0;,)
249 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
250             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n*(n))],a_badval,E,a_badval_isnan)) continue; flag = 1;,)
251 0           sum += PDL_IF_GENTYPE_REAL(
252             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
253             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
254             );
255             }} /* Close n */
256             PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
257 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
258 0 0         }PDL_BROADCASTLOOP_END_magnover_readdata
    0          
259             #undef PDL_IF_GENTYPE_INTEGER
260             #undef PDL_IF_GENTYPE_REAL
261             #undef PDL_IF_GENTYPE_UNSIGNED
262 0           } break;
263 0           case PDL_CF: {
264             #define PDL_IF_GENTYPE_INTEGER(t,f) f
265             #define PDL_IF_GENTYPE_REAL(t,f) f
266             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
267 0 0         PDL_DECLARE_PARAMS_magnover_1(PDL_CFloat,G,PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
268 0 0         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    0          
    0          
    0          
    0          
    0          
    0          
269             PDL_IF_BAD(int flag = 0;,)
270 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
271             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n*(n))],a_badval,G,a_badval_isnan)) continue; flag = 1;,)
272 0           sum += PDL_IF_GENTYPE_REAL(
273             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
274             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
275             );
276             }} /* Close n */
277             PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
278 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
279 0 0         }PDL_BROADCASTLOOP_END_magnover_readdata
    0          
280             #undef PDL_IF_GENTYPE_INTEGER
281             #undef PDL_IF_GENTYPE_REAL
282             #undef PDL_IF_GENTYPE_UNSIGNED
283 0           } break;
284 1           case PDL_CD: {
285             #define PDL_IF_GENTYPE_INTEGER(t,f) f
286             #define PDL_IF_GENTYPE_REAL(t,f) f
287             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
288 1 50         PDL_DECLARE_PARAMS_magnover_1(PDL_CDouble,C,PDL_Double,D)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
289 4 50         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    50          
    50          
    50          
    50          
    100          
    100          
290             PDL_IF_BAD(int flag = 0;,)
291 5 100         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
292             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n*(n))],a_badval,C,a_badval_isnan)) continue; flag = 1;,)
293 4           sum += PDL_IF_GENTYPE_REAL(
294             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
295             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
296             );
297             }} /* Close n */
298             PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
299 1 50         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
300 1 50         }PDL_BROADCASTLOOP_END_magnover_readdata
    50          
301             #undef PDL_IF_GENTYPE_INTEGER
302             #undef PDL_IF_GENTYPE_REAL
303             #undef PDL_IF_GENTYPE_UNSIGNED
304 1           } break;
305 0           case PDL_CLD: {
306             #define PDL_IF_GENTYPE_INTEGER(t,f) f
307             #define PDL_IF_GENTYPE_REAL(t,f) f
308             #define PDL_IF_GENTYPE_UNSIGNED(t,f) f
309 0 0         PDL_DECLARE_PARAMS_magnover_1(PDL_CLDouble,H,PDL_LDouble,E)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
310 0 0         PDL_BROADCASTLOOP_START_magnover_readdata {long double sum=0;
    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,H,a_badval_isnan)) continue; flag = 1;,)
314 0           sum += PDL_IF_GENTYPE_REAL(
315             (a_datap)[0+(__inc_a_n*(n))]*(a_datap)[0+(__inc_a_n*(n))],
316             creall((a_datap)[0+(__inc_a_n*(n))])*creall((a_datap)[0+(__inc_a_n*(n))]) + cimagl((a_datap)[0+(__inc_a_n*(n))])*cimagl((a_datap)[0+(__inc_a_n*(n))])
317             );
318             }} /* Close n */
319             PDL_IF_BAD(if ( !flag ) { (b_datap)[0]=b_badval; continue; },)
320 0 0         (b_datap)[0] = sum == 0 ? 0 : sqrtl(sum);
321 0 0         }PDL_BROADCASTLOOP_END_magnover_readdata
    0          
322             #undef PDL_IF_GENTYPE_INTEGER
323             #undef PDL_IF_GENTYPE_REAL
324             #undef PDL_IF_GENTYPE_UNSIGNED
325 0           } break;
326 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in magnover: unhandled datatype(%d), only handles (FDEGCH)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
327             }
328             #undef PDL_IF_BAD
329             }
330 12           return PDL_err;
331             }
332              
333             static pdl_datatypes pdl_magnover_vtable_gentypes[] = { PDL_D, PDL_LD, PDL_CF, PDL_CD, PDL_CLD, PDL_F, -1 };
334             static PDL_Indx pdl_magnover_vtable_realdims[] = { 1, 0 };
335             static char *pdl_magnover_vtable_parnames[] = { "a","b" };
336             static short pdl_magnover_vtable_parflags[] = {
337             0,
338             PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISREAL|PDL_PARAM_ISWRITE
339             };
340             static pdl_datatypes pdl_magnover_vtable_partypes[] = { -1, -1 };
341             static PDL_Indx pdl_magnover_vtable_realdims_starts[] = { 0, 1 };
342             static PDL_Indx pdl_magnover_vtable_realdims_ind_ids[] = { 0 };
343             static char *pdl_magnover_vtable_indnames[] = { "n" };
344             pdl_transvtable pdl_magnover_vtable = {
345             PDL_TRANS_DO_BROADCAST|PDL_TRANS_BADPROCESS, 0, pdl_magnover_vtable_gentypes, 1, 2, NULL /*CORE21*/,
346             pdl_magnover_vtable_realdims, pdl_magnover_vtable_parnames,
347             pdl_magnover_vtable_parflags, pdl_magnover_vtable_partypes,
348             pdl_magnover_vtable_realdims_starts, pdl_magnover_vtable_realdims_ind_ids, 1,
349             1, pdl_magnover_vtable_indnames,
350             NULL, pdl_magnover_readdata, NULL,
351             NULL,
352             0,"PDL::Ufunc::magnover"
353             };
354              
355              
356 12           pdl_error pdl_run_magnover(pdl *a,pdl *b) {
357 12           pdl_error PDL_err = {0, NULL, 0};
358 12 50         if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
359 12           pdl_trans *__privtrans = PDL->create_trans(&pdl_magnover_vtable);
360 12 50         if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
361 12           __privtrans->pdls[0] = a;
362 12           __privtrans->pdls[1] = b;
363 12 50         PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
364 12 50         PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
365 12           return PDL_err;
366             }