Revision b973b6f4ef165a3688527a98348f35d6819a9081 authored by Christoph Lehmann on 28 March 2023, 15:33:42 UTC, committed by Christoph Lehmann on 28 March 2023, 15:33:42 UTC
Draft: Debug output for local and global matrices See merge request ogs/ogs!4546
MeshNodesOnPoint.cpp
/**
* \file
* \copyright
* Copyright (c) 2012-2023, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*/
#include "MeshNodesOnPoint.h"
#include "MeshLib/Mesh.h"
namespace MeshGeoToolsLib
{
MeshNodesOnPoint::MeshNodesOnPoint(MeshLib::Mesh const& mesh,
GeoLib::Grid<MeshLib::Node> const& mesh_grid,
GeoLib::Point const& pnt,
double epsilon_radius,
SearchAllNodes search_all_nodes)
: _mesh(mesh), _pnt(pnt)
{
std::vector<std::size_t> vec_ids(
mesh_grid.getPointsInEpsilonEnvironment(pnt, epsilon_radius));
if (search_all_nodes == SearchAllNodes::Yes)
{
_msh_node_ids = vec_ids;
}
else
{
for (auto id : vec_ids)
{
if (MeshLib::isBaseNode(*_mesh.getNode(id),
_mesh.getElementsConnectedToNode(id)))
{
_msh_node_ids.push_back(id);
}
}
}
}
} // end namespace MeshGeoToolsLib
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...