/**
 * Stone of Orthanc
 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
 * Department, University Hospital of Liege, Belgium
 * Copyright (C) 2017-2022 Osimis S.A., Belgium
 * Copyright (C) 2021-2022 Sebastien Jodogne, ICTEAM UCLouvain, Belgium
 *
 * This program is free software: you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public License
 * as published by the Free Software Foundation, either version 3 of
 * the License, or (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this program. If not, see
 * <http://www.gnu.org/licenses/>.
 **/

#pragma once

#include <vector>
#include <utility>

#include "../Toolbox/LinearAlgebra.h"

namespace OrthancStone
{
#if 0
  struct Point3D
  {
    Point3D(double x, double y, double z) : x(x), y(y), z(z) {}
    Point3D() : x(0), y(0), z(0) {}
    double x, y, z;
  };

  struct Vector3D
  {
    Vector3D(double x, double y, double z) : x(x), y(y), z(z) {}
    Vector3D() : x(0), y(0), z(0) {}
    double x, y, z;
  };
#else
  typedef Vector Vector3D;
  typedef Vector Point3D;
#endif

  struct Point2D
  {
    Point2D(double x, double y) : x(x), y(y) {}
    Point2D() : x(0), y(0) {}
    double x, y;
  };


  /** Internal */
  struct RtStructRectangleInSlab
  {
    double xmin, xmax, ymin, ymax;
  };
  typedef std::vector<RtStructRectangleInSlab> RtStructRectanglesInSlab;

  enum RectangleBoundaryKind
  {
    RectangleBoundaryKind_Start,
    RectangleBoundaryKind_End
  };

#if 0
  /** Internal */
  void PartitionRectangleList(std::vector< std::vector<size_t> > & sets, const std::vector<RtStructRectanglesInSlab>);
#endif

  /** Internal */
  void ConvertListOfSlabsToSegments(std::vector< std::pair<Point2D, Point2D> >& segments, const std::vector<RtStructRectanglesInSlab>& slabCuts, const size_t totalRectCount);

  /** Internal */
  void AddSlabBoundaries(std::vector<std::pair<double, RectangleBoundaryKind> >& boundaries, const std::vector<RtStructRectanglesInSlab>& slabCuts, size_t iSlab);

  /** Internal */
  void ProcessBoundaryList(std::vector< std::pair<Point2D, Point2D> >& segments, const std::vector<std::pair<double, RectangleBoundaryKind> >& boundaries, double y);

}