Revision 47c3824f6c23f658b1f1bc227791082514f7d9f9 authored by Lars Bilke on 18 February 2021, 10:04:01 UTC, committed by Lars Bilke on 19 February 2021, 07:13:01 UTC
1 parent b39d6c5
PointRule1.h
/**
* \file
* \copyright
* Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*
*/
#pragma once
#include "MeshLib/MeshEnums.h"
#include "Element.h"
#include "VertexRule.h"
#include "EdgeReturn.h"
namespace MeshLib
{
/// A 0d point element.
class PointRule1 : public VertexRule
{
public:
/// Constant: The number of base nodes for this element
static const unsigned n_base_nodes = 1u;
/// Constant: The number of all nodes for this element
static const unsigned n_all_nodes = 1u;
/// Constant: The geometric type of the element
static const MeshElemType mesh_elem_type = MeshElemType::POINT;
/// Constant: The FEM type of the element
static const CellType cell_type = CellType::POINT1;
/// Constant: The number of neighbors
static const unsigned n_neighbors = 2;
/// Constant: Local node index table for edge
static const unsigned edge_nodes[1][1];
/// Edge rule
using EdgeReturn = MeshLib::NoEdgeReturn;
/// \copydoc MeshLib::Element::isPntInElement()
///
/// Specifically this function tests if the squared euclidean distance
/// between the points \c nodes and \c pnt is less then \c eps.
///
/// @param nodes the nodes of the element.
static bool isPntInElement(Node const* const* nodes,
MathLib::Point3d const& pnt, double eps);
/// Tests if the element is geometrically valid.
static ElementErrorCode validate(const Element* e);
/// Returns the ID of a face given an array of nodes.
static unsigned identifyFace(Node const* const* /*_nodes*/, Node* nodes[1]);
/// Calculates the length of a line
static double computeVolume(Node const* const* _nodes);
};
} // namespace MeshLib
Computing file changes ...