1use HTMesh;
2printf "%5s %5s %8s %8s %8s %8s %8s %8s\n",
3    qw/size level trixels stars iter/,  "us/circ", "time(5)", "time(1)";
4for my $size ( 1, 5, 10, 45, 90) {
5    printf "%s\n", "-" x 65;
6    for my $level( 3, 4, 5, 6) {
7        my $mesh = new HTMesh($level, $level);
8        my $num_triangles = $mesh->total_triangles();
9        my ($ra, $dec) = (6.75, -16.72);
10        my @circle = ($ra, $dec, $size);
11        $mesh->intersect_circle(@circle);
12        my $trixels = $mesh->result_size();
13        my $iter = 5000/ sqrt($trixels);
14        my $time = $mesh->time_circle($iter, @circle);
15        my $stars = (126000 / $num_triangles) * $trixels;
16        my $us_per_circ = 1000 * 1000 * $time/$iter;
17        my $tt1 = $us_per_circ + 5 * $stars;
18        my $tt2 = $us_per_circ + 1 * $stars;
19        printf "%5d %5d %8d %8d %8d %8d %8d %8d\n",
20            $size, $level, $trixels, $stars, $iter, $us_per_circ, $tt1, $tt2;
21    }
22    #printf "%s\n", "-" x 72;
23}
24
25
26
27
28
29
30