From e11dd6ca1dcdd73265db49bd7fc745d4745b4160 Mon Sep 17 00:00:00 2001
From: Massimiliano Corsini maxcorsini <massimiliano.corsini@isti.cnr.it>
Date: Tue, 31 Mar 2009 21:07:27 +0000
Subject: [PATCH] ASE added

---
 .../filter_perceptualmetric/perceptualmetrics.h        | 10 +++++++++-
 1 file changed, 9 insertions(+), 1 deletion(-)

diff --git a/src/meshlabplugins/filter_perceptualmetric/perceptualmetrics.h b/src/meshlabplugins/filter_perceptualmetric/perceptualmetrics.h
index 39a167b28..e6c68aa16 100644
--- a/src/meshlabplugins/filter_perceptualmetric/perceptualmetrics.h
+++ b/src/meshlabplugins/filter_perceptualmetric/perceptualmetrics.h
@@ -33,6 +33,7 @@ $Log: perceptualmetrics.h,v $
 
 #include <meshlab/meshmodel.h>
 #include <meshlab/interfaces.h>
+#include <vcg/complex/trimesh/stat.h>
 
 template <class MeshType>
 class PerceptualMetrics
@@ -144,6 +145,9 @@ private:
 		double AAb = (yEp-yFp) * xDp + (yFp-yDp) * xEp + (yDp-yEp) * xFp;  // AA = 2A
 
 		area = 0.5 * AA;
+
+		if (area < 0.0)
+			area = -area;
 	}
 
 // public methods
@@ -226,8 +230,12 @@ public:
 			Wdistortion = G * epsilon_ij_prime_epsilon_ij_prime * Sdelta;
 			W += Warea + Wdistortion;
 		}
+
+		double area1 = vcg::tri::Stat<MeshType>::ComputeMeshArea(refmesh);
+		double area2 = vcg::tri::Stat<MeshType>::ComputeMeshArea(mesh);
 		
-		return W;
+		// Average Strain Energy (ASE)
+		return W / (area1+area2);
 	}
 
 };