EMM_2PackedSymmetricTensors.h

Go to the documentation of this file.

#ifndef EMM_2PackedSymmetricTensors_H
#define EMM_2PackedSymmetricTensors_H

#include "EMM_2Tensors.h"

DEFINE_SPTR(EMM_2PackedSymmetricTensors);
class EMM_2PackedSymmetricTensors : public  EMM_2Tensors {
  
  
    SP_OBJECT(EMM_2PackedSymmetricTensors);
  
    // ATTRIBUTES
  
  
public:
  
    
private:

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

  
public:
   
   
    // CONSTRUCTORS
    inline static SP::EMM_2PackedSymmetricTensors New() {
        SP::EMM_2PackedSymmetricTensors ret(new EMM_2PackedSymmetricTensors(),
                                            EMM_2PackedSymmetricTensors::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& i,const tUSInt& j) const {
        return (i>=j)?(*this)[index][j+i*(i+1)/2]:(*this)[index][i+j*(j+1)/2];
    }

    virtual tReal& operator()(const tUIndex& index,const tUSInt& i,const tUSInt& j) {
        return (i>=j)?(*this)[index][j+i*(i+1)/2]:(*this)[index][i+j*(j+1)/2];
    }
    
    //SET methods
    
protected:
   

public:

    /* \brief vector product H= sigma T . M
     * compute the product of the tensor with a field
     * @param ncells: number of cells
     * @param dim: dim of the vector at each point of the field
     * @param incS : increment of cell's weight sigma 
     * @param sigma: weight per cell
     * @param M : value of the dim-vector per cell
     * @param H: value of the return vector
     */
    void product(const tUIndex& nCells,const tUSInt& dim,
                 const tBoolean& incS,const tReal* sigmas,
                 const tReal *M,
                 tReal *H) const;
    
    
    
    // GET methods 
public:
   
    
    
  
    virtual tUSInt getTensorSize() const {
        const tUSInt& dim=getDimension();
        return (dim*(dim+1)/2);
    }

  
   
    
   
};

#endif