DIntervalBase.h

Go to the documentation of this file.

// --------------------------------------------------------------------------
//                   OpenMS -- Open-Source Mass Spectrometry
// --------------------------------------------------------------------------
// Copyright The OpenMS Team -- Eberhard Karls University Tuebingen,
// ETH Zurich, and Freie Universitaet Berlin 2002-2022.
//
// This software is released under a three-clause BSD license:
//  * Redistributions of source code must retain the above copyright
//    notice, this list of conditions and the following disclaimer.
//  * Redistributions in binary form must reproduce the above copyright
//    notice, this list of conditions and the following disclaimer in the
//    documentation and/or other materials provided with the distribution.
//  * Neither the name of any author or any participating institution
//    may be used to endorse or promote products derived from this software
//    without specific prior written permission.
// For a full list of authors, refer to the file AUTHORS.
// --------------------------------------------------------------------------
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
// AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
// ARE DISCLAIMED. IN NO EVENT SHALL ANY OF THE AUTHORS OR THE CONTRIBUTING
// INSTITUTIONS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
// OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
// WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
// OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
// ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
//
// --------------------------------------------------------------------------
// $Maintainer: Timo Sachsenberg $
// $Authors: Clemens Groepl, Marc Sturm $
// --------------------------------------------------------------------------

#pragma once

#include <OpenMS/DATASTRUCTURES/DPosition.h>
#include <OpenMS/CONCEPT/Types.h>

#include <utility>

namespace OpenMS
{
  namespace Internal
  {
    template <UInt D>
    class DIntervalBase
    {
public:

      enum {DIMENSION = D};
      typedef DPosition<D>                            PositionType;
      typedef typename PositionType::CoordinateType CoordinateType;


      DIntervalBase() :
        min_(PositionType::maxPositive()),
        max_(PositionType::minNegative())
      {
      }

      DIntervalBase(const DIntervalBase& rhs) :
        min_(rhs.min_),
        max_(rhs.max_)
      {
      }

      DIntervalBase(DIntervalBase&&) noexcept = default;

      DIntervalBase& operator=(const DIntervalBase& rhs)
      {
        min_ = rhs.min_;
        max_ = rhs.max_;
        return *this;
      }

      ~DIntervalBase()
      {
      }

      DIntervalBase(PositionType const& minimum, PositionType const& maximum) :
        min_(minimum),
        max_(maximum)
      {
        normalize_();
      }



      inline PositionType const& minPosition() const
      {
        return min_;
      }

      inline PositionType const& maxPosition() const
      {
        return max_;
      }

      void setMin(PositionType const& position)
      {
        min_ = position;
        for (UInt i = 0; i < DIMENSION; ++i)
        {
          if (min_[i] > max_[i]) max_[i] = min_[i];
        }
      }

      void setMax(PositionType const& position)
      {
        max_ = position;
        for (UInt i = 0; i < DIMENSION; ++i)
        {
          if (min_[i] > max_[i]) min_[i] = max_[i];
        }
      }

      void setMinMax(PositionType const& min, PositionType const& max)
      {
        min_ = min;
        max_ = max;
        normalize_();
      }

      template <UInt D2>
      void assign(const DIntervalBase<D2> rhs)
      {
        for (UInt i = 0; i < std::min(D, D2); ++i)
        {
          min_[i] = rhs.minPosition()[i];
          max_[i] = rhs.maxPosition()[i];
        }
      }

      //}@

      bool operator==(const DIntervalBase& rhs) const
      {
        return (min_ == rhs.min_) && (max_ == rhs.max_);
      }

      bool operator!=(const DIntervalBase& rhs) const
      {
        return !(operator==(rhs));
      }

      DIntervalBase operator+(const PositionType& point) const
      {
        DIntervalBase result(*this);
        result += point;
        return result;
      }

      DIntervalBase& operator+=(const PositionType& point)
      {
        this->min_ += point;
        this->max_ += point;
        return *this;
      }

      DIntervalBase operator-(const PositionType& point) const
      {
        DIntervalBase result(*this);
        result -= point;
        return result;
      }

      DIntervalBase& operator-=(const PositionType& point)
      {
        this->min_ -= point;
        this->max_ -= point;
        return *this;
      }

      inline void clear()
      {
        *this = empty;
      }
      
      bool isEmpty() const
      {
        return *this == empty;
      }

      bool isEmpty(UInt dim) const
      {
        return DIntervalBase<1>(std::make_pair(min_[dim], max_[dim])) == DIntervalBase<1>::empty;
      }

      void setDimMinMax(UInt dim, const DIntervalBase<1>& min_max)
      {
        min_[dim] = min_max.min_[0];
        max_[dim] = min_max.max_[0];
      }

      // make all other templates friends to allow accessing min_ and max_ in setDimMinMax 
      template<UInt> friend class DIntervalBase; 



      PositionType center() const
      {
        PositionType center(min_);
        center += max_;
        center /= 2;
        return center;
      }

      PositionType diagonal() const
      {
        return max_ - min_;
      }

      static DIntervalBase const empty;
      static DIntervalBase const zero;

      //}@


      inline CoordinateType minX() const
      {
        return min_[0];
      }

      inline CoordinateType minY() const
      {
        return min_[1];
      }

      inline CoordinateType maxX() const
      {
        return max_[0];
      }

      inline CoordinateType maxY() const
      {
        return max_[1];
      }

      void setMinX(CoordinateType const c)
      {
        min_[0] = c;
        if (min_[0] > max_[0]) max_[0] = min_[0];
      }

      void setMinY(CoordinateType const c)
      {
        min_[1] = c;
        if (min_[1] > max_[1]) max_[1] = min_[1];
      }

      void setMaxX(CoordinateType const c)
      {
        max_[0] = c;
        if (min_[0] > max_[0]) min_[0] = max_[0];
      }

      void setMaxY(CoordinateType const c)
      {
        max_[1] = c;
        if (min_[1] > max_[1]) min_[1] = max_[1];
      }

      inline CoordinateType width() const
      {
        return max_[0] - min_[0];
      }

      inline CoordinateType height() const
      {
        return max_[1] - min_[1];
      }


protected:

      PositionType min_;

      PositionType max_;

      void normalize_()
      {
        for (UInt i = 0; i < DIMENSION; ++i)
        {
          if (min_[i] > max_[i])
          {
            std::swap(min_[i], max_[i]);
          }
        }
      }

      DIntervalBase(const std::pair<PositionType, PositionType>& pair) :
        min_(pair.first),
        max_(pair.second)
      {

      }

    };

    template <UInt D>
    DIntervalBase<D> const DIntervalBase<D>::zero
      = DIntervalBase<D>(DIntervalBase<D>::PositionType::zero(), DIntervalBase<D>::PositionType::zero());

    template <UInt D>
    DIntervalBase<D> const DIntervalBase<D>::empty
      = DIntervalBase<D>(std::make_pair(DIntervalBase<D>::PositionType::maxPositive(), DIntervalBase<D>::PositionType::minNegative()));

    template <UInt D>
    std::ostream& operator<<(std::ostream& os, const DIntervalBase<D>& rhs)
    {
      os << "--DIntervalBase BEGIN--" << std::endl;
      os << "MIN --> " << rhs.minPosition() << std::endl;
      os << "MAX --> " << rhs.maxPosition() << std::endl;
      os << "--DIntervalBase END--" << std::endl;
      return os;
    }

  } // namespace Internal

} // namespace OpenMS