CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Bugs

Change in triangle::intersection breaks stitchMesh in 1.5 and 1.6

Register Blogs Community New Posts Updated Threads Search

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   August 11, 2009, 16:47
Default Change in triangle::intersection breaks stitchMesh in 1.5 and 1.6
  #1
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51
gschaider will become famous soon enoughgschaider will become famous soon enough
This problem seems to have been introduced between 1.4.1 and 1.5: in earlier versions of OF even when the triangle was not hit the pointHit returned by triangle::intersection got some value. In 1.5 and 1.6 vector::zero is returned. Nevertheless this point value is then used in triangle::ray to calculate the distance to this point. I guess usually this distance is not used, but in slidingInterface:rojectPoints() (approx line 303 of the source file) this distance is used to discard points. Consequently stitchMesh fails (Usually complaining about "Duplicate point found in cut face.") for some meshes

After applying the attached patches to 1.5.x and 1.6 and recompiling stitchMesh happily stitched meshes on which previously it failed.

I don't know whether these patches affect other utilities as well

If you want a mesh to verify the bug/fix, I can provide you with one

Bernhard
Attached Files
File Type: gz triangleIntersect.patch.gz (519 Bytes, 25 views)
File Type: gz triangleIntersect_1.5.patch.gz (583 Bytes, 13 views)
gschaider is offline   Reply With Quote

 


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFX4.3 -build analysis form Chie Min CFX 5 July 12, 2001 23:19


All times are GMT -4. The time now is 17:37.