19 |
|
int infmt = 'a'; /* input format ('a','f','d') */ |
20 |
|
double pctcull = 95.; /* target culling percentile */ |
21 |
|
|
22 |
< |
#define dval3(ix,ox,oy) datarr[((((ix)<<log2g)+(ox))<<log2g)+(oy)] |
23 |
< |
#define dval4(ix,iy,ox,oy) datarr[((((((ix)<<log2g)+(iy))<<log2g)+(ox))<<log2g)+(oy)] |
22 |
> |
#define dval3(ix,ox,oy) datarr[((size_t)(((ix)<<log2g)+(ox))<<log2g)+(oy)] |
23 |
> |
#define dval4(ix,iy,ox,oy) datarr[((((size_t)(((ix)<<log2g)+(iy))<<log2g)+(ox))<<log2g)+(oy)] |
24 |
|
|
25 |
|
/* Tensor tree node */ |
26 |
|
typedef struct ttree_s { |
30 |
|
} TNODE; |
31 |
|
|
32 |
|
#define HISTLEN 300 /* histogram resolution */ |
33 |
< |
#define HISTMAX 10. /* maximum recorded measure in histogram */ |
33 |
> |
#define HISTMAX 4. /* maximum recorded measure in histogram */ |
34 |
|
|
35 |
|
int histo[HISTLEN]; /* histogram freq. of variance measure */ |
36 |
|
|
37 |
|
double tthresh; /* acceptance threshold (TBD) */ |
38 |
|
|
39 |
< |
#define var_measure(tp) ( ((tp)->vmax - (tp)->vmin) / \ |
39 |
> |
#define var_measure(tp) sqrt( ((tp)->vmax - (tp)->vmin) / \ |
40 |
|
(sqrt((tp)->vavg) + .03) ) |
41 |
|
#define above_threshold(tp) (var_measure(tp) > tthresh) |
42 |
|
|