File Coverage

lib/PDL/Stats/Basic-pp-corr.c

Criterion	Covered	Total	%
statement	70	114	61.4
branch	76	250	30.4
condition			n/a
subroutine			n/a
pod			n/a
total	146	364	40.1

line	stmt	bran	code
1
2			#line 453 "lib/PDL/PP.pm"
3			/*
4			* THIS FILE WAS GENERATED BY PDL::PP from lib/PDL/Stats/Basic.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_Stats_Basic
21			extern Core* PDL; /* Structure hold core C functions */
22			#line 23 "lib/PDL/Stats/Basic-pp-corr.c"
23
24			#include
25
26			#line 1857 "lib/PDL/PP.pm"
27			pdl_error pdl_corr_redodims(pdl_trans *__privtrans) {
28			pdl_error PDL_err = {0, NULL, 0};
29			#line 30 "lib/PDL/Stats/Basic-pp-corr.c"
30			#ifndef PDL_DECLARE_PARAMS_corr_0
31			#define PDL_DECLARE_PARAMS_corr_0(PDL_TYPE_OP,PDL_PPSYM_OP) \
32			PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, a, (__privtrans->pdls[0]), 0, PDL_PPSYM_OP) \
33			PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, b, (__privtrans->pdls[1]), 0, PDL_PPSYM_OP) \
34			PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, c, (__privtrans->pdls[2]), 0, PDL_PPSYM_OP)
35			#endif
36			#define PDL_IF_BAD(t,f) f
37	7		switch (__privtrans->__datatype) { /* Start generic switch */
38	0		case PDL_F: {
39	0	0	PDL_DECLARE_PARAMS_corr_0(PDL_Float,F)
		0
		0
		0
		0
		0
40	0	0	{if (__privtrans->ind_sizes[0] < 2) return PDL->make_error(PDL_EUSERERROR, "Error in corr:" "too few N");}
41	0		} break;
42	7		case PDL_D: {
43	7	50	PDL_DECLARE_PARAMS_corr_0(PDL_Double,D)
		50
		100
		50
		50
		50
44	7	50	{if (__privtrans->ind_sizes[0] < 2) return PDL->make_error(PDL_EUSERERROR, "Error in corr:" "too few N");}
45	7		} break;
46	0		default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in corr: unhandled datatype(%d), only handles (FD)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
47			}
48			#undef PDL_IF_BAD
49
50	7	50	PDL_RETERROR(PDL_err, PDL->redodims_default(__privtrans));
51	7		return PDL_err;
52			}
53
54
55			#line 1857 "lib/PDL/PP.pm"
56			pdl_error pdl_corr_readdata(pdl_trans *__privtrans) {
57			pdl_error PDL_err = {0, NULL, 0};
58			#line 59 "lib/PDL/Stats/Basic-pp-corr.c"
59	7		register PDL_Indx __n_size = __privtrans->ind_sizes[0];
60	7	50	if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in corr:" "broadcast.incs NULL");
61			/* broadcastloop declarations */
62			int __brcloopval;
63			register PDL_Indx __tind0,__tind1; /* counters along dim */
64	7		register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
65			/* dims here are how many steps along those dims */
66	7		register PDL_Indx __tinc0_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
67	7		register PDL_Indx __tinc0_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
68	7		register PDL_Indx __tinc0_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,0);
69	7		register PDL_Indx __tinc1_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
70	7		register PDL_Indx __tinc1_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
71	7		register PDL_Indx __tinc1_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,1);
72			#define PDL_BROADCASTLOOP_START_corr_readdata PDL_BROADCASTLOOP_START( \
73			readdata, \
74			__privtrans->broadcast, \
75			__privtrans->vtable, \
76			a_datap += __offsp[0]; \
77			b_datap += __offsp[1]; \
78			c_datap += __offsp[2]; \
79			, \
80			( ,a_datap += __tinc1_a - __tinc0_a * __tdims0 \
81			,b_datap += __tinc1_b - __tinc0_b * __tdims0 \
82			,c_datap += __tinc1_c - __tinc0_c * __tdims0 \
83			), \
84			( ,a_datap += __tinc0_a \
85			,b_datap += __tinc0_b \
86			,c_datap += __tinc0_c \
87			) \
88			)
89			#define PDL_BROADCASTLOOP_END_corr_readdata PDL_BROADCASTLOOP_END( \
90			__privtrans->broadcast, \
91			a_datap -= __tinc1_a * __tdims1 + __offsp[0]; \
92			b_datap -= __tinc1_b * __tdims1 + __offsp[1]; \
93			c_datap -= __tinc1_c * __tdims1 + __offsp[2]; \
94			)
95	7		register PDL_Indx __inc_a_n = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,0,0)]; (void)__inc_a_n;
96	7		register PDL_Indx __inc_b_n = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,1,0)]; (void)__inc_b_n;
97			#ifndef PDL_DECLARE_PARAMS_corr_1
98			#define PDL_DECLARE_PARAMS_corr_1(PDL_TYPE_OP,PDL_PPSYM_OP) \
99			PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, a, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
100			PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, b, (__privtrans->pdls[1]), 1, PDL_PPSYM_OP) \
101			PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, c, (__privtrans->pdls[2]), 1, PDL_PPSYM_OP)
102			#endif
103	7	100	if ( __privtrans->bvalflag ) { /* do 'bad' Code */
104			#define PDL_BAD_CODE
105			#define PDL_IF_BAD(t,f) t
106	3		switch (__privtrans->__datatype) { /* Start generic switch */
107	0		case PDL_F: {
108	0	0	PDL_DECLARE_PARAMS_corr_1(PDL_Float,F)
		0
		0
		0
		0
		0
		0
		0
		0
		0
		0
		0
109	0	0	PDL_BROADCASTLOOP_START_corr_readdata {
		0
		0
		0
		0
		0
		0
110			PDL_Float ab, sa, sb, a2, b2, cov, va, vb;
111	0		ab=0; sa=0; sb=0; a2=0; b2=0;
112	0		PDL_Indx N = PDL_IF_BAD(0,__privtrans->ind_sizes[0]);
113	0	0	{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
114	0	0	PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n(n))],a_badval,F,a_badval_isnan) \|\| PDL_ISBAD2((b_datap)[0+(__inc_b_n(n))],b_badval,F,b_badval_isnan)) continue; N++;,)
		0
		0
		0
		0
		0
115	0		ab += (a_datap)[0+(__inc_a_n(n))] (b_datap)[0+(__inc_b_n*(n))];
116	0		sa += (a_datap)[0+(__inc_a_n*(n))];
117	0		sb += (b_datap)[0+(__inc_b_n*(n))];
118	0		a2 += (a_datap)[0+(__inc_a_n(n))] (a_datap)[0+(__inc_a_n*(n))];
119	0		b2 += (b_datap)[0+(__inc_b_n(n))] (b_datap)[0+(__inc_b_n*(n))];
120			}} /* Close n */
121	0	0	if (N < 2) { (c_datap)[0]=c_badval; continue; }
122	0		cov = ab - (sa * sb) / N;
123	0		va = a2 - sa*sa / N;
124	0		vb = b2 - sb*sb / N;
125	0		(c_datap)[0] = cov / sqrt( va * vb );
126	0	0	}PDL_BROADCASTLOOP_END_corr_readdata
		0
127	0		} break;
128	3		case PDL_D: {
129	3	50	PDL_DECLARE_PARAMS_corr_1(PDL_Double,D)
		50
		50
		50
		100
		50
		50
		50
		50
		50
		50
		50
130	22	50	PDL_BROADCASTLOOP_START_corr_readdata {
		50
		50
		50
		50
		100
		100
131			PDL_Double ab, sa, sb, a2, b2, cov, va, vb;
132	11		ab=0; sa=0; sb=0; a2=0; b2=0;
133	11		PDL_Indx N = PDL_IF_BAD(0,__privtrans->ind_sizes[0]);
134	113	100	{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
135	102	50	PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n(n))],a_badval,D,a_badval_isnan) \|\| PDL_ISBAD2((b_datap)[0+(__inc_b_n(n))],b_badval,D,b_badval_isnan)) continue; N++;,)
		0
		100
		50
		0
		100
136	93		ab += (a_datap)[0+(__inc_a_n(n))] (b_datap)[0+(__inc_b_n*(n))];
137	93		sa += (a_datap)[0+(__inc_a_n*(n))];
138	93		sb += (b_datap)[0+(__inc_b_n*(n))];
139	93		a2 += (a_datap)[0+(__inc_a_n(n))] (a_datap)[0+(__inc_a_n*(n))];
140	93		b2 += (b_datap)[0+(__inc_b_n(n))] (b_datap)[0+(__inc_b_n*(n))];
141			}} /* Close n */
142	11	50	if (N < 2) { (c_datap)[0]=c_badval; continue; }
143	11		cov = ab - (sa * sb) / N;
144	11		va = a2 - sa*sa / N;
145	11		vb = b2 - sb*sb / N;
146	11		(c_datap)[0] = cov / sqrt( va * vb );
147	3	50	}PDL_BROADCASTLOOP_END_corr_readdata
		50
148	3		} break;
149	0		default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in corr: unhandled datatype(%d), only handles (FD)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
150			}
151			#undef PDL_BAD_CODE
152			#undef PDL_IF_BAD
153			} else { /* else do 'good' Code */
154			#define PDL_IF_BAD(t,f) f
155	4		switch (__privtrans->__datatype) { /* Start generic switch */
156	0		case PDL_F: {
157	0	0	PDL_DECLARE_PARAMS_corr_1(PDL_Float,F)
		0
		0
		0
		0
		0
		0
		0
		0
		0
		0
		0
158	0	0	PDL_BROADCASTLOOP_START_corr_readdata {
		0
		0
		0
		0
		0
		0
159			PDL_Float ab, sa, sb, a2, b2, cov, va, vb;
160	0		ab=0; sa=0; sb=0; a2=0; b2=0;
161	0		PDL_Indx N = PDL_IF_BAD(0,__privtrans->ind_sizes[0]);
162	0	0	{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
163			PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n(n))],a_badval,F,a_badval_isnan) \|\| PDL_ISBAD2((b_datap)[0+(__inc_b_n(n))],b_badval,F,b_badval_isnan)) continue; N++;,)
164	0		ab += (a_datap)[0+(__inc_a_n(n))] (b_datap)[0+(__inc_b_n*(n))];
165	0		sa += (a_datap)[0+(__inc_a_n*(n))];
166	0		sb += (b_datap)[0+(__inc_b_n*(n))];
167	0		a2 += (a_datap)[0+(__inc_a_n(n))] (a_datap)[0+(__inc_a_n*(n))];
168	0		b2 += (b_datap)[0+(__inc_b_n(n))] (b_datap)[0+(__inc_b_n*(n))];
169			}} /* Close n */
170	0	0	if (N < 2) { (c_datap)[0]=c_badval; continue; }
171	0		cov = ab - (sa * sb) / N;
172	0		va = a2 - sa*sa / N;
173	0		vb = b2 - sb*sb / N;
174	0		(c_datap)[0] = cov / sqrt( va * vb );
175	0	0	}PDL_BROADCASTLOOP_END_corr_readdata
		0
176	0		} break;
177	4		case PDL_D: {
178	4	50	PDL_DECLARE_PARAMS_corr_1(PDL_Double,D)
		50
		50
		50
		100
		50
		50
		50
		50
		50
		50
		50
179	36	50	PDL_BROADCASTLOOP_START_corr_readdata {
		50
		50
		50
		50
		100
		100
180			PDL_Double ab, sa, sb, a2, b2, cov, va, vb;
181	20		ab=0; sa=0; sb=0; a2=0; b2=0;
182	20		PDL_Indx N = PDL_IF_BAD(0,__privtrans->ind_sizes[0]);
183	165	100	{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
184			PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_n(n))],a_badval,D,a_badval_isnan) \|\| PDL_ISBAD2((b_datap)[0+(__inc_b_n(n))],b_badval,D,b_badval_isnan)) continue; N++;,)
185	145		ab += (a_datap)[0+(__inc_a_n(n))] (b_datap)[0+(__inc_b_n*(n))];
186	145		sa += (a_datap)[0+(__inc_a_n*(n))];
187	145		sb += (b_datap)[0+(__inc_b_n*(n))];
188	145		a2 += (a_datap)[0+(__inc_a_n(n))] (a_datap)[0+(__inc_a_n*(n))];
189	145		b2 += (b_datap)[0+(__inc_b_n(n))] (b_datap)[0+(__inc_b_n*(n))];
190			}} /* Close n */
191	20	50	if (N < 2) { (c_datap)[0]=c_badval; continue; }
192	20		cov = ab - (sa * sb) / N;
193	20		va = a2 - sa*sa / N;
194	20		vb = b2 - sb*sb / N;
195	20		(c_datap)[0] = cov / sqrt( va * vb );
196	4	50	}PDL_BROADCASTLOOP_END_corr_readdata
		50
197	4		} break;
198	0		default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in corr: unhandled datatype(%d), only handles (FD)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
199			}
200			#undef PDL_IF_BAD
201			}
202	7		return PDL_err;
203			}
204
205			static pdl_datatypes pdl_corr_vtable_gentypes[] = { PDL_F, PDL_D, -1 };
206			static PDL_Indx pdl_corr_vtable_realdims[] = { 1, 1, 0 };
207			static char *pdl_corr_vtable_parnames[] = { "a","b","c" };
208			static short pdl_corr_vtable_parflags[] = {
209			0,
210			0,
211			PDL_PARAM_ISCREAT\|PDL_PARAM_ISOUT\|PDL_PARAM_ISWRITE
212			};
213			static pdl_datatypes pdl_corr_vtable_partypes[] = { -1, -1, -1 };
214			static PDL_Indx pdl_corr_vtable_realdims_starts[] = { 0, 1, 2 };
215			static PDL_Indx pdl_corr_vtable_realdims_ind_ids[] = { 0, 0 };
216			static char *pdl_corr_vtable_indnames[] = { "n" };
217			pdl_transvtable pdl_corr_vtable = {
218			PDL_TRANS_DO_BROADCAST\|PDL_TRANS_BADPROCESS, 0, pdl_corr_vtable_gentypes, 2, 3, NULL /CORE21/,
219			pdl_corr_vtable_realdims, pdl_corr_vtable_parnames,
220			pdl_corr_vtable_parflags, pdl_corr_vtable_partypes,
221			pdl_corr_vtable_realdims_starts, pdl_corr_vtable_realdims_ind_ids, 2,
222			1, pdl_corr_vtable_indnames,
223			pdl_corr_redodims, pdl_corr_readdata, NULL,
224			NULL,
225			0,"PDL::Stats::Basic::corr"
226			};
227
228
229	7		pdl_error pdl_run_corr(pdl a,pdl b,pdl *c) {
230	7		pdl_error PDL_err = {0, NULL, 0};
231	7	50	if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
232	7		pdl_trans *__privtrans = PDL->create_trans(&pdl_corr_vtable);
233	7	50	if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
234	7		__privtrans->pdls[0] = a;
235	7		__privtrans->pdls[1] = b;
236	7		__privtrans->pdls[2] = c;
237	7	50	PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
238	7	50	PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
239	7		return PDL_err;
240			}