| 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 |
|
|
| 402 |
|
if (recipavg) |
| 403 |
|
do_reciprocity(); |
| 404 |
|
if (doheader) { |
| 405 |
< |
for (i = 0; i < argc; i++) { |
| 406 |
< |
fputs(argv[i], stdout); |
| 407 |
< |
fputc(i < argc-1 ? ' ' : '\n', stdout); |
| 408 |
< |
} |
| 405 |
> |
newheader("RADIANCE", stdout); |
| 406 |
> |
printargs(argc, argv, stdout); |
| 407 |
> |
fputnow(stdout); |
| 408 |
> |
fputformat("ascii", stdout); |
| 409 |
|
fputc('\n', stdout); |
| 410 |
|
} |
| 411 |
|
gtree.kid = NULL; /* create our tree */ |
