From 7b85798562f00abef44a457a10aa689a5e8f5581 Mon Sep 17 00:00:00 2001
From: cignoni <paolo.cignoni@isti.cnr.it>
Date: Tue, 29 May 2007 14:01:24 +0000
Subject: [PATCH] added bestmachting normal

---
 vcg/math/gen_normal.h | 25 ++++++++++++++++++++++++-
 1 file changed, 24 insertions(+), 1 deletion(-)

diff --git a/vcg/math/gen_normal.h b/vcg/math/gen_normal.h
index caa08d58..7fb7bddf 100644
--- a/vcg/math/gen_normal.h
+++ b/vcg/math/gen_normal.h
@@ -110,6 +110,29 @@ static void Perturb(std::vector<Point3<ScalarType > > &NN)
 
 }
 
+/*
+Trova la normale piu vicina a quella data.
+Assume che tutte normale in ingresso sia normalizzata;
+*/
+static int BestMatchingNormal(const Point3x &n, std::vector<Point3x> &nv)
+{
+	int ret=-1;
+	ScalarType bestang=-1;
+	ScalarType cosang;
+	typename std::vector<Point3x>::iterator ni;
+	for(ni=nv.begin();ni!=nv.end();++ni)
+		{
+			cosang=(*ni)*n;
+			if(cosang>bestang) {
+				bestang=cosang;
+				ret=ni-nv.begin();
+			}
+		}
+	assert(ret>=0 && ret <nv.size());
+	return ret;
+}
+
+
 private :
 class OctaLevel
   {
@@ -118,7 +141,7 @@ class OctaLevel
     int level;
     int sz;
 
-    Point3f &Val(int i, int j) {
+    Point3x &Val(int i, int j) {
       assert(i>=0 && i<sz);
       assert(j>=0 && j<sz);
       return v[i+j*sz];