From f3f291c3c4411f9e73b59bfff94af0be60b43934 Mon Sep 17 00:00:00 2001
From: David Kastrup <dak@gnu.org>
Date: Mon, 29 Jun 2015 10:48:29 +0200
Subject: [PATCH] Robustify Offset::direction and Offset::length

---
 flower/offset.cc | 14 ++++++++++++--
 1 file changed, 12 insertions(+), 2 deletions(-)

diff --git a/flower/offset.cc b/flower/offset.cc
index 0554e6b2e9..7cb2677481 100644
--- a/flower/offset.cc
+++ b/flower/offset.cc
@@ -90,8 +90,7 @@ Offset::angle_degrees () const
 Real
 Offset::length () const
 {
-  return sqrt (sqr (coordinate_a_[X_AXIS])
-               + sqr (coordinate_a_[Y_AXIS]));
+  return hypot (coordinate_a_[X_AXIS], coordinate_a_[Y_AXIS]);
 }
 
 bool
@@ -107,6 +106,17 @@ Offset
 Offset::direction () const
 {
   Offset d = *this;
+  if (isinf (d[X_AXIS]))
+    {
+      if (!isinf (d[Y_AXIS]))
+        return Offset ((d[X_AXIS] > 0.0 ? 1.0 : -1.0), 0.0);
+    }
+  else if (isinf (d[Y_AXIS]))
+    return Offset (0.0, (d[Y_AXIS] > 0.0 ? 1.0 : -1.0));
+  else if (d[X_AXIS] == 0.0 && d[Y_AXIS] == 0.0)
+    return d;
+  // The other cases propagate or produce NaN as appropriate.
+
   d /= length ();
   return d;
 }
-- 
2.39.5