IntegrationMethodProvider.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "IntegrationMethodRegistry.h"
#include "MeshLib/Elements/Quad.h"
 
namespace NumLib
{
template <typename T>
concept IntegrationMethodProvider = requires(T t)
{
    {
        // using Quad as a representative mesh element type
        t.template getIntegrationMethod<MeshLib::Quad>(
            std::declval<MeshLib::Element>())
        } -> std::same_as<GenericIntegrationMethod const&>;
};

class DefaultIntegrationMethodProvider
{
public:
    explicit DefaultIntegrationMethodProvider(
        IntegrationOrder const integration_order)
        : integration_order_{integration_order}
    {
    }
 
    template <typename MeshElement>
    GenericIntegrationMethod const& getIntegrationMethod(
        MeshLib::Element const& /*e*/) const
    {
        return IntegrationMethodRegistry::template getIntegrationMethod<
            MeshElement>(integration_order_);
    }
 
private:
    IntegrationOrder integration_order_;
};

inline DefaultIntegrationMethodProvider getIntegrationMethodProvider(
    NumLib::IntegrationOrder const integration_order)
{
    return DefaultIntegrationMethodProvider{integration_order};
}
 
}  // namespace NumLib