EMM_4Tensors.h

Go to the documentation of this file.

#ifndef EMM_4Tensors_H
#define EMM_4Tensors_H

#include "EMM_Tensors.h"

#include "EMM_2PackedSymmetricTensors.h"

DEFINE_SPTR(EMM_4Tensors);
class EMM_4Tensors : public  EMM_Tensors {
  
  
    SP_OBJECT(EMM_4Tensors);
  
    // ATTRIBUTES
  
  
public:
  
    
private:

   
    
public:
    // METHODS
  
    // CONSTRUCTORS 
     
    EMM_4Tensors();
    
   
   
   
    
    
    // DESTRUCTORS 
    
    
    virtual ~EMM_4Tensors(void);
   
public:
    
    // CONSTRUCTORS
  
   

  
public:
   
   
    // CONSTRUCTORS
    inline static SP::EMM_4Tensors New() {
        SP::EMM_4Tensors ret(new EMM_4Tensors(),
                             EMM_4Tensors::Delete()); 
        ret->setThis(ret);
        return ret;
    };
    
    virtual SP::EMM_Tensors NewInstance() const {
        return New();
    }
   
public:
    //OPERATORS
    virtual const tReal& operator()(const tUIndex& index,const tUSInt& p,const tUSInt& q,const tUSInt& r,const tUSInt& s) const {
        return (*this)[index][getIndex(p,q,r,s,getDimension())];
    }

    virtual tReal& operator()(const tUIndex& index,const tUSInt& p,const tUSInt& q,const tUSInt& r,const tUSInt& s) {
        return (*this)[index][getIndex(p,q,r,s,getDimension())];
    }
    
    //SET methods
    
protected:
   

public:

    inline void setValue(const tUIndex& index,
                         const tUSInt& i,
                         const tUSInt& j,
                         const tUSInt& k,
                         const tUSInt& l,
                         const tReal& v) {
        (*this)(index,i,j,k,l)=v;
    }
    
    
    // GET methods 
public:
    inline const tReal& getValue(const tUIndex& index,
                                 const tUSInt& i,
                                 const tUSInt& j,
                                 const tUSInt& k,
                                 const tUSInt& l) const {
        return (*this)(index,i,j,k,l);
    }
    
    
  
    virtual tUSInt getTensorSize() const {
        tUSInt n=1;
        const tUSInt& dim=getDimension();
        for (tUSInt k=0;k<4;k++) n*=dim;
        return n;
    }

    inline  virtual tUSInt getIndex(const tUSInt& p,const tUSInt& q,const tUSInt& r,const tUSInt& s,const tUSInt& dim) const {
        return s+dim*(r+dim*(q+dim*p));
    }
   
    //ALGEBRIC Methods
    //=================


    //vector methods
    //===============
    
    virtual void doubleDotCrossProduct(const tUIndex& n,
                                       const tUSInt& dim,
                                       const tReal& beta,
                                       const tUSInt& m,
                                       const tBoolean& incS,
                                       const tReal* S,
                                       const tBoolean& incU,
                                       const tReal *U,
                                       const tBoolean& incV,
                                       const tReal *V,
                                       const tBoolean& incW,
                                       const tReal *W,
                                       const tReal& alpha,
                                       tReal *Z) const;

    virtual void doubleDotProduct(const tUIndex& n,
                                  const tUSInt& dim,
                                  const tReal& beta,
                                  const tUSInt& m,
                                  const tBoolean& incS,
                                  const tReal* S,
                                  const EMM_2PackedSymmetricTensors& phi,
                                  const tBoolean& incU,
                                  const tReal *U,
                                  const tReal& alpha,
                                  tReal *Z) const;

    //scalar methods
    //==============
    
    virtual tReal squaredDoubleDotScalar(const EMM_2PackedSymmetricTensors& psi,
                                         const EMM_2PackedSymmetricTensors& phi) const;

    virtual tReal doubleDotCrossDoubleDotScalar(const tUIndex& n,
                                                const tUSInt& dim,
                                                const tUSInt& p,
                                                const tBoolean& incS,
                                                const tReal* S,
                                                const tBoolean& incU,
                                                const tReal *U,
                                                const tBoolean& incV,
                                                const tReal *V,
                                                const EMM_2PackedSymmetricTensors& phi) const;

    virtual tReal squaredDoubleDotCrossScalar(const tUIndex& n,
                                              const tUSInt& dim,
                                              const tUSInt& p,
                                              const tBoolean& incS,
                                              const tReal* S,
                                              const EMM_2PackedSymmetricTensors& phi,
                                              const tBoolean& incU,
                                              const tReal *U,
                                              const tBoolean& incV,
                                              const tReal *V) const;
    
    virtual tReal doubleDotCrossSquaredScalar(const tUIndex& n,
                                              const tUSInt& dim,
                                              const tUSInt& p,
                                              const tBoolean& incS,
                                              const tReal* S,
                                              const tBoolean& incU,
                                              const tReal *U,
                                              const tBoolean& incV,
                                              const tReal *V,
                                              const tBoolean& incW,
                                              const tReal *W,
                                              const tBoolean& incZ,
                                              const tReal *Z) const;
    
   
};

#endif