pztransientmat.h

Go to the documentation of this file.

#ifndef TRANSIENTMATH
#define TRANSIENTMATH

#include "TPZMaterial.h"

template<class TBASEMAT>
class TPZTransientMaterial : public TBASEMAT {
                
public:
    public:
int ClassId() const override;

                TPZTransientMaterial(int nummat, int dim, REAL TimeStep);
                
                ~TPZTransientMaterial();
                
                TPZTransientMaterial(const TPZTransientMaterial &cp);
                
                void SetExplicit();
                
                void SetImplicit();
                
                virtual void Contribute(TPZMaterialData &data,
                            REAL weight,
                            TPZFMatrix<STATE> &ek,
                            TPZFMatrix<STATE> &ef) override;
                
                virtual void ContributeBC(TPZMaterialData &data,
                                                                                                                  REAL weight,
                                                                                                                  TPZFMatrix<STATE> &ek,
                                                                                                                  TPZFMatrix<STATE> &ef,
                                                                                                                  TPZBndCond &bc) override;
                
                virtual void ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright,
                                     REAL weight,
                                     TPZFMatrix<STATE> &ek,
                                     TPZFMatrix<STATE> &ef) override;
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft,
                                       REAL weight,
                                       TPZFMatrix<STATE> &ek,
                                       TPZFMatrix<STATE> &ef,
                                       TPZBndCond &bc) override;
                
                void SetLastState();
                
                void SetCurrentState();
                
                void SetMassMatrix();
                
                void SetFluxOnly();
                
                void SetTimeStep(REAL TimeStep);
                
                REAL TimeStep();
                
                virtual bool NeedsSolutionToContribute(){
                                return true;
                }
                
                virtual bool NeedsXCoord(){
                                return (this->fStep != ELast);
                }
                
protected:
                
                enum ETemporalScheme{EImplicit = 1, EExplicit = 2};
                
                ETemporalScheme fTemporalIntegrator;
                
                enum STEPS{ENone = -1, ELast = 0, ECurrent = 1, EMassMatrix = 2, EFluxOnly = 3};
                
                STEPS fStep;
                
                REAL fTimeStep;
                
                virtual void ContributeSolutionRhs(TPZVec<STATE> &sol, TPZFMatrix<REAL> &phi, REAL weight, TPZFMatrix<STATE> &ef);
                
                virtual void ContributeTangent(TPZVec<STATE> &sol, TPZFMatrix<REAL> &phi, REAL weight, TPZFMatrix<STATE> &ek);
};

template<class TBASEMAT>
inline void TPZTransientMaterial< TBASEMAT >::SetLastState(){
                this->fStep = ELast;
}

template<class TBASEMAT>
inline void TPZTransientMaterial< TBASEMAT >::SetCurrentState(){
                this->fStep = ECurrent;
}

template<class TBASEMAT>
inline void TPZTransientMaterial< TBASEMAT >::SetMassMatrix(){
                this->fStep = EMassMatrix;
}

template<class TBASEMAT>
inline void TPZTransientMaterial< TBASEMAT >::SetFluxOnly(){
                this->fStep = EFluxOnly;
}

template<class TBASEMAT>
inline void TPZTransientMaterial< TBASEMAT >::SetTimeStep(REAL TimeStep){
                this->fTimeStep = TimeStep;
}

template<class TBASEMAT>
inline REAL TPZTransientMaterial< TBASEMAT >::TimeStep(){
                return this->fTimeStep;
}

#endif