pzcoupledtransportdarcyBC.cpp

Go to the documentation of this file.

#include "pzcoupledtransportdarcyBC.h"
#include "pzadmchunk.h"
#include "pzcmesh.h"

TPZCoupledTransportDarcyBC::TPZCoupledTransportDarcyBC(TPZCoupledTransportDarcy *material, int id) : 
TPZRegisterClassId(&TPZCoupledTransportDarcyBC::ClassId),TPZBndCond(){
                this->SetId(id);
                this->fMaterial = material;
                this->fMaterials[0] = NULL;
                this->fMaterials[1] = NULL;
}

TPZCoupledTransportDarcyBC::~TPZCoupledTransportDarcyBC(){}

void TPZCoupledTransportDarcyBC::Contribute(TPZMaterialData &data,
                                            REAL weight,
                                            TPZFMatrix<STATE> &ek,
                                            TPZFMatrix<STATE> &ef){
                
                
                TPZBndCond * bc = this->GetCurrentMaterial();
                if (!bc) return;
    int numbersol = data.dsol.size();
    if (numbersol != 1) {
        DebugStop();
    }
                UpdateConvectionDir(data.dsol[0]);
                bc->Contribute(data, weight, ek, ef);
}


void TPZCoupledTransportDarcyBC::ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright,
                                                     REAL weight,
                                                     TPZFMatrix<STATE> &ef) {
                
    int numbersol = dataleft.dsol.size();
    if (numbersol != 1) {
        DebugStop();
    }

                
                TPZBndCond * bc = this->GetCurrentMaterial();
                if (!bc) return;
                this->UpdateConvectionDirInterface(dataleft.dsol[0], dataright.dsol[0], dataleft.phi, dataright.phi);
                bc->ContributeInterface(data,dataleft, dataright, weight, ef);
}

void TPZCoupledTransportDarcyBC::
ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright,
                                                                                REAL weight,
                                                                                TPZFMatrix<STATE> &ek,
                                                                                TPZFMatrix<STATE> &ef){
                
    int numbersol = dataleft.dsol.size();
    if (numbersol != 1) {
        DebugStop();
    }

                // TPZFMatrix<REAL> &dphi = data.dphix;
                // TPZFMatrix<REAL> &dphiL = data.dphixl;
                // TPZFMatrix<REAL> &dphiR = data.dphixr;
                // TPZFMatrix<REAL> &phi = data.phi;
                TPZFMatrix<REAL> &phiL = dataleft.phi;
                TPZFMatrix<REAL> &phiR = dataright.phi;
                // TPZManVector<REAL,3> &normal = data.normal;
                // TPZManVector<REAL,3> &x = data.x;
                // int &POrder=data.p;
                // int &LeftPOrder=data.leftp;
                // int &RightPOrder=data.rightp;
                // TPZVec<REAL> &sol=data.sol;
                // TPZVec<REAL> &solL=data.soll;
                // TPZVec<REAL> &solR=data.solr;
                // TPZFMatrix<REAL> &dsol=data.dsol;
                TPZFMatrix<STATE> &dsolL=dataleft.dsol[0];
                TPZFMatrix<STATE> &dsolR=dataright.dsol[0];
                // REAL &faceSize=data.HSize;
                // TPZFMatrix<REAL> &daxesdksi=data.daxesdksi;
                // TPZFMatrix<REAL> &axes=data.axes;
                
                
                TPZBndCond * bc = this->GetCurrentMaterial();
                if (!bc) return;
                this->UpdateConvectionDirInterface(dsolL, dsolR, phiL, phiR);
                bc->ContributeInterface(data, dataleft, dataright, weight, ek, ef);
}

void TPZCoupledTransportDarcyBC::UpdateConvectionDir(TPZFMatrix<STATE> &dsol){
                TPZCoupledTransportDarcy * mat = dynamic_cast<TPZCoupledTransportDarcy*>(this->Material());
                if (!mat){
                                PZError << __PRETTY_FUNCTION__ << " FATAL ERROR" << std::endl;
                                exit(-1);
                }
                mat->UpdateConvectionDir(dsol);
}

void TPZCoupledTransportDarcyBC::UpdateConvectionDirInterface(TPZFMatrix<STATE> &dsolL, TPZFMatrix<STATE> &dsolR,
                                                              TPZFMatrix<REAL> &phiL, TPZFMatrix<REAL> &phiR){
                TPZCoupledTransportDarcy * mat = dynamic_cast<TPZCoupledTransportDarcy*>(this->Material());
                if (!mat){
                                PZError << __PRETTY_FUNCTION__ << " FATAL ERROR" << std::endl;
                                exit(-1);
                }
                if (phiL.Rows()) mat->UpdateConvectionDir(dsolL);
                if (phiR.Rows()) mat->UpdateConvectionDir(dsolR);
}

int TPZCoupledTransportDarcyBC::ClassId() const{
    return Hash("TPZCoupledTransportDarcyBC") ^ TPZBndCond::ClassId() << 1;
}