File Coverage

lib/PDL/Math-pp-asin.c
Criterion Covered Total %
statement 37 94 39.3
branch 67 504 13.2
condition n/a
subroutine n/a
pod n/a
total 104 598 17.3


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/Math.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_Math
21             extern Core* PDL; /* Structure hold core C functions */
22             #line 23 "lib/PDL/Math-pp-asin.c"
23              
24             #include
25             #include "protos.h"
26             #include "cpoly.h"
27              
28             #line 1857 "lib/PDL/PP.pm"
29             pdl_error pdl_asin_readdata(pdl_trans *__privtrans) {
30             pdl_error PDL_err = {0, NULL, 0};
31             #line 32 "lib/PDL/Math-pp-asin.c"
32 8 50         if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in asin:" "broadcast.incs NULL");
33             /* broadcastloop declarations */
34             int __brcloopval;
35             register PDL_Indx __tind0,__tind1; /* counters along dim */
36 8           register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
37             /* dims here are how many steps along those dims */
38 8           register PDL_Indx __tinc0_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
39 8           register PDL_Indx __tinc0_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
40 8           register PDL_Indx __tinc1_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
41 8           register PDL_Indx __tinc1_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
42             #define PDL_BROADCASTLOOP_START_asin_readdata PDL_BROADCASTLOOP_START( \
43             readdata, \
44             __privtrans->broadcast, \
45             __privtrans->vtable, \
46             a_datap += __offsp[0]; \
47             b_datap += __offsp[1]; \
48             , \
49             ( ,a_datap += __tinc1_a - __tinc0_a * __tdims0 \
50             ,b_datap += __tinc1_b - __tinc0_b * __tdims0 \
51             ), \
52             ( ,a_datap += __tinc0_a \
53             ,b_datap += __tinc0_b \
54             ) \
55             )
56             #define PDL_BROADCASTLOOP_END_asin_readdata PDL_BROADCASTLOOP_END( \
57             __privtrans->broadcast, \
58             a_datap -= __tinc1_a * __tdims1 + __offsp[0]; \
59             b_datap -= __tinc1_b * __tdims1 + __offsp[1]; \
60             )
61             #ifndef PDL_DECLARE_PARAMS_asin_1
62             #define PDL_DECLARE_PARAMS_asin_1(PDL_TYPE_OP,PDL_PPSYM_OP) \
63             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, a, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
64             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, b, (__privtrans->pdls[1]), 1, PDL_PPSYM_OP)
65             #endif
66 8 50         if ( __privtrans->bvalflag ) { /* ** do 'bad' Code ** */
67             #define PDL_BAD_CODE
68             #define PDL_IF_BAD(t,f) t
69 0           switch (__privtrans->__datatype) { /* Start generic switch */
70 0           case PDL_F: {
71 0 0         PDL_DECLARE_PARAMS_asin_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
72 0 0         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,F,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
73 0           (b_datap)[0] = asin((a_datap)[0]);
74 0 0         }PDL_BROADCASTLOOP_END_asin_readdata
    0          
75 0           } break;
76 0           case PDL_D: {
77 0 0         PDL_DECLARE_PARAMS_asin_1(PDL_Double,D)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
78 0 0         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,D,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
79 0           (b_datap)[0] = asin((a_datap)[0]);
80 0 0         }PDL_BROADCASTLOOP_END_asin_readdata
    0          
81 0           } break;
82 0           case PDL_LD: {
83 0 0         PDL_DECLARE_PARAMS_asin_1(PDL_LDouble,E)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
84 0 0         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,E,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
85 0           (b_datap)[0] = asin((a_datap)[0]);
86 0 0         }PDL_BROADCASTLOOP_END_asin_readdata
    0          
87 0           } break;
88 0           case PDL_CF: {
89 0 0         PDL_DECLARE_PARAMS_asin_1(PDL_CFloat,G)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
90 0 0         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,G,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
91 0           (b_datap)[0] = asin((a_datap)[0]);
92 0 0         }PDL_BROADCASTLOOP_END_asin_readdata
    0          
93 0           } break;
94 0           case PDL_CD: {
95 0 0         PDL_DECLARE_PARAMS_asin_1(PDL_CDouble,C)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
96 0 0         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,C,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
97 0           (b_datap)[0] = asin((a_datap)[0]);
98 0 0         }PDL_BROADCASTLOOP_END_asin_readdata
    0          
99 0           } break;
100 0           case PDL_CLD: {
101 0 0         PDL_DECLARE_PARAMS_asin_1(PDL_CLDouble,H)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
102 0 0         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,H,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
103 0           (b_datap)[0] = asin((a_datap)[0]);
104 0 0         }PDL_BROADCASTLOOP_END_asin_readdata
    0          
105 0           } break;
106 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in asin: unhandled datatype(%d), only handles (FDEGCH)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
107             }
108             #undef PDL_BAD_CODE
109             #undef PDL_IF_BAD
110             } else { /* ** else do 'good' Code ** */
111             #define PDL_IF_BAD(t,f) f
112 8           switch (__privtrans->__datatype) { /* Start generic switch */
113 0           case PDL_F: {
114 0 0         PDL_DECLARE_PARAMS_asin_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
115 0 0         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,F,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    0          
    0          
    0          
    0          
    0          
    0          
116 0           (b_datap)[0] = asin((a_datap)[0]);
117 0 0         }PDL_BROADCASTLOOP_END_asin_readdata
    0          
118 0           } break;
119 5           case PDL_D: {
120 5 50         PDL_DECLARE_PARAMS_asin_1(PDL_Double,D)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
121 15 50         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,D,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    50          
    50          
    50          
    50          
    100          
    100          
122 5           (b_datap)[0] = asin((a_datap)[0]);
123 5 50         }PDL_BROADCASTLOOP_END_asin_readdata
    50          
124 5           } break;
125 2           case PDL_LD: {
126 2 50         PDL_DECLARE_PARAMS_asin_1(PDL_LDouble,E)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
127 6 50         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,E,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    50          
    50          
    50          
    50          
    100          
    100          
128 2           (b_datap)[0] = asin((a_datap)[0]);
129 2 50         }PDL_BROADCASTLOOP_END_asin_readdata
    50          
130 2           } break;
131 0           case PDL_CF: {
132 0 0         PDL_DECLARE_PARAMS_asin_1(PDL_CFloat,G)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
133 0 0         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,G,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    0          
    0          
    0          
    0          
    0          
    0          
134 0           (b_datap)[0] = asin((a_datap)[0]);
135 0 0         }PDL_BROADCASTLOOP_END_asin_readdata
    0          
136 0           } break;
137 1           case PDL_CD: {
138 1 50         PDL_DECLARE_PARAMS_asin_1(PDL_CDouble,C)
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
    50          
139 3 50         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,C,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    50          
    50          
    50          
    50          
    100          
    100          
140 1           (b_datap)[0] = asin((a_datap)[0]);
141 1 50         }PDL_BROADCASTLOOP_END_asin_readdata
    50          
142 1           } break;
143 0           case PDL_CLD: {
144 0 0         PDL_DECLARE_PARAMS_asin_1(PDL_CLDouble,H)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
145 0 0         PDL_BROADCASTLOOP_START_asin_readdata {PDL_IF_BAD(if ( PDL_ISBAD2((a_datap)[0],a_badval,H,a_badval_isnan) ) { (b_datap)[0]=b_badval; } else,)
    0          
    0          
    0          
    0          
    0          
    0          
146 0           (b_datap)[0] = asin((a_datap)[0]);
147 0 0         }PDL_BROADCASTLOOP_END_asin_readdata
    0          
148 0           } break;
149 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in asin: unhandled datatype(%d), only handles (FDEGCH)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
150             }
151             #undef PDL_IF_BAD
152             }
153 8           return PDL_err;
154             }
155              
156             static pdl_datatypes pdl_asin_vtable_gentypes[] = { PDL_CF, PDL_CD, PDL_CLD, PDL_F, PDL_D, PDL_LD, -1 };
157             static PDL_Indx pdl_asin_vtable_realdims[] = { 0, 0 };
158             static char *pdl_asin_vtable_parnames[] = { "a","b" };
159             static short pdl_asin_vtable_parflags[] = {
160             0,
161             PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
162             };
163             static pdl_datatypes pdl_asin_vtable_partypes[] = { -1, -1 };
164             static PDL_Indx pdl_asin_vtable_realdims_starts[] = { 0, 0 };
165             static PDL_Indx pdl_asin_vtable_realdims_ind_ids[] = { 0 };
166             static char *pdl_asin_vtable_indnames[] = { "" };
167             pdl_transvtable pdl_asin_vtable = {
168             PDL_TRANS_DO_BROADCAST|PDL_TRANS_BADPROCESS, 0, pdl_asin_vtable_gentypes, 1, 2, NULL /*CORE21*/,
169             pdl_asin_vtable_realdims, pdl_asin_vtable_parnames,
170             pdl_asin_vtable_parflags, pdl_asin_vtable_partypes,
171             pdl_asin_vtable_realdims_starts, pdl_asin_vtable_realdims_ind_ids, 0,
172             0, pdl_asin_vtable_indnames,
173             NULL, pdl_asin_readdata, NULL,
174             NULL,
175             0,"PDL::Math::asin"
176             };
177              
178              
179 8           pdl_error pdl_run_asin(pdl *a,pdl *b) {
180 8           pdl_error PDL_err = {0, NULL, 0};
181 8 50         if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
182 8           pdl_trans *__privtrans = PDL->create_trans(&pdl_asin_vtable);
183 8 50         if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
184 8           __privtrans->pdls[0] = a;
185 8           __privtrans->pdls[1] = b;
186 8 50         PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
187 8 50         PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
188 8           return PDL_err;
189             }