get_points_of.cc

// © 2021 ETH Zurich, Mechanics and Materials Lab
//
// ae108 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 any
// later version.
//
// ae108 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 ae108. If not, see <https://www.gnu.org/licenses/>.

#include "ae108/meshing/cppgmsh/get_points_of.h"
#include <gmsh.h>
#include <set>

namespace ae108 {
namespace meshing {
namespace cppgmsh {

std::vector<int> get_points_of(const std::pair<int, int> &entity) noexcept {

  std::set<int> tags = {entity.second};
  for (int dim = entity.first; dim > 0; dim--) {
    std::set<int> dim_downward;
    for (const auto &tag : tags) {
      std::vector<int> upward, downward;
      gmsh::model::getAdjacencies(dim, tag, upward, downward);
      dim_downward.insert(downward.begin(), downward.end());
    }
    tags = dim_downward;
  }

  return std::vector(tags.begin(), tags.end());
}

} // namespace cppgmsh
} // namespace meshing
} // namespace ae108