Subversion Repositories gelsvn

Rev

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

Rev Author Line No. Line
39 bj 1 
#include <iostream>
2 
 
3 
#include "Vertex.h"
4 
#include "Manifold.h"
5 
#include "VertexCirculator.h"
6 
 
7 
namespace HMesh
8 
{
9 
 
129 jab 10 
		using namespace std;
11 
		using namespace CGLA;
39 bj 12 
 
129 jab 13 
		template<class R>
14 
		Vertex_template<R>::Vertex_template(const CGLA::Vec3f& _pos):
15 
				he(NULL_HALFEDGE_ITER),
16 
				pos(_pos),
17 
				touched(0)
18 
		{}
39 bj 19 
 
129 jab 20 
		bool is_boundary(VertexIter v)
21 
		{
22 
				return is_boundary(v->he);
23 
		}
39 bj 24 
 
129 jab 25 
		void check_boundary_consistency(VertexIter v)
26 
		{
27 
				HalfEdgeIter h = v->he;
28 
				HalfEdgeIter h0 = h;
183 jab 29 
 
129 jab 30 
				while(h->face != NULL_FACE_ITER)
31 
				{
32 
						h = h->prev->opp;
33 
						if(h==h0) return;
34 
				}
39 bj 35 
 
129 jab 36 
				v->he = h;
183 jab 37 
				HalfEdgeIter h_opp = h->opp;
129 jab 38 
				while(h_opp->face != NULL_FACE_ITER)
39 
						h_opp = h_opp->next->opp;
183 jab 40 
				h_opp->vert = v;
41 
				link(h_opp,h);
129 jab 42 
		}
39 bj 43 
 
129 jab 44 
		int valency(VertexIter v)
45 
		{
46 
				VertexCirculator vc(v);
47 
				for(;!vc.end();++vc);
48 
				return vc.no_steps();
49 
		}
39 bj 50 
 
129 jab 51 
		const Vec3f normal(VertexIter v)
52 
		{
53 
				Vec3f p0 = v->pos;
54 
				vector<Vec3f> one_ring;
55 
				for(VertexCirculator vc(v);!vc.end();++vc)
56 
				{
57 
						Vec3f p = vc.get_vertex()->pos;
58 
						Vec3f edge = p-p0;
59 
						float len = length(edge);
60 
						if(len>0.0f)
61 
								one_ring.push_back(edge/len);
62 
				}
39 bj 63 
 
129 jab 64 
				Vec3f norm(0);
65 
				int N= one_ring.size();
66 
				for(int i=0;i<N;++i)
67 
				{
68 
						Vec3f e1 = one_ring[i];
69 
						Vec3f e0 = one_ring[(i+1)%N];
70 
						Vec3f n = normalize(cross(e0,e1));
71 
						norm += n*acos(max(-1.0f,min(1.0f,dot(e0,e1))));
72 
				}
73 
				float l = sqr_length(norm);
74 
				if(l>0.0f)
75 
						return norm/sqrt(l);
39 bj 76 
 
129 jab 77 
				return norm;
78 
		}
39 bj 79 
 
129 jab 80 
		bool is_connected(VertexIter v1, VertexIter v2)
81 
		{
82 
				vector<VertexIter> vertices;
83 
				VertexCirculator vc(v1);
84 
				for(;!vc.end();++vc)
85 
						vertices.push_back(vc.get_vertex());
86 
				if(find(vertices.begin(),vertices.end(),v2) != vertices.end())
87 
						return true;
88 
				return false;
89 
		}
90 
 
91 
		template struct Vertex_template<Iters>;
92 
 
39 bj 93 
}