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
|
#include <iostream>
#include "HTMesh.h"
#include "MeshIterator.h"
int main()
{
int level = 5;
printf("level = %d\n", level);
HTMesh *mesh = new HTMesh(level, level);
mesh->setDebug(17);
double ra = 6.75 * 15.;
double dec = -16.72;
//Lookup the triangle containing (ra,dec)
long id = mesh->index(ra, dec);
const char *name = mesh->indexToName(id);
printf("(%8.4f %8.4f): %s\n", ra, dec, name);
double vr1, vd1, vr2, vd2, vr3, vd3;
mesh->vertices(id, &vr1, &vd1, &vr2, &vd2, &vr3, &vd3);
printf("\nThe three vertices of %s are:\n", name);
printf(" (%6.2f, %6.2f)\n", vr1, vd1);
printf(" (%6.2f, %6.2f)\n", vr2, vd2);
printf(" (%6.2f, %6.2f)\n", vr3, vd3);
printf("\n");
//Loop over triangles within one degree of (ra,dec)
double radius = 2.0;
//printf("Leak test! ...\n");
for (int ii = 0; ii < 1; ii++)
{
mesh->intersect(ra, dec, ra - radius, dec, ra - radius, dec + radius);
mesh->intersect(ra, dec, ra - radius, dec, ra - radius, dec + radius, ra, dec + radius);
mesh->intersect(ra, dec, radius);
//continue;
MeshIterator iterator(mesh);
printf("Number of trixels = %d\n\n", mesh->intersectSize());
printf("Triangles within %5.2f degrees of (%6.2f, %6.2f)\n", radius, ra, dec);
while (iterator.hasNext())
{
printf("%s\n", mesh->indexToName(iterator.next()));
}
}
double ra1, dec1, ra2, dec2;
ra1 = 275.874939;
ra2 = 274.882451;
dec1 = -82.470360;
dec2 = -82.482475;
printf("Testing Line intersection: (%f, %f) -->\n", ra1, dec1);
printf(" (%f, %f)\n", ra2, dec2);
mesh->intersect(ra1, dec1, ra2, dec2);
printf("found %d trixels\n", mesh->intersectSize());
return 0;
}
|
