blockmatrixmodule.h

Go to the documentation of this file.

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

class BlockMatrixModule : public DynamicBase
{
  private:
    SparsityPattern sp;
    SparseMatrix<double> M;

    Matrix m_MSi,m_MHomSi;
    OrthoMesh &_mesh;
    unsigned _m;     // total number of raw vertices : _mesh.numRawVertices()
    unsigned n_fracture_cells;
  
    VecDouble _BcPcNew,_BcPcOld,_SwIntegral, _vPcInitial,rhs_Si, rhs_Sihom, Sol, Sol_hom, Sol_old, Sol_hom_old, sat_old;
    VecDouble Solution1, Solution2;   // solution vector contain both Si_c and Si_w
    double _dt, _dx, _dy, _dz;
    double _Bx;
    VecDouble _vPor;
    VecDouble _vK;
    Function1D &_fdpc, &_fpcinv, &_fmobt, &_fmobrW; 
    
  protected:
    //create basis function for a single matrix block
    void create_bases_functions(double bc_pc, const VecDouble &previous_sw, const VecDouble &previous_si, double por, double k, VecDouble &newSi, VecDouble &newHomSi);
    void buildSparsityPattern();
    void assembleSystem(const VecDouble &Sol_old, const VecDouble &sat_old, double _k, double _por, double _pc_bc, double dt);
    void Sol_Si();
    void Sol_Sihom();

  public:
    BlockMatrixModule(OrthoMesh &mesh, unsigned n_fracture_cell, VecDouble &vPcInitial, VecDouble &vPor, VecDouble &k, Function1D &dfpc ,Function1D &fpcinv, Function1D &fmobt, Function1D &fmobrW);

    //Create the bases functions for all the matrix blocks
    void construct_bases(VecDouble &bc_pc, double dt); 
  
    double get_source(unsigned blockId, double new_bc_pc);
  
    ~BlockMatrixModule(){}
    virtual void iterate(TransportBase &trans);
    virtual void printOutput(){}

};

#endif

---