From e56eefbd2cec1decf1b506835fd40cec64b53a1d Mon Sep 17 00:00:00 2001
From: Robert Smallshire <robert@smallshire.org.uk>
Date: Fri, 30 Jan 2015 14:36:35 +0100
Subject: [PATCH] Return double-precision IEEE float values when reading
 single-precision IBM floats to prevent loss of precision.

---
 segpy/toolkit.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/segpy/toolkit.py b/segpy/toolkit.py
index a3223ff..8b56eec 100644
--- a/segpy/toolkit.py
+++ b/segpy/toolkit.py
@@ -454,7 +454,7 @@ def unpack_ibm_floats(data, count):
     Returns:
         A sequence of floats.
     """
-    return array('f', (ibm2ieee(data[i: i+4]) for i in range(0, count * 4, 4)))
+    return array('d', (ibm2ieee(data[i: i+4]) for i in range(0, count * 4, 4)))
 
 
 def unpack_values(buf, count, item_size, fmt, endian='>'):