File Coverage

lib/PDL/Stats/Basic-pp-t_corr.c
Criterion Covered Total %
statement 36 51 70.5
branch 58 208 27.8
condition n/a
subroutine n/a
pod n/a
total 94 259 36.2


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/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-t_corr.c"
23              
24             #include
25              
26             #line 1857 "lib/PDL/PP.pm"
27             pdl_error pdl_t_corr_readdata(pdl_trans *__privtrans) {
28             pdl_error PDL_err = {0, NULL, 0};
29             #line 30 "lib/PDL/Stats/Basic-pp-t_corr.c"
30 2 50         if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in t_corr:" "broadcast.incs NULL");
31             /* broadcastloop declarations */
32             int __brcloopval;
33             register PDL_Indx __tind0,__tind1; /* counters along dim */
34 2           register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
35             /* dims here are how many steps along those dims */
36 2           register PDL_Indx __tinc0_r = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
37 2           register PDL_Indx __tinc0_n = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
38 2           register PDL_Indx __tinc0_t = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,0);
39 2           register PDL_Indx __tinc1_r = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
40 2           register PDL_Indx __tinc1_n = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
41 2           register PDL_Indx __tinc1_t = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,1);
42             #define PDL_BROADCASTLOOP_START_t_corr_readdata PDL_BROADCASTLOOP_START( \
43             readdata, \
44             __privtrans->broadcast, \
45             __privtrans->vtable, \
46             r_datap += __offsp[0]; \
47             n_datap += __offsp[1]; \
48             t_datap += __offsp[2]; \
49             , \
50             ( ,r_datap += __tinc1_r - __tinc0_r * __tdims0 \
51             ,n_datap += __tinc1_n - __tinc0_n * __tdims0 \
52             ,t_datap += __tinc1_t - __tinc0_t * __tdims0 \
53             ), \
54             ( ,r_datap += __tinc0_r \
55             ,n_datap += __tinc0_n \
56             ,t_datap += __tinc0_t \
57             ) \
58             )
59             #define PDL_BROADCASTLOOP_END_t_corr_readdata PDL_BROADCASTLOOP_END( \
60             __privtrans->broadcast, \
61             r_datap -= __tinc1_r * __tdims1 + __offsp[0]; \
62             n_datap -= __tinc1_n * __tdims1 + __offsp[1]; \
63             t_datap -= __tinc1_t * __tdims1 + __offsp[2]; \
64             )
65             #ifndef PDL_DECLARE_PARAMS_t_corr_1
66             #define PDL_DECLARE_PARAMS_t_corr_1(PDL_TYPE_OP,PDL_PPSYM_OP) \
67             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, r, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
68             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, n, (__privtrans->pdls[1]), 1, PDL_PPSYM_OP) \
69             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, t, (__privtrans->pdls[2]), 1, PDL_PPSYM_OP)
70             #endif
71 2 100         if ( __privtrans->bvalflag ) { /* ** do 'bad' Code ** */
72             #define PDL_BAD_CODE
73             #define PDL_IF_BAD(t,f) t
74 1           switch (__privtrans->__datatype) { /* Start generic switch */
75 0           case PDL_F: {
76 0 0         PDL_DECLARE_PARAMS_t_corr_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
77 0 0         PDL_BROADCASTLOOP_START_t_corr_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
78 0 0         PDL_IF_BAD(
    0          
    0          
    0          
    0          
    0          
    0          
79             if (PDL_ISBAD2((r_datap)[0],r_badval,F,r_badval_isnan) || PDL_ISBAD2((n_datap)[0],n_badval,F,n_badval_isnan) ) { (t_datap)[0]=t_badval; continue; }
80             if ((n_datap)[0] <= 2) { (t_datap)[0]=t_badval; continue; }
81             ,)
82 0           (t_datap)[0] = (r_datap)[0] / pow( (1 - (r_datap)[0]*(r_datap)[0]) / ((n_datap)[0] - 2) , .5);
83 0 0         }PDL_BROADCASTLOOP_END_t_corr_readdata
    0          
84 0           } break;
85 1           case PDL_D: {
86 1 50         PDL_DECLARE_PARAMS_t_corr_1(PDL_Double,D)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
87 4 50         PDL_BROADCASTLOOP_START_t_corr_readdata {
    50          
    50          
    50          
    50          
    100          
    100          
88 1 50         PDL_IF_BAD(
    0          
    50          
    50          
    0          
    50          
    50          
89             if (PDL_ISBAD2((r_datap)[0],r_badval,D,r_badval_isnan) || PDL_ISBAD2((n_datap)[0],n_badval,D,n_badval_isnan) ) { (t_datap)[0]=t_badval; continue; }
90             if ((n_datap)[0] <= 2) { (t_datap)[0]=t_badval; continue; }
91             ,)
92 1           (t_datap)[0] = (r_datap)[0] / pow( (1 - (r_datap)[0]*(r_datap)[0]) / ((n_datap)[0] - 2) , .5);
93 1 50         }PDL_BROADCASTLOOP_END_t_corr_readdata
    50          
94 1           } break;
95 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in t_corr: unhandled datatype(%d), only handles (FD)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
96             }
97             #undef PDL_BAD_CODE
98             #undef PDL_IF_BAD
99             } else { /* ** else do 'good' Code ** */
100             #define PDL_IF_BAD(t,f) f
101 1           switch (__privtrans->__datatype) { /* Start generic switch */
102 0           case PDL_F: {
103 0 0         PDL_DECLARE_PARAMS_t_corr_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
104 0 0         PDL_BROADCASTLOOP_START_t_corr_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
105             PDL_IF_BAD(
106             if (PDL_ISBAD2((r_datap)[0],r_badval,F,r_badval_isnan) || PDL_ISBAD2((n_datap)[0],n_badval,F,n_badval_isnan) ) { (t_datap)[0]=t_badval; continue; }
107             if ((n_datap)[0] <= 2) { (t_datap)[0]=t_badval; continue; }
108             ,)
109 0           (t_datap)[0] = (r_datap)[0] / pow( (1 - (r_datap)[0]*(r_datap)[0]) / ((n_datap)[0] - 2) , .5);
110 0 0         }PDL_BROADCASTLOOP_END_t_corr_readdata
    0          
111 0           } break;
112 1           case PDL_D: {
113 1 50         PDL_DECLARE_PARAMS_t_corr_1(PDL_Double,D)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
114 4 50         PDL_BROADCASTLOOP_START_t_corr_readdata {
    50          
    50          
    50          
    50          
    100          
    100          
115             PDL_IF_BAD(
116             if (PDL_ISBAD2((r_datap)[0],r_badval,D,r_badval_isnan) || PDL_ISBAD2((n_datap)[0],n_badval,D,n_badval_isnan) ) { (t_datap)[0]=t_badval; continue; }
117             if ((n_datap)[0] <= 2) { (t_datap)[0]=t_badval; continue; }
118             ,)
119 1           (t_datap)[0] = (r_datap)[0] / pow( (1 - (r_datap)[0]*(r_datap)[0]) / ((n_datap)[0] - 2) , .5);
120 1 50         }PDL_BROADCASTLOOP_END_t_corr_readdata
    50          
121 1           } break;
122 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in t_corr: unhandled datatype(%d), only handles (FD)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
123             }
124             #undef PDL_IF_BAD
125             }
126 2           return PDL_err;
127             }
128              
129             static pdl_datatypes pdl_t_corr_vtable_gentypes[] = { PDL_F, PDL_D, -1 };
130             static PDL_Indx pdl_t_corr_vtable_realdims[] = { 0, 0, 0 };
131             static char *pdl_t_corr_vtable_parnames[] = { "r","n","t" };
132             static short pdl_t_corr_vtable_parflags[] = {
133             0,
134             0,
135             PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
136             };
137             static pdl_datatypes pdl_t_corr_vtable_partypes[] = { -1, -1, -1 };
138             static PDL_Indx pdl_t_corr_vtable_realdims_starts[] = { 0, 0, 0 };
139             static PDL_Indx pdl_t_corr_vtable_realdims_ind_ids[] = { 0 };
140             static char *pdl_t_corr_vtable_indnames[] = { "" };
141             pdl_transvtable pdl_t_corr_vtable = {
142             PDL_TRANS_DO_BROADCAST|PDL_TRANS_BADPROCESS, 0, pdl_t_corr_vtable_gentypes, 2, 3, NULL /*CORE21*/,
143             pdl_t_corr_vtable_realdims, pdl_t_corr_vtable_parnames,
144             pdl_t_corr_vtable_parflags, pdl_t_corr_vtable_partypes,
145             pdl_t_corr_vtable_realdims_starts, pdl_t_corr_vtable_realdims_ind_ids, 0,
146             0, pdl_t_corr_vtable_indnames,
147             NULL, pdl_t_corr_readdata, NULL,
148             NULL,
149             0,"PDL::Stats::Basic::t_corr"
150             };
151              
152              
153 2           pdl_error pdl_run_t_corr(pdl *r,pdl *n,pdl *t) {
154 2           pdl_error PDL_err = {0, NULL, 0};
155 2 50         if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
156 2           pdl_trans *__privtrans = PDL->create_trans(&pdl_t_corr_vtable);
157 2 50         if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
158 2           __privtrans->pdls[0] = r;
159 2           __privtrans->pdls[1] = n;
160 2           __privtrans->pdls[2] = t;
161 2 50         PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
162 2 50         PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
163 2           return PDL_err;
164             }