<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>