DBoundingBox.h

Go to the documentation of this file.

// Copyright (c) 2002-present, OpenMS Inc. -- EKU Tuebingen, ETH Zurich, and FU Berlin
// SPDX-License-Identifier: BSD-3-Clause
//
// --------------------------------------------------------------------------
// $Maintainer: Timo Sachsenberg $
// $Authors: $
// --------------------------------------------------------------------------
 
#pragma once
 
#include <OpenMS/DATASTRUCTURES/DIntervalBase.h>
#include <OpenMS/CONCEPT/HashUtils.h>
#include <OpenMS/CONCEPT/Types.h>
 
#include <functional>
 
namespace OpenMS
{
 
  template <UInt D>
  class DBoundingBox :
    public Internal::DIntervalBase<D>
  {
 
public:
 
    enum {DIMENSION = D};
    typedef Internal::DIntervalBase<D> Base;
    typedef typename Base::PositionType PositionType;
    typedef typename Base::CoordinateType CoordinateType;
 
 
    // for convenience
    using Base::min_;
    using Base::max_;
 
    DBoundingBox() :
      Base()
    {
    }
 
    DBoundingBox(const DBoundingBox& rhs) :
      Base(rhs)
    {
    }
 
    DBoundingBox& operator=(const DBoundingBox& rhs)
    {
      Base::operator=(rhs);
      return *this;
    }
 
    DBoundingBox& operator=(const Base& rhs)
    {
      Base::operator=(rhs);
      return *this;
    }
 
    ~DBoundingBox()
    {
    }
 
    DBoundingBox(const PositionType& minimum, const PositionType& maximum) :
      Base(minimum, maximum)
    {
    }
 
 
 
    void enlarge(const PositionType& p)
    {
      for (UInt i = 0; i < DIMENSION; ++i)
      {
        if (p[i] < min_[i]) min_[i] = p[i];
        if (p[i] > max_[i]) max_[i] = p[i];
      }
    }
 
    void enlarge(CoordinateType x, CoordinateType y)
    {
      enlarge(PositionType(x, y));
    }
 
    //}@
 
 
    bool operator==(const DBoundingBox& rhs) const
    {
      return Base::operator==(rhs);
    }
 
    bool operator==(const Base& rhs) const
    {
      return Base::operator==(rhs);
    }
 
    bool encloses(const PositionType& position) const
    {
      for (UInt i = 0; i < DIMENSION; ++i)
      {
        if (position[i] < min_[i] || position[i] > max_[i])
        {
          return false;
        }
      }
      return true;
    }
 
    bool encloses(CoordinateType x, CoordinateType y) const
    {
      return encloses(PositionType(x, y));
    }
 
    bool intersects(const DBoundingBox& bounding_box) const
    {
      for (UInt i = 0; i < DIMENSION; ++i)
      {
        if (bounding_box.min_[i] > max_[i]) return false;
 
        if (bounding_box.max_[i] <  min_[i]) return false;
      }
      return true;
    }
 
    bool isEmpty() const
    {
      for (UInt i = 0; i != D; i++)
      {
        if (max_[i] <= min_[i])
        {
          return true;
        }
      }
      return false;
    }
 
 
 
  };
 
  template <UInt D>
  std::ostream& operator<<(std::ostream& os, const DBoundingBox<D>& bounding_box)
  {
    os << "--DBOUNDINGBOX BEGIN--" << std::endl;
    os << "MIN --> " << bounding_box.minPosition() << std::endl;
    os << "MAX --> " << bounding_box.maxPosition() << std::endl;
    os << "--DBOUNDINGBOX END--" << std::endl;
    return os;
  }
 
} // namespace OpenMS
 
// Hash function specialization for DBoundingBox
namespace std
{
  template<OpenMS::UInt D>
  struct hash<OpenMS::DBoundingBox<D>>
  {
    std::size_t operator()(const OpenMS::DBoundingBox<D>& bb) const noexcept
    {
      // Hash both min_ and max_ positions (fields used in operator==)
      std::size_t seed = 0;
      // Hash minPosition
      const auto& min_pos = bb.minPosition();
      for (OpenMS::UInt i = 0; i < D; ++i)
      {
        OpenMS::hash_combine(seed, OpenMS::hash_float(min_pos[i]));
      }
      // Hash maxPosition
      const auto& max_pos = bb.maxPosition();
      for (OpenMS::UInt i = 0; i < D; ++i)
      {
        OpenMS::hash_combine(seed, OpenMS::hash_float(max_pos[i]));
      }
      return seed;
    }
  };
} // namespace std