AdaptiveSimulatorTimer.hpp

/*
  Copyright 2014 IRIS AS

  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 3 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/>.
*/
#ifndef OPM_ADAPTIVESIMULATORTIMER_HEADER_INCLUDED
#define OPM_ADAPTIVESIMULATORTIMER_HEADER_INCLUDED

#include <cassert>
#include <iostream>
#include <vector>

#include <algorithm>
#include <numeric>

#include <opm/simulators/timestepping/SimulatorTimerInterface.hpp>

namespace Opm
{

    class AdaptiveSimulatorTimer : public SimulatorTimerInterface
    {
    public:
        AdaptiveSimulatorTimer( const SimulatorTimerInterface& timer,
                                const double lastStepTaken,
                                const double maxTimeStep = std::numeric_limits<double>::max() );

        AdaptiveSimulatorTimer& operator++ ();

        void advance() { this->operator++ (); }

        void provideTimeStepEstimate( const double dt_estimate );

        bool initialStep () const;

        int currentStepNum () const;

        int reportStepNum() const;

        double currentStepLength () const;

        double totalTime() const;

        double simulationTimeElapsed() const;

        bool done () const;

        double averageStepLength() const;

        double maxStepLength () const;

        double minStepLength () const;

        double stepLengthTaken () const;

        void report(std::ostream& os) const;

        boost::posix_time::ptime startDateTime() const;

        bool lastStepFailed() const {return lastStepFailed_;}

        void setLastStepFailed(bool lastStepFailed) {lastStepFailed_ = lastStepFailed;}

        virtual std::unique_ptr< SimulatorTimerInterface > clone() const;

    protected:
        const boost::posix_time::ptime start_date_time_;
        const double start_time_;
        const double total_time_;
        const int report_step_;
        const double max_time_step_;

        double current_time_;
        double dt_;
        int current_step_;

        std::vector< double > steps_;
        bool lastStepFailed_;

    };

} // namespace Opm

#endif // OPM_SIMULATORTIMER_HEADER_INCLUDED