cvfe/interpolationpointdata.hh

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
//
// SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
// SPDX-License-Identifier: GPL-3.0-or-later
//
#ifndef DUMUX_CVFE_IP_DATA_HH
#define DUMUX_CVFE_IP_DATA_HH
 
#include <type_traits>
 
namespace Dumux::CVFE {
 
template<class LocalPos, class GlobalPos>
class InterpolationPointData
{
public:
    using LocalPosition = LocalPos;
    using GlobalPosition = GlobalPos;
 
    InterpolationPointData(LocalPosition&& localPos, GlobalPosition&& pos)
    : local_(std::move(localPos)), global_(std::move(pos)) {}
    InterpolationPointData(const LocalPosition& localPos, const GlobalPosition& pos)
    :  local_(localPos), global_(pos) {}
 
    const GlobalPosition& global() const
    { return global_; }
 
    const LocalPosition& local() const
    { return local_; }
 
private:
    LocalPosition local_;
    GlobalPosition global_;
};
 
template<class LocalPosition, class GlobalPosition, class LocalIndex>
class LocalDofInterpolationPointData : public InterpolationPointData<LocalPosition, GlobalPosition>
{
    using ParentType = InterpolationPointData<LocalPosition, GlobalPosition>;
public:
    LocalDofInterpolationPointData(LocalPosition&& localPos, GlobalPosition&& pos, LocalIndex index)
    : ParentType(localPos, pos), localDofIndex_(index) {}
    LocalDofInterpolationPointData(const LocalPosition& localPos, const GlobalPosition& pos, LocalIndex index)
    : ParentType(localPos, pos), localDofIndex_(index) {}
 
    LocalIndex localDofIndex() const
    { return localDofIndex_; }
 
private:
    LocalIndex localDofIndex_;
};
 
template<class LocalMapping, class GlobalPos>
class InterpolationPointDataLocalMapping
{
public:
    using LocalPosition = std::invoke_result_t<LocalMapping, const GlobalPos&>;
    using GlobalPosition = GlobalPos;
 
    InterpolationPointDataLocalMapping(LocalMapping&& mapping, GlobalPosition&& pos) : localMapping_(std::move(mapping)), global_(std::move(pos)) {}
    InterpolationPointDataLocalMapping(LocalMapping&& mapping, const GlobalPosition& pos) : localMapping_(std::move(mapping)), global_(pos) {}
 
    const GlobalPosition& global() const
    { return global_; }
 
    const LocalPosition local() const
    { return localMapping_(global_); }
 
private:
    LocalMapping localMapping_;
    GlobalPosition global_;
};
 
template<class BaseClass, class LocalIndex>
class FaceInterpolationPointData : public BaseClass
{
public:
    using GlobalPosition = std::remove_cvref_t<decltype(std::declval<BaseClass>().global())>;
    using LocalPosition = std::remove_cvref_t<decltype(std::declval<BaseClass>().local())>;
 
    template<class... Args>
    FaceInterpolationPointData(GlobalPosition&& n, LocalIndex index, Args&&... args)
    : BaseClass(std::forward<Args>(args)...), normal_(std::move(n)), scvfIndex_(index) {}
 
    template<class... Args>
    FaceInterpolationPointData(const GlobalPosition& n, LocalIndex index, Args&&... args)
    : BaseClass(std::forward<Args>(args)...), normal_(n), scvfIndex_(index) {}
 
    const GlobalPosition& unitOuterNormal() const
    { return normal_; }
 
    LocalIndex scvfIndex() const
    { return scvfIndex_; }
 
private:
    GlobalPosition normal_;
    LocalIndex scvfIndex_;
};
 
template<class IpData>
class IndexedQuadratureInterpolationPointData : public IpData
{
public:
    template<class... Args>
    IndexedQuadratureInterpolationPointData(std::size_t qpIdx, Args&&... args)
    : IpData(std::forward<Args>(args)...), qpIndex_(qpIdx) {}
 
    std::size_t qpIndex() const
    { return qpIndex_; }
 
private:
    std::size_t qpIndex_;
};
 
} // end namespace Dumux::CVFE
 
#endif