SM_BeamRandomJRange.h

#ifndef SM_BeamRandomJRange_H
#define SM_BeamRandomJRange_H
 
//inherits header
#include "SM_BeamCycle.h"
 
 
class SM_BeamRandomJRange: public SM_BeamCycle<SM_BeamRandomJRange> {
   
    //attributes
private :
 
    //class type
    typedef  SM_BeamRandomJRange SelfClass;
    typedef  SM_BeamCycle<SelfClass> SuperClass;
 
    //range [eps0,eps1,deps]
    std::array<tReal,3> mRandomJRange;
    
    tString mJFileName;//J file
    tIndex mBeamSize;
    
public:
    // CONSTRUCTORS 
    SM_BeamRandomJRange(void) {
       
    }
  
    // DESTRUCTORS 
    virtual ~SM_BeamRandomJRange(void) {
    }
    
 
 
public : 
    // CREATE class
   
 
    //SET & GET methods
 
   
 
    
      //MEMORY
    
    virtual tMemSize getMemorySize() const {
        return sizeof(*this)+getContentsMemorySize();
    }
    virtual tMemSize getContentsMemorySize() const {
        tMemSize mem=SuperClass::getContentsMemorySize();
        mem+=3*SM_Constants::DIM*sizeof(tReal);//mRandomJRange
        mem+=mJFileName.length()*sizeof(tChar);
        return mem;
    }
    
    //New method
    static CORE_UniquePointer<SelfClass> New() {
        return CORE_UniquePointer<SelfClass>(new SelfClass(),
                                             SelfClass::Delete()); 
        
    }
 
public:
  
public:
    
    inline void setStochasticRange(const std::array<tReal,3>& range) {
        const tReal *iRange=range.data();
        for(auto& Tk:mRandomJRange) {
            Tk=(*iRange);
            iRange++;
        }
       
    }
 
 
    inline void setJFileName(const tString& Jfile) {
        mJFileName=Jfile;
       
    }
    
public:
 
    
    inline tBoolean open(SM_Beam& beam,SM_StochasticOutputComponent& stochasticOutput,
                         const tInteger& outputId,tIndex& iCycleStep) {
        iCycleStep=0;
  
        //system
        SM_System& system=beam.getSystem();
        
 
        //set the temperature of the system
        system.setStochasticNoise(mRandomJRange[0]);
 
        //open  file
        if (outputId==0) {
            std::ofstream file(mJFileName.c_str(),std::ios::out);
            if (!file) {
                std::cout<<"impossible to open J Range file "<<mJFileName<<"\n";
                return false;
            }
            
            
            //print header in temperature file
            printHeader(file,beam,stochasticOutput.getOutputDescription());
            
            //close the file
            file.close();
        }
      
        if (mRandomJRange[0]==0) {
            //store beam size
            mBeamSize=beam.getBeamSize();
            //only one simulation is needed
            beam.setBeamSize(1);
        }
        
        return true;
    }
    inline void  nextCycleStep(const tIndex& iCycleStep,SM_Beam& beam,const tInteger& outputId,const SM_StochasticOutputComponent& output) {
       
        //system
        SM_System& system=beam.getSystem();
 
        //S
        //const SM_RealField& S=system.getMagneticMomentDirections();
 
 
        //save the data
        if (outputId==0) {
            std::ofstream file(mJFileName.c_str(),std::ios::app);
            const std::valarray<tReal>& outputs=output.getOutputValues();
            for(const auto& r:outputs) {
                file<<r<<"\t";
            }
            file<<"\n";
            file.close();
 
        }
        
        //set the next epsilon
        tReal eps=beam.getSystem().getStochasticNoise();
        if (eps==0) {
            //restore beam size
            beam.setBeamSize(mBeamSize);
        }
        
        eps+=mRandomJRange[2];
        
        system.setStochasticNoise(eps);
        
 
    }
    
    inline tBoolean  isCycleFinished(const tIndex& iCycle,const SM_Beam& beam) {
        return (beam.getSystem().getStochasticNoise()>=mRandomJRange[1]);
    }
    inline tBoolean close(const tIndex& iCycle,const SM_Beam& beam) {
        return true;
    }
    
    virtual tString toString() const override {
        std::stringstream ret;
        ret<<SuperClass::toString()<<"\n";
        ret<<"\t Cycle temperature range:"<<functions_array::toString(mRandomJRange)<<"\n";
        return ret.str();
    }
private:
    inline void printHeader(std::ofstream& file,const SM_Beam& beam,const std::vector<tString>& outputDescription) const {
        
        file<<"# J stochastic Range cycle Beam generated by soft:"<<CORE_Run::GetSoftName()<<" Version:"<<CORE_Run::GetVersion()<<"\n";
        tString str="J stochastic  Range :"+this->toString();
        functions_string::replaceAll("\n","\n#",str);
        file<<"#"<<str<<"\n";
        
        str="beam:"+beam.toString();
        functions_string::replaceAll("\n","\n#",str);
        file<<"#"<<str<<"\n";
 
        file<<"#";
        for(const auto& d:outputDescription) file<<"\t"<<d;
        file<<"\n";
    }
 
   
 
};
 
 
#endif