EMM_DisplacementFVM_STEGROperator.h

Go to the documentation of this file.

#ifndef EMM_DisplacementFVM_STEGROperator_H
#define EMM_DisplacementFVM_STEGROperator_H

#include "EMM_DisplacementFVM_VTEGROperator.h"

DEFINE_SPTR(EMM_DisplacementFVM_STEGROperator);
class EMM_DisplacementFVM_STEGROperator : public EMM_DisplacementFVM_VTEGROperator  {
    
    SP_OBJECT(EMM_DisplacementFVM_STEGROperator);
    // ATTRIBUTES
    
    
private:
    
    
    
public:
    
    
    // ASSOCIATION
    
private:

  
   
    
protected: 
    
    // METHODS
    
    // CONSTRUCTORS 
    
    EMM_DisplacementFVM_STEGROperator(void);
    
    // DESTRUCTORS 
    
    virtual ~EMM_DisplacementFVM_STEGROperator(void);   
    
    
public:
    
    
    //SET Methods
    //============
    
   
    
    //grad U reconstruction
    //======================
    virtual void computeGradUAtFace(const tBoolean& withConstraints,
                                    const tReal *h,
                                    const tUInteger& xCell,const tUInteger& yCell,const tUInteger& zCell,
                                    const tUCInt& f,const tUIndex& nextCell,
                                    const tUInteger& Px,const tUInteger& Py,const tUInteger& Pz,
                                    const tBoolean *periodicity,
                                    const CORE_UIndexMorseArray& neighborsIndices,
                                    const tDimension& dim,
                                    const tReal * Ucells,
                                    const tReal* Ui,
                                    const tLimitCondition* lc,const tBoolean& incU0, const tReal* U0,
                                    tReal gradU[9]) const {
        
        
        computeGradUAtFaceByTaylorExpansion(withConstraints,
                                            h,
                                            xCell,yCell,zCell,f,
                                            periodicity,
                                            Px,Py,Pz,
                                            neighborsIndices,
                                            neighborsIndices[xCell+(Px-1)*(yCell+(Py-1)*zCell)],
                                            dim,
                                            Ucells,Ui,
                                            lc,incU0,U0,
                                            gradU);
    }
    

protected:
    virtual void computeGradAlmostNullUAtFace(const tReal *h,
                                              const tUInteger& xCell,const tUInteger& yCell,const tUInteger& zCell,
                                              const tUCInt& f,
                                              const tUIndex& nextCell,
                                              const tUInteger& Px,const tUInteger& Py,const tUInteger& Pz,
                                              const tBoolean *periodicity,
                                              const CORE_UIndexMorseArray& neighborsIndices,
                                              const tLimitCondition* lc,
                                              const tDimension& dim,
                                              const tReal* Ui,
                                              tReal gradU[9]) const {
        computeGradAlmostNullUAtFaceByTaylorExpansion(h,
                                                      xCell,yCell,zCell,f,nextCell,
                                                      Px,Py,Pz,
                                                      periodicity,
                                                      neighborsIndices,
                                                      dim,Ui,gradU);
    }
    
    
  
private:
    void computeGradUAtFaceByTaylorExpansion(const tBoolean& withConstraints,
                                             const tReal *h,
                                             const tUInteger& xCell, const tUInteger& yCell,const tUInteger& zCell,
                                             const tUCInt& f,
                                             const tBoolean * periodicity,
                                             const tUInteger& Px,const tUInteger& Py,const tUInteger& Pz,
                                             const CORE_UIndexMorseArray &neighbors,const tUIndex* Ni,
                                             const tDimension& dim,
                                             const tReal *Ucells,const tReal* Ui,
                                             const tLimitCondition* lc,const tBoolean& incU0, const tReal* U0,
                                             tReal gradU[9]) const;
    

    void computeGradAlmostNullUAtFaceByTaylorExpansion(const tReal *h,
                                                       const tUInteger& xCell, const tUInteger& yCell,const tUInteger& zCell,
                                                       const tUCInt& f,
                                                       const tUIndex& nextCell,
                                                       const tUInteger& Px,const tUInteger& Py,const tUInteger& Pz,
                                                       const tBoolean* periodicity,
                                                       const CORE_UIndexMorseArray &neighbors,
                                                       const tDimension& dim,
                                                       const tReal* Ui,
                                                       tReal gradU[9]) const;
    //STRING representation
    //=======================
                                   
    virtual tString toString() const {
        tString ret=EMM_DisplacementFVMOperator::toString()+"\n";
        return ret;
    }


   
};

#endif