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
|