gistfile1.txt

import numpy as np
import time as tt

def main():
    timing = np.array([])
    log2N = np.arange(6) + 8
    for N in 2 ** log2N:
        uu = np.random.normal(size=(N, N+2))
        matrix = np.dot(uu, uu.T)
        print matrix.shape
        vector = np.random.normal(size=N)
        print "start %d" % N
        t1 = tt.time()
        np.linalg.solve(matrix, vector)
        timing = np.append(timing, [tt.time() - t1,])
        print "end %d" % N
    print timing
    print np.polyfit(log2N, np.log(timing)/np.log(2.), 1)

if __name__ == "__main__":
    main()