From 3791dc189267c714a56f256b10af58336eb04504 Mon Sep 17 00:00:00 2001
From: Tristan Carrier Baudouin <tristan.carrier@uclouvain.be>
Date: Mon, 2 Dec 2013 15:00:34 +0000
Subject: [PATCH] always the same directions on edges
---
Mesh/directions3D.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/Mesh/directions3D.cpp b/Mesh/directions3D.cpp
index 343988497a..80a479299e 100644
--- a/Mesh/directions3D.cpp
+++ b/Mesh/directions3D.cpp
@@ -109,7 +109,7 @@ STensor3 Frame_field::search(double x,double y,double z){
ANNidxArray indices;
ANNdistArray distances;
- if(field.size()==0){
+ if(field.size()<=1){
return STensor3(1.0);
}
@@ -133,7 +133,7 @@ STensor3 Frame_field::search(double x,double y,double z){
delete[] distances;
#endif
- if(fabs(distance2-distance1)<e2){
+ if(fabs(sqrt(distance2)-sqrt(distance1))<e2){
if(labels[index2]<labels[index1]){
return field[index2].second;
}
--
GitLab