From 296a644550aec58bcab0a7c55e486da19eef78f1 Mon Sep 17 00:00:00 2001
From: Pablo Garzon <gapablo2001@gmail.com>
Date: Fri, 28 Jul 2023 11:18:14 +0200
Subject: [PATCH] #3501 Geodesic Deformation

---
 packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx b/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx
index 9bcf5fe..9855950 100644
--- a/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx
+++ b/packages/itkvtk/src/bbitkvtkGeodesicMeshDeformation.cxx
@@ -32,6 +32,7 @@ void GeodesicMeshDeformation::Process()
     using MeshType      = itk::QuadEdgeMesh<double, 3>;
     std::vector<double> deformInfo;
 	bbSetOutputOut(deformInfo);
+	//std::vector<double> displacementVector;
     
     //Set up QuadEdge and filter every time polydata changes
     if ((bbGetInputIn() != polydata) && (bbGetInputActive()==true) && (bbGetInputIn() != NULL))
@@ -215,8 +216,10 @@ void GeodesicMeshDeformation::Process()
 						points->SetPoint(BegProcessedIt.Value().GetNode(), np);
 					}
 					++BegProcessedIt;
-				}				
-				std::vector<double> info{pb[0],pb[1],pb[2], (double) bbGetInputEdgeId(), s};
+				}
+				//double directionMoved[3] = {lstCenter[0]-displcement_x, lstCenter[1]-displcement_y, lstCenter[2]-displcement_z};
+				//vtkMath::Normalize(directionMoved);
+				std::vector<double> info{lstCenter[0],lstCenter[1],lstCenter[2], bbGetInputDirection()[0] , bbGetInputDirection()[1], bbGetInputDirection()[2],(double) bbGetInputEdgeId(), s};
 				bbSetOutputOut(info);
 				cout << "info updated" << endl;
 				points->Modified();
-- 
2.47.1