File Coverage

lib/PDL/Ufunc-pp-diff2.c
Criterion Covered Total %
statement 43 283 15.1
branch 62 1358 4.5
condition n/a
subroutine n/a
pod n/a
total 105 1641 6.4


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