TPZMixedDarcyFlow.h

Go to the documentation of this file.

//
//  TPZMixedDarcyFlow.h
//  HDiv
//
//  Created by Omar Durán on 3/29/19.
//

#ifndef TPZMixedDarcyFlow_h
#define TPZMixedDarcyFlow_h

#include <stdio.h>
#include "TPZMaterial.h"
#include "pzbndcond.h"
#include "pzfmatrix.h"
#include "pzaxestools.h"

class TPZMixedDarcyFlow : public TPZMaterial {
    
private:
    
    TPZFNMatrix<9,REAL> m_kappa;
    
    TPZFNMatrix<9,REAL> m_kappa_inv;
    
    REAL m_d;
    
    int m_dim;
    
public:
    
    TPZMixedDarcyFlow();
    
    ~TPZMixedDarcyFlow();
    
    TPZMixedDarcyFlow(int mat_id, int dim);
    
    TPZMixedDarcyFlow(const TPZMixedDarcyFlow &other);
    
    TPZMaterial * NewMaterial() override;
    
    TPZMixedDarcyFlow &operator=(const TPZMixedDarcyFlow &other);
    
    int Dimension() const override;
    
    void SetDimension(int dim) { m_dim = dim; }
    
    virtual int NStateVariables() const override;
    
    void Print(std::ostream & out) override;
    
    std::string Name() override;
    
    void FillDataRequirements(TPZVec<TPZMaterialData> &datavec) override;
    
    void FillBoundaryConditionDataRequirement(int type, TPZVec<TPZMaterialData> &datavec) override;
    
    virtual int ClassId() const override;
    
    virtual void Write(TPZStream &buf, int withclassid) const override;
    
    void Read(TPZStream &buf, void *context) override;
    
    void ContributeBC(TPZMaterialData &data,REAL weight,TPZFMatrix<STATE> &ef,TPZBndCond &bc) override {
        DebugStop();
    }
    
    void ContributeBC(TPZMaterialData &data,REAL weight,TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef,TPZBndCond &bc) override {
        DebugStop();
    }
    
    void Contribute(TPZMaterialData &data,REAL weight,TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef) override {
        DebugStop();
    }
    
    void Contribute(TPZMaterialData &data,REAL weight,TPZFMatrix<STATE> &ef) override {
        DebugStop();
    }
    
    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 ContributeBC(TPZVec<TPZMaterialData> &datavec,REAL weight,TPZFMatrix<STATE> &ef,TPZBndCond &bc) override;
    
    void ContributeBC(TPZVec<TPZMaterialData> &datavec,REAL weight,TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef,TPZBndCond &bc) override;
    
    int VariableIndex(const std::string &name) override;
    
    int NSolutionVariables(int var) override;
    
    void Solution(TPZVec<TPZMaterialData> &datavec, int var, TPZVec<STATE> &Solout) override;
    
    void SetPermeability(REAL kappa){
        int n_data = m_kappa.Rows();
        m_kappa.Zero();
        m_kappa_inv.Zero();
        for (int i = 0; i < 3; i++) {
            m_kappa(i,i) = kappa;
            m_kappa_inv(i,i) = 1.0/kappa;
        }
    }
    
    void SetPermeability(TPZFNMatrix<9,REAL> & kappa){
        
        if (kappa.Rows() != 3 && kappa.Cols() != 3) {
            DebugStop();
        }
        
        m_kappa = kappa;
        m_kappa.Inverse(m_kappa_inv, ELU);
    }
    
    void SetDimensionalFactor(REAL d){
        m_d = d;
    }
    
};

#endif /* TPZMixedDarcyFlow_h */