Some correction.

diff --git a/src/VTKViewer/VTKViewer_ConvexTool.cxx b/src/VTKViewer/VTKViewer_ConvexTool.cxx
index f6e19392a74b96b771ad524ef7465a29bee56089..3e22b42a6963b773550ef109bb9d32405b311375 100644 (file)
--- a/src/VTKViewer/VTKViewer_ConvexTool.cxx
+++ b/src/VTKViewer/VTKViewer_ConvexTool.cxx
@@ -192,10 +192,12 @@ bool IsConnectedFacesOnOnePlane( vtkUnstructuredGrid* theGrid,
     float b1 = vtkMath::Norm(vec_b1);
 
     float b2 = vtkMath::Norm(vec_b2);
+    float aCos = vtkMath::Dot(vec_b1,vec_b2)/(b1*b2);
     
-    float angle=180*acosf(vtkMath::Dot(vec_b1,vec_b2)/(b1*b2))/vtkMath::Pi();
+    float angle=90.0;
+    angle = aCos>=1.0 ? 0.0 : 180*acosf(aCos)/vtkMath::Pi();
     
-    if( angle <= FACE_ANGLE_TOLERANCE )
+    if( angle <= FACE_ANGLE_TOLERANCE)
       status = true;
     if (MYDEBUG){
       for(int k=0;k<4;k++){
@@ -205,7 +207,7 @@ bool IsConnectedFacesOnOnePlane( vtkUnstructuredGrid* theGrid,
        }
        cout << p[k][2] << ")   ";
       }
-      cout << "angle="<<angle<<endl;
+      cout << "angle="<<angle<<" status="<<status<<endl;
     }
     
   } else if (common_ids.size() >2){