manolas
/
vcglib
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Modified IntersectionPlaneSegment so that it returns always the same intersection independently from the segment orientation 

Return false if segment is parallel with the plane.
This commit is contained in:
Paolo Cignoni 2011-05-11 09:56:23 +00:00
parent 3bb7e7bdd9
commit db5320cc64
1 changed files with 10 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
vcg/space/intersection3.h
View File

@ -236,8 +236,16 @@ namespace vcg {
inline bool IntersectionPlaneSegment( const Plane3<T> & pl, const Segment3<T> & s, Point3<T> & p0){
T p1_proj = s.P1()*pl.Direction()-pl.Offset();
T p0_proj = s.P0()*pl.Direction()-pl.Offset();
if ( (p1_proj>0)-(p0_proj<0)) return false;
p0 = s.P0() + (s.P1()-s.P0()) * fabs(p0_proj/(p1_proj-p0_proj));
if ( (p1_proj>0)-(p0_proj<0)) return false;
if(p0_proj == p1_proj) return false;
// check that we perform the computation in a way that is independent with v0 v1 swaps
if(p0_proj < p1_proj)
p0 = s.P0() + (s.P1()-s.P0()) * fabs(p0_proj/(p1_proj-p0_proj));
if(p0_proj > p1_proj)
p0 = s.P1() + (s.P0()-s.P1()) * fabs(p1_proj/(p0_proj-p1_proj));
return true;
}