NBINS
int x[NBINS + 1], y[NBINS + 1], t[NBINS + 1];
int b1[NBINS + 1], b2[NBINS + 1], xtilde[NBINS + 1];
for (i = 0; i < NBINS; i++) {
rxsum = (lo[i + NBINS] & 0xffff) << 5 |
rosum = (hi[i + NBINS] & 0x7ff) << 16 | hi[i + NBINS] >> 16;