<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: times new roman,new york,times,serif; font-size: 12pt; color: #000000'>I see no error<br>Your fifth slab is at about 49.5; <br>(first at 9.5, ..., fifth at 49.5)<br><br>Iulian<br><br><hr id="zwchr"><blockquote style="border-left:2px solid #1010FF;margin-left:5px;padding-left:5px;color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;">I think I may have found a bug in kdtree.<br><br>My test problem consists of a tet mesh that is composed of some 4000 <br>tets. The mesh looks like that in mesh.png, a ray starts at position <br>(1,0,0) and fires along (1,0,0) and we get the following results back <br>from my test program swp.cpp<br><br>Compile the minimum working example as<br><br>g++ swp.cpp -I/<moab include dir> -L<moab lib dir> -lMOAB<br><br><br>This is as expected,<br><br># of intersections : 20<br>intersection distances are on 8.5 8.74072 9.25928 9.5 18.5 18.7407 <br>19.2593 19.5 28.5 28.7407 29.2593 29.5 38.5 38.7407 39.2593 39.5 48.5 <br>48.7407 49.2593 49.5 of ray length 900<br><br>8.5 9.5 18.5 19.5 etc are the boundaries of the faces of the mesh.<br><br>However, if we recompile and change the starting position to (1,4,0) we get<br><br># of intersections : 17<br>intersection distances are on 8.5 8.63513 9.06064 9.5 18.5 18.6351 <br>19.0606 19.5 28.5 28.6351 29.0606 29.5 38.5 38.6351 39.0606 39.5 48.5 of <br>ray length 900<br><br>As far this ray is concerned there is no mesh in the 5th slab.<br><br>What have I done wrong?<br><br>Thanks<br><br>Andy<br></blockquote><br></div></body></html>