Method manage_cell

src/lasquadtree.cpp:751–789

create or finalize the cell (in the spatial hierarchy)

// create or finalize the cell (in the spatial hierarchy) 
BOOL LASquadtree::manage_cell(const U32 cell_index, const BOOL finalize)
{
  U32 adaptive_pos = cell_index/32;
  U32 adaptive_bit = ((U32)1) << (cell_index%32);
  if (adaptive_pos >= adaptive_alloc)
  {
    if (adaptive)
    {
      size_t n_pos = (size_t)adaptive_pos * 2;
      adaptive = (U32*)realloc_las(adaptive, n_pos * sizeof(U32));
      for (size_t i = adaptive_alloc; i < n_pos; i++) adaptive[i] = 0;
      adaptive_alloc = n_pos;
    }
    else
    {
#pragma warning(push)
#pragma warning(disable : 6011)
      adaptive = (U32*)malloc_las(((size_t)adaptive_pos + 1) * sizeof(U32));
      for (size_t i = adaptive_alloc; i <= adaptive_pos; i++) adaptive[i] = 0;
      adaptive_alloc = (size_t)adaptive_pos + 1;
#pragma warning(pop)
    }
  }
  adaptive[adaptive_pos] &= ~adaptive_bit;
  U32 index;
  U32 level = get_level(cell_index);
  U32 level_index = get_level_index(cell_index, level);
  while (level)
  {
    level--;
    level_index = level_index >> 2;
    index = get_cell_index(level_index, level);
    adaptive_pos = index/32;
    adaptive_bit = ((U32)1) << (index%32);
    if (adaptive[adaptive_pos] & adaptive_bit) break;
    adaptive[adaptive_pos] |= adaptive_bit;
  }
  return TRUE;
}

// check whether the x & y coordinates fall into the tiling
BOOL LASquadtree::inside(const F64 x, const F64 y) const