You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
56 lines
3.1 KiB
56 lines
3.1 KiB
#include "cdt.h"
|
|
#include "../bounding_box.h"
|
|
#include "../triangle/triangulate.h"
|
|
#include "../remove_duplicate_vertices.h"
|
|
#include "../remove_unreferenced.h"
|
|
|
|
template <
|
|
typename DerivedV,
|
|
typename DerivedE,
|
|
typename DerivedWV,
|
|
typename DerivedWF,
|
|
typename DerivedWE,
|
|
typename DerivedJ>
|
|
IGL_INLINE void igl::triangle::cdt(
|
|
const Eigen::MatrixBase<DerivedV> & V,
|
|
const Eigen::MatrixBase<DerivedE> & E,
|
|
const std::string & flags,
|
|
Eigen::PlainObjectBase<DerivedWV> & WV,
|
|
Eigen::PlainObjectBase<DerivedWF> & WF,
|
|
Eigen::PlainObjectBase<DerivedWE> & WE,
|
|
Eigen::PlainObjectBase<DerivedJ> & J)
|
|
{
|
|
assert(V.cols() == 2);
|
|
assert(E.cols() == 2);
|
|
//MatrixX2S BV;
|
|
//Eigen::MatrixXi BE;
|
|
//igl::bounding_box(V,BV,BE);
|
|
//WV.resize(V.rows()+BV.rows(),2);
|
|
//WV<<V,BV;
|
|
//WE.resize(E.rows()+BE.rows(),2);
|
|
//WE<<E,(BE.array()+V.rows());
|
|
WV = V;
|
|
WE = E;
|
|
Eigen::VectorXi _;
|
|
igl::remove_duplicate_vertices(DerivedWV(WV),DerivedWE(WE),1e-10,WV,_,J,WE);
|
|
// Remove degenerate edges
|
|
const Eigen::Array<bool,Eigen::Dynamic,1> keep = (WE.array().col(0) != WE.array().col(1));
|
|
WE = WE(keep,Eigen::all).eval();
|
|
// c flag must be present
|
|
igl::triangle::triangulate(DerivedWV(WV),WE,DerivedWV(),flags,WV,WF);
|
|
Eigen::VectorXi UJ;
|
|
igl::remove_unreferenced(DerivedV(WV),Eigen::MatrixXi(WF),WV,WF,UJ);
|
|
for(int i=0;i<WE.rows();i++) for(int j=0;j<WE.cols();j++) WE(i,j)=UJ(WE(i,j));
|
|
// Remove edges from box
|
|
//WE.conservativeResize(WE.rows()-BE.rows(),2);
|
|
for(int i=0;i<J.size();i++) J(i)=UJ(J(i));
|
|
//J.conservativeResize(J.size()-BV.rows());
|
|
}
|
|
|
|
#ifdef IGL_STATIC_LIBRARY
|
|
// Explicit template instantiation
|
|
// generated by autoexplicit.sh
|
|
template void igl::triangle::cdt<Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, -1, 1, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, 1, 0, -1, 1> >(Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> > const&, std::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, Eigen::PlainObjectBase<Eigen::Matrix<double, -1, -1, 1, -1, -1> >&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> >&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> >&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, 1, 0, -1, 1> >&);
|
|
// generated by autoexplicit.sh
|
|
template void igl::triangle::cdt<Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, -1, 0, -1, -1>, Eigen::Matrix<int, -1, 1, 0, -1, 1> >(Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> > const&, std::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, Eigen::PlainObjectBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> >&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> >&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, -1, 0, -1, -1> >&, Eigen::PlainObjectBase<Eigen::Matrix<int, -1, 1, 0, -1, 1> >&);
|
|
#endif
|
|
|