mixedhybridbase.h

Go to the documentation of this file.

#ifndef  _MY_MixedHybridBase_
#define _MY_MixedHybridBase_
#include "dynamicbase.h"
#include <grid/tria.h>
#include <lac/sparse_matrix.h>
#include <lac/vector.h>
#include "orthomesh.h"
#include <lac/sparse_direct.h>
#include <set>
#include "sfunctions.h"


class MixedHybridBase 
{
 protected:
  

  struct _ScalarBC
  {
    OrthoMesh::Face_It face;
    double value;
    _ScalarBC(){value=NAN;}
  };
  
  typedef std::vector<_ScalarBC> ScalarBC;
  

  enum BCType {Dirichlet,Neumann};

  struct _HybridFaceBC
  {
    OrthoMesh::Face_It face; //the index of face at boundary
    BCType bcType; //the boundary condition type (Neumman or Dirichlet)
    double value; //value of the bc




    
    //construction
    _HybridFaceBC()
    {
      value=NAN;
    }

    _HybridFaceBC(OrthoMesh::Face_It &pface)
    :face(pface)
    {
      face = pface;
      value=NAN;
    }
    _HybridFaceBC(OrthoMesh::Face_It &pface,BCType ptype,double pvalue)
    {
      face=pface;
      value=pvalue;
      bcType=ptype;
    }
    void print()
    {
      if (bcType == Dirichlet)
        printf("P = %g at %d",value,face->index());
      else
        printf("Vel = %g at %d",value,face->index());
    }

  };
  typedef std::vector<_HybridFaceBC> HybridFaceBC;

 protected:

  

 protected:
  void setBoundaryCondition(HybridFaceBC &vBC,OrthoMesh &mesh,Function3D &fInwardVel,Function3D &fP,bool bCheckNoBC); 
  void addPcTerm(VecDouble &B,OrthoMesh &mesh, ScalarBC &bc,const VecDouble &KPc, const VecDouble &vPc); 
  void assemblySystem(SparseMatrix<double> &A, VecDouble &B,OrthoMesh &mesh, HybridFaceBC &bc,const VecDouble &K, const VecDouble &vMobT, const VecDouble &vGravTerm); 
  void assemblySystem(SparseMatrix<double> &A, VecDouble &B,OrthoMesh &mesh, HybridFaceBC &bc,const VecDouble &K); 


  void assemblySystem(SparseMatrix<double> &A, VecDouble &B,OrthoMesh &mesh, ScalarBC &bc,const VecDouble &vCoeff); 

  
  void getNormalVelocitiesFromPressure(VecDouble &Vq,OrthoMesh &mesh,HybridFaceBC &bc,const VecDouble &K,const VecDouble &vP,const VecDouble &vMobT, const VecDouble &vGravTerm); 
  void getNormalVelocitiesFromPressure(VecDouble &Vq,OrthoMesh &mesh,HybridFaceBC &bc,const VecDouble &K,const VecDouble &vP); 


  void getNormalVelocitiesFromPressure(VecDouble &Vq,OrthoMesh &mesh,HybridFaceBC &bc,const VecDouble &K,const VecDouble &vP,const VecDouble &vMobT, const VecDouble &vGravTerm, ScalarBC &vPcBC,const VecDouble &KPc, const VecDouble &vPc); 
  void buildSparsityPattern(OrthoMesh &mesh,SparsityPattern &sparsePattern); 


  MixedHybridBase(){}
  ~MixedHybridBase(){}

 

  


};

#endif

---