File Coverage

lib/PDL/Ufunc-pp-intover.c

Criterion	Covered	Total	%
statement	60	501	11.9
branch	85	1546	5.5
condition			n/a
subroutine			n/a
pod			n/a
total	145	2047	7.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/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_intover_readdata(pdl_trans *__privtrans) {
25			pdl_error PDL_err = {0, NULL, 0};
26			#line 27 "lib/PDL/Ufunc-pp-intover.c"
27	7		register PDL_Indx __n_size = __privtrans->ind_sizes[0];
28	7	50	if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in intover:" "broadcast.incs NULL");
29			/* broadcastloop declarations */
30			int __brcloopval;
31			register PDL_Indx __tind0,__tind1; /* counters along dim */
32	7		register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
33			/* dims here are how many steps along those dims */
34	7		register PDL_Indx __tinc0_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
35	7		register PDL_Indx __tinc0_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
36	7		register PDL_Indx __tinc1_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
37	7		register PDL_Indx __tinc1_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
38			#define PDL_BROADCASTLOOP_START_intover_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_intover_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	7		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_intover_1
59			#define PDL_DECLARE_PARAMS_intover_1(PDL_TYPE_OP,PDL_PPSYM_OP,PDL_TYPE_PARAM_b,PDL_PPSYM_PARAM_b) \
60			PDL_DECLARE_PARAMETER(PDL_TYPE_OP, a, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
61			PDL_DECLARE_PARAMETER(PDL_TYPE_PARAM_b, b, (__privtrans->pdls[1]), 1, PDL_PPSYM_PARAM_b)
62			#endif
63			#define PDL_IF_BAD(t,f) f
64	7		switch (__privtrans->__datatype) { /* Start generic switch */
65	0		case PDL_SB: {
66	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_SByte,A,PDL_Float,F)
		0
		0
		0
		0
		0
67			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
68	0		switch (__privtrans->ind_sizes[0]) {
69	0		case 1:
70	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
71	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
72	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
73	0		break;
74	0		case 2:
75	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
76	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
77	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
78	0		break;
79	0		case 3:
80	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
81	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
82	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
83	0		break;
84	0		case 4:
85	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
86	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
87	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
88	0		break;
89	0		case 5:
90	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
91	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
92	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
93	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
94	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
95	0		break;
96	0		default:
97	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
98	0		PDL_Float tmp = 0;
99	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
100	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
101	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
102	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
103	0		(b_datap)[0] = tmp;
104	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
105			}
106			}
107	0		} break;
108	0		case PDL_B: {
109	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_Byte,B,PDL_Float,F)
		0
		0
		0
		0
		0
110			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
111	0		switch (__privtrans->ind_sizes[0]) {
112	0		case 1:
113	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
114	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
115	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
116	0		break;
117	0		case 2:
118	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
119	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
120	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
121	0		break;
122	0		case 3:
123	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
124	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
125	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
126	0		break;
127	0		case 4:
128	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
129	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
130	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
131	0		break;
132	0		case 5:
133	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
134	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
135	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
136	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
137	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
138	0		break;
139	0		default:
140	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
141	0		PDL_Float tmp = 0;
142	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
143	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
144	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
145	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
146	0		(b_datap)[0] = tmp;
147	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
148			}
149			}
150	0		} break;
151	0		case PDL_S: {
152	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_Short,S,PDL_Float,F)
		0
		0
		0
		0
		0
153			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
154	0		switch (__privtrans->ind_sizes[0]) {
155	0		case 1:
156	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
157	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
158	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
159	0		break;
160	0		case 2:
161	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
162	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
163	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
164	0		break;
165	0		case 3:
166	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
167	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
168	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
169	0		break;
170	0		case 4:
171	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
172	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
173	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
174	0		break;
175	0		case 5:
176	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
177	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
178	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
179	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
180	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
181	0		break;
182	0		default:
183	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
184	0		PDL_Float tmp = 0;
185	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
186	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
187	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
188	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
189	0		(b_datap)[0] = tmp;
190	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
191			}
192			}
193	0		} break;
194	0		case PDL_US: {
195	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_Ushort,U,PDL_Float,F)
		0
		0
		0
		0
		0
196			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
197	0		switch (__privtrans->ind_sizes[0]) {
198	0		case 1:
199	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
200	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
201	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
202	0		break;
203	0		case 2:
204	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
205	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
206	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
207	0		break;
208	0		case 3:
209	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
210	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
211	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
212	0		break;
213	0		case 4:
214	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
215	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
216	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
217	0		break;
218	0		case 5:
219	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
220	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
221	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
222	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
223	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
224	0		break;
225	0		default:
226	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
227	0		PDL_Float tmp = 0;
228	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
229	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
230	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
231	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
232	0		(b_datap)[0] = tmp;
233	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
234			}
235			}
236	0		} break;
237	0		case PDL_L: {
238	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_Long,L,PDL_Float,F)
		0
		0
		0
		0
		0
239			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
240	0		switch (__privtrans->ind_sizes[0]) {
241	0		case 1:
242	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
243	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
244	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
245	0		break;
246	0		case 2:
247	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
248	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
249	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
250	0		break;
251	0		case 3:
252	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
253	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
254	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
255	0		break;
256	0		case 4:
257	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
258	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
259	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
260	0		break;
261	0		case 5:
262	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
263	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
264	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
265	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
266	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
267	0		break;
268	0		default:
269	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
270	0		PDL_Float tmp = 0;
271	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
272	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
273	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
274	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
275	0		(b_datap)[0] = tmp;
276	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
277			}
278			}
279	0		} break;
280	0		case PDL_UL: {
281	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_ULong,K,PDL_Float,F)
		0
		0
		0
		0
		0
282			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
283	0		switch (__privtrans->ind_sizes[0]) {
284	0		case 1:
285	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
286	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
287	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
288	0		break;
289	0		case 2:
290	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
291	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
292	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
293	0		break;
294	0		case 3:
295	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
296	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
297	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
298	0		break;
299	0		case 4:
300	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
301	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
302	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
303	0		break;
304	0		case 5:
305	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
306	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
307	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
308	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
309	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
310	0		break;
311	0		default:
312	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
313	0		PDL_Float tmp = 0;
314	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
315	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
316	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
317	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
318	0		(b_datap)[0] = tmp;
319	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
320			}
321			}
322	0		} break;
323	0		case PDL_IND: {
324	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_Indx,N,PDL_Float,F)
		0
		0
		0
		0
		0
325			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
326	0		switch (__privtrans->ind_sizes[0]) {
327	0		case 1:
328	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
329	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
330	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
331	0		break;
332	0		case 2:
333	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
334	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
335	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
336	0		break;
337	0		case 3:
338	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
339	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
340	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
341	0		break;
342	0		case 4:
343	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
344	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
345	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
346	0		break;
347	0		case 5:
348	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
349	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
350	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
351	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
352	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
353	0		break;
354	0		default:
355	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
356	0		PDL_Float tmp = 0;
357	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
358	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
359	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
360	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
361	0		(b_datap)[0] = tmp;
362	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
363			}
364			}
365	0		} break;
366	0		case PDL_ULL: {
367	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_ULongLong,P,PDL_Float,F)
		0
		0
		0
		0
		0
368			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
369	0		switch (__privtrans->ind_sizes[0]) {
370	0		case 1:
371	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
372	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
373	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
374	0		break;
375	0		case 2:
376	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
377	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
378	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
379	0		break;
380	0		case 3:
381	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
382	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
383	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
384	0		break;
385	0		case 4:
386	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
387	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
388	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
389	0		break;
390	0		case 5:
391	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
392	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
393	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
394	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
395	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
396	0		break;
397	0		default:
398	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
399	0		PDL_Float tmp = 0;
400	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
401	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
402	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
403	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
404	0		(b_datap)[0] = tmp;
405	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
406			}
407			}
408	0		} break;
409	0		case PDL_LL: {
410	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_LongLong,Q,PDL_Float,F)
		0
		0
		0
		0
		0
411			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
412	0		switch (__privtrans->ind_sizes[0]) {
413	0		case 1:
414	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
415	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
416	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
417	0		break;
418	0		case 2:
419	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
420	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
421	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
422	0		break;
423	0		case 3:
424	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
425	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
426	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
427	0		break;
428	0		case 4:
429	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
430	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
431	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
432	0		break;
433	0		case 5:
434	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
435	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
436	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
437	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
438	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
439	0		break;
440	0		default:
441	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
442	0		PDL_Float tmp = 0;
443	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
444	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
445	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
446	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
447	0		(b_datap)[0] = tmp;
448	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
449			}
450			}
451	0		} break;
452	0		case PDL_F: {
453	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_Float,F,PDL_Float,F)
		0
		0
		0
		0
		0
454			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
455	0		switch (__privtrans->ind_sizes[0]) {
456	0		case 1:
457	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
458	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
459	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
460	0		break;
461	0		case 2:
462	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
463	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
464	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
465	0		break;
466	0		case 3:
467	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
468	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
469	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
470	0		break;
471	0		case 4:
472	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
473	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
474	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
475	0		break;
476	0		case 5:
477	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
478	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
479	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
480	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
481	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
482	0		break;
483	0		default:
484	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
485	0		PDL_Float tmp = 0;
486	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
487	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
488	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
489	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
490	0		(b_datap)[0] = tmp;
491	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
492			}
493			}
494	0		} break;
495	7		case PDL_D: {
496	7	50	PDL_DECLARE_PARAMS_intover_1(PDL_Double,D,PDL_Double,D)
		50
		50
		50
		50
		50
497			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
498	7		switch (__privtrans->ind_sizes[0]) {
499	1		case 1:
500	4	50	PDL_BROADCASTLOOP_START_intover_readdata
		50
		50
		50
		50
		100
		100
501	1		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
502	1	50	PDL_BROADCASTLOOP_END_intover_readdata
		50
503	1		break;
504	1		case 2:
505	4	50	PDL_BROADCASTLOOP_START_intover_readdata
		50
		50
		50
		50
		100
		100
506	1		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
507	1	50	PDL_BROADCASTLOOP_END_intover_readdata
		50
508	1		break;
509	1		case 3:
510	4	50	PDL_BROADCASTLOOP_START_intover_readdata
		50
		50
		50
		50
		100
		100
511	1		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
512	1	50	PDL_BROADCASTLOOP_END_intover_readdata
		50
513	1		break;
514	1		case 4:
515	4	50	PDL_BROADCASTLOOP_START_intover_readdata
		50
		50
		50
		50
		100
		100
516	1		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
517	1	50	PDL_BROADCASTLOOP_END_intover_readdata
		50
518	1		break;
519	1		case 5:
520	4	50	PDL_BROADCASTLOOP_START_intover_readdata
		50
		50
		50
		50
		100
		100
521	1		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
522	1		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
523	1		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
524	1	50	PDL_BROADCASTLOOP_END_intover_readdata
		50
525	1		break;
526	2		default:
527	8	50	PDL_BROADCASTLOOP_START_intover_readdata
		50
		50
		50
		50
		100
		100
528	2		PDL_Double tmp = 0;
529	3	100	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
530	4	100	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
531	4	100	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
532	4	100	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
533	2		(b_datap)[0] = tmp;
534	2	50	PDL_BROADCASTLOOP_END_intover_readdata
		50
535			}
536			}
537	7		} break;
538	0		case PDL_LD: {
539	0	0	PDL_DECLARE_PARAMS_intover_1(PDL_LDouble,E,PDL_LDouble,E)
		0
		0
		0
		0
		0
540			{/* Integration formulae from Press et al 2nd Ed S 4.1 */
541	0		switch (__privtrans->ind_sizes[0]) {
542	0		case 1:
543	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
544	0		(b_datap)[0] = 0.; /* not a(n=>0); as interval has zero width */
545	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
546	0		break;
547	0		case 2:
548	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
549	0		(b_datap)[0] = 0.5((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(1))]);
550	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
551	0		break;
552	0		case 3:
553	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
554	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+4(a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))])/3.;
555	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
556	0		break;
557	0		case 4:
558	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
559	0		(b_datap)[0] = ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(3))]+3.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(2))]))0.375;
560	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
561	0		break;
562	0		case 5:
563	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
564	0		(b_datap)[0] = (14.((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n*(4))])
565	0		+64.((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n*(3))])
566	0		+24.(a_datap)[0+(__inc_a_n(2))])/45.;
567	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
568	0		break;
569	0		default:
570	0	0	PDL_BROADCASTLOOP_START_intover_readdata
		0
		0
		0
		0
		0
		0
571	0		PDL_LDouble tmp = 0;
572	0	0	{/* Open n=3:-3 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((3),0), __n_stop=(__n_size-3)); for(; n<__n_stop; n+=1) { tmp += (a_datap)[0+(__inc_a_n(n))]; }} /* Close n=3:-3 */
573	0	0	{/* Open n=-3:-2 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-3),0), __n_stop=(__n_size-2)); for(; n<__n_stop; n+=1) { tmp += (23./24.)((a_datap)[0+(__inc_a_n(2))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-3:-2 */
574	0	0	{/* Open n=-2:-1 / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-2),0), __n_stop=(__n_size-1)); for(; n<__n_stop; n+=1) { tmp += (7./6.) ((a_datap)[0+(__inc_a_n(1))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-2:-1 */
575	0	0	{/* Open n=-1: / PDL_EXPAND2(register PDL_Indx n=PDLMAX((__n_size-1),0), __n_stop=(__n_size)); for(; n<__n_stop; n+=1) { tmp += (3./8.) ((a_datap)[0+(__inc_a_n(0))]+(a_datap)[0+(__inc_a_n(n))]); }} /* Close n=-1: */
576	0		(b_datap)[0] = tmp;
577	0	0	PDL_BROADCASTLOOP_END_intover_readdata
		0
578			}
579			}
580	0		} break;
581	0		default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in intover: unhandled datatype(%d), only handles (ABSULKNPQFDE)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
582			}
583			#undef PDL_IF_BAD
584	7		return PDL_err;
585			}
586
587			static pdl_datatypes pdl_intover_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 };
588			static PDL_Indx pdl_intover_vtable_realdims[] = { 1, 0 };
589			static char *pdl_intover_vtable_parnames[] = { "a","b" };
590			static short pdl_intover_vtable_parflags[] = {
591			0,
592			PDL_PARAM_ISCREAT\|PDL_PARAM_ISOUT\|PDL_PARAM_ISTPLUS\|PDL_PARAM_ISTYPED\|PDL_PARAM_ISWRITE
593			};
594			static pdl_datatypes pdl_intover_vtable_partypes[] = { -1, PDL_F };
595			static PDL_Indx pdl_intover_vtable_realdims_starts[] = { 0, 1 };
596			static PDL_Indx pdl_intover_vtable_realdims_ind_ids[] = { 0 };
597			static char *pdl_intover_vtable_indnames[] = { "n" };
598			pdl_transvtable pdl_intover_vtable = {
599			PDL_TRANS_DO_BROADCAST\|PDL_TRANS_BADIGNORE, 0, pdl_intover_vtable_gentypes, 1, 2, NULL /CORE21/,
600			pdl_intover_vtable_realdims, pdl_intover_vtable_parnames,
601			pdl_intover_vtable_parflags, pdl_intover_vtable_partypes,
602			pdl_intover_vtable_realdims_starts, pdl_intover_vtable_realdims_ind_ids, 1,
603			1, pdl_intover_vtable_indnames,
604			NULL, pdl_intover_readdata, NULL,
605			NULL,
606			0,"PDL::Ufunc::intover"
607			};
608
609
610	7		pdl_error pdl_run_intover(pdl a,pdl b) {
611	7		pdl_error PDL_err = {0, NULL, 0};
612	7	50	if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
613	7		pdl_trans *__privtrans = PDL->create_trans(&pdl_intover_vtable);
614	7	50	if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
615	7		__privtrans->pdls[0] = a;
616	7		__privtrans->pdls[1] = b;
617	7	50	PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
618	7	50	PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
619	7		return PDL_err;
620			}