simulator.hpp

Go to the documentation of this file.

#ifndef OPM_VERTEQ_SIMULATOR_HPP_INCLUDED
#define OPM_VERTEQ_SIMULATOR_HPP_INCLUDED

// Copyright (C) 2013 Uni Research AS
// This file is licensed under the GNU General Public License v3.0

#include <memory> // unique_ptr
#include <vector>

#ifndef OPM_VERTEQ_VISIBILITY_HPP_INCLUDED
#include <opm/verteq/visibility.hpp>
#endif /* OPM_VERTEQ_VISIBILITY_HPP_INCLUDED */

#ifndef OPM_SIMULATORREPORT_HEADER_INCLUDED
#include <opm/core/simulator/SimulatorReport.hpp>
#endif /* OPM_SIMULATORREPORT_HEADER_INCLUDED */

#ifndef OPM_VERTEQ_OPMFWD_HPP_INCLUDED
#include <opm/verteq/opmfwd.hpp>
#endif /* OPM_VERTEQ_OPMFWD_HPP_INCLUDED */

namespace Opm {

struct OPM_VERTEQ_PUBLIC Simulator {
    virtual void init (
        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) = 0;

    virtual ~Simulator () {}

  virtual SimulatorReport run(
        SimulatorTimer& timer,
        TwophaseState& state,
        WellState& well_state) = 0;

    virtual Event& timestep_completed () = 0;

    virtual void sync () = 0;
};

template <typename T>
struct OPM_VERTEQ_PUBLIC SimulatorAdapter : public Simulator {
    virtual void init (
        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) {
        t_ = std::unique_ptr <T> (new T (param, grid, props, rock_comp_props,
                    wells_manager, src, bcs, linsolver,
                    gravity));
    }

    // forward all method to the underlaying instance
    virtual SimulatorReport run(
        SimulatorTimer& timer,
        TwophaseState& state,
        WellState& well_state) {
        return t_->run (timer, state, well_state);
    }
    virtual Event& timestep_completed () {
        return t_->timestep_completed ();
    }
    virtual void sync () {
        t_->sync ();
    }

private:
    std::unique_ptr <T> t_;
};

} /* namespace Opm */

#endif /* OPM_VERTEQ_SIMULATOR_HPP_INCLUDED */