CreateLocalAssemblers.h

Go to the documentation of this file.

 
#pragma once
 
#include <vector>
 
#include "BaseLib/Logging.h"
#include "LocalAssemblerFactoryPython.h"
#include "NumLib/DOF/LocalToGlobalIndexMap.h"
 
namespace ProcessLib
{
namespace BoundaryConditionAndSourceTerm
{
namespace detail
{
template <int GlobalDim,
          template <typename /* shp */, typename /* lower order shp */,
                    typename /* int meth */, int /* global dim */>
          class LocalAssemblerImplementation,
          typename LocalAssemblerInterface, typename... ExtraCtorArgs>
void createLocalAssemblersPython(
    NumLib::LocalToGlobalIndexMap const& dof_table,
    std::vector<MeshLib::Element*> const& mesh_elements,
    std::vector<std::unique_ptr<LocalAssemblerInterface>>& local_assemblers,
    ExtraCtorArgs&&... extra_ctor_args)
{
    static_assert(
        GlobalDim == 1 || GlobalDim == 2 || GlobalDim == 3,
        "Meshes with dimension greater than three are not supported.");
 
    using LocAsmFactory =
        LocalAssemblerFactoryPython<LocalAssemblerInterface,
                                    LocalAssemblerImplementation, GlobalDim,
                                    ExtraCtorArgs...>;
 
    DBUG("Create local assemblers.");
 
    LocAsmFactory factory(dof_table);
    local_assemblers.resize(mesh_elements.size());
 
    DBUG("Calling local assembler builder for all mesh elements.");
    GlobalExecutor::transformDereferenced(
        factory, mesh_elements, local_assemblers,
        std::forward<ExtraCtorArgs>(extra_ctor_args)...);
}
}  // namespace detail
template <template <typename /* shp */, typename /* lower order shp */,
                    typename /* int meth */, int /* global dim */>
          class LocalAssemblerImplementation,
          typename LocalAssemblerInterface, typename... ExtraCtorArgs>
void createLocalAssemblersPython(
    const unsigned dimension,
    std::vector<MeshLib::Element*> const& mesh_elements,
    NumLib::LocalToGlobalIndexMap const& dof_table,
    std::vector<std::unique_ptr<LocalAssemblerInterface>>& local_assemblers,
    ExtraCtorArgs&&... extra_ctor_args)
{
    DBUG("Create local assemblers.");
 
    switch (dimension)
    {
        case 1:
            detail::createLocalAssemblersPython<1,
                                                LocalAssemblerImplementation>(
                dof_table, mesh_elements, local_assemblers,
                std::forward<ExtraCtorArgs>(extra_ctor_args)...);
            break;
        case 2:
            detail::createLocalAssemblersPython<2,
                                                LocalAssemblerImplementation>(
                dof_table, mesh_elements, local_assemblers,
                std::forward<ExtraCtorArgs>(extra_ctor_args)...);
            break;
        case 3:
            detail::createLocalAssemblersPython<3,
                                                LocalAssemblerImplementation>(
                dof_table, mesh_elements, local_assemblers,
                std::forward<ExtraCtorArgs>(extra_ctor_args)...);
            break;
        default:
            OGS_FATAL(
                "Meshes with dimension greater than three are not supported.");
    }
}
}  // namespace BoundaryConditionAndSourceTerm
}  // namespace ProcessLib