[gschaider]

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

[mattijs]

Thanks. I would like to leave the fast intersection() routine intact since it is used all over so want to recalculate the intersection with plane inside the ray() routine. Can you try attached triangleI.H and/or post the testcase?

[mattijs]

Missing attachment.

[gschaider]

Originally Posted by mattijs

Thanks. I would like to leave the fast intersection() routine intact since it is used all over so want to recalculate the intersection with plane inside the ray() routine. Can you try attached triangleI.H and/or post the testcase?

Yep. That fixes it. I'll mail you about the testcase

[mattijs]

I've pushed the fix to 1.6.x and 1.5.x.

[CYMa]

Originally Posted by gschaider

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

Hallo,
Bernhard,,may you send me your recompiled stitchMesh.C,because i really don´t know where this patch should be put in the cold file..
Thank you very much! My e-mail is xgli@iwt.uni-bremen.de

[gschaider]

Originally Posted by CYMa

Hallo,
Bernhard,,may you send me your recompiled stitchMesh.C,because i really don´t know where this patch should be put in the cold file..
Thank you very much! My e-mail is xgli@iwt.uni-bremen.de

The patch discussed above goes somewhere (don't know exactly where, but the find-command is your friend) in $FOAM_SRC. You will have to recompile your OF-installation. As far as I remember no changes to stitchMesh itself were necessary.
The relevant patches should be in 1.5.x and 1.6.x so if you got one of those, everything's OK

Bernhard

[CYMa]

Thank you very much,,,I will have a try....

Originally Posted by gschaider

The patch discussed above goes somewhere (don't know exactly where, but the find-command is your friend) in $FOAM_SRC. You will have to recompile your OF-installation. As far as I remember no changes to stitchMesh itself were necessary.
The relevant patches should be in 1.5.x and 1.6.x so if you got one of those, everything's OK

Bernhard