SM_LandauLifschitzBeam.hpp

#ifndef SM_LandauLifschitzBeam_HPP
#define SM_LandauLifschitzBeam_HPP
 
template<class StochOutputImplement>
tBoolean SM_LandauLifschitzBeam::stochasticJRun(SM_StochasticOutput<StochOutputImplement>& outputSData,
                                                const std::function<void(const tReal&,const tReal& ,const tReal&,tReal& )>& F) {
    
    tString routineName="SM_Beam::stochasticJRun<O>,F";
    CORE_Profiler::StartCallRoutine(routineName);
    
    tBoolean ok=true;
    
    //get the system
    SM_LandauLifschitzSystem& system=*dynamic_cast<SM_LandauLifschitzSystem*>(&getSystem());
    
    //open the stochastic at the beginning of the simulations running of the beam
    if (!outputSData.open(*this)) {
        std::cout<<"stochastic output "<<outputSData.getIdentityString()<<" not compatible with the system \n";
        return false;
    }
    
    
    tIndex s0,s1;
    getBeamSize(s0,s1);
    
    //get the stochastic function
    SM_StochasticFunctionsInterface& sf=getStochasticFunctions()(0);
    
    for (tIndex s=s0;s<s1;s++) {
        
       
        
        //St=S0
        if (!system.initializeMagneticMomentDirections()) {
            throw CORE_Exception("stochmagnet/landauLifschitz",
                                 "SM_LandauLifschitzBeam::stochasticJRun(SM_StochasticOutput<StochOutputImplement>& outputSData,F)",
                                 "initialization of S(t=0) failed");
        }
        //open the stochastic output the simulation  0 of the system
        outputSData.open(s,system);
        
        //running the system until the equilibrium
        ok=system.stochasticJRun(getPreconditioningStepsNumber(),sf,F);
        
        //running the system   during loop steppers number
        ok= ok && system.stochasticJRun(getStepsNumber(),sf,outputSData,F);
        
        //close the outputSData for simulation s at the end of simulations
        outputSData.close(s,system,ok);
    }  
    
     
    //close the stochastic output at the end of the runnings  of the simulations
    outputSData.close(*this);
    
    CORE_Profiler::EndCallRoutine(routineName);
    
    return true;
}
 
#endif