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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
|
/////////////////////////////////////////////////////////////
// //
// 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 __TRIANGLE_H
#define __TRIANGLE_H
//-- Project includes --
#include "Foundation/vec3.h"
#include "Foundation/Matrix3.h"
//-- STL includes --
#include <utility>
using std::pair;
using std::make_pair;
//-- IO includes --
#include <iostream>
using std::ostream;
//! exception class for Triangle
class TriangleError
{
public:
TriangleError(){};
};
/*!
\class Triangle
\brief Class representing a Triangle
\author Steffen Abe
$Revision$
$$Date$
*/
class Triangle
{
public: // types
typedef Vec3 (Triangle::* VectorFieldFunction)() const;
typedef double (Triangle::* ScalarFieldFunction)() const;
private:
Matrix3 m_invtrans;
Matrix3 m_trans;
Vec3 m_p0,m_p1,m_p2;
Vec3 m_normal;
Vec3 m_force;
int m_id0,m_id1,m_id2;
int m_tri_id,m_tag;
double EdgeSep(const Vec3&, const Vec3& ,const Vec3& ) const;
public:
Triangle(int,int,int,const Vec3&,const Vec3&,const Vec3&,int,int);
double sep(const Vec3&) const;
pair<bool,double> dist(const Vec3&) const ; // signed separation according to direction of the normal
Vec3 getBoundingBoxMin() const;
Vec3 getBoundingBoxMax() const;
Vec3 getNormal() const {return m_normal;};
Vec3 toGlobal(const Vec3&);
Vec3 toLocal(const Vec3&);
bool containsEdge(const Vec3&,const Vec3&) const;
void moveNode(int,const Vec3&);
void move(const Vec3&);
inline int getID(){return m_tri_id;};
inline void applyForce(const Vec3& f){m_force+=f;};
inline void zeroForce(){m_force=Vec3(0.0,0.0,0.0);};
// 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_p0+m_p1);};
pair<int,Vec3> getP2()const{return make_pair(m_id2,m_p0+m_p2);};
// access functions
static VectorFieldFunction getVectorFieldFunction(const string&);
static ScalarFieldFunction getScalarFieldFunction(const string&);
Vec3 getForce()const {return m_force;};
double getPressure() const;
//! output for debugging purposes
friend ostream& operator<<(ostream&,const Triangle&);
};
#endif //__TRIANGLE_H
|
