Point2D.hpp

Go to the documentation of this file.

/*
  Copyright 2015 Statoil ASA.
 
  This file is part of the Open Porous Media project (OPM).
 
  OPM is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
 
  OPM 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 General Public License for more details.
 
  You should have received a copy of the GNU General Public License
  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#ifndef OPM_POINT2D_HEADER_INCLUDED
#define OPM_POINT2D_HEADER_INCLUDED
 
namespace Opm {
 
 
 
    namespace detail {
 
 
        class Point2D
        {
            public:
                Point2D(const double xi, const double yi)
                    : x_(xi),
                      y_(yi) {
 
                      }
 
                Point2D()
                {
                }
 
                double getX() const
                {
                    return x_;
                }
 
                double getY() const
                {
                    return y_;
                }
 
                void setX(const double x)
                {
                    x_ = x;
                }
 
                void setY(const double y)
                {
                    y_ = y;
                }
 
                static bool findIntersection(Point2D line_segment1[2], Point2D line2[2], Point2D& intersection_point)
                {
 
                    const double x1 = line_segment1[0].getX();
                    const double y1 = line_segment1[0].getY();
                    const double x2 = line_segment1[1].getX();
                    const double y2 = line_segment1[1].getY();
 
                    const double x3 = line2[0].getX();
                    const double y3 = line2[0].getY();
                    const double x4 = line2[1].getX();
                    const double y4 = line2[1].getY();
 
                    const double d = (x1 - x2) * (y3 - y4) - (y1 - y2) * (x3 - x4);
 
                    if (d == 0.) {
                        return false;
                    }
 
                    const double x = ((x3 - x4) * (x1 * y2 - y1 * x2) - (x1 - x2) * (x3 * y4 - y3 * x4)) / d;
                    const double y = ((y3 - y4) * (x1 * y2 - y1 * x2) - (y1 - y2) * (x3 * y4 - y3 * x4)) / d;
 
                    if (x >= std::min(x1,x2) && x <= std::max(x1,x2)) {
                        intersection_point.setX(x);
                        intersection_point.setY(y);
                        return true;
                    } else {
                       return false;
                    }
                }
 
            private:
                double x_;
                double y_;
 
        }; // class Point2D
 
 
    } // namespace detail
 
 
 
} // namespace Opm
 
#endif // OPM_POINT2D_HEADER_INCLUDED