File Coverage

lib/PDL/Primitive-pp-inner2t.c
Criterion Covered Total %
statement 0 549 0.0
branch 0 2628 0.0
condition n/a
subroutine n/a
pod n/a
total 0 3177 0.0


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/Primitive.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_Primitive
21             extern Core* PDL; /* Structure hold core C functions */
22              
23             #line 1857 "lib/PDL/PP.pm"
24             pdl_error pdl_inner2t_readdata(pdl_trans *__privtrans) {
25             pdl_error PDL_err = {0, NULL, 0};
26             #line 27 "lib/PDL/Primitive-pp-inner2t.c"
27 0           register PDL_Indx __j_size = __privtrans->ind_sizes[0];
28 0           register PDL_Indx __k_size = __privtrans->ind_sizes[1];
29 0           register PDL_Indx __m_size = __privtrans->ind_sizes[2];
30 0           register PDL_Indx __n_size = __privtrans->ind_sizes[3];
31 0 0         if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in inner2t:" "broadcast.incs NULL");
32             /* broadcastloop declarations */
33             int __brcloopval;
34             register PDL_Indx __tind0,__tind1; /* counters along dim */
35 0           register PDL_Indx __tnpdls = __privtrans->broadcast.npdls;
36             /* dims here are how many steps along those dims */
37 0           register PDL_Indx __tinc0_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,0);
38 0           register PDL_Indx __tinc0_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,0);
39 0           register PDL_Indx __tinc0_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,0);
40 0           register PDL_Indx __tinc0_tmp = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,3,0);
41 0           register PDL_Indx __tinc0_d = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,4,0);
42 0           register PDL_Indx __tinc1_a = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,0,1);
43 0           register PDL_Indx __tinc1_b = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,1,1);
44 0           register PDL_Indx __tinc1_c = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,2,1);
45 0           register PDL_Indx __tinc1_tmp = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,3,1);
46 0           register PDL_Indx __tinc1_d = PDL_BRC_INC(__privtrans->broadcast.incs,__tnpdls,4,1);
47             #define PDL_BROADCASTLOOP_START_inner2t_readdata PDL_BROADCASTLOOP_START( \
48             readdata, \
49             __privtrans->broadcast, \
50             __privtrans->vtable, \
51             a_datap += __offsp[0]; \
52             b_datap += __offsp[1]; \
53             c_datap += __offsp[2]; \
54             tmp_datap += __offsp[3]; \
55             d_datap += __offsp[4]; \
56             , \
57             ( ,a_datap += __tinc1_a - __tinc0_a * __tdims0 \
58             ,b_datap += __tinc1_b - __tinc0_b * __tdims0 \
59             ,c_datap += __tinc1_c - __tinc0_c * __tdims0 \
60             ,tmp_datap += __tinc1_tmp - __tinc0_tmp * __tdims0 \
61             ,d_datap += __tinc1_d - __tinc0_d * __tdims0 \
62             ), \
63             ( ,a_datap += __tinc0_a \
64             ,b_datap += __tinc0_b \
65             ,c_datap += __tinc0_c \
66             ,tmp_datap += __tinc0_tmp \
67             ,d_datap += __tinc0_d \
68             ) \
69             )
70             #define PDL_BROADCASTLOOP_END_inner2t_readdata PDL_BROADCASTLOOP_END( \
71             __privtrans->broadcast, \
72             a_datap -= __tinc1_a * __tdims1 + __offsp[0]; \
73             b_datap -= __tinc1_b * __tdims1 + __offsp[1]; \
74             c_datap -= __tinc1_c * __tdims1 + __offsp[2]; \
75             tmp_datap -= __tinc1_tmp * __tdims1 + __offsp[3]; \
76             d_datap -= __tinc1_d * __tdims1 + __offsp[4]; \
77             )
78 0           register PDL_Indx __inc_a_j = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,0,0)]; (void)__inc_a_j;register PDL_Indx __inc_a_n = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,0,1)]; (void)__inc_a_n;
79 0           register PDL_Indx __inc_b_n = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,1,0)]; (void)__inc_b_n;register PDL_Indx __inc_b_m = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,1,1)]; (void)__inc_b_m;
80 0           register PDL_Indx __inc_c_m = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,2,0)]; (void)__inc_c_m;register PDL_Indx __inc_c_k = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,2,1)]; (void)__inc_c_k;
81 0           register PDL_Indx __inc_d_j = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,4,0)]; (void)__inc_d_j;register PDL_Indx __inc_d_k = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,4,1)]; (void)__inc_d_k;
82 0           register PDL_Indx __inc_tmp_n = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,3,0)]; (void)__inc_tmp_n;register PDL_Indx __inc_tmp_k = __privtrans->inc_sizes[PDL_INC_ID(__privtrans->vtable,3,1)]; (void)__inc_tmp_k;
83             #ifndef PDL_DECLARE_PARAMS_inner2t_1
84             #define PDL_DECLARE_PARAMS_inner2t_1(PDL_TYPE_OP,PDL_PPSYM_OP) \
85             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, a, (__privtrans->pdls[0]), 1, PDL_PPSYM_OP) \
86             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, b, (__privtrans->pdls[1]), 1, PDL_PPSYM_OP) \
87             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, c, (__privtrans->pdls[2]), 1, PDL_PPSYM_OP) \
88             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, tmp, (__privtrans->pdls[3]), 1, PDL_PPSYM_OP) \
89             PDL_DECLARE_PARAMETER_BADVAL(PDL_TYPE_OP, d, (__privtrans->pdls[4]), 1, PDL_PPSYM_OP)
90             #endif
91 0 0         if ( __privtrans->bvalflag ) { /* ** do 'bad' Code ** */
92             #define PDL_BAD_CODE
93             #define PDL_IF_BAD(t,f) t
94 0           switch (__privtrans->__datatype) { /* Start generic switch */
95 0           case PDL_SB: {
96 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_SByte,A)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
97 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
98 0           complex long double tmp0 = 0;
99 0           PDL_IF_BAD(int flag = 0;,)
100 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
101 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,A,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,A,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
102 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
103             }} /* Close m */
104 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
105 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
106             }} /* Close n */}} /* Close k */
107 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
108 0           complex long double tmp1 = 0;
109 0           PDL_IF_BAD(int flag = 0;,)
110 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
111 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,A,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,A,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
112 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
113             }} /* Close n */
114 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
115 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
116             }} /* Close j */}} /* Close k */
117 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
118 0           } break;
119 0           case PDL_B: {
120 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Byte,B)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
121 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
122 0           complex long double tmp0 = 0;
123 0           PDL_IF_BAD(int flag = 0;,)
124 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
125 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,B,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,B,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
126 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
127             }} /* Close m */
128 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
129 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
130             }} /* Close n */}} /* Close k */
131 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
132 0           complex long double tmp1 = 0;
133 0           PDL_IF_BAD(int flag = 0;,)
134 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
135 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,B,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,B,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
136 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
137             }} /* Close n */
138 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
139 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
140             }} /* Close j */}} /* Close k */
141 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
142 0           } break;
143 0           case PDL_S: {
144 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Short,S)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
145 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
146 0           complex long double tmp0 = 0;
147 0           PDL_IF_BAD(int flag = 0;,)
148 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
149 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,S,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,S,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
150 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
151             }} /* Close m */
152 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
153 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
154             }} /* Close n */}} /* Close k */
155 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
156 0           complex long double tmp1 = 0;
157 0           PDL_IF_BAD(int flag = 0;,)
158 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
159 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,S,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,S,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
160 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
161             }} /* Close n */
162 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
163 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
164             }} /* Close j */}} /* Close k */
165 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
166 0           } break;
167 0           case PDL_US: {
168 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Ushort,U)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
169 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
170 0           complex long double tmp0 = 0;
171 0           PDL_IF_BAD(int flag = 0;,)
172 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
173 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,U,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,U,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
174 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
175             }} /* Close m */
176 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
177 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
178             }} /* Close n */}} /* Close k */
179 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
180 0           complex long double tmp1 = 0;
181 0           PDL_IF_BAD(int flag = 0;,)
182 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
183 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,U,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,U,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
184 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
185             }} /* Close n */
186 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
187 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
188             }} /* Close j */}} /* Close k */
189 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
190 0           } break;
191 0           case PDL_L: {
192 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
193 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
194 0           complex long double tmp0 = 0;
195 0           PDL_IF_BAD(int flag = 0;,)
196 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
197 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,L,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,L,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
198 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
199             }} /* Close m */
200 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
201 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
202             }} /* Close n */}} /* Close k */
203 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
204 0           complex long double tmp1 = 0;
205 0           PDL_IF_BAD(int flag = 0;,)
206 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
207 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,L,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,L,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
208 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
209             }} /* Close n */
210 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
211 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
212             }} /* Close j */}} /* Close k */
213 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
214 0           } break;
215 0           case PDL_UL: {
216 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_ULong,K)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
217 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
218 0           complex long double tmp0 = 0;
219 0           PDL_IF_BAD(int flag = 0;,)
220 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
221 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,K,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,K,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
222 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
223             }} /* Close m */
224 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
225 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
226             }} /* Close n */}} /* Close k */
227 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
228 0           complex long double tmp1 = 0;
229 0           PDL_IF_BAD(int flag = 0;,)
230 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
231 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,K,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,K,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
232 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
233             }} /* Close n */
234 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
235 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
236             }} /* Close j */}} /* Close k */
237 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
238 0           } break;
239 0           case PDL_IND: {
240 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
241 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
242 0           complex long double tmp0 = 0;
243 0           PDL_IF_BAD(int flag = 0;,)
244 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
245 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,N,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,N,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
246 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
247             }} /* Close m */
248 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
249 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
250             }} /* Close n */}} /* Close k */
251 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
252 0           complex long double tmp1 = 0;
253 0           PDL_IF_BAD(int flag = 0;,)
254 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
255 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,N,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,N,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
256 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
257             }} /* Close n */
258 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
259 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
260             }} /* Close j */}} /* Close k */
261 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
262 0           } break;
263 0           case PDL_ULL: {
264 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_ULongLong,P)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
265 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
266 0           complex long double tmp0 = 0;
267 0           PDL_IF_BAD(int flag = 0;,)
268 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
269 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,P,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,P,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
270 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
271             }} /* Close m */
272 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
273 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
274             }} /* Close n */}} /* Close k */
275 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
276 0           complex long double tmp1 = 0;
277 0           PDL_IF_BAD(int flag = 0;,)
278 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
279 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,P,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,P,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
280 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
281             }} /* Close n */
282 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
283 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
284             }} /* Close j */}} /* Close k */
285 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
286 0           } break;
287 0           case PDL_LL: {
288 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_LongLong,Q)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
289 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
290 0           complex long double tmp0 = 0;
291 0           PDL_IF_BAD(int flag = 0;,)
292 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
293 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,Q,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,Q,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
294 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
295             }} /* Close m */
296 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
297 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
298             }} /* Close n */}} /* Close k */
299 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
300 0           complex long double tmp1 = 0;
301 0           PDL_IF_BAD(int flag = 0;,)
302 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
303 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,Q,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,Q,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
304 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
305             }} /* Close n */
306 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
307 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
308             }} /* Close j */}} /* Close k */
309 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
310 0           } break;
311 0           case PDL_F: {
312 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
313 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
314 0           complex long double tmp0 = 0;
315 0           PDL_IF_BAD(int flag = 0;,)
316 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
317 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,F,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,F,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
318 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
319             }} /* Close m */
320 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
321 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
322             }} /* Close n */}} /* Close k */
323 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
324 0           complex long double tmp1 = 0;
325 0           PDL_IF_BAD(int flag = 0;,)
326 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
327 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,F,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,F,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
328 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
329             }} /* Close n */
330 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
331 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
332             }} /* Close j */}} /* Close k */
333 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
334 0           } break;
335 0           case PDL_D: {
336 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Double,D)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
337 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
338 0           complex long double tmp0 = 0;
339 0           PDL_IF_BAD(int flag = 0;,)
340 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
341 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,D,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,D,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
342 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
343             }} /* Close m */
344 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
345 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
346             }} /* Close n */}} /* Close k */
347 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
348 0           complex long double tmp1 = 0;
349 0           PDL_IF_BAD(int flag = 0;,)
350 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
351 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,D,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,D,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
352 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
353             }} /* Close n */
354 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
355 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
356             }} /* Close j */}} /* Close k */
357 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
358 0           } break;
359 0           case PDL_LD: {
360 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_LDouble,E)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
361 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
362 0           complex long double tmp0 = 0;
363 0           PDL_IF_BAD(int flag = 0;,)
364 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
365 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,E,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,E,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
366 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
367             }} /* Close m */
368 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
369 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
370             }} /* Close n */}} /* Close k */
371 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
372 0           complex long double tmp1 = 0;
373 0           PDL_IF_BAD(int flag = 0;,)
374 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
375 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,E,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,E,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
376 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
377             }} /* Close n */
378 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
379 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
380             }} /* Close j */}} /* Close k */
381 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
382 0           } break;
383 0           case PDL_CF: {
384 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_CFloat,G)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
385 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
386 0           complex long double tmp0 = 0;
387 0           PDL_IF_BAD(int flag = 0;,)
388 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
389 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,G,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,G,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
390 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
391             }} /* Close m */
392 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
393 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
394             }} /* Close n */}} /* Close k */
395 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
396 0           complex long double tmp1 = 0;
397 0           PDL_IF_BAD(int flag = 0;,)
398 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
399 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,G,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,G,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
400 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
401             }} /* Close n */
402 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
403 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
404             }} /* Close j */}} /* Close k */
405 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
406 0           } break;
407 0           case PDL_CD: {
408 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_CDouble,C)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
409 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
410 0           complex long double tmp0 = 0;
411 0           PDL_IF_BAD(int flag = 0;,)
412 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
413 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,C,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,C,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
414 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
415             }} /* Close m */
416 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
417 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
418             }} /* Close n */}} /* Close k */
419 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
420 0           complex long double tmp1 = 0;
421 0           PDL_IF_BAD(int flag = 0;,)
422 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
423 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,C,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,C,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
424 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
425             }} /* Close n */
426 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
427 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
428             }} /* Close j */}} /* Close k */
429 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
430 0           } break;
431 0           case PDL_CLD: {
432 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_CLDouble,H)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
433 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
434 0           complex long double tmp0 = 0;
435 0           PDL_IF_BAD(int flag = 0;,)
436 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
437 0 0         PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,H,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,H,c_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
438 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
439             }} /* Close m */
440 0 0         PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
441 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
442             }} /* Close n */}} /* Close k */
443 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
444 0           complex long double tmp1 = 0;
445 0           PDL_IF_BAD(int flag = 0;,)
446 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
447 0 0         PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,H,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,H,tmp_badval_isnan)) continue;flag = 1;,)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
448 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
449             }} /* Close n */
450 0 0         PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
451 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
452             }} /* Close j */}} /* Close k */
453 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
454 0           } break;
455 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in inner2t: unhandled datatype(%d), only handles (ABSULKNPQFDEGCH)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
456             }
457             #undef PDL_BAD_CODE
458             #undef PDL_IF_BAD
459             } else { /* ** else do 'good' Code ** */
460             #define PDL_IF_BAD(t,f) f
461 0           switch (__privtrans->__datatype) { /* Start generic switch */
462 0           case PDL_SB: {
463 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_SByte,A)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
464 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
465 0           complex long double tmp0 = 0;
466             PDL_IF_BAD(int flag = 0;,)
467 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
468             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,A,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,A,c_badval_isnan)) continue;flag = 1;,)
469 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
470             }} /* Close m */
471             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
472 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
473             }} /* Close n */}} /* Close k */
474 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
475 0           complex long double tmp1 = 0;
476             PDL_IF_BAD(int flag = 0;,)
477 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
478             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,A,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,A,tmp_badval_isnan)) continue;flag = 1;,)
479 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
480             }} /* Close n */
481             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
482 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
483             }} /* Close j */}} /* Close k */
484 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
485 0           } break;
486 0           case PDL_B: {
487 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Byte,B)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
488 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
489 0           complex long double tmp0 = 0;
490             PDL_IF_BAD(int flag = 0;,)
491 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
492             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,B,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,B,c_badval_isnan)) continue;flag = 1;,)
493 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
494             }} /* Close m */
495             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
496 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
497             }} /* Close n */}} /* Close k */
498 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
499 0           complex long double tmp1 = 0;
500             PDL_IF_BAD(int flag = 0;,)
501 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
502             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,B,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,B,tmp_badval_isnan)) continue;flag = 1;,)
503 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
504             }} /* Close n */
505             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
506 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
507             }} /* Close j */}} /* Close k */
508 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
509 0           } break;
510 0           case PDL_S: {
511 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Short,S)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
512 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
513 0           complex long double tmp0 = 0;
514             PDL_IF_BAD(int flag = 0;,)
515 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
516             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,S,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,S,c_badval_isnan)) continue;flag = 1;,)
517 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
518             }} /* Close m */
519             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
520 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
521             }} /* Close n */}} /* Close k */
522 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
523 0           complex long double tmp1 = 0;
524             PDL_IF_BAD(int flag = 0;,)
525 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
526             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,S,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,S,tmp_badval_isnan)) continue;flag = 1;,)
527 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
528             }} /* Close n */
529             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
530 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
531             }} /* Close j */}} /* Close k */
532 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
533 0           } break;
534 0           case PDL_US: {
535 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Ushort,U)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
536 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
537 0           complex long double tmp0 = 0;
538             PDL_IF_BAD(int flag = 0;,)
539 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
540             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,U,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,U,c_badval_isnan)) continue;flag = 1;,)
541 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
542             }} /* Close m */
543             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
544 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
545             }} /* Close n */}} /* Close k */
546 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
547 0           complex long double tmp1 = 0;
548             PDL_IF_BAD(int flag = 0;,)
549 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
550             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,U,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,U,tmp_badval_isnan)) continue;flag = 1;,)
551 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
552             }} /* Close n */
553             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
554 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
555             }} /* Close j */}} /* Close k */
556 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
557 0           } break;
558 0           case PDL_L: {
559 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Long,L)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
560 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
561 0           complex long double tmp0 = 0;
562             PDL_IF_BAD(int flag = 0;,)
563 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
564             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,L,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,L,c_badval_isnan)) continue;flag = 1;,)
565 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
566             }} /* Close m */
567             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
568 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
569             }} /* Close n */}} /* Close k */
570 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
571 0           complex long double tmp1 = 0;
572             PDL_IF_BAD(int flag = 0;,)
573 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
574             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,L,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,L,tmp_badval_isnan)) continue;flag = 1;,)
575 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
576             }} /* Close n */
577             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
578 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
579             }} /* Close j */}} /* Close k */
580 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
581 0           } break;
582 0           case PDL_UL: {
583 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_ULong,K)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
584 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
585 0           complex long double tmp0 = 0;
586             PDL_IF_BAD(int flag = 0;,)
587 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
588             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,K,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,K,c_badval_isnan)) continue;flag = 1;,)
589 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
590             }} /* Close m */
591             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
592 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
593             }} /* Close n */}} /* Close k */
594 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
595 0           complex long double tmp1 = 0;
596             PDL_IF_BAD(int flag = 0;,)
597 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
598             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,K,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,K,tmp_badval_isnan)) continue;flag = 1;,)
599 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
600             }} /* Close n */
601             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
602 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
603             }} /* Close j */}} /* Close k */
604 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
605 0           } break;
606 0           case PDL_IND: {
607 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Indx,N)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
608 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
609 0           complex long double tmp0 = 0;
610             PDL_IF_BAD(int flag = 0;,)
611 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
612             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,N,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,N,c_badval_isnan)) continue;flag = 1;,)
613 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
614             }} /* Close m */
615             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
616 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
617             }} /* Close n */}} /* Close k */
618 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
619 0           complex long double tmp1 = 0;
620             PDL_IF_BAD(int flag = 0;,)
621 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
622             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,N,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,N,tmp_badval_isnan)) continue;flag = 1;,)
623 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
624             }} /* Close n */
625             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
626 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
627             }} /* Close j */}} /* Close k */
628 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
629 0           } break;
630 0           case PDL_ULL: {
631 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_ULongLong,P)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
632 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
633 0           complex long double tmp0 = 0;
634             PDL_IF_BAD(int flag = 0;,)
635 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
636             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,P,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,P,c_badval_isnan)) continue;flag = 1;,)
637 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
638             }} /* Close m */
639             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
640 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
641             }} /* Close n */}} /* Close k */
642 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
643 0           complex long double tmp1 = 0;
644             PDL_IF_BAD(int flag = 0;,)
645 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
646             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,P,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,P,tmp_badval_isnan)) continue;flag = 1;,)
647 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
648             }} /* Close n */
649             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
650 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
651             }} /* Close j */}} /* Close k */
652 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
653 0           } break;
654 0           case PDL_LL: {
655 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_LongLong,Q)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
656 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
657 0           complex long double tmp0 = 0;
658             PDL_IF_BAD(int flag = 0;,)
659 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
660             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,Q,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,Q,c_badval_isnan)) continue;flag = 1;,)
661 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
662             }} /* Close m */
663             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
664 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
665             }} /* Close n */}} /* Close k */
666 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
667 0           complex long double tmp1 = 0;
668             PDL_IF_BAD(int flag = 0;,)
669 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
670             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,Q,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,Q,tmp_badval_isnan)) continue;flag = 1;,)
671 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
672             }} /* Close n */
673             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
674 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
675             }} /* Close j */}} /* Close k */
676 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
677 0           } break;
678 0           case PDL_F: {
679 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Float,F)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
680 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
681 0           complex long double tmp0 = 0;
682             PDL_IF_BAD(int flag = 0;,)
683 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
684             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,F,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,F,c_badval_isnan)) continue;flag = 1;,)
685 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
686             }} /* Close m */
687             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
688 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
689             }} /* Close n */}} /* Close k */
690 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
691 0           complex long double tmp1 = 0;
692             PDL_IF_BAD(int flag = 0;,)
693 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
694             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,F,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,F,tmp_badval_isnan)) continue;flag = 1;,)
695 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
696             }} /* Close n */
697             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
698 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
699             }} /* Close j */}} /* Close k */
700 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
701 0           } break;
702 0           case PDL_D: {
703 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_Double,D)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
704 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
705 0           complex long double tmp0 = 0;
706             PDL_IF_BAD(int flag = 0;,)
707 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
708             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,D,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,D,c_badval_isnan)) continue;flag = 1;,)
709 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
710             }} /* Close m */
711             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
712 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
713             }} /* Close n */}} /* Close k */
714 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
715 0           complex long double tmp1 = 0;
716             PDL_IF_BAD(int flag = 0;,)
717 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
718             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,D,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,D,tmp_badval_isnan)) continue;flag = 1;,)
719 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
720             }} /* Close n */
721             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
722 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
723             }} /* Close j */}} /* Close k */
724 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
725 0           } break;
726 0           case PDL_LD: {
727 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_LDouble,E)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
728 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
729 0           complex long double tmp0 = 0;
730             PDL_IF_BAD(int flag = 0;,)
731 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
732             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,E,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,E,c_badval_isnan)) continue;flag = 1;,)
733 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
734             }} /* Close m */
735             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
736 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
737             }} /* Close n */}} /* Close k */
738 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
739 0           complex long double tmp1 = 0;
740             PDL_IF_BAD(int flag = 0;,)
741 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
742             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,E,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,E,tmp_badval_isnan)) continue;flag = 1;,)
743 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
744             }} /* Close n */
745             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
746 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
747             }} /* Close j */}} /* Close k */
748 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
749 0           } break;
750 0           case PDL_CF: {
751 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_CFloat,G)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
752 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
753 0           complex long double tmp0 = 0;
754             PDL_IF_BAD(int flag = 0;,)
755 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
756             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,G,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,G,c_badval_isnan)) continue;flag = 1;,)
757 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
758             }} /* Close m */
759             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
760 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
761             }} /* Close n */}} /* Close k */
762 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
763 0           complex long double tmp1 = 0;
764             PDL_IF_BAD(int flag = 0;,)
765 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
766             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,G,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,G,tmp_badval_isnan)) continue;flag = 1;,)
767 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
768             }} /* Close n */
769             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
770 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
771             }} /* Close j */}} /* Close k */
772 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
773 0           } break;
774 0           case PDL_CD: {
775 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_CDouble,C)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
776 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
777 0           complex long double tmp0 = 0;
778             PDL_IF_BAD(int flag = 0;,)
779 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
780             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,C,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,C,c_badval_isnan)) continue;flag = 1;,)
781 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
782             }} /* Close m */
783             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
784 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
785             }} /* Close n */}} /* Close k */
786 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
787 0           complex long double tmp1 = 0;
788             PDL_IF_BAD(int flag = 0;,)
789 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
790             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,C,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,C,tmp_badval_isnan)) continue;flag = 1;,)
791 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
792             }} /* Close n */
793             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
794 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
795             }} /* Close j */}} /* Close k */
796 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
797 0           } break;
798 0           case PDL_CLD: {
799 0 0         PDL_DECLARE_PARAMS_inner2t_1(PDL_CLDouble,H)
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
800 0 0         PDL_BROADCASTLOOP_START_inner2t_readdata {{/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
    0          
    0          
    0          
    0          
    0          
    0          
    0          
801 0           complex long double tmp0 = 0;
802             PDL_IF_BAD(int flag = 0;,)
803 0 0         {/* Open m */ PDL_EXPAND2(register PDL_Indx m=0, __m_stop=(__m_size)); for(; m<__m_stop; m+=1) {
804             PDL_IF_BAD(if (PDL_ISBAD2((b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))],b_badval,H,b_badval_isnan) || PDL_ISBAD2((c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))],c_badval,H,c_badval_isnan)) continue;flag = 1;,)
805 0           tmp0 += (b_datap)[0+(__inc_b_n*(n))+(__inc_b_m*(m))] * (c_datap)[0+(__inc_c_m*(m))+(__inc_c_k*(k))];
806             }} /* Close m */
807             PDL_IF_BAD(if (!flag) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))]=tmp_badval; }
808 0           else,) { (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))] = tmp0; }
809             }} /* Close n */}} /* Close k */
810 0 0         {/* Open j */ PDL_EXPAND2(register PDL_Indx j=0, __j_stop=(__j_size)); for(; j<__j_stop; j+=1) {{/* Open k */ PDL_EXPAND2(register PDL_Indx k=0, __k_stop=(__k_size)); for(; k<__k_stop; k+=1) {
    0          
811 0           complex long double tmp1 = 0;
812             PDL_IF_BAD(int flag = 0;,)
813 0 0         {/* Open n */ PDL_EXPAND2(register PDL_Indx n=0, __n_stop=(__n_size)); for(; n<__n_stop; n+=1) {
814             PDL_IF_BAD(if (PDL_ISBAD2((a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))],a_badval,H,a_badval_isnan) || PDL_ISBAD2((tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))],tmp_badval,H,tmp_badval_isnan)) continue;flag = 1;,)
815 0           tmp1 += (a_datap)[0+(__inc_a_j*(j))+(__inc_a_n*(n))] * (tmp_datap)[0+(__inc_tmp_n*(n))+(__inc_tmp_k*(k))];
816             }} /* Close n */
817             PDL_IF_BAD(if (!flag) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))]=d_badval; }
818 0           else,) { (d_datap)[0+(__inc_d_j*(j))+(__inc_d_k*(k))] = tmp1; }
819             }} /* Close j */}} /* Close k */
820 0 0         }PDL_BROADCASTLOOP_END_inner2t_readdata
    0          
821 0           } break;
822 0           default: return PDL->make_error(PDL_EUSERERROR, "PP INTERNAL ERROR in inner2t: unhandled datatype(%d), only handles (ABSULKNPQFDEGCH)! PLEASE MAKE A BUG REPORT\n", __privtrans->__datatype);
823             }
824             #undef PDL_IF_BAD
825             }
826 0           return PDL_err;
827             }
828              
829             static pdl_datatypes pdl_inner2t_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, PDL_CF, PDL_CD, PDL_CLD, -1 };
830             static PDL_Indx pdl_inner2t_vtable_realdims[] = { 2, 2, 2, 2, 2 };
831             static char *pdl_inner2t_vtable_parnames[] = { "a","b","c","tmp","d" };
832             static short pdl_inner2t_vtable_parflags[] = {
833             0,
834             0,
835             0,
836             PDL_PARAM_ISCREAT|PDL_PARAM_ISTEMP|PDL_PARAM_ISWRITE,
837             PDL_PARAM_ISCREAT|PDL_PARAM_ISOUT|PDL_PARAM_ISWRITE
838             };
839             static pdl_datatypes pdl_inner2t_vtable_partypes[] = { -1, -1, -1, -1, -1 };
840             static PDL_Indx pdl_inner2t_vtable_realdims_starts[] = { 0, 2, 4, 6, 8 };
841             static PDL_Indx pdl_inner2t_vtable_realdims_ind_ids[] = { 0, 3, 3, 2, 2, 1, 3, 1, 0, 1 };
842             static char *pdl_inner2t_vtable_indnames[] = { "j","k","m","n" };
843             pdl_transvtable pdl_inner2t_vtable = {
844             PDL_TRANS_DO_BROADCAST|PDL_TRANS_BADPROCESS, 0, pdl_inner2t_vtable_gentypes, 3, 5, NULL /*CORE21*/,
845             pdl_inner2t_vtable_realdims, pdl_inner2t_vtable_parnames,
846             pdl_inner2t_vtable_parflags, pdl_inner2t_vtable_partypes,
847             pdl_inner2t_vtable_realdims_starts, pdl_inner2t_vtable_realdims_ind_ids, 10,
848             4, pdl_inner2t_vtable_indnames,
849             NULL, pdl_inner2t_readdata, NULL,
850             NULL,
851             0,"PDL::Primitive::inner2t"
852             };
853              
854              
855 0           pdl_error pdl_run_inner2t(pdl *a,pdl *b,pdl *c,pdl *d) {
856 0           pdl_error PDL_err = {0, NULL, 0};
857 0 0         if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
858 0           pdl_trans *__privtrans = PDL->create_trans(&pdl_inner2t_vtable);
859 0 0         if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
860 0           __privtrans->pdls[0] = a;
861 0           __privtrans->pdls[1] = b;
862 0           __privtrans->pdls[2] = c;
863 0           __privtrans->pdls[4] = d;
864 0 0         PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
865 0 0         PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));
866 0           return PDL_err;
867             }