pzgradientflux.h

Go to the documentation of this file.

#ifndef TPZGRADIENTFLUX_H
#define TPZGRADIENTFLUX_H

#include "pzreal.h"  // for STATE, REAL
#include "pzvec.h"   // for TPZVec
class TPZMaterialData;

class TPZGradientFlux{
                
public:
                TPZGradientFlux();
                TPZGradientFlux(const TPZGradientFlux &cp);
                ~TPZGradientFlux();
                
                void ComputeFlux(TPZVec<STATE> &solL, TPZVec<STATE> &solR, const TPZVec<REAL> &normal, TPZVec<STATE> & F);
                
                void ApplyLimiter(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright);
                
private:
                
                void ApplyVanAlbadaLimiter(STATE &soll, STATE &solr,
                                                                                                                   const TPZVec<STATE>& gradL, const TPZVec<STATE> &gradR, 
                                                                                                                   const TPZVec<STATE> &normal, 
                                                                                                                   const TPZVec<STATE> &dL, const TPZVec<STATE> & dR);
                
                void ApplyMinModLimiter(STATE &soll, STATE &solr,
                                                                                                                const TPZVec<STATE>& gradL, const TPZVec<STATE> &gradR, 
                                                                                                                const TPZVec<STATE> &normal, 
                                                                                                                const TPZVec<STATE> &dL, const TPZVec<STATE> & dR);

                STATE Dot(const TPZVec<STATE> &A, const TPZVec<STATE> &B) {
                                double result = 0.;
                                int n = A.NElements();
                                for(int i = 0; i < n; i++) result += A[i]*B[i];
                                return result;
                }
                
};

#endif