Subversion Repositories gelsvn

Rev

Rev 119 | Rev 121 | 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 
{
13 
	template<class R>
14 
	Face_template<R>::Face_template():
15 
		last(NULL_HALFEDGE_ITER),
16 
		touched(0)
17 
	{}
18 
 
19 
	int no_edges(FaceIter f)
20 
	{
21 
		FaceCirculator fc(f);
22 
		while(!fc.end()) ++fc;
23 
		return fc.no_steps();
24 
	}
25 
 
62 jab 26 
	Vec3f normal(FaceIter f)
39 bj 27 
	{
28 
		vector<Vec3f> v;
29 
		FaceCirculator fc(f);
30 
		int k;
31 
		for(k=0;!fc.end();++fc,++k)
32 
			{
62 jab 33 
				Vec3f p = fc.get_vertex()->pos;
39 bj 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 
	}
48 
 
62 jab 49 
	float area(FaceIter f)
39 bj 50 
	{
51 
		FaceCirculator fc(f);
52 
		int k;
53 
 
54 
		// M is a matrix that projects a vector onto the orthogonal
55 
		// complement of the face normal
62 jab 56 
		Vec3f n = normal(f);
54 jab 57 
		Vec3f a,b;
58 
		orthogonal(n, a, b);
59 
		Mat3x3f M(a,b,n);
39 bj 60 
 
61 
		// Get all projected vertices
62 
		vector<Vec2f> v;
63 
		for(k=0;!fc.end();++fc,++k)
64 
			{
62 jab 65 
				Vec3f p = M * fc.get_vertex()->pos;
39 bj 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 
	}
75 
 
62 jab 76 
	Vec3f centre(FaceIter f)
39 bj 77 
	{
78 
		Vec3f c(0);
79 
		FaceCirculator fc(f);
80 
		while(!fc.end())
81 
			{
62 jab 82 
				c += fc.get_vertex()->pos;
39 bj 83 
				++fc;
84 
			}
85 
		c /= fc.no_steps();
86 
		return c;
87 
	}
88 
 
89 
}