/Users/jack/Code/basso_dev/inc/DiffusionFormulation.h

Go to the documentation of this file.

#ifndef _DIFFUSION_FORUMULATION_BASSO_H_
#define _DIFFUSION_FORUMULATION_BASSO_H_

// std includes
#include <list>

// Basso includes
#include "basso.h"
#include "Formulation.h"

namespace basso {


class DiffusionFormulation : public Formulation {

public:
        
        // PUBLIC DEFS
        
        // CONSTRUCTORS 
        DiffusionFormulation( list< Element > &elem, const Array< Node > &nodes, const Material *material ) 
                : Formulation( elem, nodes ) { fMat=material; Initialize(); }
        
        virtual ~DiffusionFormulation() {  }
        
        // ACCESSORS
        virtual int NumLocalDofs() const { return 1; }
        
protected:
        virtual void GetElementStiffnessMatrix( nMatrix &ke, list<Element>::const_iterator eItr  );
        
        virtual void Initialize(  );
        
        virtual bool SetElementParameters( list<Element>::const_iterator eItr, int qorder=2 );

        virtual Numeric MassFactor() const { return fMat->Density(); }
        
        virtual Numeric StiffnessFactor() const { return fMat->Diffusivity(); }
        
protected:
        
        const Material *fMat;  
};

void DiffusionFormulation::ActivateDofs( DofMap &dofmap )  
{
        fDofMap=&dofmap;
        list< Element >::const_iterator eItr;
        Element::const_iterator nItr;
        for ( eItr=this->fElement->begin(); eItr!=this->fElement->end(); ++eItr ) 
                for ( nItr=eItr->begin(); nItr!=eItr->end(); ++nItr ) 
                        for ( int s=0; s<NumLocalDofs(); ++s )
                                fDofMap->AddNodalDof( *nItr, this->activeLocalDofs[s] );
        fDofMap->RenumberGlobalDofs();          
}



void DiffusionFormulation::Initialize(  ) 
{ 
        resize( activeLocalDofs, NumLocalDofs() );
        activeLocalDofs[0]=kCONCENTRATION; 
}


} // end namespace

#endif

---