CreateEmbeddedAnchor.cpp

Go to the documentation of this file.

 
#include "CreateEmbeddedAnchor.h"
 
#include "BaseLib/ConfigTree.h"
#include "BaseLib/Logging.h"
#include "EmbeddedAnchor.h"
#include "ParameterLib/Utils.h"
 
namespace ProcessLib
{
template <int GlobalDim>
std::unique_ptr<SourceTermBase> createEmbeddedAnchor(
    BaseLib::ConfigTree const& config, MeshLib::Mesh const& st_mesh,
    MeshLib::Mesh const& bulk_mesh,
    NumLib::LocalToGlobalIndexMap const& dof_table_bulk,
    std::size_t const source_term_mesh_id, const int variable_id)
{
    DBUG("Constructing EmbeddedAnchor from config.");
    // ...special is necessary, because EmbeddedAnchor has no configuration via
    // child tags, hence it wouldn't be handled properly by the documentation
    // scripts otherwise.
    config.checkConfigParameter("type", "EmbeddedAnchor");
 
    for (MeshLib::Element const* const element : st_mesh.getElements())
    {
        if (element->getNumberOfNodes() != 2)
        {
            OGS_FATAL(
                "Every anchor element needs to have precisely two nodes.");
        }
    }
#ifdef USE_PETSC
    OGS_FATAL(
        "The EmbeddedAnchor source term has not been tested with PETSc yet.");
#endif
 
    return std::make_unique<EmbeddedAnchor<GlobalDim>>(
        bulk_mesh, dof_table_bulk, source_term_mesh_id, st_mesh, variable_id);
}
 
template std::unique_ptr<SourceTermBase> createEmbeddedAnchor<2>(
    BaseLib::ConfigTree const& config, MeshLib::Mesh const& st_mesh,
    MeshLib::Mesh const& bulk_mesh,
    NumLib::LocalToGlobalIndexMap const& dof_table_bulk,
    std::size_t const source_term_mesh_id, const int variable_id);
 
template std::unique_ptr<SourceTermBase> createEmbeddedAnchor<3>(
    BaseLib::ConfigTree const& config, MeshLib::Mesh const& st_mesh,
    MeshLib::Mesh const& bulk_mesh,
    NumLib::LocalToGlobalIndexMap const& dof_table_bulk,
    std::size_t const source_term_mesh_id, const int variable_id);
}  // namespace ProcessLib