From 9691602218739df966d15498c3fe9f1f4c5ac18c Mon Sep 17 00:00:00 2001
From: nico <nico.pietroni@gmail.com>
Date: Wed, 21 Nov 2018 23:18:17 +1100
Subject: [PATCH] made the LaplacianReproject independent wrt scalar type used
 for the tri mesh

---
 vcg/complex/algorithms/polygonal_algorithms.h | 29 ++++++++++++-------
 1 file changed, 18 insertions(+), 11 deletions(-)

diff --git a/vcg/complex/algorithms/polygonal_algorithms.h b/vcg/complex/algorithms/polygonal_algorithms.h
index 972d4af0..d27c6b57 100644
--- a/vcg/complex/algorithms/polygonal_algorithms.h
+++ b/vcg/complex/algorithms/polygonal_algorithms.h
@@ -598,14 +598,16 @@ public:
                                    int nstep=100,
                                    ScalarType Damp=0.5)
     {
-        typedef typename TriMeshType::FaceType FaceType;
-        typedef vcg::GridStaticPtr<FaceType, ScalarType> TriMeshGrid;
+        typedef typename TriMeshType::FaceType TriFaceType;
+        typedef typename TriMeshType::ScalarType TriScalarType;
+        typedef typename TriMeshType::CoordType TriCoordType;
+        typedef vcg::GridStaticPtr<TriFaceType, TriScalarType> TriMeshGrid;
         TriMeshGrid grid;
 
         //initialize the grid
         grid.Set(tri_mesh.face.begin(),tri_mesh.face.end());
 
-        ScalarType MaxD=tri_mesh.bbox.Diag();
+        TriScalarType MaxD=tri_mesh.bbox.Diag();
 
         for (int s=0;s<nstep;s++)
         {
@@ -622,15 +624,20 @@ public:
 
             for (size_t i=0;i<poly_m.vert.size();i++)
             {
-                CoordType testPos=poly_m.vert[i].P();
-                CoordType closestPt;
-                ScalarType minDist;
-                FaceType *f=NULL;
-                CoordType norm,ip;
+                TriCoordType testPos;
+                testPos.Import(poly_m.vert[i].P());
+                TriCoordType closestPt;
+                TriScalarType minDist;
+                TriFaceType *f=NULL;
+                TriCoordType norm,ip;
                 f=vcg::tri::GetClosestFaceBase(tri_mesh,grid,testPos,MaxD,minDist,closestPt,norm,ip);
+                CoordType closestImp;
+                closestImp.Import(closestPt);
                 poly_m.vert[i].P()=poly_m.vert[i].P()*Damp+
-                        closestPt*(1-Damp);
-                poly_m.vert[i].N()=norm;
+                                   closestImp*(1-Damp);
+                CoordType normalImp;
+                normalImp.Import(norm);
+                poly_m.vert[i].N()=normalImp;
             }
         }
 
@@ -648,7 +655,7 @@ public:
         vcg::tri::UpdateNormal<TempMesh>::PerVertexNormalizedPerFace(GuideSurf);
         vcg::tri::UpdateTopology<TempMesh>::FaceFace(GuideSurf);
         vcg::tri::UpdateFlags<TempMesh>::FaceBorderFromFF(GuideSurf);
-        LaplacianReproject(poly_m,GuideSurf,nstep,Damp=0.5);
+        LaplacianReproject<TempMesh>(poly_m,GuideSurf,nstep,Damp=0.5);
     }
 
     static void Laplacian(PolyMeshType &poly_m,