__author__ = 'chunk'

import struct

import numpy as np

import math

from  common import *

def float2bits(f, fmt='bin'):
    if fmt == 'hex':
        return hex(struct.unpack('!l', struct.pack('!f', f))[0])
    return bin(struct.unpack('!l', struct.pack('!f', f))[0])


def double2bits(d, fmt='bin'):
    if fmt == 'hex':
        return hex(struct.unpack('!q', struct.pack('!d', d))[0])
    return bin(struct.unpack('!q', struct.pack('!d', d))[0])


def float2bin(num):
    # http://stackoverflow.com/questions/16444726/binary-representation-of-float-in-python-bits-not-hex
    bits = [bin(ord(c)).replace('0b', '').rjust(8, '0') for c in struct.pack('!f', num)]
    print bits
    return ''.join(bits)


def double2bin(num):

    bits = [bin(ord(c)).replace('0b', '').rjust(8, '0') for c in struct.pack('!d', num)]
    print bits
    return ''.join(bits)

def infinite_limit():

    s = np.float32(0)
    tmp = np.float32(-1)

    i = 0

    while s != tmp:

        i += 1
        tmp = s

        s += np.float32(1.0 / i)

    print i, s, tmp, float2bin(s)

def infinite_float(n):
    s = np.float32(0)
    for i in range(1, n + 1):
        s += np.float32(1.0 / i)
    return s

def infinite_double(n):
    s = np.float64(0)
    for i in range(1, n + 1):
        s += np.float64(1.0 / i)
    return s


def testn0():
    for i in range(100):
        print i, np.log(i)
        if np.log(i) >= 8:
            print i

if __name__ == '__main__':

    # testn0()

    infinite_limit()

    for i in range(1,11):
        print infinite_float(i)
    for i in range(1,11):
        print infinite_double(i)

    # i = 1
    # while infinite_double(i) < 16:
    #     i += 1
    # print i, infinite_double(i)

    timer = Timer()

    for i in range(1,8):
        n = math.pow(10,i)
        timer.mark()
        infinite_double(int(n)) # 1125899906842624
        timer.report()

    pass