bbox.hpp

Go to the documentation of this file.

// --------------------------------------------------------------------------
//
// Copyright (C) 2025-present by Pablo Antolin
//
// This file is part of the QUGaR library.
//
// SPDX-License-Identifier:    MIT
//
// --------------------------------------------------------------------------
 
#ifndef QUGAR_LIBRARY_BBOX_HPP
 
#define QUGAR_LIBRARY_BBOX_HPP

 
 
#include <qugar/concepts.hpp>
#include <qugar/point.hpp>
#include <qugar/types.hpp>
 
#include <algoim/hyperrectangle.hpp>
 
#include <array>
#include <cstddef>
 
namespace qugar {

template<int dim> class BoundBox
{
 
public:


  BoundBox();

  BoundBox(const std::array<real, static_cast<std::size_t>(dim)> &min,
    const std::array<real, static_cast<std::size_t>(dim)> &max);

  BoundBox(const real min, const real max);

  BoundBox(const Point<dim> &min, const Point<dim> &max);

  explicit BoundBox(const ::algoim::HyperRectangle<real, dim> &rectangle);

 
private:
  Point<dim> min_;

  Point<dim> max_;
 
public:
  void set(const Point<dim> &min, const Point<dim> &max);

  void extend(const Point<dim> &point);



  [[nodiscard]] real min(int dir) const;

  [[nodiscard]] real max(int dir) const;

  template<int dim_aux = dim>
    requires Is1D<dim_aux, dim>
  [[nodiscard]] real min() const;

  template<int dim_aux = dim>
    requires Is1D<dim_aux, dim>
  [[nodiscard]] real max() const;

  [[nodiscard]] const Point<dim> &min_corner() const;

  [[nodiscard]] const Point<dim> &max_corner() const;

  [[nodiscard]] BoundBox<dim> extend(real delta) const;

  ::algoim::HyperRectangle<real, dim> to_hyperrectangle() const;

  [[nodiscard]] Point<dim> get_lengths() const;

  [[nodiscard]] real length(int dir) const;

  template<int dim_aux = dim>
    requires Is1D<dim_aux, dim>
  [[nodiscard]] real length() const;

  [[nodiscard]] real volume() const;

  [[nodiscard]] Point<dim> mid_point() const;

  [[nodiscard]] Point<dim> scale_to_new_domain(const BoundBox<dim> &new_domain, const Point<dim> &point) const;

  [[nodiscard]] Point<dim> scale_to_0_1(const Point<dim> &point) const;

  [[nodiscard]] Point<dim> scale_from_0_1(const Point<dim> &point_01) const;

  template<int dim_aux = dim>
    requires(dim_aux == dim && dim > 1)
  [[nodiscard]] BoundBox<dim - 1> slice(const int const_dir) const
  {
    Point<dim - 1> _min_corner;
    Point<dim - 1> _max_corner;
    for (int dir = 0, dir2 = 0; dir < dim; ++dir) {
      if (dir != const_dir) {
        _min_corner(dir2) = this->min(dir);
        _max_corner(dir2) = this->max(dir);
        ++dir2;
      }
    }
    return BoundBox<dim - 1>(_min_corner, _max_corner);
  }
 

};

using Interval = BoundBox<1>;
 
}// namespace qugar
 
#endif// QUGAR_LIBRARY_BBOX_HPP