[MOAB-dev] Bug in kdtree?

[Iulian Grindeanu]

OK, obb tree works fine with the same model, default setting for everything 
(you do no have to put tetras in the build, only the faces) 

So there might be an issue with the kdtree. 
I will look at it. 

Iulian 
----- Original Message -----

| The problem is on the second run, moving the s tarting point in the
| to y=4, then we get no intersections
| bey ond 48.5.

| On 07/24/2013 10:50 AM, Iulian Grindeanu wrote:

| | I see no error
| 
| | Your fifth slab is at about 49.5;
| 
| | (first at 9.5, ..., fifth at 49.5)
| 

| | Iulian
| 

| | ----- Original Message -----
| 

| | | I think I may have found a bug in kdtree.
| | 
| 

| | | My test problem consists of a tet mesh that is composed of some
| | | 4000
| | 
| 
| | | tets. The mesh looks like that in mesh.png, a ray starts at
| | | position
| | 
| 
| | | (1,0,0) and fires along (1,0,0) and we get the following results
| | | back
| | 
| 
| | | from my test program swp.cpp
| | 
| 

| | | Compile the minimum working example as
| | 
| 

| | | g++ swp.cpp -I/<moab include dir> -L<moab lib dir> -lMOAB
| | 
| 

| | | This is as expected,
| | 
| 

| | | # of intersections : 20
| | 
| 
| | | intersection distances are on 8.5 8.74072 9.25928 9.5 18.5
| | | 18.7407
| | 
| 
| | | 19.2593 19.5 28.5 28.7407 29.2593 29.5 38.5 38.7407 39.2593 39.5
| | | 48.5
| | 
| 
| | | 48.7407 49.2593 49.5 of ray length 900
| | 
| 

| | | 8.5 9.5 18.5 19.5 etc are the boundaries of the faces of the
| | | mesh.
| | 
| 

| | | However, if we recompile and change the starting position to
| | | (1,4,0)
| | | we get
| | 
| 

| | | # of intersections : 17
| | 
| 
| | | intersection distances are on 8.5 8.63513 9.06064 9.5 18.5
| | | 18.6351
| | 
| 
| | | 19.0606 19.5 28.5 28.6351 29.0606 29.5 38.5 38.6351 39.0606 39.5
| | | 48.5
| | | of
| | 
| 
| | | ray length 900
| | 
| 

| | | As far this ray is concerned there is no mesh in the 5th slab.
| | 
| 

| | | What have I done wrong?
| | 
| 

| | | Thanks
| | 
| 

| | | Andy
| | 
| 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20130724/f75f86c9/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: obb.cpp
Type: text/x-c++src
Size: 4587 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20130724/f75f86c9/attachment.cpp>