QuadraticMeshGenerator.cpp File Reference

Detailed Description

Copyright

Copyright (c) 2012-2024, OpenGeoSys Community (http://www.opengeosys.org) Distributed under a Modified BSD License. See accompanying file LICENSE.txt or http://www.opengeosys.org/project/license

Definition in file QuadraticMeshGenerator.cpp.

#include "QuadraticMeshGenerator.h"
#include <set>
#include "MeshLib/Elements/Element.h"
#include "MeshLib/Elements/Hex.h"
#include "MeshLib/Elements/Line.h"
#include "MeshLib/Elements/Prism.h"
#include "MeshLib/Elements/Pyramid.h"
#include "MeshLib/Elements/Quad.h"
#include "MeshLib/Elements/Tet.h"
#include "MeshLib/Elements/Tri.h"
#include "MeshLib/Mesh.h"
#include "MeshLib/Node.h"
#include "MeshLib/Utils/DuplicateMeshComponents.h"

Include dependency graph for QuadraticMeshGenerator.cpp:

Go to the source code of this file.

Classes
struct	nodeByCoordinatesComparator

Namespaces
namespace	MeshToolsLib

Functions
template<typename QuadraticElement >
std::unique_ptr< QuadraticElement >	convertLinearToQuadratic (MeshLib::Element const &e)
template<>
std::unique_ptr< MeshLib::Quad9 >	convertLinearToQuadratic< MeshLib::Quad9 > (MeshLib::Element const &e)
	Special case for Quad-9 adding a centre node too.
template<>
std::unique_ptr< MeshLib::Prism15 >	convertLinearToQuadratic< MeshLib::Prism15 > (MeshLib::Element const &e)
	Special case for Prism15.
std::unique_ptr< MeshLib::Element >	createQuadraticElement (MeshLib::Element const &e, bool const add_centre_node)
	Return a new quadratic element corresponding to the linear element's type.
std::unique_ptr< MeshLib::Mesh >	MeshToolsLib::createQuadraticOrderMesh (MeshLib::Mesh const &linear_mesh, bool const add_centre_node)

Function Documentation

convertLinearToQuadratic()

template<typename QuadraticElement >

std::unique_ptr< QuadraticElement > convertLinearToQuadratic

(

MeshLib::Element const &

e

)

Given an (linear) element divide all its edges by inserting a point in the middle and return a new element.

Definition at line 29 of file QuadraticMeshGenerator.cpp.

{
    int const n_all_nodes = QuadraticElement::n_all_nodes;
    int const n_base_nodes = QuadraticElement::n_base_nodes;
    assert(n_base_nodes == e.getNumberOfBaseNodes());
 
    // Copy base nodes of element to the quadratic element new nodes'.
    std::array<MeshLib::Node*, n_all_nodes> nodes{};
    for (int i = 0; i < n_base_nodes; i++)
    {
        nodes[i] = const_cast<MeshLib::Node*>(e.getNode(i));
    }
 
    // For each edge create a middle node.
    int const number_of_edges = e.getNumberOfEdges();
    for (int i = 0; i < number_of_edges; i++)
    {
        auto const& a = *e.getEdgeNode(i, 0);
        auto const& b = *e.getEdgeNode(i, 1);
 
        nodes[n_base_nodes + i] = new MeshLib::Node(
            (a[0] + b[0]) / 2, (a[1] + b[1]) / 2, (a[2] + b[2]) / 2);
    }
 
    return std::make_unique<QuadraticElement>(nodes, e.getID());
}

References MeshLib::Element::getEdgeNode(), MeshLib::Element::getID(), MeshLib::Element::getNode(), MeshLib::Element::getNumberOfBaseNodes(), MeshLib::Element::getNumberOfEdges(), and MeshLib::Node.

Referenced by createQuadraticElement().

convertLinearToQuadratic< MeshLib::Prism15 >()

template<>

std::unique_ptr< MeshLib::Prism15 > convertLinearToQuadratic< MeshLib::Prism15 >

(

MeshLib::Element const &

e

)

Special case for Prism15.

Definition at line 29 of file QuadraticMeshGenerator.cpp.

{
    int const n_all_nodes = MeshLib::Prism15::n_all_nodes;
    int const n_base_nodes = MeshLib::Prism15::n_base_nodes;
    assert(n_base_nodes == e.getNumberOfBaseNodes());
 
    // Copy base nodes of element to the quadratic element new nodes'.
    std::array<MeshLib::Node*, n_all_nodes> nodes{};
    for (int i = 0; i < n_base_nodes; i++)
    {
        nodes[i] = const_cast<MeshLib::Node*>(e.getNode(i));
    }
 
    // For each edge create a middle node.
    int const number_of_edges = e.getNumberOfEdges();
    for (int i = 0; i < 3; i++)
    {
        auto const& a = *e.getEdgeNode(i, 0);
        auto const& b = *e.getEdgeNode(i, 1);
 
        nodes[n_base_nodes + i] = new MeshLib::Node(
            (a[0] + b[0]) / 2, (a[1] + b[1]) / 2, (a[2] + b[2]) / 2);
    }
    for (int i = 3; i < 6; i++)
    {
        auto const& a = *e.getEdgeNode(i + 3, 0);
        auto const& b = *e.getEdgeNode(i + 3, 1);
 
        nodes[n_base_nodes + i] = new MeshLib::Node(
            (a[0] + b[0]) / 2, (a[1] + b[1]) / 2, (a[2] + b[2]) / 2);
    }
    for (int i = 6; i < number_of_edges; i++)
    {
        auto const& a = *e.getEdgeNode(i - 3, 0);
        auto const& b = *e.getEdgeNode(i - 3, 1);
 
        nodes[n_base_nodes + i] = new MeshLib::Node(
            (a[0] + b[0]) / 2, (a[1] + b[1]) / 2, (a[2] + b[2]) / 2);
    }
 
    return std::make_unique<MeshLib::Prism15>(nodes, e.getID());
}

Referenced by createQuadraticElement().

convertLinearToQuadratic< MeshLib::Quad9 >()

template<>

std::unique_ptr< MeshLib::Quad9 > convertLinearToQuadratic< MeshLib::Quad9 >

(

MeshLib::Element const &

e

)

Special case for Quad-9 adding a centre node too.

Definition at line 29 of file QuadraticMeshGenerator.cpp.

{
    int const n_all_nodes = MeshLib::Quad9::n_all_nodes;
    int const n_base_nodes = MeshLib::Quad9::n_base_nodes;
    assert(n_base_nodes == e.getNumberOfBaseNodes());
 
    // Copy base nodes of element to the quadratic element new nodes'.
    std::array<MeshLib::Node*, n_all_nodes> nodes{};
    for (int i = 0; i < n_base_nodes; i++)
    {
        nodes[i] = const_cast<MeshLib::Node*>(e.getNode(i));
    }
 
    // For each edge create a middle node.
    int const number_of_edges = e.getNumberOfEdges();
    for (int i = 0; i < number_of_edges; i++)
    {
        auto const& a = *e.getEdgeNode(i, 0);
        auto const& b = *e.getEdgeNode(i, 1);
 
        nodes[n_base_nodes + i] = new MeshLib::Node(
            (a[0] + b[0]) / 2, (a[1] + b[1]) / 2, (a[2] + b[2]) / 2);
    }
 
    // Compute the centre point coordinates.
    auto* centre_node = new MeshLib::Node(0, 0, 0);
    for (int i = 0; i < n_base_nodes; i++)
    {
        for (int d = 0; d < 3; d++)
        {
            (*centre_node)[d] += (*nodes[i])[d] / n_base_nodes;
        }
    }
    nodes[n_all_nodes - 1] = centre_node;
 
    return std::make_unique<MeshLib::Quad9>(nodes, e.getID());
}

Referenced by createQuadraticElement().

createQuadraticElement()

std::unique_ptr< MeshLib::Element > createQuadraticElement	(	MeshLib::Element const &	e,
		bool const	add_centre_node )

Return a new quadratic element corresponding to the linear element's type.

Definition at line 145 of file QuadraticMeshGenerator.cpp.

{
    if (e.getCellType() == MeshLib::CellType::LINE2)
    {
        return convertLinearToQuadratic<MeshLib::Line3>(e);
    }
    if (e.getCellType() == MeshLib::CellType::TRI3)
    {
        return convertLinearToQuadratic<MeshLib::Tri6>(e);
    }
    if (e.getCellType() == MeshLib::CellType::TET4)
    {
        return convertLinearToQuadratic<MeshLib::Tet10>(e);
    }
    if (e.getCellType() == MeshLib::CellType::QUAD4)
    {
        if (add_centre_node)
        {
            return convertLinearToQuadratic<MeshLib::Quad9>(e);
        }
        return convertLinearToQuadratic<MeshLib::Quad8>(e);
    }
    if (e.getCellType() == MeshLib::CellType::HEX8)
    {
        return convertLinearToQuadratic<MeshLib::Hex20>(e);
    }
    if (e.getCellType() == MeshLib::CellType::PRISM6)
    {
        return convertLinearToQuadratic<MeshLib::Prism15>(e);
    }
    if (e.getCellType() == MeshLib::CellType::PYRAMID5)
    {
        return convertLinearToQuadratic<MeshLib::Pyramid13>(e);
    }
 
    OGS_FATAL("Mesh element type {:s} is not supported",
              MeshLib::CellType2String(e.getCellType()));
}

References MeshLib::CellType2String(), convertLinearToQuadratic(), convertLinearToQuadratic< MeshLib::Prism15 >(), convertLinearToQuadratic< MeshLib::Quad9 >(), MeshLib::Element::getCellType(), MeshLib::HEX8, MeshLib::LINE2, OGS_FATAL, MeshLib::PRISM6, MeshLib::PYRAMID5, MeshLib::QUAD4, MeshLib::TET4, and MeshLib::TRI3.

Referenced by MeshToolsLib::createQuadraticOrderMesh().