EMM_DisplacementFVM_VIGROperator.h

Go to the documentation of this file.

#ifndef EMM_DisplacementFVM_VIGROperator_H
#define EMM_DisplacementFVM_VIGROperator_H

#include "EMM_DisplacementFVM_VGROperator.h"

DEFINE_SPTR(EMM_DisplacementFVM_VIGROperator);
class EMM_DisplacementFVM_VIGROperator : public EMM_DisplacementFVM_VGROperator  {
    
    SP_OBJECT(EMM_DisplacementFVM_VIGROperator);
    // ATTRIBUTES
    
    
private:
    
    
    
public:
    
    
    // ASSOCIATION
    
private:

    //gradient U at cells
    mutable SP::EMM_RealField mUp;
   
   
    
protected: 
    
    // METHODS
    
    // CONSTRUCTORS 
    
    EMM_DisplacementFVM_VIGROperator(void);
    
    // DESTRUCTORS 
    
    virtual ~EMM_DisplacementFVM_VIGROperator(void);
    
    
public:

    //Get methods

   
    
    inline const EMM_RealField& getUVertexInterpolation() const {
        return *mUp.get();
    }
    
    inline EMM_RealField& getUVertexInterpolation()  {
        return *mUp.get();
    }
    

   


    //Fields methods
    //=============
    
    virtual void computeElasticTensor(const EMM_RealField& U, EMM_2PackedSymmetricTensors& eTensor) {
        
        //interpolate U at vertices from the displacement fields
        interpolateUAtVertices(true,
                               getSegmentsNumber(),
                               isDirectionPeriodic(),
                               getNeighborsIndices(),
                               U,
                               getLimitConditionOnPoints(),
                               getDisplacementOnDirichletBoundary(),
                               *mUp.get());
        
        EMM_DisplacementFVMOperator::computeElasticTensor(U,eTensor);
    
    }
    virtual void initializeEquilibriumSolver( const EMM_RealField& U0) {
        //Z=0;
        EMM_RealField &Z=getAccelerator();
        Z.setDimension(getDisplacement().getDimension());
        Z.setSize(getDisplacement().getSize());
        Z.initField(0.);
        
        //interpolate U at vertices only from displacement field
        interpolateUAtVertices(true,
                               getSegmentsNumber(),
                               isDirectionPeriodic(),
                               getNeighborsIndices(),
                               Z,
                               getLimitConditionOnPoints(),
                               getDisplacementOnDirichletBoundary(),
                               *mUp.get());
        
        EMM_DisplacementFVMOperator::initializeEquilibriumSolver(U0);
    }
    
    virtual void computeElasticStressMatrixProduct(const tUIndex& nData,const tDimension& dim, const tReal* U,tReal* D) const {
        
        //interpolate U at vertices from the displacement fields
        interpolateUAtVertices(false,
                               getSegmentsNumber(),
                               isDirectionPeriodic(),
                               getNeighborsIndices(),
                               nData,dim,U,
                               getLimitConditionOnPoints(),getDisplacementOnDirichletBoundary(),
                               *mUp.get());
        
        EMM_DisplacementFVMOperator::computeElasticStress(false,-1,nData,dim,U,D);
    }

   
};

#endif