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
|
#include "morpher/quadrilateral.h"
#include <assert.h>
#include <stdlib.h>
static inline bool neighbors_equals(TriangleMap *neighbors[],
TriangleMap *n1, TriangleMap *n2, TriangleMap *n3) {
return neighbors[0] == n1 && neighbors[1] == n2 && neighbors[2] == n3;
}
static inline bool vertices_equals(CartesianMapping maps[],
CartesianMapping m1, CartesianMapping m2, CartesianMapping m3) {
return mappings_equals(maps[0], m1) && mappings_equals(maps[1], m2) && mappings_equals(maps[2], m3);
}
/*
* -
* / \ t
* ---
* r |/| l
* ---
*/
static void test_quadrilateral_flip_diagonal() {
CartesianMapping A = m(50, 0), B = m(0, 100), C = m(100, 100), D = m(0, 200), E = m(100, 200);
TriangleMap *t = trianglemap_create(A, B, C);
TriangleMap *r = trianglemap_create(C, B, E);
TriangleMap *l = trianglemap_create(B, D, E);
trianglemap_set_neighbors(t, NULL, r, NULL, r);
trianglemap_set_neighbors(r, t, l, NULL, l);
trianglemap_set_neighbors(l, NULL, NULL, r, NULL);
quadrilateral_flip_diagonal(r, l);
assert(vertices_equals(t->vertices, A, B, C));
assert(vertices_equals(r->vertices, B, D, C));
assert(vertices_equals(l->vertices, D, E, C));
assert(neighbors_equals(t->neighbors, NULL, r, NULL));
assert(neighbors_equals(r->neighbors, NULL, l, t));
assert(neighbors_equals(l->neighbors, NULL, NULL, r));
while (t != NULL) t = trianglemap_destroy(t);
}
static void test_quadrilateral_is_delaunay(CartesianMapping A, CartesianMapping B, CartesianMapping C,
CartesianMapping D, bool expected) {
TriangleMap *l = trianglemap_create(A, B, C);
TriangleMap *r = trianglemap_create(A, C, D);
trianglemap_set_neighbors(l, NULL, NULL, r, r);
trianglemap_set_neighbors(r, l, NULL, NULL, NULL);
assert(quadrilateral_is_delaunay(l, r) == expected);
trianglemap_destroy(l);
trianglemap_destroy(r);
}
int main(int argc, char **argv) {
test_quadrilateral_flip_diagonal();
test_quadrilateral_is_delaunay(m(0, 0), m(0, 3), m(3, 3), m(2, 1), false);
test_quadrilateral_is_delaunay(m(0, 0), m(0, 3), m(3, 3), m(4, -1), true);
test_quadrilateral_is_delaunay(m(0, 0), m(0, 1), m(1, 1), m(1, 0), true);
return 0;
}
|