"""Utilities for handling IEEE 754 floating point special values

This python module implements constants and functions for working with
IEEE754 double-precision special values.  It provides constants for
Not-a-Number (NaN), Positive Infinity (PosInf), and Negative Infinity
(NegInf), as well as functions to test for these values.

The code is implemented in pure python by taking advantage of the
'struct' standard module. Care has been taken to generate proper
results on both big-endian and little-endian machines. Some efficiency
could be gained by translating the core routines into C.

See <http://babbage.cs.qc.edu/courses/cs341/IEEE-754references.html>
for reference material on the IEEE 754 floating point standard.

Further information on this package is available at
<http://software.biostat.washington.edu/statsoft/snake/fpconst>.

Author:    Gregory R. Warnes <gregory_r_warnes@groton.pfizer.com>
Date::     2003-04-08
Copyright: (c) 2003, Pfizer, Inc.
"""

__version__ = "0.6.0"
ident = "$Id: fpconst.py,v 1.8 2003/05/12 15:14:00 warnes Exp $"

import struct

# check endianess
_big_endian = struct.pack('i',1)[0] != '\x01'

# and define appropriate constants
if(_big_endian): 
    _HW = 0
    _LW = 1
    
    NaN = struct.unpack('d', '\x7F\xFF\xFF\xFF\xFF\xFF\xFF\xFF')[0]
    PosInf = struct.unpack('d', '\x7F\xF0\x00\x00\x00\x00\x00\x00')[0]
    NegInf = -PosInf
    
else:
    _HW = 1
    _LW = 0

    NaN = struct.unpack('d', '\x00\x00\x00\x00\x00\x00\xf8\xff')[0]
    PosInf = struct.unpack('d', '\x00\x00\x00\x00\x00\x00\xf0\x7f')[0]
    NegInf = -PosInf

def _double_as_longs(dval):
    "Use struct.unpack to decode a double precision float into two longs"
    tmp = struct.unpack('ll',struct.pack('d', dval))
    return (tmp[_HW], tmp[_LW])


##
## Functions to extract components of the IEEE 754 floating point format
##

def _sign(dval):
    "Extract the sign bit from a double-precision floating point value"
    ll = _double_as_longs(dval)
    return ll[0] >> 31 & 0x01 

def _exponent(dval):
    """Extract the exponentent bits from a double-precision floating
    point value.

    Note that for normalized values, the exponentdent bits have an offset
    of 1023. As a consequence, the actual exponentent is obtained
    by subtracting 1023 for the value returned by this function
    """
    ll = _double_as_longs(dval)
    return (ll[0] >> 20) & 0x7ff

def _mantissa(dval):
    """Extract the _mantissa bits from a double-precision floating
    point value."""

    ll = _double_as_longs(dval)
    mantissa0 = (ll[0] & 0x000fffffL) << 32
    mantissa1 = ll[1] 
    return mantissa0 + mantissa1

##
## Functions to test for IEEE 754 special values
##

def isNaN(value):
    "Determine if the argument is a IEEE 754 NaN (Not a Number) value."
    return (_exponent(value)==0x7ff and _mantissa(value)!=0) 

def isInf(value):
    """Determine if the argument is an infinite IEEE 754 value (positive
    or negative inifinity)"""
    return (_exponent(value)==0x7ff and _mantissa(value)== 0)

def isFinite(value):
    """Determine if the argument is an finite IEEE 754 value (i.e., is
    not NaN, positive or negative inifinity)"""
    return (_exponent(value)!=0x7ff)


def isPosInf(value):
    "Determine if the argument is a IEEE 754 positive infinity value"
    return (_sign(value)==0 and _exponent(value)==0x7ff and \
            _mantissa(value)== 0) 

def isNegInf(value):
    "Determine if the argument is a IEEE 754 negative infinity value"
    return (_sign(value)==1 and _exponent(value)==0x7ff and \
            _mantissa(value)== 0) 
 23:55:12 +0000'>2009-01-24</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=369df82d012294a89abb51293855c11b5a8f54eb'>tcltk: Fix the native packages after last tcltk change</a></td><td>Tim Ellis</td><td>2</td></tr>
<tr><td><span title='2009-01-24 21:23:09 +0100'>2009-01-24</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=2d441c0f1b596a90dc4b8ecac238b435132de4c1'>tcltk: 8.4.11 -&gt; 8.4.19; try harder to get proper library packaging</a></td><td>Michael 'Mickey' Lauer</td><td>5</td></tr>
<tr><td><span title='2009-01-01 23:35:15 +0000'>2009-01-01</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=3129426514932a31dc4820d72fbe21364b8055fa'>tk-native: Reordered lines to fix dependencies.</a></td><td>Stanislav Brabec</td><td>1</td></tr>
<tr><td><span title='2008-12-16 15:21:11 +0000'>2008-12-16</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=6c3031554f0ffa98fef3032af1c906c19d3b538f'>tcl-native, tk-native: New packages. Needed for compilation of geda/pcb.</a></td><td>Stanislav Brabec</td><td>2</td></tr>
<tr><td><span title='2008-10-18 18:02:07 +0200'>2008-10-18</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=c0b6dd9b41279ff00c5a90e49108945a0fee9cdd'>bb files: Revert the introduction of FILE_PR</a></td><td>Holger Hans Peter Freyther</td><td>2</td></tr>
<tr><td><span title='2008-10-15 21:41:43 +0200'>2008-10-15</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=aa03004ace2fabb135a6208ef8c2d2b312aa7b0b'>[PR] Change PR to FILE_PR and intro PR as FILE_PR + DISTRO_PR</a></td><td>Holger Hans Peter Freyther</td><td>2</td></tr>
<tr><td><span title='2008-10-15 15:27:49 +0200'>2008-10-15</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=869ba9e1fb2237b75560feed13545a35f3c3569d'>mtn2git: remove .mtn2git_empty files in non-empty directories</a></td><td>Jan Luebbe</td><td>4</td></tr>
<tr><td><span title='2008-02-21 19:38:56 +0000'>2008-02-21</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=6f9b918ba4c22fec0628869ab323f43d7cc9ce5e'>tcl: drop -L${libdir} in tclConfig.sh.  Fixes 3771.  Thanks, zecke.</a></td><td>Rolf Leggewie</td><td>1</td></tr>
<tr><td><span title='2007-01-06 05:34:00 +0000'>2007-01-06</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=f90462217bda8121fd05da79f862c0ec93c8dde5'>tcl: make sure ${STAGING_BINDIR_CROSS} exists before trying to install tclCon...</a></td><td>Rod Whitby</td><td>1</td></tr>
<tr><td><span title='2006-11-29 11:34:04 +0000'>2006-11-29</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=081432e434818a616a523f9f2bb18e707e4e68e3'>Introduce STAGING_BINDIR_CROSS and STAGING_BINDIR_NATIVE as discussed on the ...</a></td><td>Richard Purdie</td><td>2</td></tr>
<tr><td><span title='2006-09-29 20:32:52 +0000'>2006-09-29</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=c62509d5242fcec2865b0115aae60435239db294'>tk: fix dirs in tkConfig.sh so tk related stuff will be buildable after clean...</a></td><td>Marcin Juszkiewicz</td><td>1</td></tr>
<tr><td><span title='2006-09-29 20:32:18 +0000'>2006-09-29</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=0d4e5824544fd99dff585af2a0075f0be6ca709f'>tcl: fix dirs in tclConfig.sh so tcl related stuff will be buildable after cl...</a></td><td>Marcin Juszkiewicz</td><td>1</td></tr>
<tr><td><span title='2006-07-29 01:56:36 +0000'>2006-07-29</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=1e2a6ff6ffac1ac6b3005f716bb5c11395396434'>x11: change virtual/x11 to virtual/libx11</a></td><td>Justin Patrin</td><td>1</td></tr>
<tr><td><span title='2006-07-26 17:07:51 +0000'>2006-07-26</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=175e2780753ef215331518d0c555657af95a48fa'>x11: switch to virtual/x11 so that diet or full x11 can be selected</a></td><td>Justin Patrin</td><td>1</td></tr>
<tr><td><span title='2006-07-24 16:21:42 +0000'>2006-07-24</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=25b369f82c032ee41b33d52a55f204fcce0b11f0'>tcl - Fix for a crazy configure script in TCL.</a></td><td>Raymond Danks</td><td>1</td></tr>
<tr><td><span title='2006-05-13 12:38:32 +0000'>2006-05-13</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=743d7e938e0e052bf169e30ab422bd4fc8acbe68'>tk: depend on libxt - close #966</a></td><td>Marcin Juszkiewicz</td><td>1</td></tr>
<tr><td><span title='2006-04-16 14:57:33 +0000'>2006-04-16</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=71b52e6a24560bb80742794c15a2495b1400c86e'>tcl/tk: add pnum=2 to make the latest patch apply</a></td><td>Koen Kooi</td><td>2</td></tr>
<tr><td><span title='2006-04-16 14:35:51 +0000'>2006-04-16</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=256e671e1944961ca4e19c13684cec34da1966bf'>tcltk: apply patch to fix a bogus line in the configure scripts</a></td><td>Michael Lauer</td><td>6</td></tr>
<tr><td><span title='2006-03-01 18:21:32 +0000'>2006-03-01</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=2fb5f01e23b04e171d4dea6e5d1af5640386efac'>all over the place: adapt DEPENDS to new xlib names</a></td><td>Philipp Zabel</td><td>1</td></tr>
<tr><td><span title='2005-09-26 22:55:18 +0000'>2005-09-26</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=8682e21eca56ef918222942b53694e8333842406'>disapproval of revision 05dd4ce9d2ee0995688e45e60a21534114fbc45f</a></td><td>Michael Lauer</td><td>1</td></tr>
<tr><td><span title='2005-09-26 17:05:54 +0000'>2005-09-26</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=2888db41c577bcce515ca6f0767ddcb3604d265a'>tk: DEPEND on libx11, not x11</a></td><td>Michael Lauer</td><td>1</td></tr>
<tr><td><span title='2005-09-25 14:28:40 +0000'>2005-09-25</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=c5ad04d9323ce0152ec164ec4f38af0004ed1b0d'>tcl/tk: add SONAME patches, ship .so in $PN, not $PN-dev</a></td><td>Michael Lauer</td><td>4</td></tr>
<tr><td><span title='2005-09-25 00:45:57 +0000'>2005-09-25</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=5f4f99dd916ac86e40b58decb2a0163067edce58'>tcl/tk: install header files into $STAGING_DIR</a></td><td>Michael Lauer</td><td>2</td></tr>
<tr><td><span title='2005-09-24 22:25:03 +0000'>2005-09-24</span></td><td><a href='/cgit/multitech-oe.git/commit/packages/tcltk?id=68d6601e8f93243347b58324ce3d4f02eb3a84a9'>add tcl (tool command language) and the tk (toolkit) extension</a></td><td>Michael Lauer</td><td>5</td></tr>