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