SM_LinearSliceOperator.h

#ifndef SM_LinearSliceOperator_H
#define SM_LinearSliceOperator_H
 
//super class header
#include "SM_Object.h"
 
class SM_LinearSliceOperator : public virtual SM_Object {
    
    //attributes
private :
    
   
protected:
    // CONSTRUCTORS 
    SM_LinearSliceOperator(void) {
      
    }
    
    // DESTRUCTORS 
    virtual ~SM_LinearSliceOperator(void) {
    }
    
public :
 
    virtual tMemSize getMemorySize() const {
        return sizeof(*this)+this->getContentsMemorySize();
    }
 
    inline static  tReal ComputeEnergy(const SM_Network& network,
                                       const SM_Material& material,
                                       const tIndex& startIndex,
                                       const tIndex& endIndex,
                                       const tReal *S,
                                       const tInteger& nH,
                                       const tReal *H)  {
        
        tReal E=0;
        if ((S==null) || (H==null) || (nH==0)) return E;
        
        //number of particles of the slice
        tIndex s=endIndex-startIndex;
    
        //H values
        const tReal *iH=H;
        if (nH>s) iH+=startIndex*SM_Constants::DIM;
        
        //begin iterator on S
        const tReal *iS=S;iS+=startIndex*SM_Constants::DIM;
        
        //end iterator on S
        const tReal *eS=S;eS+=endIndex*SM_Constants::DIM;
    
    
        while (iS!=eS) {
            E+=functions_array::scalarProduct<tReal,SM_Constants::DIM>(iS,iH);
            iS+=SM_Constants::DIM;
            iH+=SM_Constants::DIM;
            
            
        }
          
            
        E*=-0.5;
        
        //adimensionized spin moment
        E*=material.getAtomicSpinMoment();
        
        //material.getAdimensionizedDerivativeEnergyFactor()=
        // - 1/Mu_B for dimensionized variable
        // - 1 for adimensionized variable
        //E*=Mu_B or 1
        E/=material.getAdimensionizedDerivativeEnergyFactor();
 
        return E;
    }
    
    
};
 
 
#endif