diff --git a/range.cpp b/range.cpp
index 83e51d74..4155a323 100644
--- a/range.cpp
+++ b/range.cpp
@@ -29,7 +29,7 @@
 #include <cstdint>
 
 
-//-----------------------------------------------------------------------------
+///////////////////////////////////////////////////////////////////////////////
 template <typename T>
 util::range<T>::range (T _min, T _max):
         min (_min),
@@ -39,7 +39,16 @@ util::range<T>::range (T _min, T _max):
 }
 
 
-//-----------------------------------------------------------------------------
+///////////////////////////////////////////////////////////////////////////////
+template <typename T>
+T
+util::range<T>::magnitude (void) const
+{
+    return max - min;
+}
+
+
+///////////////////////////////////////////////////////////////////////////////
 template <typename T>
 bool
 util::range<T>::contains (T val) const
@@ -57,28 +66,6 @@ util::range<T>::contains (const range <T> &r) const
 }
 
 
-//-----------------------------------------------------------------------------
-template <typename T>
-void
-util::range<T>::sanity (void) const
-{
-    CHECK (min <= max);
-}
-
-
-//-----------------------------------------------------------------------------
-namespace util {
-    template <>
-    void
-    range<double>::sanity (void) const
-    {
-        if (std::isnan (min) || std::isnan (max))
-            return;
-        CHECK (min <= max);
-    }
-}
-
-
 //-----------------------------------------------------------------------------
 template <typename T>
 T
@@ -100,7 +87,7 @@ util::range<T>::expand (T val)
 }
 
 
-//-----------------------------------------------------------------------------
+///////////////////////////////////////////////////////////////////////////////
 template <typename T>
 util::range<T>&
 util::range<T>::operator*= (T val)
@@ -121,7 +108,7 @@ util::range<T>::operator* (T val) const
 }
 
 
-//-----------------------------------------------------------------------------
+///////////////////////////////////////////////////////////////////////////////
 namespace util {
     template <>
     double
@@ -150,7 +137,7 @@ util::range<T>::random (void) const
 }
 
 
-//-----------------------------------------------------------------------------
+///////////////////////////////////////////////////////////////////////////////
 namespace util {
     template <>
     bool
@@ -180,7 +167,7 @@ util::range<T>::operator ==(const util::range<T> &rhs) const
 }
 
 
-//-----------------------------------------------------------------------------
+///////////////////////////////////////////////////////////////////////////////
 template <typename T> 
 const util::range<T>
 util::range<T>::UNLIMITED (std::numeric_limits<T>::has_infinity ? -std::numeric_limits<T>::infinity () :
@@ -201,6 +188,28 @@ const util::range<T>
 util::range<T>::UNIT (0.0, 1.0);
 
 
+///////////////////////////////////////////////////////////////////////////////
+template <typename T>
+void
+util::range<T>::sanity (void) const
+{
+    CHECK (min <= max);
+}
+
+
+//-----------------------------------------------------------------------------
+namespace util {
+    template <>
+    void
+    range<double>::sanity (void) const
+    {
+        if (std::isnan (min) || std::isnan (max))
+            return;
+        CHECK (min <= max);
+    }
+}
+
+
 //-----------------------------------------------------------------------------
 namespace util {
     template struct range<double>;
diff --git a/range.hpp b/range.hpp
index 62082698..d8910182 100644
--- a/range.hpp
+++ b/range.hpp
@@ -36,6 +36,8 @@ namespace util {
 
         range (T _min, T _max);
 
+        T magnitude (void) const;
+
         /// Check whether value falls within this range (inclusive)
         bool contains (T val) const;
         /// Check whether a range falls completely within (inclusive) this range