compdiffusivestep.h

Go to the documentation of this file.

#ifndef  _MY_CompDiffusiveStep_
#define _MY_CompDiffusiveStep_
#include "mixedhybridbase.h"
#include "vecincaccess.h"
#include "linearsolver.h"
#include "flashcompositional.h"
#include "diffusivestep.h"

class CompDiffusiveStep : public DiffusiveStep, public MixedHybridBase 
{
 private:
  struct FlashBC
  {
    double Saq;
    //FlashData data;

    //FlashBC(unsigned i,unsigned j)
    //:data(i,j){data.allocateOwnMemory();}
  };
  VecIncAccess<FlashBC> m_flashBC;
  
 protected:

  VecDouble m_vSaq,m_cV,m_vV;
  VecDouble m_vSaqPrev;
  OrthoMesh &m_mesh;
  ScalarBC m_swBC;
  SparsityPattern m_sparsePattern;
  SparseMatrix<double> m_A;
  VecDouble m_sol;
  VecDouble m_B;  
  VecDouble m_vCoeff,m_vCoeff2;
  LinearSolver& m_solver;
  FlashCompositional& m_flash;
  ScalarBC m_bc;
  const VecDouble &m_vK,&m_vPor;
  Function1D &m_fMobPc,&m_DfPc;
 public:
  CompDiffusiveStep(OrthoMesh &mesh, Function3D &fPrescribedPressure, Function3D &fTransportBC,const VecDouble &vK, const VecDouble &vPor,Function1D &fMobPc,  Function1D &DfPc,FlashCompositional &flash,LinearSolver &solver);
  virtual ~CompDiffusiveStep();
  void setSaqBC(VecIncAccess<FlashBC> &flashbc,ScalarBC &swbc,const Function3D &fPressureBC, const Function3D &fTransBC,FlashCompositional &flash); 
  
  virtual void iterate(double dt); 
  virtual void printOutput(){}
  
};

#endif

---