Revision 3c652339ac466c5fcc10155a752fb17d5a8ca197 authored by Tobias Meisel on 18 December 2020, 19:17:01 UTC, committed by Tobias Meisel on 18 December 2020, 19:17:01 UTC
[MeL/IO] Parallel writer based on HDF5 and MPI IO See merge request ogs/ogs!3307
Node.cpp
/**
* \file
* \author Karsten Rink
* \date 2012-05-02
* \brief Implementation of the Node class.
*
* \copyright
* Copyright (c) 2012-2020, 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 "MeshLib/Node.h"
#include "Elements/Element.h"
namespace MeshLib {
Node::Node(const double coords[3], std::size_t id)
: MathLib::Point3dWithID(
std::array<double,3>{{coords[0], coords[1], coords[2]}}, id)
{
}
Node::Node(std::array<double, 3> const& coords, std::size_t id)
: MathLib::Point3dWithID(coords, id)
{
}
Node::Node(double x, double y, double z, std::size_t id)
: MathLib::Point3dWithID(std::array<double,3>({{x, y, z}}), id)
{
}
Node::Node(const Node &node)
: MathLib::Point3dWithID(node._x, node.getID())
{
}
void Node::updateCoordinates(double x, double y, double z)
{
_x[0] = x;
_x[1] = y;
_x[2] = z;
const std::size_t nElements (this->_elements.size());
for (std::size_t i = 0; i < nElements; i++)
{
_elements[i]->computeVolume();
}
}
bool isBaseNode(Node const& node)
{
// Check if node is connected.
if (node.getNumberOfElements() == 0)
{
return true;
}
// In a mesh a node always belongs to at least one element.
auto const e = node.getElement(0);
auto const n_base_nodes = e->getNumberOfBaseNodes();
auto const local_index = e->getNodeIDinElement(&node);
return local_index < n_base_nodes;
}
} // namespace MeshLib
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...