<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font size="+1">The problem is <font size="+1">on the second<font
size="+1"> run, moving the s<font size="+1">tarting point in
the <font size="+1">to y=4, then <font size="+1">we <font
size="+1">get <font size="+1">no intersections<br>
<font size="+1">bey<font size="+1">ond 48.5.</font></font><br>
</font></font></font></font></font></font></font></font>
<div class="moz-cite-prefix">On 07/24/2013 10:50 AM, Iulian
Grindeanu wrote:<br>
</div>
<blockquote
cite="mid:195029797.14314638.1374681019383.JavaMail.root@mcs.anl.gov"
type="cite">
<style type="text/css">p { margin: 0; }</style>
<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>
</blockquote>
<br>
</body>
</html>