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) |