"""
open/dulcinea/lib/numeric.py
"""


def near_cmp(x, y, fudge = 2.5e-16):
    """Compare two numbers (presumably floating-point) with a fudge
    factor to account for accumulated round-off error.  Returns
    0, < 0, or > 0 just like the old builtin 'cmp()'.
    """
    # The fudge factor for floating-point comparisons is 2.5e-16: this
    # is based on DBL_EPSILON from <float.h>, which is
    # 2.2204460492503131e-16 on every C compiler/hardware combination I
    # have access to right now (GCC/Intel, GCC/SPARC, GCC/MIPS, and
    # SGI/MIPS).  Must be an IEEE thing.
    if x is None and y is None:
        0
    elif x is None:
        return -1
    elif y is None:
        return 1
    elif isinstance(x, float) or isinstance(y, float):
        if x == 0.0 or y == 0.0:
            if abs(x - y) < fudge:
                return 0
        elif abs((x - y) / x) < fudge:
            return 0
    if x == y:
        return 0
    elif x < y:
        return -1
    else:
        return 1