EMM_DisplacementFVM_SSGROperator.h

Go to the documentation of this file.

#ifndef EMM_DisplacementFVM_SSGROperator_H
#define EMM_DisplacementFVM_SSGROperator_H

#include "EMM_DisplacementFVM_VOGGROperator.h"

DEFINE_SPTR(EMM_DisplacementFVM_SSGROperator);
class EMM_DisplacementFVM_SSGROperator : public EMM_DisplacementFVM_VOGGROperator  {
    
    SP_OBJECT(EMM_DisplacementFVM_SSGROperator);
    // ATTRIBUTES
    
    
private:
    
    
    
public:
    
    
    // ASSOCIATION
    
private:

  
   
    
protected: 
    
    // METHODS
    
    // CONSTRUCTORS 
    
    EMM_DisplacementFVM_SSGROperator(void);
    
    // DESTRUCTORS 
    
    virtual ~EMM_DisplacementFVM_SSGROperator(void);   
    
    
public:
    
    
    //SET Methods
    //============
    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 {
        //get the values of U on all points
        const tReal *Upoints=null;
        tUIndex n;
        getUVertexInterpolation().getValues(n,Upoints);
        
        computeGradUAtFaceByStokesIntegration(withConstraints,
                                              h,
                                              xCell,yCell,zCell,f,nextCell,
                                              Px,Py,Pz,
                                              periodicity,
                                              dim,Ucells,Ui,
                                              lc,Upoints,
                                              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 {
        
        computeGradAlmostNullUAtFaceByStokesIntegration(h,xCell,yCell,zCell,f,nextCell,
                                                        Px,Py,Pz,
                                                        periodicity,neighborsIndices,lc,dim,Ui,gradU);
    }
    
    
private:
    
    //computing grad U at face reconstruction
    //=======================================

    void computeGradUAtFaceByStokesIntegration(const tBoolean& withConstraints,
                                               const tReal *h,
                                               const tUInteger& xCell,const tUInteger& yCell,const tUInteger& zCell,
                                               const tUCInt& f,
                                               const tUIndex& nextCell,
                                               const tUInteger& nPx,const tUInteger& nPy,const tUInteger& nPz,
                                               const tBoolean *periodicity,
                                               const tDimension& dim,
                                               const tReal * Ucells,
                                               const tReal* Ui,
                                               const tLimitCondition* lc,
                                               const tReal *Upoints,
                                               tReal gradU[9]) const;
    

   
   void  computeGradAlmostNullUAtFaceByStokesIntegration(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 tLimitCondition* lc,
                                                         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