setup_trailer = '''
import random, math
width,height = 512,512
values = range(256)
pixels = [ tuple( random.choice(values) for _ in xrange(3)) 
    for _ in xrange(width*height) ]
'''


distance_1_setup = '''
def distance(x,y,z,x1,y1,z1):
    dist = math.sqrt((x-x1)*(x-x1) + (y-y1)*(y-y1) + (z-z1)*(z-z1))
    return dist

''' + setup_trailer

distance_2_setup = '''
def distance(x,y,z,x1,y1,z1):
    dist = (x-x1)*(x-x1) + (y-y1)*(y-y1) + (z-z1)*(z-z1)
    return dist

''' + setup_trailer

distance_3_setup = '''
def distance(x,y,z,x1,y1,z1):
    return (
        (x1-x)**2 +
        (y1-y)**2 +
        (z1-z)**2 )

''' + setup_trailer



statement = '''
index = 0
for x in xrange(width):
    for y in xrange(height):
        pixel = pixels[index]
        distance(pixel[0],pixel[1],pixel[2],127,127,127)
        index += 1
'''

import timeit

print timeit.timeit(statement,distance_1_setup,number=20)
print timeit.timeit(statement,distance_2_setup,number=20)
print timeit.timeit(statement,distance_3_setup,number=20)

tripwire fucked around with this message at 04:32 on Jun 22, 2011

Computer viking fucked around with this message at 19:04 on Jun 22, 2011

>>> import pyximport; pyximport.install()
>>> distance_4_setup="from distance_test import distance;" + setup_trailer
>>> print timeit.timeit(statement,distance_4_setup,number=20)
1.72920703888

def distance(int x, int y, int z, int x1, int y1, int z1):
    dist = (x-x1)*(x-x1) + (y-y1)*(y-y1) + (z-z1)*(z-z1)
    return dist

Computer viking fucked around with this message at 19:34 on Jun 22, 2011

import numpy as np

infile = open("distance.txt", "r")
header = infile.readline().strip().split("\t")
header_len = len(header)
data = []
rownames = []
for line in infile:
	parts = line.strip().split("\t")
	if len(parts) < header_len:
		break
	rownames.append(parts[0])
	values = [int(p) for p in parts[1:] ]
	data.append(values)

data_array = np.array(data)

import operator
sum_of_distances = map(sum,data_array)
name_with_distance = zip(header,sum_of_distances)
nwd_sorted = sorted(name_with_distance, key=operator.itemgetter(1))

>>> header
['a', 'b', 'c', 'd', 'e']
>>> data_array
array([[0, 1, 1, 5, 4],
       [1, 0, 2, 5, 4],
       [1, 2, 0, 3, 3],
       [5, 5, 3, 0, 5],
       [4, 4, 3, 5, 0]])
>>> name_with_distance
[('a', 11), ('b', 12), ('c', 9), ('d', 18), ('e', 16)]
>>> nwd_sorted
[('c', 9), ('a', 11), ('b', 12), ('e', 16), ('d', 18)]

Computer viking fucked around with this message at 20:17 on Jun 24, 2011

word_items[words] = word_counts.items()

from collections import defaultdict
word_tally = defaultdict(int)
#a default dict takes a factory function as an argument. whenever you lookup a key
#which isn't in the dictionary, it uses that factory function to make a value for
#the key. In this case, int() is used to return the integer value 0.

with open('gettysburg.txt', 'r') as myfile:
    for line in myfile:
        for word in line.strip().split(' '):
            word_tally[word] += 1

for word in sorted( word_tally.keys() ):
    print word, word_tally[word]

tripwire fucked around with this message at 17:48 on Jun 26, 2011