ultimhex/tag_face.cpp at master · lbinria/ultimhex · GitHub

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
#include "tag_face.h"
#include <ultimaille/all.h>

using namespace UM;

namespace algo {

	int tag(Triangle3 &t) {
		vec3 n = t.normal();
		for (int d = 0; d < 3; d++) {

			if (std::abs(n[d]) < std::abs(n[(d + 1) % 3]) ||
				std::abs(n[d]) < std::abs(n[(d + 2) % 3]))
				continue;

			if (n[d] < 0)
				return d;
			else
				return d+3;
		}

		return -1;
	}

	int tag(Triangle3 &t, bool constraints[3]) {
		vec3 n = t.normal();

		for (int d = 0; d < 3; d++) {

			if (constraints[d] &&
				(!constraints[(d + 1) % 3] || constraints[(d + 1) % 3] && std::abs(n[d]) > std::abs(n[(d + 1) % 3])) &&
				(!constraints[(d + 2) % 3] || constraints[(d + 2) % 3] && std::abs(n[d]) > std::abs(n[(d + 2) % 3]))) {

				if (n[d] < 0)
					return d;
				else
					return d+3;

			}
		}

		return -1;
	}

	void naive_tag(Tetrahedra &m, CellFacetAttribute<int> &tag_attr) {
		if (!m.connected())
			m.connect();

		for (auto f : m.iter_facets()) {
			if (!f.on_boundary())
				continue;

			Triangle3 t = f;
			tag_attr[f] = tag(t);
		}
	}

	void naive_tag(Triangles &m, FacetAttribute<int> &tag_attr) {
		for (auto f : m.iter_facets()) {
			Triangle3 t = f;
			tag_attr[f] = tag(t);
		}
	}

	void naive_tag(std::vector<Surface::Facet> facets, FacetAttribute<int> &tag_attr, bool constraints[3]) {
		for (auto f : facets) {
			Triangle3 t = f;
			tag_attr[f] = tag(t, constraints);
		}
	}


	void naive_tag(std::vector<Surface::Facet> facets, GEO::Attribute<GEO::signed_index_t> &flag_attr, bool constraints[3]) {
		for (auto f : facets) {
			Triangle3 t = f;
			flag_attr[f] = tag(t, constraints);
		}
	}

}