libigl/include/igl/predicates/point_inside_convex_polygon.cpp

// This file is part of libigl, a simple c++ geometry processing library.
//
// Copyright (C) 2019 Hanxiao Shen <[email protected]>
//
// This Source Code Form is subject to the terms of the Mozilla Public License
// v. 2.0. If a copy of the MPL was not distributed with this file, You can
// obtain one at http://mozilla.org/MPL/2.0/.

#include "point_inside_convex_polygon.h"
#include "orient2d.h"

template <typename DerivedP, typename DerivedQ>
IGL_INLINE bool igl::predicates::point_inside_convex_polygon(
  const Eigen::MatrixBase<DerivedP>& P,
  const Eigen::MatrixBase<DerivedQ>& q
){
  EIGEN_STATIC_ASSERT_MATRIX_SPECIFIC_SIZE(DerivedP, Eigen::Dynamic, 2);
  EIGEN_STATIC_ASSERT_MATRIX_SPECIFIC_SIZE(DerivedQ, 1, 2);
  typedef typename DerivedP::Scalar Scalar;
  for(int i=0;i<P.rows();i++){
    int i_1 = (i+1) % P.rows();
    Eigen::Matrix<Scalar,1,2> a = P.row(i);
    Eigen::Matrix<Scalar,1,2> b = P.row(i_1);
    auto r = igl::predicates::orient2d(a,b,q);
    if(r == igl::predicates::Orientation::COLLINEAR || 
       r == igl::predicates::Orientation::NEGATIVE)
      return false;
  }
  return true;
}

#ifdef IGL_STATIC_LIBRARY
template bool igl::predicates::point_inside_convex_polygon<Eigen::Matrix<double, -1, 2, 0, -1, 2>, Eigen::Matrix<double, 1, 2, 1, 1, 2> >(Eigen::MatrixBase<Eigen::Matrix<double, -1, 2, 0, -1, 2> > const&, Eigen::MatrixBase<Eigen::Matrix<double, 1, 2, 1, 1, 2> > const&);
#endif