doubleporositydiff.h

Go to the documentation of this file.

#ifndef  _MY__DoublePorosityDiff_
#define _MY_DoublePorosityDiff_
#include "diffusivestep.h"
#include "sfunctions.h"
#include "orthomesh.h"
#include "blockmatrixmodule.h"

class DoublePorosityDiff : public DiffusiveStep
{
 private:
  OrthoMesh &mesh;
  const VecDouble &_por;
  const VecDouble &_K;
  Function1D &_fmob;
  Function1D &_dPC;
  Function1D &_PC;
  VecDouble S_new;
  VecDouble _S_oldTime;
  VecDouble _vPot;
  BlockMatrixModule &_blockMatrix;
  double q_mf, dq_mf;
  unsigned _numTimeStep;

 protected:

 public:
  DoublePorosityDiff(OrthoMesh &mesh, const VecDouble &K,const VecDouble &por, Function1D &fmob, Function1D & dPC,Function1D &_PC,BlockMatrixModule& blockMatrix);
  virtual ~DoublePorosityDiff(){printf("I HAVE BEEN CALLED \n");}
  virtual void iterate(double dt);
  void iterate_linear(double dt);
  double fracturePc(double sat);
  VecDouble& fracturePot(const VecDouble &sat,VecDouble &C);
  virtual void printOutput();
};

class block_mat
{
 private:
  double _potBar; // in reality it should be of VecDouble type
  double _potHat; // in reality it should be of VecDouble type
 public:
  block_mat()
  {
    _potBar = 0;
    _potHat = 0;
  }
  ~block_mat(){}
  void basis_fcts(double potAt_n)
  {
    _potBar = potAt_n;
    _potHat = potAt_n;
  }
  double getSourceTerm(double potAt_n1)
  {
    return 0.0001;//_potBar+potAt_n1*_potHat;
  }
};


#endif

---