File Coverage

lib/PDL/Stats/GLM-pp-dvrs.c
Criterion Covered Total %
statement 38 55 69.0
branch 62 212 29.2
condition n/a
subroutine n/a
pod n/a
total 100 267 37.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/Stats/GLM.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_GLM
21             extern Core* PDL; /* Structure hold core C functions */
22             #line 23 "lib/PDL/Stats/GLM-pp-dvrs.c"
23              
24             #include
25             #include
26             #include
27              
28              
29             #line 1857 "lib/PDL/PP.pm"
30             pdl_error pdl_dvrs_readdata(pdl_trans *__privtrans) {
31             pdl_error PDL_err = {0, NULL, 0};
32             #line 33 "lib/PDL/Stats/GLM-pp-dvrs.c"
33 2 50         if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in dvrs:" "broadcast.incs NULL");
34             /* broadcastloop declarations */
35             int __brcloopval;
36             register PDL_Indx __tind0,__tind1; /* counters along dim */
37 2           register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
38             /* dims here are how many steps along those dims */
39 2           register PDL_Indx __tinc0_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
40 2           register PDL_Indx __tinc0_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
41 2           register PDL_Indx __tinc0_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,0);
42 2           register PDL_Indx __tinc1_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
43 2           register PDL_Indx __tinc1_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
44 2           register PDL_Indx __tinc1_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,1);
45             #define PDL_BROADCASTLOOP_START_dvrs_readdata PDL_BROADCASTLOOP_START( \
46             readdata, \
47             __privtrans->broadcast, \
48             __privtrans->vtable, \
49             a_datap += __offsp[0]; \
50             b_datap += __offsp[1]; \
51             c_datap += __offsp[2]; \
52             , \
53             ( ,a_datap += __tinc1_a - __tinc0_a * __tdims0 \
54             ,b_datap += __tinc1_b - __tinc0_b * __tdims0 \
55             ,c_datap += __tinc1_c - __tinc0_c * __tdims0 \
56             ), \
57             ( ,a_datap += __tinc0_a \
58             ,b_datap += __tinc0_b \
59             ,c_datap += __tinc0_c \
60             ) \
61             )
62             #define PDL_BROADCASTLOOP_END_dvrs_readdata PDL_BROADCASTLOOP_END( \
63             __privtrans->broadcast, \
64             a_datap -= __tinc1_a * __tdims1 + __offsp[0]; \
65             b_datap -= __tinc1_b * __tdims1 + __offsp[1]; \
66             c_datap -= __tinc1_c * __tdims1 + __offsp[2]; \
67             )
68             #ifndef PDL_DECLARE_PARAMS_dvrs_1
69             #define PDL_DECLARE_PARAMS_dvrs_1(PDL_TYPE_OP,PDL_PPSYM_OP) \
70             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, a, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
71             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, b, (__privtrans->pdls[1]), 1, PDL_PPSYM_OP) \
72             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, c, (__privtrans->pdls[2]), 1, PDL_PPSYM_OP)
73             #endif
74 2 100         if ( __privtrans->bvalflag ) { /* ** do 'bad' Code ** */
75             #define PDL_BAD_CODE
76             #define PDL_IF_BAD(t,f) t
77 1           switch (__privtrans->__datatype) { /* Start generic switch */
78 0           case PDL_F: {
79 0 0         PDL_DECLARE_PARAMS_dvrs_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
80 0 0         PDL_BROADCASTLOOP_START_dvrs_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
81 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0],a_badval,F,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,F,b_badval_isnan)) { (c_datap)[0]=c_badval; continue; },)
    0          
    0          
    0          
    0          
    0          
82 0           (c_datap)[0] = (a_datap)[0]? sqrt( -2 * log((b_datap)[0]) )
83 0 0         : -1 * sqrt( -2 * log(1-(b_datap)[0]) )
84             ;
85 0 0         }PDL_BROADCASTLOOP_END_dvrs_readdata
    0          
86 0           } break;
87 1           case PDL_D: {
88 1 50         PDL_DECLARE_PARAMS_dvrs_1(PDL_Double,D)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
89 9 50         PDL_BROADCASTLOOP_START_dvrs_readdata {
    50          
    50          
    50          
    50          
    100          
    100          
90 6 50         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0],a_badval,D,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,D,b_badval_isnan)) { (c_datap)[0]=c_badval; continue; },)
    0          
    100          
    50          
    0          
    50          
91 5           (c_datap)[0] = (a_datap)[0]? sqrt( -2 * log((b_datap)[0]) )
92 5 100         : -1 * sqrt( -2 * log(1-(b_datap)[0]) )
93             ;
94 1 50         }PDL_BROADCASTLOOP_END_dvrs_readdata
    50          
95 1           } break;
96 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in dvrs: unhandled datatype(%d), only handles (FD)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
97             }
98             #undef PDL_BAD_CODE
99             #undef PDL_IF_BAD
100             } else { /* ** else do 'good' Code ** */
101             #define PDL_IF_BAD(t,f) f
102 1           switch (__privtrans->__datatype) { /* Start generic switch */
103 0           case PDL_F: {
104 0 0         PDL_DECLARE_PARAMS_dvrs_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
105 0 0         PDL_BROADCASTLOOP_START_dvrs_readdata {
    0          
    0          
    0          
    0          
    0          
    0          
106             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0],a_badval,F,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,F,b_badval_isnan)) { (c_datap)[0]=c_badval; continue; },)
107 0           (c_datap)[0] = (a_datap)[0]? sqrt( -2 * log((b_datap)[0]) )
108 0 0         : -1 * sqrt( -2 * log(1-(b_datap)[0]) )
109             ;
110 0 0         }PDL_BROADCASTLOOP_END_dvrs_readdata
    0          
111 0           } break;
112 1           case PDL_D: {
113 1 50         PDL_DECLARE_PARAMS_dvrs_1(PDL_Double,D)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
114 8 50         PDL_BROADCASTLOOP_START_dvrs_readdata {
    50          
    50          
    50          
    50          
    100          
    100          
115             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0],a_badval,D,a_badval_isnan) || PDL_ISBAD2((b_datap)[0],b_badval,D,b_badval_isnan)) { (c_datap)[0]=c_badval; continue; },)
116 5           (c_datap)[0] = (a_datap)[0]? sqrt( -2 * log((b_datap)[0]) )
117 5 100         : -1 * sqrt( -2 * log(1-(b_datap)[0]) )
118             ;
119 1 50         }PDL_BROADCASTLOOP_END_dvrs_readdata
    50          
120 1           } break;
121 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in dvrs: unhandled datatype(%d), only handles (FD)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
122             }
123             #undef PDL_IF_BAD
124             }
125 2           return PDL_err;
126             }
127              
128             static pdl_datatypes pdl_dvrs_vtable_gentypes[] = { PDL_F, PDL_D, -1 };
129             static PDL_Indx pdl_dvrs_vtable_realdims[] = { 0, 0, 0 };
130             static char *pdl_dvrs_vtable_parnames[] = { "a","b","c" };
131             static short pdl_dvrs_vtable_parflags[] = {
132             0,
133             0,
134             PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
135             };
136             static pdl_datatypes pdl_dvrs_vtable_partypes[] = { -1, -1, -1 };
137             static PDL_Indx pdl_dvrs_vtable_realdims_starts[] = { 0, 0, 0 };
138             static PDL_Indx pdl_dvrs_vtable_realdims_ind_ids[] = { 0 };
139             static char *pdl_dvrs_vtable_indnames[] = { "" };
140             pdl_transvtable pdl_dvrs_vtable = {
141             PDL_TRANS_DO_BROADCAST|PDL_TRANS_BADPROCESS, 0, pdl_dvrs_vtable_gentypes, 2, 3, NULL /*CORE21*/,
142             pdl_dvrs_vtable_realdims, pdl_dvrs_vtable_parnames,
143             pdl_dvrs_vtable_parflags, pdl_dvrs_vtable_partypes,
144             pdl_dvrs_vtable_realdims_starts, pdl_dvrs_vtable_realdims_ind_ids, 0,
145             0, pdl_dvrs_vtable_indnames,
146             NULL, pdl_dvrs_readdata, NULL,
147             NULL,
148             0,"PDL::Stats::GLM::dvrs"
149             };
150              
151              
152 2           pdl_error pdl_run_dvrs(pdl *a,pdl *b,pdl *c) {
153 2           pdl_error PDL_err = {0, NULL, 0};
154 2 50         if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
155 2           pdl_trans *__privtrans = PDL->create_trans(&pdl_dvrs_vtable);
156 2 50         if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
157 2           __privtrans->pdls[0] = a;
158 2           __privtrans->pdls[1] = b;
159 2           __privtrans->pdls[2] = c;
160 2 50         PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
161 2 50         PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
162 2           return PDL_err;
163             }