CSG Coincidence Speckle
==========================
FIXED : Very Thin Cylinder Speckles in CSG difference holes
---------------------------------------------------------------
* ~/opticks_refs/speckle_lvTopESR.png
* FIXED using new primitive CSG_DISC
Complex box cuts from hollowed cylinder speckles
----------------------------------------------------
* ~/opticks_refs/edge_speckle.png
* ~/opticks_refs/speckle_lvSstTopCirRibBase
* how to generalize nudges to avoid subtractions that lead to ghost speckle coincidence ???
* have implemeted bileaf uncoincidence detection and nudge fixing for box3 - box3,
but thats not directly applicable as the box cuts from an already
composite: cy - cy (for cylinder with inner radius)
* could reorder this inner subtraction ?? ie cut out the boxes from the cylinder endcap
* suspect this situation is rather rare, so fixing it via metadata
labelled nudge hints seems appropriate
* could add inner radius param handling to cylinder primitive,
but thats significant work, will probably eventually do this but not a priority
::
tboolean-;tboolean-rip sc
In [3]: target.lv.solid
Out[3]:
[259] Subtraction SstTopCirRibBase0xc264f78
l:[257] Subtraction SstTopCirRibPri-ChildForSstTopCirRibBase0xc264e78
l:[255] Subtraction SstTopCirRibPri-ChildForSstTopCirRibBase0xbf755c0
l:[253] Subtraction SstTopCirRibPri-ChildForSstTopCirRibBase0xc354ef8
l:[251] Tube SstTopCirRibPri0xc2648b8 mm rmin 1200.0 rmax 1220.0 x 0.0 y 0.0 z 231.89
r:[252] Box Cutbox10xc264960 mm rmin 0.0 rmax 0.0 x 2460.0 y 20.0 z 231.89
r:[254] Box Cutbox20xc265a38 mm rmin 0.0 rmax 0.0 x 2460.0 y 100.0 z 20.0
r:[256] Box Cutbox30xc265b98 mm rmin 0.0 rmax 0.0 x 2460.0 y 20.0 z 231.89
r:[258] Box Cutbox40xc265cf8 mm rmin 0.0 rmax 0.0 x 2460.0 y 100.0 z 20.0
In [6]: obj.dump(detailed=True)
[2017-06-16 20:34:57,767] p44626 {/Users/blyth/opticks/analytic/csg.py:712} INFO - CSG.dump name:SstTopCirRibBase0xc264f78_balanced
in(in(in(cy,!cy),in(!bo,!bo)),in(!bo,!bo)) height:3 totnodes:15
intersection;SstTopCirRibBase0xc264f78_balanced : None None
intersection;treebuilder_midop : None None
intersection;treebuilder_bileaf : None None
cylinder;SstTopCirRibPri0xc2648b8_outer : array([ 0., 0., 0., 1220.], dtype=float32) array([-115.945, 115.945, 0. , 0. ], dtype=float32)
cylinder;SstTopCirRibPri0xc2648b8_inner : array([ 0., 0., 0., 1200.], dtype=float32) array([-117.1044, 117.1044, 0. , 0. ], dtype=float32)
intersection;treebuilder_bileaf : None None
box3;Cutbox10xc264960 : array([ 2460. , 20. , 231.89, 0. ], dtype=float32) array([ 0., 0., 0., 0.], dtype=float32)
box3;Cutbox20xc265a38 : array([ 2460., 100., 20., 0.], dtype=float32) array([ 0., 0., 0., 0.], dtype=float32)
intersection;treebuilder_bileaf : None None
box3;Cutbox30xc265b98 : array([ 2460. , 20. , 231.89, 0. ], dtype=float32) array([ 0., 0., 0., 0.], dtype=float32)
box3;Cutbox40xc265cf8 : array([ 2460., 100., 20., 0.], dtype=float32) array([ 0., 0., 0., 0.], dtype=float32)
in
in in
in in !bo !bo
cy !cy !bo !bo
::
In [7]: 231.89/2.
Out[7]: 115.945
In [9]: 1.01*231.89/2.
Out[9]: 117.10445 # z size Tube subtraction inner is nudged up by 1%
In [10]: 1220*2. # 2460 is more than the diameter, each subtracted does cuts on both sides
Out[10]: 2440.0
::
[2017-06-16 20:58:34,978] p45431 {/Users/blyth/opticks/analytic/csg.py:712} INFO - BALANCED name:SstTopCirRibBase0xc264f78_balanced
in(in(in(cy,!cy),in(!bo,!bo)),in(!bo,!bo)) height:3 totnodes:15
intersection;SstTopCirRibBase0xc264f78_balanced : None None
None
intersection;treebuilder_midop : None None
None
intersection;treebuilder_bileaf : None None
None
cylinder;SstTopCirRibPri0xc2648b8_outer : array([ 0., 0., 0., 1220.], dtype=float32) array([-115.945, 115.945, 0. , 0. ], dtype=float32)
None
cylinder;SstTopCirRibPri0xc2648b8_inner : array([ 0., 0., 0., 1200.], dtype=float32) array([-117.1044, 117.1044, 0. , 0. ], dtype=float32)
None
intersection;treebuilder_bileaf : None None
None
box3;Cutbox10xc264960 : array([ 2460. , 20. , 231.89, 0. ], dtype=float32) array([ 0., 0., 0., 0.], dtype=float32)
array([[ 1., 0., 0., 0.],
[ 0., 1., 0., 0.],
[ 0., 0., 1., 0.],
[ 0., 0., 0., 1.]], dtype=float32)
box3;Cutbox20xc265a38 : array([ 2460., 100., 20., 0.], dtype=float32) array([ 0., 0., 0., 0.], dtype=float32)
array([[ 1. , 0. , 0. , 0. ],
[ 0. , 1. , 0. , 0. ],
[ 0. , 0. , 1. , 0. ],
[ 0. , 0. , -105.945, 1. ]], dtype=float32)
intersection;treebuilder_bileaf : None None
None
box3;Cutbox30xc265b98 : array([ 2460. , 20. , 231.89, 0. ], dtype=float32) array([ 0., 0., 0., 0.], dtype=float32)
array([[ 0.7071, -0.7071, 0. , 0. ],
[ 0.7071, 0.7071, 0. , 0. ],
[ 0. , 0. , 1. , 0. ],
[ 0. , 0. , 0. , 1. ]], dtype=float32)
box3;Cutbox40xc265cf8 : array([ 2460., 100., 20., 0.], dtype=float32) array([ 0., 0., 0., 0.], dtype=float32)
array([[ 0.7071, -0.7071, 0. , 0. ],
[ 0.7071, 0.7071, 0. , 0. ],
[ 0. , 0. , 1. , 0. ],
[ 0. , 0. , -105.945 , 1. ]], dtype=float32)
in
in in
in in !bo !bo
cy !cy !bo !bo
[2017-06-16 20:58:34,983] p45431 {/Users/blyth/opticks/analytic/csg.py:321} INFO - CSG.Serialize : writing 2 trees to directory /tmp/blyth/opticks/tboolean-sc--
analytic=1_csgpath=/tmp/blyth/opticks/tboolean-sc--_name=tboolean-sc--_mode=PyCsgInBox