181 |
|
norm_radians(p) : if(-p - PI/4, p + 2*PI, p); |
182 |
|
in_disk_r = .999995*sqrt(Dx*Dx + Dy*Dy); |
183 |
|
in_disk_phi = norm_radians(atan2(Dy, Dx)); |
184 |
< |
in_disk_rgn = floor((in_disk_phi + PI/4)/(PI/2)) + 1; |
184 |
> |
in_disk_rgn = floor((.999995*in_disk_phi + PI/4)/(PI/2)) + 1; |
185 |
|
out_square_a = select(in_disk_rgn, |
186 |
|
in_disk_r, |
187 |
|
(PI/2 - in_disk_phi)*in_disk_r/(PI/4), |
278 |
|
if ($pctcull >= 0) { |
279 |
|
$cmd .= "-of $td/" . ($bmodnm,$fmodnm)[$forw] . ".flt " . |
280 |
|
"| rttree_reduce -h -ff -t $pctcull -r $tensortree -g $ttlog2"; |
281 |
+ |
$cmd .= " -a" if ($tensortree == 3); |
282 |
|
system "$cmd" || die "Failure running rttree_reduce"; |
283 |
|
} else { |
284 |
|
$cmd .= "$td/" . ($bmodnm,$fmodnm)[$forw] . ".flt"; |