wrapper.hpp

Go to the documentation of this file.

#ifndef OPM_VERTEQ_WRAPPER_HPP_INCLUDED
#define OPM_VERTEQ_WRAPPER_HPP_INCLUDED
 
// Copyright (C) 2013 Uni Research AS
// This file is licensed under the GNU General Public License v3.0
 
#include <vector>
#include <memory> // unique_ptr
 
#ifndef OPM_VERTEQ_VISIBILITY_HPP_INCLUDED
#include <opm/verteq/visibility.hpp>
#endif /* OPM_VERTEQ_VISIBILITY_HPP_INCLUDED */
 
#ifndef OPM_VERTEQ_OPMFWD_HPP_INCLUDED
#include <opm/verteq/opmfwd.hpp>
#endif /* OPM_VERTEQ_OPMFWD_HPP_INCLUDED */
 
#ifndef OPM_VERTEQ_SIMULATOR_HPP_INCLUDED
#include <opm/verteq/simulator.hpp>
#endif /* OPM_VERTEQ_SIMULATOR_HPP_INCLUDED */
 
namespace Opm {
 
// forward declaration
class VertEq;
 
struct OPM_VERTEQ_PUBLIC VertEqWrapperBase {
protected:
    VertEqWrapperBase (
        std::unique_ptr <Simulator> underlaying,
        const parameter::ParameterGroup& param,
        const UnstructuredGrid& grid,
        const IncompPropertiesInterface& props,
        const RockCompressibility* rock_comp_props,
        WellsManager& wells_manager,
        const std::vector<double>& src,
        const FlowBoundaryConditions* bcs,
        LinearSolverInterface& linsolver,
        const double* gravity);
public:
    virtual ~VertEqWrapperBase ();
 
  virtual SimulatorReport run(
        SimulatorTimer& timer,
        TwophaseState& state,
        WellState& well_state);
 
    Event& timestep_completed ();
 
    void sync ();
 
private:
    // underlaying simulator to use for 2D
    std::unique_ptr <Simulator> sim;
 
    // vertical equilibrium model
    VertEq* ve;
 
    // list of translated wells
    WellsManager* wells_mgr;
 
    // list of handlers for timestep notifications
    std::unique_ptr <EventSource> timestep_callbacks;
 
    // current state of any simulation we are doing
    TwophaseState* fineState;
    TwophaseState* coarseState;
 
    // flag that determines whether we have synced or not
    bool syncDone;
    void resetSyncFlag ();
};
 
template <typename SimulatorType>
struct OPM_VERTEQ_PUBLIC VertEqWrapper : public VertEqWrapperBase {
    VertEqWrapper (const parameter::ParameterGroup& param,
                   const UnstructuredGrid& grid,
                   const IncompPropertiesInterface& props,
                   const RockCompressibility* rock_comp_props,
                   WellsManager& wells_manager,
                   const std::vector<double>& src,
                   const FlowBoundaryConditions* bcs,
                   LinearSolverInterface& linsolver,
                   const double* gravity)
        : VertEqWrapperBase (std::unique_ptr <Simulator> (
                               new SimulatorAdapter <SimulatorType> ()),
                             param, grid, props, rock_comp_props, wells_manager,
                             src, bcs, linsolver, gravity) {
    }
};
 
} /* namespace Opm */
 
#endif /* OPM_VERTEQ_WRAPPER_HPP_INCLUDED */