CreateLocalAssemblersTaylorHood.h

Go to the documentation of this file.

 
#pragma once
 
#include <vector>
 
#include "BaseLib/Logging.h"
#include "LocalAssemblerFactoryTaylorHood.h"
#include "NumLib/DOF/LocalToGlobalIndexMap.h"
 
namespace ProcessLib
{
namespace detail
{
template <template <typename /*LocalAssemblerInterface*/,
                    template <typename /* shp fct */,
                              typename /* lower order shp fct */,
                              int /* global dim */>
                    class /*LocalAssemblerImplementation*/,
                    int /* global dim */, typename... /*ConstructorArgs*/>
          class LocalAssemblerFactory,
          int GlobalDim,
          template <typename /* shp fct */, typename /* lower order shp fct */,
                    int /* global dim */>
          class LocalAssemblerImplementation,
          typename LocalAssemblerInterface, typename... ExtraCtorArgs>
void createLocalAssemblersTaylorHood(
    std::vector<MeshLib::Element*> const& mesh_elements,
    NumLib::LocalToGlobalIndexMap const& dof_table,
    std::vector<std::unique_ptr<LocalAssemblerInterface>>& local_assemblers,
    NumLib::IntegrationOrder const integration_order,
    ExtraCtorArgs&&... extra_ctor_args)
{
    using LocAsmFac = LocalAssemblerFactory<LocalAssemblerInterface,
                                            LocalAssemblerImplementation,
                                            GlobalDim, ExtraCtorArgs...>;
 
    DBUG("Create local assemblers.");
 
    LocAsmFac factory(dof_table, integration_order);
    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 <int GlobalDim,
          template <typename /* shp fct */, typename /* lower order shp fct */,
                    int /* global dim */>
          class LocalAssemblerImplementation,
          typename LocalAssemblerInterface, typename... ExtraCtorArgs>
void createLocalAssemblersHM(
    std::vector<MeshLib::Element*> const& mesh_elements,
    NumLib::LocalToGlobalIndexMap const& dof_table,
    std::vector<std::unique_ptr<LocalAssemblerInterface>>& local_assemblers,
    NumLib::IntegrationOrder const integration_order,
    ExtraCtorArgs&&... extra_ctor_args)
{
    detail::createLocalAssemblersTaylorHood<LocalAssemblerFactoryHM, GlobalDim,
                                            LocalAssemblerImplementation,
                                            LocalAssemblerInterface>(
        mesh_elements, dof_table, local_assemblers, integration_order,
        std::forward<ExtraCtorArgs>(extra_ctor_args)...);
}
 
template <int GlobalDim,
          template <typename /* shp fct */, typename /* lower order shp fct */,
                    int /* global dim */>
          class LocalAssemblerImplementation,
          typename LocalAssemblerInterface, typename... ExtraCtorArgs>
void createLocalAssemblersStokes(
    std::vector<MeshLib::Element*> const& mesh_elements,
    NumLib::LocalToGlobalIndexMap const& dof_table,
    std::vector<std::unique_ptr<LocalAssemblerInterface>>& local_assemblers,
    NumLib::IntegrationOrder const integration_order,
    ExtraCtorArgs&&... extra_ctor_args)
{
    detail::createLocalAssemblersTaylorHood<
        LocalAssemblerFactoryStokes, GlobalDim, LocalAssemblerImplementation,
        LocalAssemblerInterface>(
        mesh_elements, dof_table, local_assemblers, integration_order,
        std::forward<ExtraCtorArgs>(extra_ctor_args)...);
}
 
}  // namespace ProcessLib