view run.py @ 5:d7b2784ff5a3

* Fixed some miscalculations from the way python does float division
* Made data files more precise

import cortical_model as cm
import numpy as np
import utils

t = np.array(range(44100))
sine = np.sin(2*np.pi*t*1000/44100)
sine1 = np.sin(2*np.pi*t*100/44100)
x = np.append(sine, sine1)

exc_i = cm.get_excitation_i(x, 44100, 60, verbose = True)
sl = cm.get_specific_loudness(exc_i)

#for i in range(39):
#    utils.plot_waveform(sl[i])

print "Loudness"
loudness = 2*np.sum(sl,0)

utils.plot_waveform(loudness)