LineSegmentIntersection.h

Go to the documentation of this file.

/* This file is part of MoDELib, the Mechanics Of Defects Evolution Library.
 *
 *
 * MoDELib is distributed without any warranty under the
 * GNU General Public License (GPL) v2 <http://www.gnu.org/licenses/>.
 */
 
#ifndef model_LineSegmentIntersection_H_
#define model_LineSegmentIntersection_H_
 
#include <tuple>
#include <map>
#include <stdexcept>
#include <Eigen/Dense>
 
namespace model
{
    
    template <int dim>
    class LineSegmentIntersection
    {
        
    public:
        enum IntersectionType {PARALLEL,COINCIDENT,INCIDENT,SKEW,OUTSIDE};
        
        
        typedef Eigen::Matrix<double,dim,1> VectorDimD;
        typedef std::tuple<IntersectionType,VectorDimD,VectorDimD,double,double> SolutionType;
        
    private:
        const SolutionType sol;
        
        /**********************************************************************/
        static SolutionType findIntersections(const VectorDimD& A0,
                                              const VectorDimD& D0,
                                              const VectorDimD& A1,
                                              const VectorDimD& D1);
        
    public:
        
        const IntersectionType& type;
        const VectorDimD& x0;
        const VectorDimD& x1;
        const double& u0;
        const double& u1;
        
        /**********************************************************************/
        LineSegmentIntersection(const VectorDimD& A0,
                             const VectorDimD& D0,
                             const VectorDimD& P0,
                             const VectorDimD& P1);
        
    };
    
    /******************************************************************/
    /******************************************************************/
} /* namespace model */
#endif