Method crop_to_tile_impl

tiledb/sm/subarray/subarray.cc:2743–2773

template <class T, class SubarrayT>
void Subarray::crop_to_tile_impl(const T* tile_coords, SubarrayT& ret) const {
  T new_range[2];
  bool overlaps;

  // Get tile subarray based on the input coordinates
  const auto& array_schema = array_->array_schema_latest();
  std::vector<T> tile_subarray(2 * dim_num());
  array_schema.domain().get_tile_subarray(tile_coords, &tile_subarray[0]);

  // Compute cropped subarray
  for (unsigned d = 0; d < dim_num(); ++d) {
    auto r_size{2 * array_schema.dimension_ptr(d)->coord_size()};
    uint64_t i = 0;
    for (size_t r = 0; r < range_subset_[d].num_ranges(); ++r) {
      const auto& range = range_subset_[d][r];
      rectangle::overlap(
          (const T*)range.data(),
          &tile_subarray[2 * d],
          1,
          new_range,
          &overlaps);

      if (overlaps) {
        ret.add_range_unsafe(d, Range(new_range, r_size));
        ret.original_range_idx_unsafe().resize(dim_num());
        ret.original_range_idx_unsafe()[d].resize(i + 1);
        ret.original_range_idx_unsafe()[d][i++] = r;
      }
    }
  }
}

// Explicit instantiations
template void Subarray::compute_tile_coords<int8_t>();