TPZMixedPoissonParabolic.h

Go to the documentation of this file.

//
//  mixedpoisson.h
//  PZ
//
//  Created by Agnaldo Farias on 5/28/12.
//  Copyright (c) 2012 LabMec-Unicamp. All rights reserved.
//

#ifndef TPZMIXEDPOISSONPARABOLICH
#define TPZMIXEDPOISSONPARABOLICH

#include "mixedpoisson.h"

class TPZMixedPoissonParabolic : public TPZMixedPoisson {
    
protected:
    
    // thermal capacity or compressibility
    STATE fRho;
    
    // timestep used
    STATE fDeltaT;
    
public:
    TPZMixedPoissonParabolic();
    
    TPZMixedPoissonParabolic(int matid, int dim);
    
    virtual ~TPZMixedPoissonParabolic();
    
    TPZMixedPoissonParabolic(const TPZMixedPoissonParabolic &cp);
    
    TPZMixedPoissonParabolic &operator=(const TPZMixedPoissonParabolic &copy);
    
    void SetParabolicConstants(STATE rho, STATE deltat)
    {
        fRho = rho;
        fDeltaT = deltat;
    }
    
    void SetDeltaT(STATE deltat)
    {
        fDeltaT = deltat;
    }
    
    void FillDataRequirements(TPZVec<TPZMaterialData > &datavec);
    
    virtual void FillBoundaryConditionDataRequirement(int type,TPZMaterialData &data)
    {
        data.SetAllRequirements(false);
        if (type == 3) {
            data.fNeedsNormal = true;
        }
    }

    virtual void FillBoundaryConditionDataRequirement(int type,TPZVec<TPZMaterialData > &datavec)
    {
        // default is no specific data requirements
        int nref = datavec.size();
        for (int iref = 0; iref <nref; iref++) {
            datavec[iref].SetAllRequirements(false);
        }
        datavec[1].fNeedsSol = true;
        if(type == 50)
        {
            for(int iref = 0; iref<nref; iref++){
                datavec[iref].fNeedsSol = true;
            }
        }
    }

    
    virtual void Print(std::ostream & out);
                
                virtual std::string Name() { return "TPZMixedPoissonParabolic"; }
    
    virtual void Contribute(TPZVec<TPZMaterialData> &datavec, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef);
                
    virtual void Contribute(TPZVec<TPZMaterialData> &datavec, REAL weight, TPZFMatrix<STATE> &ef);
    
    virtual void ContributeBC(TPZVec<TPZMaterialData> &datavec, REAL weight, TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef,TPZBndCond &bc);

    
};

#endif