pzhybridpoisson.cpp

Go to the documentation of this file.

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

#include "pzhybridpoisson.h"
#include "pzlog.h"
#include "pzbndcond.h"
#include "pzfmatrix.h"
#include "pzaxestools.h"


#include <iostream>

#ifdef LOG4CXX
static LoggerPtr logdata(Logger::getLogger("pz.hybridpoisson"));
#endif

TPZHybridPoisson::TPZHybridPoisson(): TPZMatPoisson3d(){
}

TPZHybridPoisson::TPZHybridPoisson(int matid, int dim): TPZMatPoisson3d(matid,dim){
  
}

TPZHybridPoisson::~TPZHybridPoisson(){
}


void TPZHybridPoisson::Print(std::ostream &out) {
                out << "name of material : " << Name() << "\n";
                out << "Base Class properties :";
                TPZMaterial::Print(out);
                out << "\n";
}

void TPZHybridPoisson:: ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef)
{
    if (dataleft.fShapeType != TPZMaterialData::EVecandShape || dataright.fShapeType != TPZMaterialData::EScalarShape) {
        DebugStop();
    }
    // loop over the flux functions
    int nvec = dataleft.fVecShapeIndex.size();
    int npressure = dataleft.numberdualfunctions;
    int nlagrange = dataright.phi.Rows();
    for (int ivec=0; ivec<nvec; ivec++) {
        int vec = dataleft.fVecShapeIndex[ivec].first;
        int ish = dataleft.fVecShapeIndex[ivec].second;
        REAL normal = 0.;
        for (int in=0; in<3; in++) {
            normal += data.normal[in]*dataleft.fNormalVec(in,vec);
        }
        REAL normalishape = normal * dataleft.phi(ish,0);
        for (int jlagr=0; jlagr < nlagrange; jlagr++) {
            ek(ivec,nvec+npressure+jlagr) += weight*normalishape*dataright.phi(jlagr,0);
            ek(nvec+npressure+jlagr,ivec) += weight*normalishape*dataright.phi(jlagr,0);
        }
    }
    // compute their normal component
    // loop over the pressure functions
    // contribute the Langrange multiplier
    
}

void TPZHybridPoisson::ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, REAL weight, TPZFMatrix<STATE> &ef)
{
    if (dataleft.fShapeType != TPZMaterialData::EVecandShape || dataright.fShapeType != TPZMaterialData::EScalarShape) {
        DebugStop();
    }
    DebugStop();
    
}