SM_Beam.hpp

Go to the documentation of this file.

#ifndef SM_Beam_HPP
#define SM_Beam_HPP
 
#include "SMS_BoostNormalFunction.h"
#include "SMS_STDNormalFunction.h"
 
template<class SystemImpl,class NoiseImpl>
tIndex SM_Beam::templatedSimulate(SystemImpl& system,const NoiseImpl& noise) {
    
    //get the size of the beam
    const tIndex & beamSize=getBeamSize();
 
    //iterator on trajectories
    tIndex s;
    
    //the size of mu is mBeamSize x mTimeStepsNumber x mParticlesNumber x mDimension
    //the value of the k-coordinates of the particle p at time t for the simulation s is at index:
    //index=k+mDimension*(p+mParticlesNumber*(t+mTimeStepsNumber*(s)))
    
    //the mu by simulation
    std::valarray<tReal>& mu=getMagneticMoment();
 
    //the energy of simulations
    std::valarray<tReal>& Es=getEnergies();
    
    //mu Size for one trajectory
    tIndex muSize=mu.size()/beamSize;
 
    //E Size for one trajectory
    tIndex ESize=Es.size()/beamSize;
        
    //get the stochastic function
    const SM_StochasticFunction & sf=getStochasticFunction();
 
    //mu at simulation 0
    tReal *mu_s=&mu[0];
 
    //Es at simulation 0
    tReal *Es_s=&Es[0];
 
    //number of simulations without errors
    tIndex nEffectiveSimulations=0;
    
    if (sf.isInstanceOf<SMS_BoostNormalFunction>()) {
        //Boost normal function
        
        for (s=0;s<beamSize;s++) {
            if (system.makeTemplatedRelaxation(*dynamic_cast<const SMS_BoostNormalFunction*>(&sf),
                                               noise,
                                               mu_s,Es_s)) {
                nEffectiveSimulations++;
                
                //next trajectory
                mu_s+=muSize;
                Es_s+=ESize;
                
            }
           
        }   
    } else  {
        // STD normal function
        for (s=0;s<beamSize;s++) {
            if (system.makeTemplatedRelaxation(*dynamic_cast<const SMS_STDNormalFunction*>(&sf),
                                               noise,
                                               mu_s,Es_s)) {
                nEffectiveSimulations++;
                
                //next trajectory
                mu_s+=muSize;
                Es_s+=ESize;
                
            }
          
        }   
    } 
    
    //main loop
    return nEffectiveSimulations;
  
}
#endif