From 46f11a0a29e9a11679d9b84f798594115429c906 Mon Sep 17 00:00:00 2001
From: Christina Vaz <cvaz@andrew.cmu.edu>
Date: Mon, 18 Jun 2018 19:11:51 -0400
Subject: [PATCH] fix for multiple sources

---
 .../include/CGAL/Heat_method_3/Heat_method_3.h     | 14 ++++++++------
 1 file changed, 8 insertions(+), 6 deletions(-)

diff --git a/Heat_method_3/include/CGAL/Heat_method_3/Heat_method_3.h b/Heat_method_3/include/CGAL/Heat_method_3/Heat_method_3.h
index 60ca5647cfb..65473b4933e 100644
--- a/Heat_method_3/include/CGAL/Heat_method_3/Heat_method_3.h
+++ b/Heat_method_3/include/CGAL/Heat_method_3/Heat_method_3.h
@@ -28,7 +28,7 @@
 #include <set>
 
 #include <CGAL/property_map.h>
-
+#include <CGAL/double.h>
 #include <Eigen/Cholesky>
 #include <Eigen/Sparse>
 
@@ -357,7 +357,7 @@ namespace Heat_method_3 {
     Eigen::VectorXd value_at_source_set(const Eigen::VectorXd& phi, int dimension) const
     {
       Eigen::VectorXd source_set_val(dimension,1);
-      if(!(sources.empty()))
+      if(sources.empty())
       {
         for(int k = 0; k<dimension; k++)
         {
@@ -369,7 +369,7 @@ namespace Heat_method_3 {
       {
         for(int i = 0; i<dimension; i++)
         {
-          int min_val = INT_MAX;
+          double min_val = INT_MAX;
           vertex_iterator current;
           Index current_Index;
           current = sources.begin();
@@ -377,9 +377,10 @@ namespace Heat_method_3 {
           for(int j = 0; j<sources.size(); j++)
           {
             current_Index = get(vertex_id_map, *current);
-            if(std::abs(phi.coeff(i,0)-phi.coeff(current_Index,0))< min_val)
+            double new_d = fabs(-phi.coeff(current_Index,0)+phi.coeff(i,0));
+            if(new_d < min_val)
             {
-              min_val = std::abs(phi.coeff(i,0)-phi.coeff(current_Index,0));
+              min_val = new_d;
             }
             current = ++current;
           }
@@ -422,7 +423,8 @@ namespace Heat_method_3 {
       build();
       BOOST_FOREACH(vertex_descriptor vd, vertices(tm)){
         Index i_d = get(vertex_id_map, vd);
-        d =solved_phi(i_d,0);
+        d = solved_phi(i_d,0);
+        std::cout<<d<<"\n";
         put(vdm,vd,d);
       }
     }