diff options
Diffstat (limited to 'forcescheme.cpp')
-rw-r--r-- | forcescheme.cpp | 3 |
1 files changed, 2 insertions, 1 deletions
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; |