discretization/cellcentered/mpfa/gridfluxvariablescache.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_DISCRETIZATION_CCMPFA_GRID_FLUXVARSCACHE_HH
#define DUMUX_DISCRETIZATION_CCMPFA_GRID_FLUXVARSCACHE_HH
 
// make the local view function available whenever we use this class
#include <dumux/discretization/localview.hh>
#include <dumux/discretization/cellcentered/mpfa/elementfluxvariablescache.hh>
 
namespace Dumux {
 
template<class ModelTraits>
struct IvDataHandlePhysicsTraits
{
    static constexpr bool enableAdvection = ModelTraits::enableAdvection();
    static constexpr bool enableMolecularDiffusion = ModelTraits::enableMolecularDiffusion();
    static constexpr bool enableHeatConduction = ModelTraits::enableEnergyBalance();
 
    static constexpr int numPhases = ModelTraits::numFluidPhases();
    static constexpr int numComponents = ModelTraits::numFluidComponents();
};
 
template<class P,
         class FVC, class FVCF,
         class PIV, class SIV,
         class PDH, class SDH>
struct CCMpfaDefaultGridFluxVariablesCacheTraits
{
    using Problem = P;
    using FluxVariablesCache = FVC;
    using FluxVariablesCacheFiller = FVCF;
 
    using PrimaryInteractionVolume = PIV;
    using SecondaryInteractionVolume = SIV;
    using PrimaryIvDataHandle = PDH;
    using SecondaryIvDataHandle = SDH;
 
    template<class GridFluxVariablesCache, bool cachingEnabled>
    using LocalView = CCMpfaElementFluxVariablesCache<GridFluxVariablesCache, cachingEnabled>;
 
    // Reserve memory (over-) estimate for interaction volumes and corresponding data.
    // The overestimate doesn't hurt as we are not in a memory-limited configuration.
    // We need to avoid reallocation because in the caches we store pointers to the data handles.
    // Default -> each facet has two neighbors (local adaption) and all scvfs belongs to different ivs.
    // If you want to use higher local differences change the parameter below.
    static constexpr std::size_t maxLocalElementLevelDifference()
    { return 2; };
};
 
template<class Traits, bool cachingEnabled>
class CCMpfaGridFluxVariablesCache;
 
template<class TheTraits>
class CCMpfaGridFluxVariablesCache<TheTraits, true>
{
    using Problem = typename TheTraits::Problem;
    using ThisType = CCMpfaGridFluxVariablesCache<TheTraits, true>;
 
    using FluxVariablesCacheFiller = typename TheTraits::FluxVariablesCacheFiller;
public:
    using Traits = TheTraits;
 
    using PrimaryInteractionVolume = typename Traits::PrimaryInteractionVolume;
    using SecondaryInteractionVolume = typename Traits::SecondaryInteractionVolume;
 
    using PrimaryIvDataHandle = typename Traits::PrimaryIvDataHandle;
    using SecondaryIvDataHandle = typename Traits::SecondaryIvDataHandle;
 
    using FluxVariablesCache = typename Traits::FluxVariablesCache;
 
    static constexpr bool cachingEnabled = true;
 
    using LocalView = typename Traits::template LocalView<ThisType, cachingEnabled>;
 
    CCMpfaGridFluxVariablesCache(const Problem& problem)
    : problemPtr_(&problem)
    {}
 
    template<class GridGeometry, class GridVolumeVariables, class SolutionVector>
    void update(const GridGeometry& gridGeometry,
                const GridVolumeVariables& gridVolVars,
                const SolutionVector& sol,
                bool forceUpdate = false)
    {
        // Update only if the filler puts solution-dependent
        // stuff into the caches or if update is enforced
        if (FluxVariablesCacheFiller::isSolDependent || forceUpdate)
        {
            // clear previous data if forced update is desired
            if (forceUpdate)
            {
                clear_();
 
                const auto& gridIvIndexSets = gridGeometry.gridInteractionVolumeIndexSets();
                const auto numPrimaryIvs = gridIvIndexSets.numPrimaryInteractionVolumes();
                const auto numSecondaryIVs = gridIvIndexSets.numSecondaryInteractionVolumes();
                ivDataStorage_.primaryInteractionVolumes.reserve(numPrimaryIvs);
                ivDataStorage_.secondaryInteractionVolumes.reserve(numSecondaryIVs);
                ivDataStorage_.primaryDataHandles.reserve(numPrimaryIvs);
                ivDataStorage_.secondaryDataHandles.reserve(numSecondaryIVs);
 
                // reserve memory estimate for caches, interaction volumes and corresponding data
                fluxVarsCache_.resize(gridGeometry.numScvf());
            }
 
            // instantiate helper class to fill the caches
            FluxVariablesCacheFiller filler(problem());
 
            // set all the caches to "outdated"
            for (auto& cache : fluxVarsCache_)
                cache.setUpdateStatus(false);
 
            auto fvGeometry = localView(gridGeometry);
            auto elemVolVars = localView(gridVolVars);
            for (const auto& element : elements(gridGeometry.gridView()))
            {
                fvGeometry.bind(element);
                elemVolVars.bind(element, fvGeometry, sol);
 
                // Prepare all caches of the scvfs inside the corresponding interaction volume. Skip
                // those ivs that are touching a boundary, we only store the data on interior ivs here.
                for (const auto& scvf : scvfs(fvGeometry))
                    if (!isEmbeddedInBoundaryIV_(scvf, gridGeometry) && !fluxVarsCache_[scvf.index()].isUpdated())
                        filler.fill(*this, fluxVarsCache_[scvf.index()], ivDataStorage_, fvGeometry, elemVolVars, scvf, forceUpdate);
            }
        }
    }
 
    template<class FVElementGeometry, class ElementVolumeVariables>
    void updateElement(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
                       const FVElementGeometry& fvGeometry,
                       const ElementVolumeVariables& elemVolVars)
    {
        // Update only if the filler puts
        // solution-dependent stuff into the caches
        if (FluxVariablesCacheFiller::isSolDependent)
        {
            const auto& gridGeometry = fvGeometry.gridGeometry();
            const auto& assemblyMapI = gridGeometry.connectivityMap()[gridGeometry.elementMapper().index(element)];
 
            // helper class to fill flux variables caches
            FluxVariablesCacheFiller filler(problem());
 
            // first, set all the caches to "outdated"
            for (const auto& scvf : scvfs(fvGeometry))
                fluxVarsCache_[scvf.index()].setUpdateStatus(false);
            for (const auto& dataJ : assemblyMapI)
                for (const auto scvfIdx : dataJ.scvfsJ)
                    fluxVarsCache_[scvfIdx].setUpdateStatus(false);
 
            // go through the caches maybe update them
            for (const auto& scvf : scvfs(fvGeometry))
            {
                auto& scvfCache = fluxVarsCache_[scvf.index()];
                if (!isEmbeddedInBoundaryIV_(scvf, gridGeometry) && !scvfCache.isUpdated())
                    filler.fill(*this, scvfCache, ivDataStorage_, fvGeometry, elemVolVars, scvf);
            }
 
            for (const auto& dataJ : assemblyMapI)
            {
                for (const auto scvfIdx : dataJ.scvfsJ)
                {
                    auto& scvfCache = fluxVarsCache_[scvfIdx];
                    const auto& scvf = fvGeometry.scvf(scvfIdx);
                    if (!isEmbeddedInBoundaryIV_(scvf, gridGeometry) && !scvfCache.isUpdated())
                        filler.fill(*this, scvfCache, ivDataStorage_, fvGeometry, elemVolVars, scvf);
                }
            }
        }
    }
 
    template<class SubControlVolumeFace>
    const FluxVariablesCache& operator [](const SubControlVolumeFace& scvf) const
    { return fluxVarsCache_[scvf.index()]; }
 
    template<class SubControlVolumeFace>
    FluxVariablesCache& operator [](const SubControlVolumeFace& scvf)
    { return fluxVarsCache_[scvf.index()]; }
 
    template<class SubControlVolumeFace>
    const PrimaryInteractionVolume& primaryInteractionVolume(const SubControlVolumeFace& scvf) const
    { return ivDataStorage_.primaryInteractionVolumes[ (*this)[scvf].ivIndexInContainer() ]; }
 
    template<class SubControlVolumeFace>
    const PrimaryIvDataHandle& primaryDataHandle(const SubControlVolumeFace& scvf) const
    { return ivDataStorage_.primaryDataHandles[ (*this)[scvf].ivIndexInContainer() ]; }
 
    template<class SubControlVolumeFace>
    const SecondaryInteractionVolume& secondaryInteractionVolume(const SubControlVolumeFace& scvf) const
    { return ivDataStorage_.secondaryInteractionVolumes[ (*this)[scvf].ivIndexInContainer() ]; }
 
    template<class SubControlVolumeFace>
    const SecondaryIvDataHandle& secondaryDataHandle(const SubControlVolumeFace& scvf) const
    { return ivDataStorage_.secondaryDataHandles[ (*this)[scvf].ivIndexInContainer() ]; }
 
    const Problem& problem() const
    { return *problemPtr_; }
 
private:
    template<class SubControlVolumeFace, class GridGeometry>
    bool isEmbeddedInBoundaryIV_(const SubControlVolumeFace& scvf, const GridGeometry& gridGeometry) const
    {
        const auto& gridIvIndexSets = gridGeometry.gridInteractionVolumeIndexSets();
        if (gridGeometry.vertexUsesSecondaryInteractionVolume(scvf.vertexIndex()))
            return gridIvIndexSets.secondaryIndexSet(scvf).nodalIndexSet().numBoundaryScvfs() > 0;
        else
            return gridIvIndexSets.primaryIndexSet(scvf).nodalIndexSet().numBoundaryScvfs() > 0;
    }
 
    void clear_()
    {
        fluxVarsCache_.clear();
        ivDataStorage_.primaryInteractionVolumes.clear();
        ivDataStorage_.secondaryInteractionVolumes.clear();
        ivDataStorage_.primaryDataHandles.clear();
        ivDataStorage_.secondaryDataHandles.clear();
    }
 
    const Problem* problemPtr_;
    std::vector<FluxVariablesCache> fluxVarsCache_;
 
    // stored interaction volumes and handles
    using IVDataStorage = InteractionVolumeDataStorage<PrimaryInteractionVolume,
                                                       PrimaryIvDataHandle,
                                                       SecondaryInteractionVolume,
                                                       SecondaryIvDataHandle>;
    IVDataStorage ivDataStorage_;
};
 
template<class TheTraits>
class CCMpfaGridFluxVariablesCache<TheTraits, false>
{
    using Problem = typename TheTraits::Problem;
    using ThisType = CCMpfaGridFluxVariablesCache<TheTraits, false>;
 
    using FluxVariablesCacheFiller = typename TheTraits::FluxVariablesCacheFiller;
public:
    using Traits = TheTraits;
 
    using PrimaryInteractionVolume = typename Traits::PrimaryInteractionVolume;
    using SecondaryInteractionVolume = typename Traits::SecondaryInteractionVolume;
 
    using PrimaryIvDataHandle = typename Traits::PrimaryIvDataHandle;
    using SecondaryIvDataHandle = typename Traits::SecondaryIvDataHandle;
 
    using FluxVariablesCache = typename Traits::FluxVariablesCache;
 
    static constexpr bool cachingEnabled = false;
 
    using LocalView = typename Traits::template LocalView<ThisType, cachingEnabled>;
 
    CCMpfaGridFluxVariablesCache(const Problem& problem) : problemPtr_(&problem) {}
 
    template<class GridGeometry, class GridVolumeVariables, class SolutionVector>
    void update(const GridGeometry& gridGeometry,
                const GridVolumeVariables& gridVolVars,
                const SolutionVector& sol,
                bool forceUpdate = false) {}
 
    template<class FVElementGeometry, class ElementVolumeVariables>
    void updateElement(const typename FVElementGeometry::GridGeometry::GridView::template Codim<0>::Entity& element,
                       const FVElementGeometry& fvGeometry,
                       const ElementVolumeVariables& elemVolVars) {}
 
    const Problem& problem() const
    { return *problemPtr_; }
 
private:
    const Problem* problemPtr_;
};
 
} // end namespace Dumux
 
#endif