File Coverage

lib/PDL/ImageND-pp-rebin.c

Criterion	Covered	Total	%
statement	0	258	0.0
branch	0	468	0.0
condition			n/a
subroutine			n/a
pod			n/a
total	0	726	0.0

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