From 27301dd71b83a65e2652b66bdc9f02597edde688 Mon Sep 17 00:00:00 2001
From: Peter Reischig <peter.reischig@esrf.fr>
Date: Wed, 20 Jun 2012 15:39:51 +0000
Subject: [PATCH] Added correction for numerical errors to get proper zeros.

git-svn-id: https://svn.code.sf.net/p/dct/code/trunk@552 4c865b51-4357-4376-afb4-474e03ccb993
---
 zUtil_Cryst/gtCrystRodriguesFundZone.m | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/zUtil_Cryst/gtCrystRodriguesFundZone.m b/zUtil_Cryst/gtCrystRodriguesFundZone.m
index a344f626..6fddca21 100644
--- a/zUtil_Cryst/gtCrystRodriguesFundZone.m
+++ b/zUtil_Cryst/gtCrystRodriguesFundZone.m
@@ -83,8 +83,12 @@ end
 
 
 % Normalise direction vectors
-ldir = sqrt(sum(faces(:,4:6).*faces(:,4:6),2));
-faces(:,4:6) = faces(:,4:6)./repmat(ldir,1,3);
+ll = sqrt(sum(faces(:,4:6).*faces(:,4:6),2));
+faces(:,4:6) = faces(:,4:6)./ll(:,[1 1 1]);
+
+
+% Correct numerical errors
+faces(abs(faces)<1e-15) = 0;
 
 
 end
\ No newline at end of file
-- 
GitLab