PressureGauge.h

#ifndef __RPK_PRESSURE_GAUGE_H__
#define __RPK_PRESSURE_GAUGE_H__

#include <boost/array.hpp>
#include <boost/circular_buffer.hpp>
#include <deque>

namespace MDSimulation
{
    const size_t PGAUGE_ORDER = 2;

    typedef std::pair<double, double> TimeSeries;

    class PressureGauge
    {
    public:

    #ifdef WINDOWED_PGAUGE
        typedef std::deque<TimeSeries> BufferType;
    #else
        typedef boost::circular_buffer<TimeSeries> BufferType;
    #endif

        PressureGauge();

        void AddSample(const double time, const double pressure);


        double PressureTime() const;

        double SampleGauge(const double radius);

        double ComputePressure(const double currentTime,
                               const double radius,
                               const double drdt) const;

        double ComputeDpdt(const double currentTime,
                           const double radius,
                           const double dpdt) const;

        void Serialize(FILE* file, bool save = true);

        template <typename Writer>
        void Serialize(Writer& pt) const;

        // void Read(const boost::property_tree::ptree& tree);

        bool IsReady() const;

        bool CanSample() const;

        void SetPressureStats(const double t, const double p, const double dpdt);

    private:

        // Circular buffer of pressure samples.
        BufferType Buffer;

        boost::array<double, PGAUGE_ORDER> Coefficients;

        double SampleRadius;
        double SampleTime;

        double LastPressure;
        double LastDpdt;
        double LastTime;
        double CurrentPressure;
        double CurrentDpdt;
        double CurrentTime;

        bool Ready;
    };

    template <typename Writer>
    void PressureGauge::Serialize(Writer& writer) const
    {
        writer.StartObject();

        writer.String("LastPressure")    , writer.Double(LastPressure);
        writer.String("LastDpdt")        , writer.Double(LastDpdt);
        writer.String("LastTime")        , writer.Double(LastTime);
        writer.String("CurrentPressure") , writer.Double(CurrentPressure);
        writer.String("CurrentDpdt")     , writer.Double(CurrentDpdt);
        writer.String("CurrentTime")     , writer.Double(CurrentTime);

        // Write the buffer as an array of pairs. Each element of the array is
        // an object with two values: "first" and "second".
        writer.String("Buffer");
        writer.BeginArray();
        std::for_each(Buffer.begin(), Buffer.end(), [&writer](const TimeSeries& ts)
                {
                    writer.BeginObject();
                    writer.String("first"),  writer.Double(ts.first);
                    writer.String("second"), writer.Double(ts.second);
                    writer.EndObject();
                });

        writer.EndArray();

        // End the pressure gauge object.
        writer.EndObject();
    }
}

#endif