| 212 |
|
newct[0] = cent[0] - newax[0]; |
| 213 |
|
newct[1] = cent[1] - newax[1]; |
| 214 |
|
newct[2] = cent[2] - newax[2]; |
| 215 |
< |
npl = cyl_partit(ro, pt, pi, mp*3/4, newct, newax, d2); |
| 215 |
> |
npl = cyl_partit(ro, pt, pi, mp/2, newct, newax, d2); |
| 216 |
|
/* upper half */ |
| 217 |
|
newct[0] = cent[0] + newax[0]; |
| 218 |
|
newct[1] = cent[1] + newax[1]; |
| 219 |
|
newct[2] = cent[2] + newax[2]; |
| 220 |
< |
npu = cyl_partit(ro, pt, pi, mp-npl, newct, newax, d2); |
| 220 |
> |
npu = cyl_partit(ro, pt, pi, mp/2, newct, newax, d2); |
| 221 |
|
/* return total */ |
| 222 |
|
return(npl + npu); |
| 223 |
|
} |
| 282 |
|
newct[0] = cent[0] - newax[0]; |
| 283 |
|
newct[1] = cent[1] - newax[1]; |
| 284 |
|
newct[2] = cent[2] - newax[2]; |
| 285 |
< |
npl = flt_partit(ro, pt, pi, mp*3/4, newct, u, v, du2, dv2); |
| 285 |
> |
npl = flt_partit(ro, pt, pi, mp/2, newct, u, v, du2, dv2); |
| 286 |
|
/* upper half */ |
| 287 |
|
newct[0] = cent[0] + newax[0]; |
| 288 |
|
newct[1] = cent[1] + newax[1]; |
| 289 |
|
newct[2] = cent[2] + newax[2]; |
| 290 |
< |
npu = flt_partit(ro, pt, pi, mp-npl, newct, u, v, du2, dv2); |
| 290 |
> |
npu = flt_partit(ro, pt, pi, mp/2, newct, u, v, du2, dv2); |
| 291 |
|
/* return total */ |
| 292 |
|
return(npl + npu); |
| 293 |
|
} |