~esys-p-dev/esys-particle/trunk

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
/////////////////////////////////////////////////////////////
//                                                         //
// 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 __WALLINTERACTION_H
#define __WALLINTERACTION_H

#include "Wall.h"
#include "Interaction.h"
#include "Particle.h"

/*! 
  \class AWallInteraction
  \brief Abstract base for all interactions between a particle and a wall
  
  \author Steffen Abe
  $Revision$
  $Date$  
*/
template <class T>
class AWallInteraction : public AInteraction
{
protected:
  T *m_p;
  CWall *m_wall;
  /*!
    flag showing if particle is in the inner area of the
    local particle array - needed for global force summation
  */
  bool m_inner_flag; 

public:

  AWallInteraction(T*,CWall*,bool iflag=true);

  virtual ~AWallInteraction(){};

  virtual bool hasTag(int,int) const;
  virtual Vec3 getPosFirst() const {return m_p->getPos();};

  inline bool isInner(){return m_inner_flag;};
  virtual void calcForces()=0;
  virtual double getStiffness(){return 0.0;};
};

#include "WallInteraction.hpp"

#endif //__WALLINTERACTION_H