Branch Coverage

Transform.xs
Criterion Covered Total %
branch 368 1766 20.8


line true false branch
53 737701 0 if (slimit < 6.0)
59 1475402 737701 for (i=0; i
60 2950804 1475402 for (j=0; j
66 3735371 360821 while (RotCount != 0 && SweepCount <= slimit)
3358491 376880 while (RotCount != 0 && SweepCount <= slimit)
70 3358491 3358491 for (j=0; j
72 3358491 3358491 for (k=j+1; k
75 6716982 3358491 for (i=0; i
81 526731 2831760 if (q >= r)
83 526731 0 if (q<=e2*Z[0] || fabs(p)<=tol*q) RotCount--;
360821 165910 if (q<=e2*Z[0] || fabs(p)<=tol*q) RotCount--;
88 663640 165910 for (i=0; i
99 99020 2732740 if (p<0) s0 = -s0;
101 11327040 2831760 for (i=0; i
109 0 3358491 while (EstColRank>=3 && Z[(EstColRank-1)]<=Z[0]*tol+tol*tol)
0 0 while (EstColRank>=3 && Z[(EstColRank-1)]<=Z[0]*tol+tol*tol)
149 1475402 737701 for( i=offset=0; i
154 1475402 737701 for( i=0; i
157 1471456 3946 char symmetric = !(bot || top);
1467510 3946 char symmetric = !(bot || top);
160 1471456 3946 ohi = ((PDLA_Double *)map->data) + ( offset + ( top ? 0 : diff ));
161 1471456 3946 olo = ((PDLA_Double *)map->data) + ( offset - ( bot ? 0 : diff ));
163 2950804 1475402 for( j=0; j
169 2935020 15784 if(symmetric)
186 1475402 737701 for (i=0;i
192 1475402 737701 for(i=0; i
194 2950804 1475402 for(j=0;j
201 1475402 737701 for(i=0;i
203 720800 754602 if(*aptr < sv_min)
205 360821 1114581 if(*aptr > max_sv )
218 1475402 737701 for(i=0;i
219 2950804 1475402 for(j=0;j
222 5901608 2950804 for(k=0;k
253 0 0 PDLA_Byte * k0_datap = ((PDLA_Byte *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Byte * k0_datap = ((PDLA_Byte *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
259 0 0 PDLA_Short * k0_datap = ((PDLA_Short *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Short * k0_datap = ((PDLA_Short *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
265 0 0 PDLA_Ushort * k0_datap = ((PDLA_Ushort *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Ushort * k0_datap = ((PDLA_Ushort *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
271 0 5 PDLA_Long * k0_datap = ((PDLA_Long *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Long * k0_datap = ((PDLA_Long *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
277 0 0 PDLA_Indx * k0_datap = ((PDLA_Indx *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Indx * k0_datap = ((PDLA_Indx *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
283 0 0 PDLA_LongLong * k0_datap = ((PDLA_LongLong *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_LongLong * k0_datap = ((PDLA_LongLong *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
289 0 0 PDLA_Float * k0_datap = ((PDLA_Float *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Float * k0_datap = ((PDLA_Float *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
295 0 7 PDLA_Double * k0_datap = ((PDLA_Double *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Double * k0_datap = ((PDLA_Double *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
323 12 0 if(!hdrp &&
0 12 if(!hdrp &&
324 0 0 __privtrans->pdls[0]->hdrsv &&
330 0 12 if (hdrp) {
331 0 0 if(hdrp == &PL_sv_undef)
339 0 0 PUSHMARK(SP) ;
340 0 0 XPUSHs( hdrp );
344 0 0 if(count != 1)
349 0 0 if(hdr_copy && hdr_copy != &PL_sv_undef) {
0 0 if(hdr_copy && hdr_copy != &PL_sv_undef) {
353 0 0 FREETMPS ;
360 0 0 if(hdr_copy != &PL_sv_undef)
384 0 0 for(i=0; i<__copy->vtable->npdls; i++)
388 0 0 if(__copy->__ddone) {
405 0 0 PDLA_Byte * k0_datap = ((PDLA_Byte *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Byte * k0_datap = ((PDLA_Byte *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
410 0 0 if ( PDLA->startthreadloop(&(__privtrans->__pdlthread),__privtrans->vtable->readdata, __tr) ) return;
419 0 0 for( __tind2 = 0 ;
425 0 0 for( __tind1 = 0 ;
464 0 0 PDLA_Byte badval = SvNV(__privtrans->bv);
509 0 0 bstr = SvPV(__privtrans->boundary,blen);
511 0 0 if(blen == 0) {
514 0 0 for (i=0;i
518 0 0 for(i=0;i
519 0 0 switch(bstr[i < blen ? i : blen-1 ]) {
550 0 0 big = labs((PDLA_Long) (SvNV(__privtrans->big)));
551 0 0 if(big <= 0)
554 0 0 blur = fabs((PDLA_Double) (SvNV(__privtrans->blur)));
555 0 0 if(blur < 0)
558 0 0 sv_min = fabs((PDLA_Double) (SvNV(__privtrans->sv_min)));
559 0 0 if(sv_min < 0)
562 0 0 flux = (SvNV(__privtrans->flux) != 0);
567 0 0 mstr = SvPV(__privtrans->method,mlen);
569 0 0 if(mlen==0)
574 0 0 if( needs_hanning_calc ) {
576 0 0 for(i=0;i
585 0 0 0 :
592 0 0 if( needs_gaussian_calc ) {
594 0 0 for(i=0;i
622 0 0 for(i=0;i
651 0 0 if(psize <= big) {
657 0 0 for (i=0;i
664 0 0 for(i=0;i
670 0 0 for(i=0; i < in->dims[ndims]; i++)
676 0 0 for(i=0;i < in->dims[ndims]; i++)
681 0 0 for(i=0;i < in->dims[ndims]; i++)
698 0 0 for(i=0;i
700 0 0 if(j<0 || j >= in->dims[i]) {
0 0 if(j<0 || j >= in->dims[i]) {
709 0 0 if(j<0)
715 0 0 if(j<0)
721 0 0 if(j<0)
724 0 0 if(j<0) {
742 0 0 for(i=0;i
765 0 0 if( !t_vio ) {
770 0 0 for(i=0; i
784 0 0 for(i=0; i
791 0 0 for(j=0;j
801 0 0 if(beta > 0) {
802 0 0 if(beta >= blur) {
808 0 0 beta -= lodex; if(lodex > HANNING_LOOKUP_SIZE)
822 0 0 for(i=0; i
825 0 0 for(j=0;j
828 0 0 if( dd > 1 ) {
841 0 0 for(i=0; i
844 0 0 for(j=0;j
848 0 0 if(sum > GAUSSIAN_MAXVAL) {
853 0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
873 0 0 for(i=0; i
876 0 0 for(j=0;j
880 0 0 if(sum > 4) /* 2 pixels -- four half-widths */
884 0 0 if(sum > GAUSSIAN_MAXVAL)
903 0 0 for( i=0; i < max; i++ ) {
904 0 0 if( (badval==0) || (*dat != badval) ) {
0 0 if( (badval==0) || (*dat != badval) ) {
917 0 0 for(i=0; i
0 0 for(i=0; i
923 0 0 if( j > 0 && j <= in->dims[i]-1 ) {
0 0 if( j > 0 && j <= in->dims[i]-1 ) {
931 0 0 if( j == 0 )
933 0 0 else if( j == in->dims[i] )
939 0 0 if((j % in->dims[i]) == 0) {
949 0 0 if( j!=0 && j!= -in->dims[i] ) {
0 0 if( j!=0 && j!= -in->dims[i] ) {
950 0 0 if(j<0)
960 0 0 if(ivec[i] <= psize) {
963 0 0 for(k=0;k
970 0 0 if(bounds[i]==1) {
972 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
976 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
984 0 0 } while(carry==0); /* end of total data accumulation loop (termination condition has carry on last dim) */
993 0 0 for(i=0;i
996 0 0 if(!flux) {
998 0 0 for(i=0; i < out->dims[ndims]; i++) {
999 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
1011 0 0 for(i=0; i < out->dims[ndims]; i++) {
1012 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
1028 0 0 for(i=0;i
1030 0 0 for(i=0;idims[ndims];i++) {
1038 0 0 for(i=0;
1039 0 0 (i
1040 0 0 (map_ptr += map->dimincs[i+1]) && /* Funky pre-test increment */
1047 0 0 } while(i
1056 0 0 } while(PDLA->iterthreadloop(&__privtrans->__pdlthread,2)); } break; case PDLA_S: {
1057 0 0 PDLA_Short * k0_datap = ((PDLA_Short *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Short * k0_datap = ((PDLA_Short *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
1062 0 0 if ( PDLA->startthreadloop(&(__privtrans->__pdlthread),__privtrans->vtable->readdata, __tr) ) return;
1071 0 0 for( __tind2 = 0 ;
1077 0 0 for( __tind1 = 0 ;
1116 0 0 PDLA_Short badval = SvNV(__privtrans->bv);
1161 0 0 bstr = SvPV(__privtrans->boundary,blen);
1163 0 0 if(blen == 0) {
1166 0 0 for (i=0;i
1170 0 0 for(i=0;i
1171 0 0 switch(bstr[i < blen ? i : blen-1 ]) {
1202 0 0 big = labs((PDLA_Long) (SvNV(__privtrans->big)));
1203 0 0 if(big <= 0)
1206 0 0 blur = fabs((PDLA_Double) (SvNV(__privtrans->blur)));
1207 0 0 if(blur < 0)
1210 0 0 sv_min = fabs((PDLA_Double) (SvNV(__privtrans->sv_min)));
1211 0 0 if(sv_min < 0)
1214 0 0 flux = (SvNV(__privtrans->flux) != 0);
1219 0 0 mstr = SvPV(__privtrans->method,mlen);
1221 0 0 if(mlen==0)
1226 0 0 if( needs_hanning_calc ) {
1228 0 0 for(i=0;i
1237 0 0 0 :
1244 0 0 if( needs_gaussian_calc ) {
1246 0 0 for(i=0;i
1274 0 0 for(i=0;i
1303 0 0 if(psize <= big) {
1309 0 0 for (i=0;i
1316 0 0 for(i=0;i
1322 0 0 for(i=0; i < in->dims[ndims]; i++)
1328 0 0 for(i=0;i < in->dims[ndims]; i++)
1333 0 0 for(i=0;i < in->dims[ndims]; i++)
1350 0 0 for(i=0;i
1352 0 0 if(j<0 || j >= in->dims[i]) {
0 0 if(j<0 || j >= in->dims[i]) {
1361 0 0 if(j<0)
1367 0 0 if(j<0)
1373 0 0 if(j<0)
1376 0 0 if(j<0) {
1394 0 0 for(i=0;i
1417 0 0 if( !t_vio ) {
1422 0 0 for(i=0; i
1436 0 0 for(i=0; i
1443 0 0 for(j=0;j
1453 0 0 if(beta > 0) {
1454 0 0 if(beta >= blur) {
1460 0 0 beta -= lodex; if(lodex > HANNING_LOOKUP_SIZE)
1474 0 0 for(i=0; i
1477 0 0 for(j=0;j
1480 0 0 if( dd > 1 ) {
1493 0 0 for(i=0; i
1496 0 0 for(j=0;j
1500 0 0 if(sum > GAUSSIAN_MAXVAL) {
1505 0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
1525 0 0 for(i=0; i
1528 0 0 for(j=0;j
1532 0 0 if(sum > 4) /* 2 pixels -- four half-widths */
1536 0 0 if(sum > GAUSSIAN_MAXVAL)
1555 0 0 for( i=0; i < max; i++ ) {
1556 0 0 if( (badval==0) || (*dat != badval) ) {
0 0 if( (badval==0) || (*dat != badval) ) {
1569 0 0 for(i=0; i
0 0 for(i=0; i
1575 0 0 if( j > 0 && j <= in->dims[i]-1 ) {
0 0 if( j > 0 && j <= in->dims[i]-1 ) {
1583 0 0 if( j == 0 )
1585 0 0 else if( j == in->dims[i] )
1591 0 0 if((j % in->dims[i]) == 0) {
1601 0 0 if( j!=0 && j!= -in->dims[i] ) {
0 0 if( j!=0 && j!= -in->dims[i] ) {
1602 0 0 if(j<0)
1612 0 0 if(ivec[i] <= psize) {
1615 0 0 for(k=0;k
1622 0 0 if(bounds[i]==1) {
1624 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
1628 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
1636 0 0 } while(carry==0); /* end of total data accumulation loop (termination condition has carry on last dim) */
1645 0 0 for(i=0;i
1648 0 0 if(!flux) {
1650 0 0 for(i=0; i < out->dims[ndims]; i++) {
1651 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
1663 0 0 for(i=0; i < out->dims[ndims]; i++) {
1664 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
1680 0 0 for(i=0;i
1682 0 0 for(i=0;idims[ndims];i++) {
1690 0 0 for(i=0;
1691 0 0 (i
1692 0 0 (map_ptr += map->dimincs[i+1]) && /* Funky pre-test increment */
1699 0 0 } while(i
1708 0 0 } while(PDLA->iterthreadloop(&__privtrans->__pdlthread,2)); } break; case PDLA_US: {
1709 0 0 PDLA_Ushort * k0_datap = ((PDLA_Ushort *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Ushort * k0_datap = ((PDLA_Ushort *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
1714 0 0 if ( PDLA->startthreadloop(&(__privtrans->__pdlthread),__privtrans->vtable->readdata, __tr) ) return;
1723 0 0 for( __tind2 = 0 ;
1729 0 0 for( __tind1 = 0 ;
1768 0 0 PDLA_Ushort badval = SvNV(__privtrans->bv);
1813 0 0 bstr = SvPV(__privtrans->boundary,blen);
1815 0 0 if(blen == 0) {
1818 0 0 for (i=0;i
1822 0 0 for(i=0;i
1823 0 0 switch(bstr[i < blen ? i : blen-1 ]) {
1854 0 0 big = labs((PDLA_Long) (SvNV(__privtrans->big)));
1855 0 0 if(big <= 0)
1858 0 0 blur = fabs((PDLA_Double) (SvNV(__privtrans->blur)));
1859 0 0 if(blur < 0)
1862 0 0 sv_min = fabs((PDLA_Double) (SvNV(__privtrans->sv_min)));
1863 0 0 if(sv_min < 0)
1866 0 0 flux = (SvNV(__privtrans->flux) != 0);
1871 0 0 mstr = SvPV(__privtrans->method,mlen);
1873 0 0 if(mlen==0)
1878 0 0 if( needs_hanning_calc ) {
1880 0 0 for(i=0;i
1889 0 0 0 :
1896 0 0 if( needs_gaussian_calc ) {
1898 0 0 for(i=0;i
1926 0 0 for(i=0;i
1955 0 0 if(psize <= big) {
1961 0 0 for (i=0;i
1968 0 0 for(i=0;i
1974 0 0 for(i=0; i < in->dims[ndims]; i++)
1980 0 0 for(i=0;i < in->dims[ndims]; i++)
1985 0 0 for(i=0;i < in->dims[ndims]; i++)
2002 0 0 for(i=0;i
2004 0 0 if(j<0 || j >= in->dims[i]) {
0 0 if(j<0 || j >= in->dims[i]) {
2013 0 0 if(j<0)
2019 0 0 if(j<0)
2025 0 0 if(j<0)
2028 0 0 if(j<0) {
2046 0 0 for(i=0;i
2069 0 0 if( !t_vio ) {
2074 0 0 for(i=0; i
2088 0 0 for(i=0; i
2095 0 0 for(j=0;j
2105 0 0 if(beta > 0) {
2106 0 0 if(beta >= blur) {
2112 0 0 beta -= lodex; if(lodex > HANNING_LOOKUP_SIZE)
2126 0 0 for(i=0; i
2129 0 0 for(j=0;j
2132 0 0 if( dd > 1 ) {
2145 0 0 for(i=0; i
2148 0 0 for(j=0;j
2152 0 0 if(sum > GAUSSIAN_MAXVAL) {
2157 0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
2177 0 0 for(i=0; i
2180 0 0 for(j=0;j
2184 0 0 if(sum > 4) /* 2 pixels -- four half-widths */
2188 0 0 if(sum > GAUSSIAN_MAXVAL)
2207 0 0 for( i=0; i < max; i++ ) {
2208 0 0 if( (badval==0) || (*dat != badval) ) {
0 0 if( (badval==0) || (*dat != badval) ) {
2221 0 0 for(i=0; i
0 0 for(i=0; i
2227 0 0 if( j > 0 && j <= in->dims[i]-1 ) {
0 0 if( j > 0 && j <= in->dims[i]-1 ) {
2235 0 0 if( j == 0 )
2237 0 0 else if( j == in->dims[i] )
2243 0 0 if((j % in->dims[i]) == 0) {
2253 0 0 if( j!=0 && j!= -in->dims[i] ) {
0 0 if( j!=0 && j!= -in->dims[i] ) {
2254 0 0 if(j<0)
2264 0 0 if(ivec[i] <= psize) {
2267 0 0 for(k=0;k
2274 0 0 if(bounds[i]==1) {
2276 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
2280 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
2288 0 0 } while(carry==0); /* end of total data accumulation loop (termination condition has carry on last dim) */
2297 0 0 for(i=0;i
2300 0 0 if(!flux) {
2302 0 0 for(i=0; i < out->dims[ndims]; i++) {
2303 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
2315 0 0 for(i=0; i < out->dims[ndims]; i++) {
2316 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
2332 0 0 for(i=0;i
2334 0 0 for(i=0;idims[ndims];i++) {
2342 0 0 for(i=0;
2343 0 0 (i
2344 0 0 (map_ptr += map->dimincs[i+1]) && /* Funky pre-test increment */
2351 0 0 } while(i
2360 0 0 } while(PDLA->iterthreadloop(&__privtrans->__pdlthread,2)); } break; case PDLA_L: {
2361 0 5 PDLA_Long * k0_datap = ((PDLA_Long *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Long * k0_datap = ((PDLA_Long *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
2366 0 5 if ( PDLA->startthreadloop(&(__privtrans->__pdlthread),__privtrans->vtable->readdata, __tr) ) return;
2375 5 5 for( __tind2 = 0 ;
2381 5 5 for( __tind1 = 0 ;
2420 0 5 PDLA_Long badval = SvNV(__privtrans->bv);
2465 5 0 bstr = SvPV(__privtrans->boundary,blen);
2467 0 5 if(blen == 0) {
2470 0 0 for (i=0;i
2474 10 5 for(i=0;i
2475 5 5 switch(bstr[i < blen ? i : blen-1 ]) {
2506 0 5 big = labs((PDLA_Long) (SvNV(__privtrans->big)));
2507 0 5 if(big <= 0)
2510 5 0 blur = fabs((PDLA_Double) (SvNV(__privtrans->blur)));
2511 0 5 if(blur < 0)
2514 5 0 sv_min = fabs((PDLA_Double) (SvNV(__privtrans->sv_min)));
2515 0 5 if(sv_min < 0)
2518 5 0 flux = (SvNV(__privtrans->flux) != 0);
2523 5 0 mstr = SvPV(__privtrans->method,mlen);
2525 0 5 if(mlen==0)
2530 1 0 if( needs_hanning_calc ) {
2532 2500 1 for(i=0;i
2541 1 0 0 :
2548 1 1 if( needs_gaussian_calc ) {
2550 4000 1 for(i=0;i
2578 10 5 for(i=0;i
2607 737280 0 if(psize <= big) {
2613 1474560 737280 for (i=0;i
2620 1474560 737280 for(i=0;i
2626 737280 737280 for(i=0; i < in->dims[ndims]; i++)
2632 737280 737280 for(i=0;i < in->dims[ndims]; i++)
2637 737280 737280 for(i=0;i < in->dims[ndims]; i++)
2654 1474560 737280 for(i=0;i
2656 356280 1118280 if(j<0 || j >= in->dims[i]) {
0 356280 if(j<0 || j >= in->dims[i]) {
2665 1118280 0 if(j<0)
2671 0 0 if(j<0)
2677 0 0 if(j<0)
2680 0 0 if(j<0) {
2698 1474560 737280 for(i=0;i
2721 3043295 9358625 if( !t_vio ) {
2726 6086590 3043295 for(i=0; i
2740 873347 608659 for(i=0; i
2747 1746694 873347 for(j=0;j
2757 865595 7752 if(beta > 0) {
2758 511183 354412 if(beta >= blur) {
2764 0 354412 beta -= lodex; if(lodex > HANNING_LOOKUP_SIZE)
2778 873355 608659 for(i=0; i
2781 1746710 873355 for(j=0;j
2784 511175 362180 if( dd > 1 ) {
2797 2360082 1217318 for(i=0; i
2800 4720164 2360082 for(j=0;j
2804 308982 2051100 if(sum > GAUSSIAN_MAXVAL) {
2809 908336 308982 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
900584 7752 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 900584 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
2829 1100068 608659 for(i=0; i
2832 2200136 1100068 for(j=0;j
2836 316983 783085 if(sum > 4) /* 2 pixels -- four half-widths */
2840 123548 485111 if(sum > GAUSSIAN_MAXVAL)
2859 3043295 3043295 for( i=0; i < max; i++ ) {
2860 0 3043295 if( (badval==0) || (*dat != badval) ) {
0 0 if( (badval==0) || (*dat != badval) ) {
2873 24803840 2932640 for(i=0; i
15334560 9469280 for(i=0; i
2879 6289185 9045375 if( j > 0 && j <= in->dims[i]-1 ) {
6289185 0 if( j > 0 && j <= in->dims[i]-1 ) {
2887 521970 8523405 if( j == 0 )
2889 0 8523405 else if( j == in->dims[i] )
2895 0 0 if((j % in->dims[i]) == 0) {
2905 0 0 if( j!=0 && j!= -in->dims[i] ) {
0 0 if( j!=0 && j!= -in->dims[i] ) {
2906 0 0 if(j<0)
2916 11664640 3669920 if(ivec[i] <= psize) {
2919 2195360 11664640 for(k=0;k
2926 3669920 0 if(bounds[i]==1) {
2928 1531800 2138120 if( j < 0 || j >= in->dims[i] )
0 1531800 if( j < 0 || j >= in->dims[i] )
2932 1009830 2660090 if( j < 0 || j >= in->dims[i] )
0 1009830 if( j < 0 || j >= in->dims[i] )
2940 11664640 737280 } while(carry==0); /* end of total data accumulation loop (termination condition has carry on last dim) */
2949 1474560 737280 for(i=0;i
2952 737280 0 if(!flux) {
2954 737280 737280 for(i=0; i < out->dims[ndims]; i++) {
2955 139110 598170 if(*wg && (*wg2 / *wg) < 1.5 ) {
138062 1048 if(*wg && (*wg2 / *wg) < 1.5 ) {
2967 0 0 for(i=0; i < out->dims[ndims]; i++) {
2968 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
2984 0 0 for(i=0;i
2986 0 0 for(i=0;idims[ndims];i++) {
2994 739200 5 for(i=0;
2995 739200 0 (i
2996 1925 737275 (map_ptr += map->dimincs[i+1]) && /* Funky pre-test increment */
3003 737275 5 } while(i
3012 0 5 } while(PDLA->iterthreadloop(&__privtrans->__pdlthread,2)); } break; case PDLA_IND: {
3013 0 0 PDLA_Indx * k0_datap = ((PDLA_Indx *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Indx * k0_datap = ((PDLA_Indx *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
3018 0 0 if ( PDLA->startthreadloop(&(__privtrans->__pdlthread),__privtrans->vtable->readdata, __tr) ) return;
3027 0 0 for( __tind2 = 0 ;
3033 0 0 for( __tind1 = 0 ;
3072 0 0 PDLA_Indx badval = SvNV(__privtrans->bv);
3117 0 0 bstr = SvPV(__privtrans->boundary,blen);
3119 0 0 if(blen == 0) {
3122 0 0 for (i=0;i
3126 0 0 for(i=0;i
3127 0 0 switch(bstr[i < blen ? i : blen-1 ]) {
3158 0 0 big = labs((PDLA_Long) (SvNV(__privtrans->big)));
3159 0 0 if(big <= 0)
3162 0 0 blur = fabs((PDLA_Double) (SvNV(__privtrans->blur)));
3163 0 0 if(blur < 0)
3166 0 0 sv_min = fabs((PDLA_Double) (SvNV(__privtrans->sv_min)));
3167 0 0 if(sv_min < 0)
3170 0 0 flux = (SvNV(__privtrans->flux) != 0);
3175 0 0 mstr = SvPV(__privtrans->method,mlen);
3177 0 0 if(mlen==0)
3182 0 0 if( needs_hanning_calc ) {
3184 0 0 for(i=0;i
3193 0 0 0 :
3200 0 0 if( needs_gaussian_calc ) {
3202 0 0 for(i=0;i
3230 0 0 for(i=0;i
3259 0 0 if(psize <= big) {
3265 0 0 for (i=0;i
3272 0 0 for(i=0;i
3278 0 0 for(i=0; i < in->dims[ndims]; i++)
3284 0 0 for(i=0;i < in->dims[ndims]; i++)
3289 0 0 for(i=0;i < in->dims[ndims]; i++)
3306 0 0 for(i=0;i
3308 0 0 if(j<0 || j >= in->dims[i]) {
0 0 if(j<0 || j >= in->dims[i]) {
3317 0 0 if(j<0)
3323 0 0 if(j<0)
3329 0 0 if(j<0)
3332 0 0 if(j<0) {
3350 0 0 for(i=0;i
3373 0 0 if( !t_vio ) {
3378 0 0 for(i=0; i
3392 0 0 for(i=0; i
3399 0 0 for(j=0;j
3409 0 0 if(beta > 0) {
3410 0 0 if(beta >= blur) {
3416 0 0 beta -= lodex; if(lodex > HANNING_LOOKUP_SIZE)
3430 0 0 for(i=0; i
3433 0 0 for(j=0;j
3436 0 0 if( dd > 1 ) {
3449 0 0 for(i=0; i
3452 0 0 for(j=0;j
3456 0 0 if(sum > GAUSSIAN_MAXVAL) {
3461 0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
3481 0 0 for(i=0; i
3484 0 0 for(j=0;j
3488 0 0 if(sum > 4) /* 2 pixels -- four half-widths */
3492 0 0 if(sum > GAUSSIAN_MAXVAL)
3511 0 0 for( i=0; i < max; i++ ) {
3512 0 0 if( (badval==0) || (*dat != badval) ) {
0 0 if( (badval==0) || (*dat != badval) ) {
3525 0 0 for(i=0; i
0 0 for(i=0; i
3531 0 0 if( j > 0 && j <= in->dims[i]-1 ) {
0 0 if( j > 0 && j <= in->dims[i]-1 ) {
3539 0 0 if( j == 0 )
3541 0 0 else if( j == in->dims[i] )
3547 0 0 if((j % in->dims[i]) == 0) {
3557 0 0 if( j!=0 && j!= -in->dims[i] ) {
0 0 if( j!=0 && j!= -in->dims[i] ) {
3558 0 0 if(j<0)
3568 0 0 if(ivec[i] <= psize) {
3571 0 0 for(k=0;k
3578 0 0 if(bounds[i]==1) {
3580 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
3584 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
3592 0 0 } while(carry==0); /* end of total data accumulation loop (termination condition has carry on last dim) */
3601 0 0 for(i=0;i
3604 0 0 if(!flux) {
3606 0 0 for(i=0; i < out->dims[ndims]; i++) {
3607 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
3619 0 0 for(i=0; i < out->dims[ndims]; i++) {
3620 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
3636 0 0 for(i=0;i
3638 0 0 for(i=0;idims[ndims];i++) {
3646 0 0 for(i=0;
3647 0 0 (i
3648 0 0 (map_ptr += map->dimincs[i+1]) && /* Funky pre-test increment */
3655 0 0 } while(i
3664 0 0 } while(PDLA->iterthreadloop(&__privtrans->__pdlthread,2)); } break; case PDLA_LL: {
3665 0 0 PDLA_LongLong * k0_datap = ((PDLA_LongLong *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_LongLong * k0_datap = ((PDLA_LongLong *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
3670 0 0 if ( PDLA->startthreadloop(&(__privtrans->__pdlthread),__privtrans->vtable->readdata, __tr) ) return;
3679 0 0 for( __tind2 = 0 ;
3685 0 0 for( __tind1 = 0 ;
3724 0 0 PDLA_LongLong badval = SvNV(__privtrans->bv);
3769 0 0 bstr = SvPV(__privtrans->boundary,blen);
3771 0 0 if(blen == 0) {
3774 0 0 for (i=0;i
3778 0 0 for(i=0;i
3779 0 0 switch(bstr[i < blen ? i : blen-1 ]) {
3810 0 0 big = labs((PDLA_Long) (SvNV(__privtrans->big)));
3811 0 0 if(big <= 0)
3814 0 0 blur = fabs((PDLA_Double) (SvNV(__privtrans->blur)));
3815 0 0 if(blur < 0)
3818 0 0 sv_min = fabs((PDLA_Double) (SvNV(__privtrans->sv_min)));
3819 0 0 if(sv_min < 0)
3822 0 0 flux = (SvNV(__privtrans->flux) != 0);
3827 0 0 mstr = SvPV(__privtrans->method,mlen);
3829 0 0 if(mlen==0)
3834 0 0 if( needs_hanning_calc ) {
3836 0 0 for(i=0;i
3845 0 0 0 :
3852 0 0 if( needs_gaussian_calc ) {
3854 0 0 for(i=0;i
3882 0 0 for(i=0;i
3911 0 0 if(psize <= big) {
3917 0 0 for (i=0;i
3924 0 0 for(i=0;i
3930 0 0 for(i=0; i < in->dims[ndims]; i++)
3936 0 0 for(i=0;i < in->dims[ndims]; i++)
3941 0 0 for(i=0;i < in->dims[ndims]; i++)
3958 0 0 for(i=0;i
3960 0 0 if(j<0 || j >= in->dims[i]) {
0 0 if(j<0 || j >= in->dims[i]) {
3969 0 0 if(j<0)
3975 0 0 if(j<0)
3981 0 0 if(j<0)
3984 0 0 if(j<0) {
4002 0 0 for(i=0;i
4025 0 0 if( !t_vio ) {
4030 0 0 for(i=0; i
4044 0 0 for(i=0; i
4051 0 0 for(j=0;j
4061 0 0 if(beta > 0) {
4062 0 0 if(beta >= blur) {
4068 0 0 beta -= lodex; if(lodex > HANNING_LOOKUP_SIZE)
4082 0 0 for(i=0; i
4085 0 0 for(j=0;j
4088 0 0 if( dd > 1 ) {
4101 0 0 for(i=0; i
4104 0 0 for(j=0;j
4108 0 0 if(sum > GAUSSIAN_MAXVAL) {
4113 0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
4133 0 0 for(i=0; i
4136 0 0 for(j=0;j
4140 0 0 if(sum > 4) /* 2 pixels -- four half-widths */
4144 0 0 if(sum > GAUSSIAN_MAXVAL)
4163 0 0 for( i=0; i < max; i++ ) {
4164 0 0 if( (badval==0) || (*dat != badval) ) {
0 0 if( (badval==0) || (*dat != badval) ) {
4177 0 0 for(i=0; i
0 0 for(i=0; i
4183 0 0 if( j > 0 && j <= in->dims[i]-1 ) {
0 0 if( j > 0 && j <= in->dims[i]-1 ) {
4191 0 0 if( j == 0 )
4193 0 0 else if( j == in->dims[i] )
4199 0 0 if((j % in->dims[i]) == 0) {
4209 0 0 if( j!=0 && j!= -in->dims[i] ) {
0 0 if( j!=0 && j!= -in->dims[i] ) {
4210 0 0 if(j<0)
4220 0 0 if(ivec[i] <= psize) {
4223 0 0 for(k=0;k
4230 0 0 if(bounds[i]==1) {
4232 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
4236 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
4244 0 0 } while(carry==0); /* end of total data accumulation loop (termination condition has carry on last dim) */
4253 0 0 for(i=0;i
4256 0 0 if(!flux) {
4258 0 0 for(i=0; i < out->dims[ndims]; i++) {
4259 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
4271 0 0 for(i=0; i < out->dims[ndims]; i++) {
4272 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
4288 0 0 for(i=0;i
4290 0 0 for(i=0;idims[ndims];i++) {
4298 0 0 for(i=0;
4299 0 0 (i
4300 0 0 (map_ptr += map->dimincs[i+1]) && /* Funky pre-test increment */
4307 0 0 } while(i
4316 0 0 } while(PDLA->iterthreadloop(&__privtrans->__pdlthread,2)); } break; case PDLA_F: {
4317 0 0 PDLA_Float * k0_datap = ((PDLA_Float *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Float * k0_datap = ((PDLA_Float *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
4322 0 0 if ( PDLA->startthreadloop(&(__privtrans->__pdlthread),__privtrans->vtable->readdata, __tr) ) return;
4331 0 0 for( __tind2 = 0 ;
4337 0 0 for( __tind1 = 0 ;
4376 0 0 PDLA_Float badval = SvNV(__privtrans->bv);
4421 0 0 bstr = SvPV(__privtrans->boundary,blen);
4423 0 0 if(blen == 0) {
4426 0 0 for (i=0;i
4430 0 0 for(i=0;i
4431 0 0 switch(bstr[i < blen ? i : blen-1 ]) {
4462 0 0 big = labs((PDLA_Long) (SvNV(__privtrans->big)));
4463 0 0 if(big <= 0)
4466 0 0 blur = fabs((PDLA_Double) (SvNV(__privtrans->blur)));
4467 0 0 if(blur < 0)
4470 0 0 sv_min = fabs((PDLA_Double) (SvNV(__privtrans->sv_min)));
4471 0 0 if(sv_min < 0)
4474 0 0 flux = (SvNV(__privtrans->flux) != 0);
4479 0 0 mstr = SvPV(__privtrans->method,mlen);
4481 0 0 if(mlen==0)
4486 0 0 if( needs_hanning_calc ) {
4488 0 0 for(i=0;i
4497 0 0 0 :
4504 0 0 if( needs_gaussian_calc ) {
4506 0 0 for(i=0;i
4534 0 0 for(i=0;i
4563 0 0 if(psize <= big) {
4569 0 0 for (i=0;i
4576 0 0 for(i=0;i
4582 0 0 for(i=0; i < in->dims[ndims]; i++)
4588 0 0 for(i=0;i < in->dims[ndims]; i++)
4593 0 0 for(i=0;i < in->dims[ndims]; i++)
4610 0 0 for(i=0;i
4612 0 0 if(j<0 || j >= in->dims[i]) {
0 0 if(j<0 || j >= in->dims[i]) {
4621 0 0 if(j<0)
4627 0 0 if(j<0)
4633 0 0 if(j<0)
4636 0 0 if(j<0) {
4654 0 0 for(i=0;i
4677 0 0 if( !t_vio ) {
4682 0 0 for(i=0; i
4696 0 0 for(i=0; i
4703 0 0 for(j=0;j
4713 0 0 if(beta > 0) {
4714 0 0 if(beta >= blur) {
4720 0 0 beta -= lodex; if(lodex > HANNING_LOOKUP_SIZE)
4734 0 0 for(i=0; i
4737 0 0 for(j=0;j
4740 0 0 if( dd > 1 ) {
4753 0 0 for(i=0; i
4756 0 0 for(j=0;j
4760 0 0 if(sum > GAUSSIAN_MAXVAL) {
4765 0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
4785 0 0 for(i=0; i
4788 0 0 for(j=0;j
4792 0 0 if(sum > 4) /* 2 pixels -- four half-widths */
4796 0 0 if(sum > GAUSSIAN_MAXVAL)
4815 0 0 for( i=0; i < max; i++ ) {
4816 0 0 if( (badval==0) || (*dat != badval) ) {
0 0 if( (badval==0) || (*dat != badval) ) {
4829 0 0 for(i=0; i
0 0 for(i=0; i
4835 0 0 if( j > 0 && j <= in->dims[i]-1 ) {
0 0 if( j > 0 && j <= in->dims[i]-1 ) {
4843 0 0 if( j == 0 )
4845 0 0 else if( j == in->dims[i] )
4851 0 0 if((j % in->dims[i]) == 0) {
4861 0 0 if( j!=0 && j!= -in->dims[i] ) {
0 0 if( j!=0 && j!= -in->dims[i] ) {
4862 0 0 if(j<0)
4872 0 0 if(ivec[i] <= psize) {
4875 0 0 for(k=0;k
4882 0 0 if(bounds[i]==1) {
4884 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
4888 0 0 if( j < 0 || j >= in->dims[i] )
0 0 if( j < 0 || j >= in->dims[i] )
4896 0 0 } while(carry==0); /* end of total data accumulation loop (termination condition has carry on last dim) */
4905 0 0 for(i=0;i
4908 0 0 if(!flux) {
4910 0 0 for(i=0; i < out->dims[ndims]; i++) {
4911 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
4923 0 0 for(i=0; i < out->dims[ndims]; i++) {
4924 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
4940 0 0 for(i=0;i
4942 0 0 for(i=0;idims[ndims];i++) {
4950 0 0 for(i=0;
4951 0 0 (i
4952 0 0 (map_ptr += map->dimincs[i+1]) && /* Funky pre-test increment */
4959 0 0 } while(i
4968 0 0 } while(PDLA->iterthreadloop(&__privtrans->__pdlthread,2)); } break; case PDLA_D: {
4969 0 7 PDLA_Double * k0_datap = ((PDLA_Double *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
0 0 PDLA_Double * k0_datap = ((PDLA_Double *)(PDLA_REPRP_TRANS((__privtrans->pdls[0]),(__privtrans->vtable->per_pdl_flags[0]))));
4974 0 7 if ( PDLA->startthreadloop(&(__privtrans->__pdlthread),__privtrans->vtable->readdata, __tr) ) return;
4983 7 7 for( __tind2 = 0 ;
4989 7 7 for( __tind1 = 0 ;
5028 1 6 PDLA_Double badval = SvNV(__privtrans->bv);
5073 7 0 bstr = SvPV(__privtrans->boundary,blen);
5075 0 7 if(blen == 0) {
5078 0 0 for (i=0;i
5082 14 7 for(i=0;i
5083 8 6 switch(bstr[i < blen ? i : blen-1 ]) {
5114 0 7 big = labs((PDLA_Long) (SvNV(__privtrans->big)));
5115 0 7 if(big <= 0)
5118 6 1 blur = fabs((PDLA_Double) (SvNV(__privtrans->blur)));
5119 0 7 if(blur < 0)
5122 7 0 sv_min = fabs((PDLA_Double) (SvNV(__privtrans->sv_min)));
5123 0 7 if(sv_min < 0)
5126 7 0 flux = (SvNV(__privtrans->flux) != 0);
5131 7 0 mstr = SvPV(__privtrans->method,mlen);
5133 0 7 if(mlen==0)
5138 1 5 if( needs_hanning_calc ) {
5140 2500 1 for(i=0;i
5149 6 0 0 :
5156 1 0 if( needs_gaussian_calc ) {
5158 4000 1 for(i=0;i
5186 14 7 for(i=0;i
5215 421 0 if(psize <= big) {
5221 842 421 for (i=0;i
5228 842 421 for(i=0;i
5234 421 421 for(i=0; i < in->dims[ndims]; i++)
5240 421 421 for(i=0;i < in->dims[ndims]; i++)
5245 421 421 for(i=0;i < in->dims[ndims]; i++)
5262 842 421 for(i=0;i
5264 594 248 if(j<0 || j >= in->dims[i]) {
120 474 if(j<0 || j >= in->dims[i]) {
5273 208 60 if(j<0)
5279 20 30 if(j<0)
5285 0 50 if(j<0)
5288 50 0 if(j<0) {
5306 842 421 for(i=0;i
5329 7182 4519 if( !t_vio ) {
5334 14364 7182 for(i=0; i
5348 8353 6341 for(i=0; i
5355 16706 8353 for(j=0;j
5365 6875 1478 if(beta > 0) {
5366 5707 1168 if(beta >= blur) {
5372 0 1168 beta -= lodex; if(lodex > HANNING_LOOKUP_SIZE)
5386 0 0 for(i=0; i
5389 0 0 for(j=0;j
5392 0 0 if( dd > 1 ) {
5405 1682 841 for(i=0; i
5408 3364 1682 for(j=0;j
5412 100 1582 if(sum > GAUSSIAN_MAXVAL) {
5417 741 100 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
741 0 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
0 741 if( sum > GAUSSIAN_MAXVAL || !isfinite(sum) || isnan(sum) ) {
5437 0 0 for(i=0; i
5440 0 0 for(j=0;j
5444 0 0 if(sum > 4) /* 2 pixels -- four half-widths */
5448 0 0 if(sum > GAUSSIAN_MAXVAL)
5467 7182 7182 for( i=0; i < max; i++ ) {
5468 306 6876 if( (badval==0) || (*dat != badval) ) {
306 0 if( (badval==0) || (*dat != badval) ) {
5481 23402 2203 for(i=0; i
13904 9498 for(i=0; i
5487 12626 1278 if( j > 0 && j <= in->dims[i]-1 ) {
7651 4975 if( j > 0 && j <= in->dims[i]-1 ) {
5495 642 3631 if( j == 0 )
5497 950 2681 else if( j == in->dims[i] )
5503 100 230 if((j % in->dims[i]) == 0) {
5513 1400 250 if( j!=0 && j!= -in->dims[i] ) {
1150 250 if( j!=0 && j!= -in->dims[i] ) {
5514 1000 150 if(j<0)
5524 11280 2624 if(ivec[i] <= psize) {
5527 1782 11280 for(k=0;k
5534 2024 600 if(bounds[i]==1) {
5536 2024 0 if( j < 0 || j >= in->dims[i] )
1130 894 if( j < 0 || j >= in->dims[i] )
5540 1382 642 if( j < 0 || j >= in->dims[i] )
180 1202 if( j < 0 || j >= in->dims[i] )
5548 11280 421 } while(carry==0); /* end of total data accumulation loop (termination condition has carry on last dim) */
5557 842 421 for(i=0;i
5560 421 0 if(!flux) {
5562 421 421 for(i=0; i < out->dims[ndims]; i++) {
5563 326 95 if(*wg && (*wg2 / *wg) < 1.5 ) {
326 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
5575 0 0 for(i=0; i < out->dims[ndims]; i++) {
5576 0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
0 0 if(*wg && (*wg2 / *wg) < 1.5 ) {
5592 0 0 for(i=0;i
5594 0 0 for(i=0;idims[ndims];i++) {
5602 474 7 for(i=0;
5603 474 0 (i
5604 60 414 (map_ptr += map->dimincs[i+1]) && /* Funky pre-test increment */
5611 414 7 } while(i
5620 0 7 } while(PDLA->iterthreadloop(&__privtrans->__pdlthread,2)); break;}
5637 12 0 if(__privtrans->__ddone) {
5712 1 11 if (badflag_cache) __privtrans->bvalflag = 1;
5713 12 0 __privtrans->__datatype = 0;if(__privtrans->__datatype < k0->datatype) {
5716 12 0 if(__privtrans->__datatype == PDLA_B) {}
5717 12 0 else if(__privtrans->__datatype == PDLA_S) {}
5718 12 0 else if(__privtrans->__datatype == PDLA_U) {}
5719 7 5 else if(__privtrans->__datatype == PDLA_L) {}
5720 7 0 else if(__privtrans->__datatype == PDLA_N) {}
5721 7 0 else if(__privtrans->__datatype == PDLA_Q) {}
5722 7 0 else if(__privtrans->__datatype == PDLA_F) {}
5723 0 7 else if(__privtrans->__datatype == PDLA_D) {}
5725 0 12 if(__privtrans->__datatype != k0->datatype) {
5744 1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 1 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 1 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
1 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 1 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 1 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
0 0 if (SvTRUE (ERRSV)) Perl_croak(aTHX_ "%s",SvPV_nolen (ERRSV));
5746 0 1 if (CoreSV==NULL)
5748 1 0 PDLA = INT2PTR(Core*,SvIV( CoreSV )); /* Core* value */
5749 0 1 if (PDLA->Version != PDLA_CORE_VERSION)