From d49f0ed1d0de677dcbf8b98140ba78df618bb7af Mon Sep 17 00:00:00 2001
From: Thomas Fischer <thomas.fischer@ufz.de>
Date: Wed, 19 Dec 2012 12:40:45 +0100
Subject: [PATCH] Set delta to a small value.
In case delta is really zero, set it to a small value to omit division by zero.
---
GeoLib/Grid.h | 3 +++
1 file changed, 3 insertions(+)
diff --git a/GeoLib/Grid.h b/GeoLib/Grid.h
index 9011e555f02..21e161b1fe0 100644
--- a/GeoLib/Grid.h
+++ b/GeoLib/Grid.h
@@ -122,6 +122,9 @@ public:
// some frequently used expressions to fill the grid vectors
for (std::size_t k(0); k < 3; k++) {
+ if (fabs(delta[k]) < std::numeric_limits<double>::epsilon()) {
+ delta[k] = std::numeric_limits<double>::epsilon();
+ }
_step_sizes[k] = delta[k] / _n_steps[k];
_inverse_step_sizes[k] = 1.0 / _step_sizes[k];
}
--
GitLab