// // IGL Lib - Simple C++ mesh library // // Copyright 2011, Daniele Panozzo. All rights reserved. #ifndef TT_H #define TT_H #include #include #include "ismanifold.h" #include namespace igl { // Preprocessing template inline void tt_preprocess(const Eigen::Matrix& V, const Eigen::MatrixXi& F, std::vector >& TTT) { for(int f=0;f v2) std::swap(v1,v2); std::vector r(4); r[0] = v1; r[1] = v2; r[2] = f; r[3] = i; TTT.push_back(r); } std::sort(TTT.begin(),TTT.end()); } // Extract the face adjacencies inline void tt_extractTT(const Eigen::MatrixXi& F, std::vector >& TTT, Eigen::MatrixXi& TT) { TT = Eigen::MatrixXi::Constant((int)(F.rows()),3,-1); for(int i=1;i& r1 = TTT[i-1]; std::vector& r2 = TTT[i]; if ((r1[0] == r2[0]) && (r1[1] == r2[1])) { TT(r1[2],r1[3]) = r2[2]; TT(r2[2],r2[3]) = r1[2]; } } } // Extract the face adjacencies indices (needed for fast traversal) inline void tt_extractTTi(const Eigen::MatrixXi& F, std::vector >& TTT, Eigen::MatrixXi& TTi) { TTi = Eigen::MatrixXi::Constant((int)(F.rows()),3,-1); for(int i=1;i& r1 = TTT[i-1]; std::vector& r2 = TTT[i]; if ((r1[0] == r2[0]) && (r1[1] == r2[1])) { TTi(r1[2],r1[3]) = r2[3]; TTi(r2[2],r2[3]) = r1[3]; } } } // Compute triangle-triangle adjacency template inline void tt(const Eigen::Matrix& V, const Eigen::MatrixXi& F, Eigen::MatrixXi& TT) { assert(igl::isManifold(V,F)); std::vector > TTT; tt_preprocess(V,F,TTT); tt_extractTT(F,TTT,TT); } // Compute triangle-triangle adjacency with indices template inline void tt(const Eigen::Matrix& V, const Eigen::MatrixXi& F, Eigen::MatrixXi& TT, Eigen::MatrixXi& TTi) { assert(igl::isManifold(V,F)); std::vector > TTT; tt_preprocess(V,F,TTT); tt_extractTT(F,TTT,TT); tt_extractTTi(F,TTT,TTi); } } #endif