[MOAB-dev] Bug in kdtree?

[Iulian Grindeanu]

OK, I see; It is a problem, indeed. 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/055ab449/attachment.html>