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 */