| 77 |
|
{ Compute oriented axis values } |
| 78 |
|
inc_dz = -Dx*rNx-Dy*rNy-Dz*rNz; |
| 79 |
|
inc_rx = -RHS*(Dx*(Uy*rNz-Uz*rNy) + Dy*(Uz*rNx-Ux*rNz) + Dz*(Ux*rNy-Uy*rNx)); |
| 80 |
< |
inc_ry = -Dx*Ux-Dy*Uy-Dz*Uz - inc_dz*(rNx*Ux+rNy*Uy+rNz*Uz); |
| 80 |
> |
inc_ry = Dx*Ux+Dy*Uy+Dz*Uz + inc_dz*(rNx*Ux+rNy*Uy+rNz*Uz); |
| 81 |
|
inc_den2 = inc_rx*inc_rx + inc_ry*inc_ry; |
| 82 |
|
inc_radf = if(inc_den2-1e-7, sqrt((1 - inc_dz*inc_dz)/inc_den2), 0); |
| 83 |
|
{ Pass to formulas in first section } |
| 84 |
|
in_disk_x = inc_rx*inc_radf; |
| 85 |
< |
in_disk_y = inc_ry*inc_radf; |
| 85 |
> |
in_disk_y = -inc_ry*inc_radf; |
| 86 |
|
{ Compute final bin (-1 if behind surface) } |
| 87 |
|
scbin = if(inc_dz, |
| 88 |
|
floor(out_square_x*SCdim)*SCdim + floor(out_square_y*SCdim), |
