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