ESyS-Particle  4.0.1
ElasticInteractionGroup.h
00001 
00002 //                                                         //
00003 // Copyright (c) 2003-2011 by The University of Queensland //
00004 // Earth Systems Science Computational Centre (ESSCC)      //
00005 // http://www.uq.edu.au/esscc                              //
00006 //                                                         //
00007 // Primary Business: Brisbane, Queensland, Australia       //
00008 // Licensed under the Open Software License version 3.0    //
00009 // http://www.opensource.org/licenses/osl-3.0.php          //
00010 //                                                         //
00012 
00013 #ifndef __ELASTICINTERACTIONGROUP_H
00014 #define __ELASTICINTERACTIONGROUP_H
00015 
00016 #include "Model/InteractionGroup.h"
00017 #include "Model/ElasticInteraction.h"
00018 #include "Model/IGParam.h"
00019 
00020 //--- IO includes ---
00021 #include <iostream>
00022 using std::ostream;
00023 using std::endl;
00024 
00025 
00029 template <class T>
00030 class CElasticInteractionGroup : public APairInteractionGroup<T>
00031 {
00032  protected:
00033   vector<CElasticInteraction> m_interactions;
00034   AParallelInteractionStorage* m_exIG; //<! if an interaction is in m_exIG, it can't be in m_interactions
00035   double m_k; //<! spring constant
00036  
00037  public:
00038   CElasticInteractionGroup();
00039   CElasticInteractionGroup(const CElasticIGP*);
00040   virtual ~CElasticInteractionGroup(){};
00041   
00042   virtual void setExIG(AParallelInteractionStorage* eg){m_exIG=eg;};
00043   void setParam(const CElasticIGP*);
00044 
00045   virtual void calcForces();
00046 
00050   virtual void setTimeStepSize(double dt)
00051   {
00052   }
00053   
00054   virtual void Update(ParallelParticleArray<T>*);
00055   friend ostream& operator<< <>(ostream&,const CElasticInteractionGroup<T>&); 
00056 };
00057 
00058 #include "Model/ElasticInteractionGroup.hpp"
00059 
00060 #endif //__ELASTICINTERACTIONGROUP_H