COPSS.h

#pragma once 

#include "Constraint.h"
#include <iostream>
#include <iomanip>
#include <boost/mpi.hpp>
#include <fstream>

namespace SSAGES
{
    class COPSS : public Constraint
    {
    private:
        std::ofstream myout_;

    public:

        COPSS(boost::mpi::communicator& comm,
                   unsigned int frequency) : 
        Constraint(frequency, comm)
        {
            myout_.open("test.out");
        }

        void PreSimulation(Snapshot*, const CVList&) override
        {
            myout_<<"In PreSimulation"<<std::endl;
        }

        void PostIntegration(Snapshot* snapshot, const CVList&) override
        {

            using std::setw;
            using std::right;
            using std::setprecision;
            

            myout_ << snapshot->GetWalkerID() 
            << " " << snapshot->GetCommunicator().rank() << std::endl;
            auto& f = snapshot->GetForces();
            f[0][0] = 1.000000*f[0][0];
        }

        void PostSimulation(Snapshot*, const CVList&) override
        {
            myout_ <<" Post simulation"<<std::endl;
        }


        void Serialize(Json::Value& /* json */) const override
        {

        }

        ~COPSS() { myout_.close(); }
    };
}