unstructuredgridvanguard.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:
/*
  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 EWOMS_UNSTRUCTURED_GRID_VANGUARD_HH
#define EWOMS_UNSTRUCTURED_GRID_VANGUARD_HH
 
#include <dune/grid/io/file/dgfparser/gridptr.hh>
 
#include <opm/models/io/basevanguard.hh>
#include <opm/models/utils/parametersystem.hpp>
#include <opm/models/utils/propertysystem.hh>
 
#ifdef HAVE_OPM_GRID
#include <opm/grid/UnstructuredGrid.h>
#include <string>
#endif
 
namespace Opm {

template <class TypeTag>
class UnstructuredGridVanguard : public BaseVanguard<TypeTag>
{
    using ParentType = BaseVanguard<TypeTag>;
    using Scalar = GetPropType<TypeTag, Properties::Scalar>;
    using Simulator = GetPropType<TypeTag, Properties::Simulator>;
    using Grid = GetPropType<TypeTag, Properties::Grid>;
 
    using GridPointer = Dune::GridPtr<Grid>;
 
public:
    static void registerParameters() {
        Parameters::Register<Parameters::GridGlobalRefinements>
            ("The number of global refinements of the grid "
             "executed after it was loaded");
        Parameters::Register<Parameters::GridFile>
            ("The file name of the file to load");
    }

    explicit UnstructuredGridVanguard(Simulator& simulator)
        : ParentType(simulator)
    {
#ifdef HAVE_OPM_GRID
        const std::string gridFileName = Parameters::Get<Parameters::GridFile>();
        const int numRefinments = Parameters::Get<Parameters::GridGlobalRefinements>();
 
        UnstructuredGrid* ugrid = read_grid(gridFileName.c_str());
        if (ugrid == nullptr) {
            throw std::runtime_error("RuntimeError: UnstructuredGridVanguard could not read grid file: " +
                                     gridFileName + ". Are you sure the filename is correct?");
        }
        ugPtr_.reset(std::move(ugrid));
        //GridPointer polygrid(new Grid(*ugPtr));
        gridPtr_ = new Grid(*ugPtr_);//std::move(polygrid);
        if (numRefinments > 0) {
            gridPtr_->globalRefine(numRefinments);
        }
        this->finalizeInit_();
#endif
    }

    Grid& grid()
    { return *gridPtr_; }

    const Grid& grid() const
    { return *gridPtr_; }
 
private:
    GridPointer gridPtr_;
    typename Grid::UnstructuredGridPtr ugPtr_;
};
 
}  // namespace Opm
 
#endif