ESphereBodyInteractionGroup.h

Go to the documentation of this file.

//                                                         //
// Copyright (c) 2003-2014 by The University of Queensland //
// Centre for Geoscience Computing                         //
// http://earth.uq.edu.au/centre-geoscience-computing      //
//                                                         //
// Primary Business: Brisbane, Queensland, Australia       //
// Licensed under the Open Software License version 3.0    //
// http://www.opensource.org/licenses/osl-3.0.php          //
//                                                         //

#ifndef __ESPHEREBODYINTERACTIONGROUP_H
#define __ESPHEREBODYINTERACTIONGROUP_H

//--- project includes ---
#include "Model/ESphereBodyInteraction.h"
#include "Model/SphereBodyIG.h"
#include "Model/ElasticInteractionGroup.h"
#include "Model/IGParam.h"

template <class T> class ParallelParticleArray;

//--- STL includes ---
#include <map>

using std::map;

class CESphereBodyIGP : public CElasticIGP
{
protected:
  std::string m_spherename;
public:

  CESphereBodyIGP(const std::string&,const std::string&,double);
  virtual void packInto(CVarMPIBuffer*) const;
  std::string getSphereBodyName() const {return m_spherename;};
  friend ostream& operator<<(ostream&,const CESphereBodyIGP&);
};

CESphereBodyIGP* extractESphereBodyIGP(AMPIBuffer*);


// --- Forward decl ---
template <class T> class CESphereBodyInteractionGroup;
template <class T> ostream& operator<< (ostream &, const CESphereBodyInteractionGroup<T> &);

template<class T>
class CESphereBodyInteractionGroup : public ASphereBodyInteractionGroup<T>
{
 protected:
  vector<CElasticSphereBodyInteraction<T> > m_interactions;
  double m_k; 
  double m_k_global; 
        double k_local;

 public:
  CESphereBodyInteractionGroup(TML_Comm*);
  CESphereBodyInteractionGroup(TML_Comm*,CSphereBody*,const CESphereBodyIGP*);
  virtual ~CESphereBodyInteractionGroup(){}

  virtual void setTimeStepSize(double dt)
  {
  }
  
  virtual void calcForces();
  virtual void applyForce(const Vec3&);
  virtual void Update(ParallelParticleArray<T>*);

  friend ostream& operator<< <>(ostream &, const CESphereBodyInteractionGroup &);
};

#include  "ESphereBodyInteractionGroup.hpp"

#endif //__ESPHEREBODYINTERACTIONGROUP_H