MUMM Class Reference

#include <mumm.h>

Inheritance diagram for MUMM:

Collaboration diagram for MUMM:

Public Member Functions
	MUMM (OrthoMesh &mesh, VecDouble &K, Function1D &fMobT)
virtual void	getNormalVelocityAtFaces (VecDouble &vFNC)
virtual void	iterate (TransportBase &trans)
virtual void	printOutput ()
virtual const VecDouble &	getPressureAtCells ()
virtual	~MUMM ()
Private Attributes
OrthoMesh &	m_mesh
VecDouble	m_P
VecDouble &	m_K
Function1D &	_fMobT
VecDouble	_vMobT

Detailed Description

Definition at line 9 of file mumm.h.

Constructor & Destructor Documentation

MUMM::MUMM	(	OrthoMesh &	mesh,
		VecDouble &	K,
		Function1D &	fMobT
	)

Definition at line 5 of file mumm.cpp.

00006   :m_mesh(mesh),m_K(K),_fMobT(fMobT)
00007 {
00008   assert(mesh.numElemZ()==1);
00009   _vMobT.reinit(mesh.numCells());
00010 }

virtual MUMM::~MUMM ( ) [inline, virtual]

Definition at line 26 of file mumm.h.

00026 {}

Member Function Documentation

void MUMM::getNormalVelocityAtFaces ( VecDouble & vFNC ) [virtual]

Reimplemented from DynamicBase.

Definition at line 24 of file mumm.cpp.

00025 {
00026   OrthoMesh::Cell_It cell = m_mesh.begin_cell();
00027   OrthoMesh::Cell_It endc = m_mesh.end_cell();
00028   for(;cell!=endc;cell++)
00029   {
00030     double QL,QR,QD,QU; 
00031 
00032     //Get the QL,QR,QD,QU .....cell->index()
00033 
00034     
00035     vFNC(cell->face_index_up()) = 0;
00036     vFNC(cell->face_index_bottom()) =0;
00037 
00038     vFNC(cell->face_index_left())=QL;
00039     vFNC(cell->face_index_right())=QR;
00040     vFNC(cell->face_index_front())=QD;
00041     vFNC(cell->face_index_back())=QU;
00042   }
00043 
00044 }

const VecDouble & MUMM::getPressureAtCells ( ) [virtual]

Reimplemented from DynamicBase.

Definition at line 46 of file mumm.cpp.

00047 {
00048   return  m_P;
00049 }

void MUMM::iterate ( TransportBase & trans ) [virtual]

Implements DynamicBase.

Definition at line 12 of file mumm.cpp.

00013 {
00014   //vMobT = Transport.SW
00015   trans.getSolutionAtCells().getComponent(_vMobT,0);
00016   _fMobT.evaluate(_vMobT);
00017   const double *d=&(_vMobT(0));
00018   
00019 }