From dc50effe61c833da9c15b886d5ebba65620b99f0 Mon Sep 17 00:00:00 2001
From: Samuel Fadel <samuelfadel@gmail.com>
Date: Tue, 2 Feb 2016 17:42:04 -0200
Subject: Using std::numeric_limits where appropriate.

---
 forcescheme.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

(limited to 'forcescheme.cpp')

diff --git a/forcescheme.cpp b/forcescheme.cpp
index 22ab293..9e3feaa 100644
--- a/forcescheme.cpp
+++ b/forcescheme.cpp
@@ -1,6 +1,7 @@
 #include "mp.h"
 
 #include <algorithm>
+#include <limits>
 
 static const double EPSILON = 1e-3;
 
@@ -18,7 +19,7 @@ arma::mat mp::forceScheme(const arma::mat &D,
         i[k] = j[k] = k;
     }
 
-    double prevDeltaSum = 1. / 0.;
+    double prevDeltaSum = std::numeric_limits<double>::infinity();
     for (size_t iter = 0; iter < maxIter; iter++) {
         double deltaSum = 0;
 
-- 
cgit v1.2.3