- node = &fNodes[tnode/2 + (tnode%2) - 1];
- if((TMath::Abs(point[node->fAxis] - node->fValue) > fkNNdist[k-1])
- &&
- ((point[node->fAxis] > node->fValue && tnode%2) ||
- (point[node->fAxis] < node->fValue && !tnode%2))) {
- //printf("\tREMOVE NODE %d\n", tnode/2 + (tnode%2) - 1);
+ Int_t inode = (tnode-1)>>1; //tnode/2 + (tnode%2) - 1;
+ ax = fAxis[inode];
+ val = fValue[inode];
+ dif = point[ax] - val;
+ if((TMath::Abs(dif) > fkNNdist[k-1]) &&
+ ((dif > 0. && (tnode&1)) || (dif < 0. && !(tnode&1)))) {
+ if(debug>=1) printf("\tremove %d\n", (tnode-1)>>1);