GenericCpGridVanguard.hpp

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:
/*
  This file is part of the Open Porous Media project (OPM).
 
  OPM is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 2 of the License, or
  (at your option) any later version.
 
  OPM is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.
 
  You should have received a copy of the GNU General Public License
  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
 
  Consult the COPYING file in the top-level source directory of this
  module for the precise wording of the license and the list of
  copyright holders.
*/
#ifndef OPM_GENERIC_CPGRID_VANGUARD_HPP
#define OPM_GENERIC_CPGRID_VANGUARD_HPP
 
#include <opm/grid/CpGrid.hpp>
 
#include <opm/simulators/flow/FlowGenericVanguard.hpp>
 
#include <functional>
#include <memory>
#include <optional>
#include <vector>
 
#if HAVE_MPI
#include <filesystem>
#endif
 
namespace Opm {
    class EclipseState;
    class Schedule;
    class Well;
    class ParallelEclipseState;
    class LgrCollection;
}
 
namespace Opm {
 
#if HAVE_MPI
namespace details {
class MPIPartitionFromFile
{
public:
    explicit MPIPartitionFromFile(const std::filesystem::path& partitionFile)
        : partitionFile_(partitionFile)
    {}
 
    std::vector<int> operator()(const Dune::CpGrid& grid) const;
 
private:
    std::filesystem::path partitionFile_{};
};
 
} // namespace Opm::details
#endif // HAVE_MPI
 
 
extern std::optional<std::function<std::vector<int> (const Dune::CpGrid&)>> externalLoadBalancer;
 
template<class ElementMapper, class GridView, class Scalar>
class GenericCpGridVanguard {
protected:
    using CartesianIndexMapper = Dune::CartesianIndexMapper<Dune::CpGrid>;
    using Element = typename GridView::template Codim<0>::Entity;
 
public:
    GenericCpGridVanguard();
 
    virtual ~GenericCpGridVanguard() = default;
 
    Dune::CpGrid& grid()
    { return *grid_; }
 
    const Dune::CpGrid& grid() const
    { return *grid_; }
 
    const Dune::CpGrid& equilGrid() const;
 
    void releaseEquilGrid();
 
    static void setExternalLoadBalancer(const std::function<std::vector<int> (const Dune::CpGrid&)>& loadBalancer)
    {
        externalLoadBalancer = loadBalancer;
    }
 
    const CartesianIndexMapper& cartesianIndexMapper() const;
 
    const CartesianIndexMapper& equilCartesianIndexMapper() const;
 
    const std::vector<int>& cellPartition() const
    {
        return this->cell_part_;
    }
 
protected:
#if HAVE_MPI
    void doLoadBalance_(const Dune::EdgeWeightMethod             edgeWeightsMethod,
                        const bool                               ownersFirst,
                        const bool                               serialPartitioning,
                        const bool                               enableDistributedWells,
                        const double                             zoltanImbalanceTol,
                        const GridView&                          gridView,
                        const Schedule&                          schedule,
                        EclipseState&                            eclState,
                        FlowGenericVanguard::ParallelWellStruct& parallelWells,
                        const int                                numJacobiBlocks);
 
    void distributeFieldProps_(EclipseState& eclState);
 
private:
    std::vector<double> extractFaceTrans(const GridView& gridView) const;
 
    void distributeGrid(const Dune::EdgeWeightMethod             edgeWeightsMethod,
                        const bool                               ownersFirst,
                        const bool                               serialPartitioning,
                        const bool                               enableDistributedWells,
                        const double                             zoltanImbalanceTol,
                        const bool                               loadBalancerSet,
                        const std::vector<double>&               faceTrans,
                        const std::vector<Well>&                 wells,
                        EclipseState&                            eclState,
                        FlowGenericVanguard::ParallelWellStruct& parallelWells);
 
    void distributeGrid(const Dune::EdgeWeightMethod             edgeWeightsMethod,
                        const bool                               ownersFirst,
                        const bool                               serialPartitioning,
                        const bool                               enableDistributedWells,
                        const double                             zoltanImbalanceTol,
                        const bool                               loadBalancerSet,
                        const std::vector<double>&               faceTrans,
                        const std::vector<Well>&                 wells,
                        ParallelEclipseState*                    eclState,
                        FlowGenericVanguard::ParallelWellStruct& parallelWells);
 
protected:
    virtual const std::string& zoltanParams() const = 0;
 
#endif  // HAVE_MPI
 
    void allocCartMapper();
 
    void doCreateGrids_(EclipseState& eclState);
    void addLgrsUpdateLeafView(const LgrCollection& lgrCollection, const int lgrsSize);
 
    virtual void allocTrans() = 0;
    virtual double getTransmissibility(unsigned I, unsigned J) const = 0;
 
    // removing some connection located in inactive grid cells
    void doFilterConnections_(Schedule& schedule);
 
    Scalar computeCellThickness(const Element& element) const;
 
    std::unique_ptr<Dune::CpGrid> grid_;
    std::unique_ptr<Dune::CpGrid> equilGrid_;
    std::unique_ptr<CartesianIndexMapper> cartesianIndexMapper_;
    std::unique_ptr<CartesianIndexMapper> equilCartesianIndexMapper_;
 
    int mpiRank;
    std::vector<int> cell_part_{};
};
 
} // namespace Opm
 
#endif // OPM_GENERIC_CPGRID_VANGUARD_HPP