TPZTracerFlow.h

Go to the documentation of this file.

//
//  TPZTracerFlow.h


#ifndef TPZTracerFlow_h
#define TPZTracerFlow_h

#include <stdio.h>
#include "pzdiscgal.h"
#include "pzbndcond.h"

class TPZTracerFlow : public TPZDiscontinuousGalerkin {
    
private:
     // TODO:: Candidate for deletion
    int m_dimension;
    
    int m_mat_id;
    
    bool m_mass_matrix_Q;
    
    REAL m_dt;
    
    REAL m_phi;
    
    REAL m_fracture_epsilon;
    
public:
    
    TPZTracerFlow();
    
    TPZTracerFlow(int matid, int dimension);
    
    TPZTracerFlow(const TPZTracerFlow &other);
    
    TPZTracerFlow &operator=(const TPZTracerFlow &other);
    
    ~TPZTracerFlow();
    
    virtual void FillDataRequirements(TPZVec<TPZMaterialData> &datavec) override ;
    
    virtual void FillBoundaryConditionDataRequirement(int type, TPZVec<TPZMaterialData> &datavec) override;
    
    virtual void FillDataRequirementsInterface(TPZMaterialData &data) override;
    
    virtual void FillDataRequirementsInterface(TPZMaterialData &data, TPZVec<TPZMaterialData > &datavec_left, TPZVec<TPZMaterialData > &datavec_right) override;
    
    virtual std::string Name() override{
        return "TPZTracerFlow";
    }
    
    int Dimension() const override {return m_dimension;}
    
    void SetDimension(int dim) { m_dimension = dim; }
    
    int NStateVariables() const override {return 1;} // Deprecated, must to be removed
    
    virtual TPZMaterial *NewMaterial() override
    {
        return new TPZTracerFlow(*this);
    }
    
    void Print(std::ostream &out = std::cout) override;
    
    int VariableIndex(const std::string &name) override;
    
    int NSolutionVariables(int var) override;
    
    void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef)  override { DebugStop();}
    void ContributeBC(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef, TPZBndCond &bc)  override{DebugStop();}
    void ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef)  override{DebugStop();}
    void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft, REAL weight, TPZFMatrix<STATE> &ef,TPZBndCond &bc) override {DebugStop();}
    void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft, REAL weight, TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef,TPZBndCond &bc)  override {DebugStop();}
    void ContributeBC(TPZVec<TPZMaterialData> &datavec, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef, TPZBndCond &bc) override {DebugStop();}
    
    
    // Contribute Methods being used
    
    void Solution(TPZVec<TPZMaterialData> &datavec, int var, TPZVec<REAL> &Solout) override;
    
    void Contribute(TPZVec<TPZMaterialData> &datavec, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef) override;
    
    void Contribute(TPZVec<TPZMaterialData> &datavec, REAL weight, TPZFMatrix<STATE> &ef) override;
    
    void ContributeInterface(TPZMaterialData &data, TPZVec<TPZMaterialData> &datavecleft, TPZVec<TPZMaterialData> &datavecright, REAL weight, TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef) override;
    
    void ContributeInterface(TPZMaterialData &data, TPZVec<TPZMaterialData> &datavecleft, TPZVec<TPZMaterialData> &datavecright, REAL weight, TPZFMatrix<STATE> &ef) override;
    
    void ContributeBCInterface(TPZMaterialData &data, TPZVec<TPZMaterialData> &datavecleft, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef, TPZBndCond &bc) override;
    
    void ContributeBCInterface(TPZMaterialData &data, TPZVec<TPZMaterialData> &datavecleft, REAL weight, TPZFMatrix<STATE> &ef, TPZBndCond &bc);
    
    
    int ClassId() const override;
    
    void Write(TPZStream &buf, int withclassid);
    
    void Read(TPZStream &buf, void *context) override;
    
    void SetMassMatrixAssembly(bool mass_matrix_Q){
        m_mass_matrix_Q = mass_matrix_Q;
    }
    
    bool GetMassMatrixAssembly(){
        return m_mass_matrix_Q;
    }
    
    void SetTimeStep(REAL dt){
        m_dt = dt;
    }
    
    REAL GetTimeStep(){
        return m_dt;
    }
    
    void SetPorosity(REAL phi){
        m_phi = phi;
    }
    
    REAL GetPorosity(){
        return m_phi;
    }
    
    void SetFractureCrossLength(REAL fracture_epsilon){
        m_fracture_epsilon = fracture_epsilon;
    }
    
    REAL GetFractureCrossLength(){
        return m_fracture_epsilon;
    }
    
    REAL FractureFactor(TPZMaterialData & data);

};

#endif /* defined(TPZTracerFlow) */