File: Edge2D.h

package info (click to toggle)
esys-particle 2.1-4
  • links: PTS, VCS
  • area: main
  • in suites: wheezy
  • size: 7,284 kB
  • sloc: cpp: 77,304; python: 5,647; makefile: 1,176; sh: 10
file content (79 lines) | stat: -rw-r--r-- 2,326 bytes parent folder | download
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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
/////////////////////////////////////////////////////////////
//                                                         //
// Copyright (c) 2003-2011 by The University of Queensland //
// Earth Systems Science Computational Centre (ESSCC)      //
// http://www.uq.edu.au/esscc                              //
//                                                         //
// Primary Business: Brisbane, Queensland, Australia       //
// Licensed under the Open Software License version 3.0    //
// http://www.opensource.org/licenses/osl-3.0.php          //
//                                                         //
/////////////////////////////////////////////////////////////

#ifndef __EDGE2D_H
#define __EDGE2D_H

//-- Project includes --
#include "Foundation/vec3.h"
#include "Foundation/Matrix3.h"
#include "Geometry/AEdge.h"

//-- STL includes --
#include <utility>
using std::pair;
using std::make_pair;

//-- IO includes --
#include <iostream>
using std::ostream;

/*!
  \class Edge2D
  \brief class for edge in 2D "mesh"

  \author Steffen Abe
  $Revision$
  $Date$
*/
class Edge2D : public AEdge
{
 public: // types
  typedef Vec3 (Edge2D::* VectorFieldFunction)() const;
  typedef double (Edge2D::* ScalarFieldFunction)() const;
  
 private:
  Vec3 m_normal;
  Vec3 m_force;
  int m_id0,m_id1; // corner ids
  int m_edge_id,m_tag;

 public:
  Edge2D(int,int,const Vec3&,const Vec3&,int,int);
  void moveNode(int,const Vec3&);

  inline int getID(){return m_edge_id;};
  inline void applyForce(const Vec3& f){m_force+=f;};
  inline void zeroForce(){m_force=Vec3(0.0,0.0,0.0);};
  Vec3 getNormal() const {return m_normal;};
  Vec3 toGlobal(const Vec3&);
  Vec3 toLocal(const Vec3&);

  // get id/pos pairs for each node -> mainly for checkpointing
  pair<int,Vec3> getP0()const{return make_pair(m_id0,m_p0);};
  pair<int,Vec3> getP1()const{return make_pair(m_id1,m_p1);};
  
  // access functions
  static VectorFieldFunction getVectorFieldFunction(const string&);
  static ScalarFieldFunction getScalarFieldFunction(const string&);

  Vec3 getForce() const {return m_force;};
  Vec3 getForceDensity() const {return m_force/((m_p1-m_p0).norm());};
  double getPressure() const;

  //! output for debugging purposes
  friend ostream& operator<<(ostream&,const Edge2D&); 
  void print();
};


#endif // __EDGE2D_H