Subversion Repositories gelsvn

Rev

Rev 103 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
39 bj 1 
#include <vector>
2 
#include "CGLA/Mat3x3f.h"
3 
#include "CGLA/Vec2f.h"
4 
#include "Face.h"
5 
#include "Manifold.h"
6 
#include "FaceCirculator.h"
7 
 
8 
using namespace CGLA;
9 
using namespace std;
10 
 
11 
namespace HMesh
12 
{
129 jab 13 
		template<class R>
14 
		Face_template<R>::Face_template():
15 
				last(NULL_HALFEDGE_ITER),
16 
				touched(0)
17 
		{}
39 bj 18 
 
129 jab 19 
		int no_edges(FaceIter f)
20 
		{
21 
				FaceCirculator fc(f);
22 
				while(!fc.end()) ++fc;
23 
				return fc.no_steps();
24 
		}
39 bj 25 
 
129 jab 26 
		Vec3f normal(FaceIter f)
27 
		{
28 
				vector<Vec3f> v;
29 
				FaceCirculator fc(f);
30 
				int k;
31 
				for(k=0;!fc.end();++fc,++k)
32 
				{
33 
						Vec3f p = fc.get_vertex()->pos;
34 
						v.push_back(p);
35 
				}
36 
				Vec3f norm(0);
37 
				for(int i=0;i<k;++i)
38 
				{
39 
						norm[0] += (v[i][1]-v[(i+1)%k][1])*(v[i][2]+v[(i+1)%k][2]);
40 
						norm[1] += (v[i][2]-v[(i+1)%k][2])*(v[i][0]+v[(i+1)%k][0]);
41 
						norm[2] += (v[i][0]-v[(i+1)%k][0])*(v[i][1]+v[(i+1)%k][1]);
42 
				}
43 
				float l = norm.length();
44 
				if(l>0.0f)
45 
						norm /= l;
46 
				return norm;
47 
		}
39 bj 48 
 
129 jab 49 
		float area(FaceIter f)
50 
		{
51 
				FaceCirculator fc(f);
52 
				int k;
39 bj 53 
 
129 jab 54 
				// M is a matrix that projects a vector onto the orthogonal
55 
				// complement of the face normal
56 
				Vec3f n = normal(f);
57 
				Vec3f a,b;
58 
				orthogonal(n, a, b);
59 
				Mat3x3f M(a,b,n);
39 bj 60 
 
129 jab 61 
				// Get all projected vertices
62 
				vector<Vec2f> v;
63 
				for(k=0;!fc.end();++fc,++k)
64 
				{
65 
						Vec3f p = M * fc.get_vertex()->pos;
66 
						v.push_back(Vec2f(p[0], p[1]));
67 
				}
68 
				float area = 0;
69 
				for(int i=0;i<k;++i)
70 
				{
71 
						area += 0.5 * cross(v[i], v[(i+1)%k]);
72 
				}
73 
				return fabs(area);
74 
		}
39 bj 75 
 
129 jab 76 
		Vec3f centre(FaceIter f)
77 
		{
78 
				Vec3f c(0);
79 
				FaceCirculator fc(f);
80 
				while(!fc.end())
81 
				{
82 
						c += fc.get_vertex()->pos;
83 
						++fc;
84 
				}
85 
				c /= fc.no_steps();
86 
				return c;
87 
		}
39 bj 88 
 
129 jab 89 
		template struct Face_template<Iters>;
90 
 
39 bj 91 
}