From acdcb30360936f160b9d253772bfd58e5dc35324 Mon Sep 17 00:00:00 2001
From: nico <nico.pietroni@gmail.com>
Date: Wed, 9 May 2018 22:02:55 +1000
Subject: [PATCH] added InitQualityFaceTorsion and InitQualityFaceBending

---
 vcg/complex/algorithms/polygonal_algorithms.h | 27 +++++++++++++++++++
 1 file changed, 27 insertions(+)

diff --git a/vcg/complex/algorithms/polygonal_algorithms.h b/vcg/complex/algorithms/polygonal_algorithms.h
index a839af28..1bf179d6 100644
--- a/vcg/complex/algorithms/polygonal_algorithms.h
+++ b/vcg/complex/algorithms/polygonal_algorithms.h
@@ -1113,6 +1113,33 @@ public:
         }
     }
 
+    static ScalarType InitQualityFaceTorsion(PolyMeshType &poly_m)
+    {
+        UpdateFaceNormalByFitting(poly_m);
+        vcg::tri::UpdateNormal<PolyMeshType>::PerVertexFromCurrentFaceNormal(poly_m);
+        ScalarType MaxA=0;
+        for (size_t i=0;i<poly_m.face.size();i++)
+        {
+            poly_m.face[i].Q()=PolygonTorsion(poly_m.face[i]);
+            MaxA=std::max(MaxA,poly_m.face[i].Q());
+        }
+        return MaxA;
+    }
+
+    static ScalarType InitQualityFaceBending(PolyMeshType &poly_m)
+    {
+        UpdateFaceNormalByFitting(poly_m);
+        vcg::tri::UpdateNormal<PolyMeshType>::PerVertexFromCurrentFaceNormal(poly_m);
+        ScalarType MaxA=0;
+        for (size_t i=0;i<poly_m.face.size();i++)
+        {
+            poly_m.face[i].Q()=PolygonBending(poly_m.face[i]);
+            MaxA=std::max(MaxA,poly_m.face[i].Q());
+        }
+        return MaxA;
+    }
+
+
     static void InitQualityVertEdgeLenght(PolyMeshType &poly_m)
     {
         for (size_t i=0;i<poly_m.vert.size();i++)